Method and system for determining target geologic body in mining area three-dimensional geologic modeling

Abstract

The invention provides a method and system for determining a target geologic body in mining area three-dimensional geological modeling, and belongs to the technical field of geological exploration, and the method comprises the following steps: S1, obtaining mining area geological mineral data; s2, based on the stratum age, the lithology combination and the structural occurrence, a preliminary geologic body is divided in combination with actually measured geological profile data; s3, collecting a physical property measurement sample of the preliminary geologic body; and S4, merging or splitting the preliminary geologic body according to the physical property parameter difference threshold value and the distinguishable degree. According to the method and system for determining the target geologic body in three-dimensional geological modeling of the mining area, samples are collected through a physical property parameter analysis module, density and other parameters are measured, a mean value, a standard deviation and the like are calculated, differences are analyzed through T inspection and a distinguishable degree matrix, and a target geologic body determination module combines or splits a primary geologic body according to a difference threshold value and distinguishable degrees; quantitative analysis is achieved, and the problem that boundaries are difficult to distinguish due to physical property parameter overlapping is solved.

Description

Technical Field

[0001] This application relates to the field of geological exploration technology, specifically to a method and system for determining target geological bodies in three-dimensional geological modeling of mining areas. Background Technology

[0002] In existing technologies, the determination of target geological bodies mainly relies on the surface characteristics of geological age and lithology, lacking quantitative difference analysis of physical property parameters such as density and magnetic susceptibility. This makes it difficult to accurately distinguish the boundaries of adjacent geological bodies due to the overlap of physical property parameters. At the same time, the lack of a multi-source data collaborative constraint mechanism results in insufficient spatial coupling between measured geological profiles, borehole cores and geophysical data. In structurally complex areas, the boundary of the target geological body is easily misjudged due to the deviation of a single data point, resulting in significant differences between the three-dimensional model and the actual geological conditions, increasing drilling verification costs and exploration risks. Summary of the Invention

[0003] In view of the technical problems existing in the background art, this application provides a method and system for determining target geological bodies in three-dimensional geological modeling of mining areas. The method involves collecting samples and measuring parameters such as density through a physical property parameter analysis module, calculating the mean, standard deviation, etc., and using T-test and distinguishability matrix to analyze differences. The target geological body determination module merges or splits preliminary geological bodies based on difference thresholds and distinguishability to achieve quantitative analysis, thereby achieving the technical effect of high division accuracy and close integration of multi-source data.

[0004] In a first aspect, embodiments of this application provide a method for determining target geological bodies in three-dimensional geological modeling of a mining area, the specific steps of which include: Obtain geological and mineral data and measured geological profile data of the mining area; Based on the geological and mineral data of the mining area and the measured geological profile data, the geological bodies are divided, and a preliminary geological body division table is obtained; Collect physical property samples of the initial geological body, process the physical performance parameters of the measured physical property samples, calculate the frequency distribution histogram of the physical performance parameters, and determine the peak interval and overlapping area; The T-test method was used to analyze and determine the differences in physical properties among geological bodies; Establish a parameter distinguishability matrix for the physical performance parameters, and determine the overlap of the physical performance parameters between geological bodies based on the parameter distinguishability matrix; Set a threshold for the difference in physical property parameters, and process the geological bodies in the preliminary geological body division table according to the difference and overlap of physical property parameters to obtain the target geological bodies for three-dimensional geological modeling of the mining area.

[0005] Furthermore, in this embodiment, the physical performance parameters include density, magnetic susceptibility, resistivity, and longitudinal wave velocity; the threshold values ​​for the differences in these physical property parameters include a density difference of 0.08-0.1 g / cm³ and a magnetic susceptibility difference of 15 × 10⁻⁶ g / cm³. -6 -20×10 -6 SI, resistivity difference 120-150 ohm-meters, wave velocity difference 400-500 meters per second, adjacent geological bodies in the preliminary geological body division table whose physical property parameters are all lower than the physical property parameter difference threshold and whose overlap is not less than 40%-50% are merged into one target geological body, and any two adjacent geological bodies in the preliminary geological body division table whose physical property parameters are higher than the physical property parameter difference threshold and whose overlap is not higher than 20% are divided into independent target geological bodies.

[0006] Furthermore, in this embodiment, the method further includes: after merging two adjacent geological bodies into one target geological body, calculating the standard deviation of the physical performance parameters of the merged target geological body and the mean of the standard deviations of the physical performance parameters of the two preliminary geological bodies corresponding to the target geological body; when the standard deviation of the target geological body is higher than 80% of the mean of the standard deviations of the physical performance parameters of the two preliminary geological bodies, canceling the merging of the target geological body.

[0007] Furthermore, in this embodiment, after dividing two adjacent geological bodies into independent target geological bodies, the boundary data of the target geological bodies and the physical field characteristic data of the mining area are obtained. The physical characteristic data is processed to obtain geophysical anomaly area data. The degree of consistency between the boundary data of the target geological bodies and the geophysical anomaly area data is calculated, and the rationality of the boundary data of the target geological bodies is determined based on the degree of consistency.

[0008] Furthermore, in this embodiment, the method for measuring the physical performance parameters includes: measuring density using the ring cutter method, measuring magnetic susceptibility using a proton precession magnetometer, measuring resistivity using the quadrupole method, and measuring longitudinal wave velocity using an ultrasonic detector.

[0009] Secondly, this application also provides a system for determining target geological bodies in three-dimensional geological modeling of mining areas, including: The data acquisition module is used to acquire geological and mineral data and measured geological profile data of the mining area; The preliminary geological body division module is used to divide geological bodies according to the geological and mineral data of the mining area and the measured geological profile data, and to obtain a preliminary geological body division table; The physical property parameter analysis module is used to collect physical property measurement samples of the initial geological body, process the physical property parameters of the measured physical property samples, calculate the frequency distribution histogram of the physical property parameters, determine the peak interval and overlapping area, use the T-test method to analyze and determine the differences in physical property parameters between geological bodies, establish the parameter distinguishability matrix of the physical property parameters, and determine the overlap of the physical property parameters between geological bodies based on the parameter distinguishability matrix. The target geological body determination module is used to set the threshold for differences in physical property parameters and process the geological bodies in the preliminary geological body division table according to the differences and overlaps in physical property parameters to obtain the target geological bodies for three-dimensional geological modeling of the mining area.

[0010] Furthermore, in this embodiment, in the target geological body determination module, the physical performance parameters include density, magnetic susceptibility, resistivity, and P-wave velocity; the threshold values ​​for the differences in physical property parameters include a density difference of 0.08-0.1 g / cm³ and a magnetic susceptibility difference of 15 × 10⁻⁶ g / cm³. -6 -20×10 -6 SI, resistivity difference 120-150 ohm-meters, wave velocity difference 400-500 meters per second. When the difference in physical property parameters between two adjacent geological bodies in the preliminary geological body division table is lower than the difference threshold of physical property parameters and the overlap is not less than 40%-50%, the adjacent preliminary geological bodies are merged into the same target geological body. When the difference in physical property parameters between any two adjacent geological bodies in the preliminary geological body division table is higher than the difference threshold of physical property parameters and the overlap is not higher than 20%, they are divided into independent target geological bodies.

[0011] Furthermore, in this embodiment, it also includes: a target geological body verification module, used to calculate the standard deviation of the physical performance parameters of the merged target geological body and the mean of the standard deviations of the physical performance parameters of the two preliminary geological bodies corresponding to the target geological body after merging two adjacent geological bodies into one target geological body; when the standard deviation of the target geological body is higher than 80% of the mean of the standard deviations of the physical performance parameters of the two preliminary geological bodies, the merging of the target geological body is cancelled.

[0012] Furthermore, in this embodiment, the target geological body verification module further includes: After dividing two adjacent geological bodies into independent target geological bodies, the boundary data of the target geological bodies and the physical field characteristic data of the mining area are obtained. The physical characteristic data is processed to obtain geophysical anomaly area data. The degree of agreement between the boundary data of the target geological bodies and the geophysical anomaly area data is calculated, and the rationality of the boundary data of the target geological bodies is determined based on the degree of agreement.

[0013] Furthermore, in this embodiment, the method for measuring the physical performance parameters includes: measuring density using the ring cutter method, measuring magnetic susceptibility using a proton precession magnetometer, measuring resistivity using the quadrupole method, and measuring longitudinal wave velocity using an ultrasonic detector.

[0014] Beneficial effects: The present application provides a method and system for determining target geological bodies in three-dimensional geological modeling of mining areas. The method collects samples and measures parameters such as density through a physical property parameter analysis module, calculates the mean and standard deviation, and analyzes the differences using T-test and distinguishability matrix. The target geological body determination module merges or splits preliminary geological bodies based on the difference threshold and distinguishability, thereby achieving quantitative analysis and solving the problem of difficult-to-distinguish boundaries due to overlapping physical property parameters.

[0015] Furthermore, the method and system for determining target geological bodies in the three-dimensional geological modeling of the mining area integrates multi-source geological data to delineate preliminary geological bodies, establishes a multi-source data collaborative constraint mechanism, improves spatial coupling, reduces boundary misjudgment caused by single data deviation, and reduces drilling costs and risks. It solves the problems of existing technologies where the delineation of target geological bodies depends on surface features, lacks analysis of physical property parameters, and has low multi-source data coupling.

[0016] The above description is only an overview of the technical solution of this application. In order to better understand the technical means of this application and to implement it in accordance with the contents of the specification, and to make the above and other objects, features and advantages of this application more obvious and understandable, the following are specific embodiments of this application. Attached Figure Description

[0017] To more clearly illustrate the technical solutions of this application, the accompanying drawings used in this application will be briefly described below. Obviously, the drawings described below are merely some embodiments of this application. For those skilled in the art, other drawings can be obtained from these drawings without any creative effort.

[0018] Figure 1 This application provides a system framework diagram for determining target geological bodies in three-dimensional geological modeling of a mining area. Figure 2 This application provides a flowchart of a method for determining target geological bodies in three-dimensional geological modeling of a mining area. Detailed Implementation

[0019] The embodiments of the technical solution of this application will now be described in detail with reference to the accompanying drawings. These embodiments are only used to more clearly illustrate the technical solution of this application and are therefore merely examples, and should not be used to limit the scope of protection of this application.

[0020] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application pertains; the terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the application; the terms “comprising” and “having”, and any variations thereof, in the specification, claims, and foregoing description of the drawings are intended to cover non-exclusive inclusion.

[0021] In the description of the embodiments of this application, technical terms such as "first" and "second" are used only to distinguish different objects and should not be construed as indicating or implying relative importance or implicitly specifying the number, specific order, or primary and secondary relationship of the indicated technical features. In the description of the embodiments of this application, "multiple" means two or more, unless otherwise explicitly defined.

[0022] In this document, the term "embodiment" means that a particular feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of this application. The appearance of this phrase in various places throughout the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.

[0023] In the description of the embodiments in this application, the term "and / or" is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A existing alone, A and B existing simultaneously, and B existing alone. Additionally, the character " / " in this document generally indicates that the preceding and following related objects have an "or" relationship.

[0024] In the description of the embodiments of this application, the term "multiple" refers to two or more (including two), similarly, "multiple sets" refers to two or more (including two sets), and "multiple pieces" refers to two or more (including two pieces).

[0025] In the description of the embodiments of this application, the technical terms "center," "longitudinal," "lateral," "length," "width," "thickness," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," "clockwise," "counterclockwise," "axial," "radial," and "circumferential" indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are only for the convenience of describing the embodiments of this application and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on the embodiments of this application.

[0026] In the description of the embodiments of this application, unless otherwise expressly specified and limited, technical terms such as "installation," "connection," "joining," and "fixing" should be interpreted broadly. For example, they can refer to a fixed connection, a detachable connection, or an integral part; they can refer to a mechanical connection or an electrical connection; they can refer to a direct connection or an indirect connection through an intermediate medium; they can refer to the internal communication of two components or the interaction between two components. For those skilled in the art, the specific meaning of the above terms in the embodiments of this application can be understood according to the specific circumstances.

[0027] In the 3D geological modeling of a mining area, the determination of the target geological body is achieved by comprehensively analyzing multi-source data such as geology and physical properties, dividing the geological units into those with clear spatial distribution and physical property differences. This provides the core foundation for constructing a 3D model that truly reflects the geological structure of the mining area. Its accuracy directly affects the reliability of subsequent work such as mineral resource reserve estimation, metallogenic regularity analysis, and mining scheme design.

[0028] In existing technologies, the determination of target geological bodies mainly relies on the surface characteristics of geological age and lithology, lacking quantitative difference analysis of physical property parameters such as density and magnetic susceptibility. This leads to difficulties in accurately distinguishing the boundaries of adjacent geological bodies due to overlapping physical property parameters. For example, traditional methods often result in incorrect delineation of strata with resistivity differences of less than 150 ohm-meters. At the same time, the lack of a multi-source data collaborative constraint mechanism results in insufficient spatial coupling between measured geological profiles, borehole cores, and geophysical data. In structurally complex areas, misjudgment of target geological body boundaries can easily occur due to deviations in a single data source, leading to significant differences between the 3D model and the actual geological conditions. This increases drilling verification costs and exploration risks. Therefore, a method and system for determining target geological bodies in 3D geological modeling of mining areas is proposed.

[0029] To address the technical challenge of providing a method and system for determining target geological bodies in 3D geological modeling of mining areas, achieving advantages such as high segmentation accuracy and tight integration of multi-source data, this application provides a method and system for determining target geological bodies in 3D geological modeling of mining areas. The system involves collecting samples and measuring parameters such as density through a physical property parameter analysis module, calculating the mean and standard deviation, and analyzing differences using T-tests and discriminability matrices. The target geological body determination module merges or splits preliminary geological bodies based on difference thresholds and discriminability, achieving quantitative analysis and thus achieving the technical effects of high segmentation accuracy and tight integration of multi-source data.

[0030] Please refer to Figure 1 , Figure 1 This application provides a system framework diagram for determining target geological bodies in three-dimensional geological modeling of a mining area.

[0031] like Figure 1As shown, the system comprises a data acquisition module, a preliminary geological body division module, a physical property parameter analysis module, and a target geological body determination module. The data acquisition module acquires geological and mineral data and measured geological profile data for the mining area. The preliminary geological body division module divides geological bodies based on the geological and mineral data and measured geological profile data, generating a preliminary geological body division table. The physical property parameter analysis module collects physical property measurement samples from the preliminary geological bodies, processes the physical property parameters of these samples, calculates the frequency distribution histogram of the physical property parameters, determines peak intervals and overlapping areas, uses the T-test to analyze and determine the differences in physical property parameters between geological bodies, establishes a parameter distinguishability matrix for physical property parameters, and determines the overlap of physical property parameters between geological bodies based on the parameter distinguishability matrix. The target geological body determination module sets thresholds for differences in physical property parameters and processes the geological bodies in the preliminary geological body division table based on the differences in physical property parameters and the overlap, obtaining the target geological bodies used for three-dimensional geological modeling of the mining area.

[0032] The physical property parameter analysis module collects samples and measures parameters such as density, calculates the mean and standard deviation, and analyzes the differences using T-tests and discriminability matrices. The target geological body determination module merges or splits preliminary geological bodies based on difference thresholds and discriminability, realizing quantitative analysis and solving the problem of difficult-to-distinguish boundaries due to overlapping physical property parameters.

[0033] For example, such as Figure 2 As shown in this embodiment, the specific execution method of a target geological body determination system in three-dimensional geological modeling of a mining area is as follows: S1. Obtain geological and mineral data and measured geological profile data of the mining area; Obtain geological and mineral data for the mining area: First, carry out data collection work, and collect complete geological reports, geological maps at scales of 1:3000-1:5000, structural outline maps, geological maps of the mining area, and mineral deposit exploration reports for the mining area. The focus is on extracting geological and mineral data such as stratigraphic age, lithological assemblage, and structural occurrence.

[0034] Obtain measured geological profile data: Using a total station with an angle measurement accuracy of ±1-±2 seconds and a distance measurement accuracy of ±(1mm+1ppm)-±(2mm+2ppm), profiles are arranged along the main exploration line of the mining area, with the profile spacing controlled at 60-100m. The outcrops of each stratigraphic profile are measured using the total station to obtain the dip angle, dip direction and thickness of the outcrops. Drilling work is carried out using drilling equipment (the core recovery rate of each hole is over 90%-95%), and the depth, lithological stratification and contact relationship of the geological bodies exposed by the drilling are recorded in detail through core logging. The spatial correspondence between surface and deep geological bodies is established by measuring the dip angle, dip direction, and thickness of strata outcrops, as well as the burial depth, lithological stratification, and contact relationships of geological bodies.

[0035] S2. Based on the geological and mineral data of the mining area and the measured geological profile data, the geological bodies are divided and a preliminary geological body division table is obtained; For example, in the preliminary geological body division stage, the strata are sorted out from oldest to youngest, and metamorphic rocks, limestone, volcanic rocks, etc. of different ages are divided into different preliminary geological bodies. Each geological body is recorded in detail with its lithological composition (e.g., the proportion of phyllite and schist in metamorphic rocks of a certain age is different), thickness (calculated by profile data and borehole interpolation, the thickness range may be tens to hundreds of meters), and formation age (combined with the regional geological timeline), generating a division table containing 10-12 preliminary geological bodies. The table clearly defines the code of each body, lithological characteristics and spatial distribution range.

[0036] S3. Collect physical property measurement samples for each geological body in the preliminary geological body division table, measure the physical performance parameters of the physical property measurement samples, and process the physical performance parameters of the physical property measurement samples to determine the differences and overlaps of physical property parameters between geological bodies. For example, specific methods include: S31. Physical property sample collection: Fresh rock samples are selected for each geological body in the preliminary geological body classification table. The number of samples for each geological body is no less than 15-20, covering different degrees of weathering and mineral assemblage types, so as to provide sufficient samples for subsequent physical property parameter measurement.

[0037] S32. Measure the physical properties of samples corresponding to each geological body, including density, magnetic susceptibility, resistivity, and longitudinal wave velocity.

[0038] Density measurement was performed using the ring cutter method. Samples with dimensions of φ40-φ50mm×40-50mm were prepared. An electronic densitometer with an accuracy of ±0.005-±0.01 g / cm³ was used. Each sample was measured 2-3 times and the average value was taken. For example, the average density of a limestone sample from a certain stratum was between 2.65-2.7 g / cm³, with a standard deviation of 0.03-0.05 g / cm³.

[0039] Before magnetic susceptibility measurement, the oxide layer on the sample surface was mechanically polished for pretreatment. A proton precession magnetometer was used to control the ambient magnetic field interference within ±3-±5 nT, achieving a measurement accuracy of ±2×10⁻⁶. -6 -±3×10 -6 SI, the average magnetic susceptibility of a certain volcanic rock sample is 300 × 10⁻⁶. -6 -350×10 -6 SI indicates moderate magnetic characteristics.

[0040] Resistivity measurements were performed using a four-electrode method, with the electrode spacing adjusted according to the thickness of the geological body (15-50 cm for this mining area). A high-precision resistance meter (measuring range 0.5 ohm-meter to 1.5 × 10⁻⁶) was used. 6 (Ohm-meter, accuracy ±3%-±5%), apply a constant voltage of 8-10V to each sample, collect current data to calculate resistivity. The resistivity of a certain metamorphic rock sample is 800-900 ohm-meters, exhibiting high resistivity characteristics.

[0041] An ultrasonic detector was used to measure the longitudinal wave velocity, with the transmission frequency set to 50-60kHz. Coupling agent was applied to both ends of the sample to ensure acoustic coupling effect. The measurement accuracy was ±1%-±2%. The average longitudinal wave velocity of the limestone sample was 4200-4500 meters per second, reflecting its dense structural characteristics.

[0042] S33. Process the physical property parameters of the measured samples to determine the differences and overlaps in physical property parameters between geological bodies.

[0043] In the physical performance parameter processing stage, the frequency distribution histogram of each physical property parameter is first drawn to determine the peak range and the overlapping area of ​​adjacent geological bodies. For example, the magnetic susceptibility overlap between a certain volcanic rock and a certain sandstone is between 30% and 40%.

[0044] The T-test was used to analyze the significance of the difference in means, with a significance level of α = 0.03-0.05. When the P-value was less than the set value, the parameter difference was determined to be significant. Then, the differences in physical parameters between geological bodies were analyzed and determined. For example, the P-values ​​of the densities of a certain limestone and a certain shale were 0.02-0.03, which were significant.

[0045] Establish a parameter distinguishability matrix for physical performance parameters, calculate the ratio of the overlapping area of ​​parameters between adjacent geological bodies to the total area (i.e., the smaller the value, the higher the distinguishability), and determine the degree of overlap of physical performance parameters between geological bodies based on the parameter distinguishability matrix. For example, the resistivity distinguishability matrix value of a certain metamorphic rock and a certain limestone is 0.2-0.25, which shows a high distinguishability.

[0046] S4. Set the threshold for the difference in physical property parameters, and process the geological bodies in the preliminary geological body division table according to the difference and overlap of physical property parameters to obtain the target geological bodies for three-dimensional geological modeling of the mining area. In this embodiment, a threshold for the difference in physical property parameters is set. The threshold for the difference in parameters includes a density difference of 0.08-0.1 g / cm³ and a magnetic susceptibility difference of 15 × 10⁻⁶ g / cm³. -6 -20×10 -6 SI, resistivity difference 120-150 ohm-meters, wave velocity difference 400-500 meters per second.

[0047] When the difference in physical property parameters between two adjacent geological bodies in the preliminary geological body division table is lower than the threshold for difference in physical property parameters and the overlap is not less than 40%-50%, the merging procedure is initiated to merge the two geological bodies with the difference in physical property parameters being lower than the threshold for difference in physical property parameters and the overlap being not less than 40%-50% into the same target geological body.

[0048] If the difference in physical property parameters between any two adjacent geological bodies in the preliminary geological body division table is higher than the threshold for difference in physical property parameters and the overlap is not higher than 15%-20%, the original geological body division is maintained, and the two geological bodies with any difference in physical property parameters higher than the threshold for difference in physical property parameters and the overlap is not higher than 15%-20% are treated as target geological bodies separately.

[0049] In some embodiments, the system further includes: a target geological body verification module, used to verify the target geological body, calculate the standard deviation of the physical performance parameters of the merged target geological body and the mean of the standard deviations of the physical performance parameters of the two preliminary geological bodies corresponding to the target geological body; when the standard deviation of the target geological body is higher than 80% of the mean of the standard deviations of the physical performance parameters of the two preliminary geological bodies, the merging of the target geological body is cancelled.

[0050] Specifically, merging adjacent preliminary geological bodies is based on their spatial continuity and geometric similarity (such as boundary fit and topological relationships). The goal is to merge spatially adjacent preliminary geological bodies with similar physical properties (belonging to the same geological unit) into a larger target geological body to simplify the model and reflect continuous geological units. However, spatial proximity is a necessary condition for merging, but not a sufficient one. Two spatially adjacent geological bodies may have vastly different physical properties (such as gold grade). In 3D geological modeling of mining areas, within the same geological body (such as a specific ore body, strata, or alteration zone), its key physical properties (especially the grade that determines economic value) should be similar. A certain degree of continuity and relative homogeneity is required (natural variations exist, but should be within a reasonable range). Therefore, after merging two adjacent geological bodies into a single target geological body, the rationality of the merger is verified by calculating the standard deviation of the physical property parameters of the merged target geological body. The standard deviation measures the dispersion or variability of the internal physical property parameters (such as density, magnetic susceptibility, resistivity, grade, etc.) of the new target geological body. When the standard deviation of the target geological body is less than 80% of the mean standard deviation of the physical property parameters of the two preliminary geological bodies, the internal properties of the merged target geological body are sufficiently homogeneous, meeting the expectations of a single geological unit, and the merger operation is confirmed. However, when the standard deviation of the target geological body is greater than 80% of the mean standard deviation of the physical property parameters of the two preliminary geological bodies, the internal properties of the merged target geological body are too different, indicating that preliminary geological bodies with different properties may have been mistakenly merged. In this case, the merger operation is canceled / rolled back, and the two preliminary geological bodies remain independent.

[0051] In some embodiments, the target geological body verification module is further configured to, after dividing two adjacent geological bodies into independent target geological bodies, acquire the boundary data of the target geological bodies and the physical field characteristic data of the mining area geology, process the physical characteristic data to obtain geophysical anomaly area data, calculate the consistency between the boundary data of the target geological bodies and the geophysical anomaly area data, and determine the rationality of the boundary data of the target geological bodies based on the consistency.

[0052] Understandably, the boundary data of a target geological body refers to the distribution boundary of the target geological body in three-dimensional space, that is, the actual contact interface between geological bodies, such as stratigraphic interfaces, fault lines, and lithological boundaries.

[0053] Geophysical anomaly data: This refers to areas whose physical field characteristics (such as density and magnetism) differ significantly from the surrounding normal physical field, determined through geophysical exploration methods such as gravity and magnetic methods, by measuring and inverting the physical field characteristics of the geological bodies in the mining area. The process of determining geophysical anomaly data involves three steps: Data acquisition: Using equipment such as gravimeters and magnetometers, collect physical field data such as gravitational acceleration and magnetic field strength in the mining area; Data processing and inversion: The measured data is processed using professional software. First, surface interference and normal field background values ​​are removed. Then, through inversion calculation, the specific distribution range of physical field anomalies is obtained. Geological interpretation: By combining the existing geological data of the mining area, we will analyze the correlation between the physical field anomaly area and the possible geological bodies (such as high-density ore bodies and strongly magnetic magmatic rocks) to clarify the geological significance of the "geophysical anomaly area".

[0054] In this embodiment, after dividing two adjacent geological bodies into independent target geological bodies, the boundary data of the target geological body and the geophysical anomaly area data are calculated based on the differences in physical properties such as density and magnetic susceptibility of the geological bodies (a boundary is considered reasonable only if the consistency reaches more than 70%), thereby verifying the rationality of the boundary in terms of physical properties from the data.

[0055] In some optional embodiments, after dividing two adjacent geological bodies into independent target geological bodies, the rationality of the boundary data of the target geological bodies can also be determined by the geological contact relationship of the target geological bodies.

[0056] Specifically, geological contacts refer to the spatial contact boundaries and interrelationships of different geological bodies (such as strata, rock masses, and structural surfaces). In this embodiment, after the target geological body is divided according to geological characteristics such as stratigraphic age, lithological assemblage, and structural occurrence, it is determined whether the spatial distribution of the target geological body is consistent with the measured geological contact relationship of the target geological body (such as stratigraphic dip angle and lithological stratification interface), so as to ensure that the boundary conforms to the geological structure law in terms of nature.

[0057] For example, in the process of determining whether the spatial distribution of the target geological body is consistent with the measured geological contact relationship of the target geological body, when the spatial distribution of the target geological body is inconsistent with the measured geological contact relationship of the target geological body, the boundary data of the target geological body can be adjusted based on the measured geological profile, the stratigraphic depth and lithological stratification data compiled from the borehole core, and the results of gravity and magnetic exploration inversion. For example, the boundary coordinates can be modified or the range can be redefined until the spatial distribution is reasonable.

[0058] For example, geophysical anomaly verification is carried out after the initial merging or splitting of the target geological body as a step to finally determine the boundary; while geological contact relationship check is checked immediately after splitting. If it is found that the split boundary and geological contact relationship do not match (for example, the boundary strike and the actual dip angle of the strata are significantly different), even if the geophysical consistency is sufficient, the boundary must be adjusted, and the geophysical consistency must be checked again until the two are consistent.

[0059] For example, Taking a volcanic rock (V1) and a contemporaneous subvolcanic clastic rock (V2), as well as a limestone (L1) and a shale (S1), as examples, the density difference between the volcanic rock (V1) and the contemporaneous subvolcanic clastic rock (V2) is 0.07-0.08 g / cm³, and the magnetic susceptibility difference is 12 × 10⁻⁶ g / cm³. -6 -15×10 -6 SI, resistivity difference 100-120 ohm-meters, wave velocity difference 350-400 meters per second, overlap 55%-65%, meet the merging conditions, and the merged body is named a certain volcanic rock assemblage (V1-V2). The standard deviation of its physical property parameters is verified. The standard deviation of density after merging is 0.05-0.06 grams per cubic centimeter, which is less than 70%-80% of the mean standard deviation of the two geological bodies before merging (0.06-0.07 grams per cubic centimeter), confirming the validity of the merging.

[0060] For the merged target geological bodies, the rationality of the boundary is verified by combining the results of gravity and magnetic exploration inversion. For example, if the magnetic anomaly boundary of a certain volcanic rock assemblage matches the boundary of the merged geological body with a degree of 70%-75% (≥70%), the boundary division is confirmed to be reasonable.

[0061] If the resistivity difference between a certain limestone (L1) and a certain shale (S1) reaches 700-800 ohm-meters (exceeding the 150 ohm-meter threshold), the wave velocity difference is 1000-1200 meters per second, and the overlap is 12%-15%, then the two rocks are separated if the difference between any two parameters exceeds the threshold and the overlap is not higher than 15%-20%. If the separation conditions are met, the original independent geological bodies are maintained.

[0062] After splitting, check whether the spatial distribution conforms to the geological contact relationship. Verify the interlayer contact interface between limestone and shale through three-dimensional modeling to ensure that the splitting boundary is consistent with the measured profile and borehole data.

[0063] The target geological bodies are integrated to generate a classification table of the target geological bodies, which clarifies the boundary coordinates, average physical property parameters and geological age information of each body, providing core data for 3D modeling.

[0064] In some embodiments, after determining the division of the target geological body, it is also necessary to perform three-dimensional geological modeling of the mining area based on the divided target geological body. In the process of three-dimensional geological modeling of the mining area, the collected geological and mineral data are first integrated, including measured profile data (25-30 lines), borehole core data (20-25 boreholes), and geophysical inversion data (8-10 gravity and magnetic survey lines each). A spatial database is constructed using a GIS platform to store parameters such as the three-dimensional coordinates of the stratigraphic interface (X, Y, and Z accuracy all reach 0.3-0.5m), the attitude of structural fractures (dip angle, dip direction, and spacing), and the extent of magmatic intrusion (boundary coordinates and thickness). A multi-source data association index is established to realize efficient management and retrieval of spatial information of geological bodies.

[0065] During the 3D raster model construction phase, the raster size is set to 8-10m×8-10m×8-10m, covering the entire mining area (1500-2000m east-west, 1200-1500m north-south, and 800-1000m deep), generating a total of 2-3 million raster units.

[0066] For each preliminary geological body, raster attributes are calculated using borehole data and profile interpolation, and information such as lithology code, formation age, and mean physical property parameters are assigned. For example, a raster unit in a certain limestone distribution area is assigned the lithology code L1, age 450-500 Ma, and density 2.6-2.7 g / cm³.

[0067] During the model construction process, structural attitude data is introduced to constrain the grid boundary, ensuring that the dip angle (average 30°-35°) and dip (NE direction) of the strata are consistent with the actual measurements, thereby improving the geometric accuracy of the model.

[0068] The final generated 3D grid model was verified by borehole data, and the burial depth error of the strata was less than 3-5m, which meets the accuracy requirements for 3D modeling of the mining area.

[0069] Outlier detection was performed on the collected physical property measurement data. The Z-score method was used to calculate the standard score of each sample parameter. A threshold of |Z|>2.5-3 was set as an outlier. For example, if the magnetic susceptibility Z value of a volcanic rock sample reached 3-3.5, it was judged as an outlier and removed. This ensured that the number of valid samples for each geological body was no less than 12-15 (the final number of valid samples was 18-25).

[0070] After outliers are removed, the mean and standard deviation are recalculated. For example, the mean density of a certain limestone is adjusted to 2.67-2.69 grams per cubic centimeter, and the standard deviation is 0.03-0.04 grams per cubic centimeter, making the data distribution more concentrated.

[0071] In the three-dimensional contour mapping stage, based on an 8-10m×8-10m×8-10m grid model, spatial interpolation was performed on parameters such as density and magnetic susceptibility. Kriging interpolation was used to generate a continuous data volume, and the contour spacing was set to 8%-10% of the parameter mean (e.g., for a density mean of 2.65 g / cm³, the spacing is 0.21-0.265 g / cm³). Three-dimensional contour maps of each physical property parameter were then drawn. For the magnetic susceptibility parameter, a high-value anomaly zone (>450×10m) was delineated in the eastern part of the mining area. -6 -500×10 -6 The SI value matches the measured magmatic intrusion range by 75%-80%; the resistivity contour map shows a deep low resistivity zone (<400-500 ohm-meters), which is consistent with the location of the mineralized alteration zone revealed by the borehole.

[0072] By using 3D visualization software (such as GOCAD) to overlay contour maps and geological body models, the spatial distribution characteristics of physical parameters can be intuitively displayed, providing visualization support for the verification of target geological body boundaries and the prediction of ore-forming target areas. Actual drilling verification shows that the anomaly area of ​​the contour lines matches the ore body by 70%-75%, effectively improving the practicality of multi-source data fusion analysis.

[0073] Secondly, embodiments of this application provide a method for determining target geological bodies in three-dimensional geological modeling of a mining area, the specific method including: Obtain geological and mineral data and measured geological profile data of the mining area; Based on the geological and mineral data of the mining area and the measured geological profile data, the geological bodies are divided, and a preliminary geological body division table is obtained; Collect physical property samples of the initial geological body, process the physical performance parameters of the measured physical property samples, calculate the frequency distribution histogram of the physical performance parameters, and determine the peak interval and overlapping area; The T-test method was used to analyze and determine the differences in physical properties among geological bodies; Establish a parameter distinguishability matrix for the physical performance parameters, and determine the overlap of the physical performance parameters between geological bodies based on the parameter distinguishability matrix; Set the threshold for the difference in physical property parameters, and process the geological bodies in the preliminary geological body division table according to the difference and overlap of physical property parameters to obtain the target geological bodies for three-dimensional geological modeling of the mining area; In some embodiments, the physical performance parameters include density, magnetic susceptibility, resistivity, and longitudinal wave velocity; the threshold values ​​for differences in physical property parameters include a density difference of 0.08-0.1 g / cm³ and a magnetic susceptibility difference of 15 × 10⁻⁶. -6 -20×10 -6 SI, resistivity difference 120-150 ohm-meters, wave velocity difference 400-500 meters per second, adjacent geological bodies in the preliminary geological body division table whose physical property parameters are all below the physical property parameter difference threshold and whose overlap is not less than 40%-50% are merged into one target geological body. If the difference in physical property parameters between any two adjacent geological bodies in the preliminary geological body division table is higher than the physical property parameter difference threshold and whose overlap is not higher than 20%, they are divided into independent target geological bodies.

[0074] For example, the threshold values ​​for the differences in physical properties such as density, magnetic susceptibility, resistivity, and longitudinal wave velocity are set to a density difference of 0.08-0.1 g / cm³, a magnetic susceptibility difference of 15 × 10⁻⁶ g / cm³, and a magnetic susceptibility difference of 15 × 10⁻⁶ g / cm³, respectively. -6 -20×10 -6 SI, resistivity difference of 120-150 ohm-meters, wave velocity difference of 400-500 meters per second. By setting thresholds for differences in physical property parameters, the quantitative standards are clarified and the accuracy of the classification is improved.

[0075] In some embodiments, the method further includes: after merging two adjacent geological bodies into a target geological body, calculating the standard deviation of the physical performance parameters of the merged target geological body and the mean of the standard deviations of the physical performance parameters of the two preliminary geological bodies corresponding to the target geological body; when the standard deviation of the target geological body is higher than 80% of the mean of the standard deviations of the physical performance parameters of the two preliminary geological bodies, canceling the merging of the target geological bodies.

[0076] In some embodiments, after dividing two adjacent geological bodies into independent target geological bodies, the boundary data of the target geological bodies and the physical field characteristic data of the mining area are obtained. The physical characteristic data is processed to obtain geophysical anomaly area data. The degree of agreement between the boundary data of the target geological bodies and the geophysical anomaly area data is calculated, and the rationality of the boundary data of the target geological bodies is determined based on the degree of agreement.

[0077] In some embodiments, the methods for measuring the physical performance parameters include: measuring density using the ring cutter method, measuring magnetic susceptibility using a proton precession magnetometer, measuring resistivity using the quadrupole method, and measuring longitudinal wave velocity using an ultrasonic detector.

[0078] In summary, the method for determining target geological bodies in the 3D geological modeling of this mining area involves collecting samples and measuring parameters such as density through a physical property parameter analysis module, calculating the mean and standard deviation, and using T-tests and discriminability matrices to analyze differences. The target geological body determination module merges or splits preliminary geological bodies based on difference thresholds and discriminability, achieving quantitative analysis and solving the problem of difficult-to-distinguish boundaries due to overlapping physical property parameters.

[0079] Furthermore, the method and system for determining target geological bodies in the three-dimensional geological modeling of the mining area integrates multi-source geological data to delineate preliminary geological bodies, establishes a multi-source data collaborative constraint mechanism, improves spatial coupling, reduces boundary misjudgment caused by single data deviation, and reduces drilling costs and risks. It solves the problems of existing technologies where the delineation of target geological bodies depends on surface features, lacks analysis of physical property parameters, and has low multi-source data coupling.

[0080] It should be noted that this application is not limited to the above-described embodiments. The above embodiments are merely examples, and any embodiments with the same structure and effect as the technical concept within the scope of this application are included in the technical scope of this application. Furthermore, various modifications that can be conceived by those skilled in the art to the embodiments, and other ways of constructing by combining some of the constituent elements of the embodiments, without departing from the spirit of this application, are also included in the scope of this application.

Claims

1. A method for determining target geological bodies in three-dimensional geological modeling of a mining area, characterized in that, The specific steps include: Obtain geological and mineral data and measured geological profile data of the mining area; Based on the geological and mineral data of the mining area and the measured geological profile data, the geological bodies are divided, and a preliminary geological body division table is obtained; Collect physical property samples of the initial geological body, process the physical performance parameters of the measured physical property samples, calculate the frequency distribution histogram of the physical performance parameters, and determine the peak interval and overlapping area; The T-test method was used to analyze and determine the differences in physical properties among geological bodies; Establish a parameter distinguishability matrix for the physical performance parameters, and determine the overlap of the physical performance parameters between geological bodies based on the parameter distinguishability matrix; Set a threshold for the difference in physical property parameters, and process the geological bodies in the preliminary geological body division table according to the difference and overlap of physical property parameters to obtain the target geological bodies for three-dimensional geological modeling of the mining area.

2. The method for determining target geological bodies in three-dimensional geological modeling of a mining area according to claim 1, characterized in that, The physical performance parameters include density, magnetic susceptibility, resistivity, and longitudinal wave velocity; the threshold values ​​for the differences in these physical property parameters include a density difference of 0.08-0.1 g / cm³ and a magnetic susceptibility difference of 15 × 10⁻⁶. -6 -20×10 -6 SI, resistivity difference 120-150 ohm-meters, wave velocity difference 400-500 meters per second, adjacent geological bodies in the preliminary geological body division table whose physical property parameters are all lower than the physical property parameter difference threshold and whose overlap is not less than 40%-50% are merged into one target geological body, and any two adjacent geological bodies in the preliminary geological body division table whose physical property parameters are higher than the physical property parameter difference threshold and whose overlap is not higher than 20% are divided into independent target geological bodies.

3. The method for determining target geological bodies in three-dimensional geological modeling of a mining area according to claim 2, characterized in that, Also includes: After merging two adjacent geological bodies into one target geological body, the standard deviation of the physical performance parameters of the merged target geological body and the mean standard deviation of the physical performance parameters of the two preliminary geological bodies corresponding to the target geological body are calculated; when the standard deviation of the target geological body is higher than 80% of the mean standard deviation of the physical performance parameters of the two preliminary geological bodies, the merging of the target geological body is cancelled.

4. The method for determining target geological bodies in three-dimensional geological modeling of a mining area according to claim 2, characterized in that, After dividing two adjacent geological bodies into independent target geological bodies, the boundary data of the target geological bodies and the physical field characteristic data of the mining area are obtained. The physical characteristic data is processed to obtain geophysical anomaly area data. The degree of agreement between the boundary data of the target geological bodies and the geophysical anomaly area data is calculated, and the rationality of the boundary data of the target geological bodies is determined based on the degree of agreement.

5. The method for determining target geological bodies in three-dimensional geological modeling of a mining area according to claim 1, characterized in that, The methods for measuring the physical performance parameters include: measuring density using the ring cutter method, measuring magnetic susceptibility using a proton precession magnetometer, measuring resistivity using the quadrupole method, and measuring longitudinal wave velocity using an ultrasonic detector.

6. A system for determining target geological bodies in three-dimensional geological modeling of mining areas, characterized in that, include: The data acquisition module is used to acquire geological and mineral data and measured geological profile data of the mining area; The preliminary geological body division module is used to divide geological bodies based on the geological and mineral data of the mining area and the measured geological profile data, and to obtain a preliminary geological body division table; The physical property parameter analysis module is used to collect physical property measurement samples of the initial geological body, process the physical property parameters of the measured physical property samples, calculate the frequency distribution histogram of the physical property parameters, determine the peak interval and overlapping area, use the T-test method to analyze and determine the differences in physical property parameters between geological bodies, establish the parameter distinguishability matrix of the physical property parameters, and determine the overlap of the physical property parameters between geological bodies based on the parameter distinguishability matrix. The target geological body determination module is used to set the threshold for differences in physical property parameters and process the geological bodies in the preliminary geological body division table according to the differences and overlaps in physical property parameters to obtain the target geological bodies for three-dimensional geological modeling of the mining area.

7. The system for determining target geological bodies in three-dimensional geological modeling of mining areas according to claim 6, characterized in that, In the target geological body identification module, the physical performance parameters include density, magnetic susceptibility, resistivity, and P-wave velocity; the threshold values ​​for the differences in these physical property parameters include a density difference of 0.08-0.1 g / cm³ and a magnetic susceptibility difference of 15 × 10⁻⁶. -6 -20×10 -6 SI, resistivity difference 120-150 ohm-meters, wave velocity difference 400-500 meters per second. When the difference in physical property parameters between two adjacent geological bodies in the preliminary geological body division table is lower than the difference threshold of physical property parameters and the overlap is not less than 40%-50%, the adjacent preliminary geological bodies are merged into the same target geological body. When the difference in physical property parameters between any two adjacent geological bodies in the preliminary geological body division table is higher than the difference threshold of physical property parameters and the overlap is not higher than 20%, they are divided into independent target geological bodies.

8. The system for determining target geological bodies in three-dimensional geological modeling of mining areas according to claim 7, characterized in that, Also includes: The target geological body verification module is used to calculate the standard deviation of the physical performance parameters of the merged target geological body and the mean standard deviation of the physical performance parameters of the two preliminary geological bodies corresponding to the target geological body after merging two adjacent geological bodies into one target geological body; when the standard deviation of the target geological body is higher than 80% of the mean standard deviation of the physical performance parameters of the two preliminary geological bodies, the merging of the target geological body is cancelled.

9. The system for determining target geological bodies in three-dimensional geological modeling of mining areas according to claim 7, characterized in that, The target geological body verification module also includes: After dividing two adjacent geological bodies into independent target geological bodies, the boundary data of the target geological bodies and the physical field characteristic data of the mining area are obtained. The physical characteristic data is processed to obtain geophysical anomaly area data. The degree of agreement between the boundary data of the target geological bodies and the geophysical anomaly area data is calculated, and the rationality of the boundary data of the target geological bodies is determined based on the degree of agreement.

10. The system for determining target geological bodies in three-dimensional geological modeling of mining areas according to claim 6, characterized in that, The methods for measuring the physical performance parameters include: measuring density using the ring cutter method, measuring magnetic susceptibility using a proton precession magnetometer, measuring resistivity using the quadrupole method, and measuring longitudinal wave velocity using an ultrasonic detector.

Images