Three-dimensional geological modeling method and system for non-ferrous metal mine based on exploration line profile

Abstract

This application provides a method and system for three-dimensional geological modeling of non-ferrous metal deposits based on exploration line profiles, relating to the field of geological exploration technology. This application identifies curvature abrupt change points at the corresponding stratigraphic interfaces by acquiring coordinate data of stratigraphic interfaces on multiple exploration line profiles. Priority identifiers are assigned to the normal vector, curvature value, and position coordinates of curvature abrupt change points for each stratigraphic interface to form a comprehensive dataset with priority identifiers. A three-dimensional model framework is constructed based on the extreme values ​​of the three-dimensional coordinates of each exploration line profile. Surface feature data corresponding to each stratigraphic grid unit is obtained using three-dimensional texture mapping technology. Structural identifier data is generated based on the model position coordinates of each curvature abrupt change point within the three-dimensional model framework. Based on the priority order of the comprehensive dataset, a three-dimensional geological model of the non-ferrous metal deposit is generated, achieving high-precision, structurally constrained, and geologically significant three-dimensional geological modeling.

Description

Technical Field

[0001] This application relates to the field of geological exploration technology, and in particular to a method and system for three-dimensional geological modeling of non-ferrous metal mines based on exploration line profiles. Background Technology

[0002] In the process of non-ferrous metal exploration, in order to support high-precision resource assessment and mining design, it is urgent to accurately reproduce the complex stratigraphic interface morphology contained in the exploration line profile into a three-dimensional geological model.

[0003] The current mainstream approach is to use a combination of multi-profile interpolation and conventional texture mapping for 3D modeling. First, the profile stratigraphic contour is constructed using borehole and geophysical data. Then, a 3D mesh surface is generated using radial basis functions or kriging interpolation algorithms. Finally, the geometric properties in the profile are mapped to the entire model surface using a uniform texture.

[0004] However, existing solutions have significant shortcomings in handling abrupt changes in local morphology of strata: their interpolation process is too smooth, making it difficult to retain key curvature anomaly information in the original profile; at the same time, texture mapping does not distinguish the physical characteristics of different tectonic regions, resulting in a lack of effective visual and semantic distinction between fault fracture zones and intact rock strata, and failing to meet the modeling requirements for high-fidelity restoration of micro-tectonic features. Summary of the Invention

[0005] The purpose of this application is to provide a method and system for three-dimensional geological modeling of non-ferrous metal mines based on exploration line profiles, so as to solve the problem that it is difficult to retain key curvature anomaly information in the original profile in the existing technology.

[0006] Firstly, this application provides a three-dimensional geological modeling method for non-ferrous metal deposits based on exploration line profiles, including:

[0007] Obtain coordinate data of stratigraphic interfaces of non-ferrous metal deposits on multiple exploration profiles;

[0008] Geometric calculations are performed on the coordinate data to obtain the normal vector and curvature value of each stratigraphic interface. Numerical variation analysis is then performed on all curvature values ​​of each stratigraphic interface to identify the curvature abrupt change points of the corresponding stratigraphic interface.

[0009] Priority labels are assigned to the normal vector, curvature value, and location coordinates of curvature abrupt change points of each stratigraphic interface to form a comprehensive dataset with priority labels;

[0010] A three-dimensional model framework is constructed based on the extreme values ​​of the three-dimensional coordinates of each exploration line profile;

[0011] The curvature value of each stratigraphic interface is bound to the corresponding stratigraphic mesh unit in the three-dimensional model framework using three-dimensional texture mapping technology to obtain the surface feature data corresponding to each stratigraphic mesh unit.

[0012] Based on the model position coordinates of each curvature change point in the 3D model framework, and combined with the preset construction type labels, construction identification data is generated.

[0013] Based on the priority order of the comprehensive dataset, the three-dimensional model framework, surface feature data under the normal vector constraints of each stratigraphic interface, structural identification data, and coordinate data are jointly constructed to generate a three-dimensional geological model of non-ferrous metal deposits.

[0014] Optionally, a three-dimensional model framework is constructed based on the extreme values ​​of the three-dimensional coordinates of each exploration line profile, including:

[0015] Based on the extreme values ​​of the three-dimensional coordinates of each exploration line profile, determine the spatial coordinate boundaries of all exploration line profiles;

[0016] Based on the preset grid size and coordinate data, the three-dimensional space within the spatial coordinate boundary is divided into multiple stratigraphic grid units corresponding to each stratigraphic interface;

[0017] Based on the stratigraphic interface, the distribution range of strata within the stratigraphic profile, and the exploration line number of each exploration line profile, the coordinate data of adjacent exploration line profiles are interpolated and extended to obtain the attribute information of each stratigraphic grid cell. The attribute information includes the regional cell number, the regional coordinate range, and the exploration line number.

[0018] Based on the attribute information of all stratigraphic grid cells, spatial coordinate matching and attribute consistency verification are performed on all stratigraphic grid cells to obtain stratigraphic grid cells. The target stratigraphic grid cells are then stitched together to obtain a three-dimensional model framework.

[0019] Optionally, a 3D texture mapping technique is used to bind the curvature value of each stratigraphic interface to the corresponding stratigraphic mesh cell in the 3D model framework to obtain the surface feature data corresponding to each stratigraphic mesh cell, including:

[0020] Based on the curvature value of each stratigraphic interface, the range of the first vertical coordinate of each stratigraphic interface is divided to obtain multiple sub-intervals;

[0021] The overlap between the sub-interval of each stratigraphic interface and the second vertical coordinate range of the corresponding stratigraphic grid cell in the 3D model frame is calculated. Sub-intervals with an overlap greater than or equal to a preset overlap threshold are associated with the corresponding stratigraphic grid cells to obtain multiple association relationships.

[0022] According to the preset texture mapping rules, the texture level corresponding to each curvature value is determined by the three-dimensional texture mapping technology, and each curvature value is converted into a texture parameter according to the texture level corresponding to each curvature value.

[0023] Based on the association relationship, the texture parameters of each stratigraphic interface are associated and bound with the corresponding stratigraphic grid cell to obtain multiple feature association results;

[0024] Add a corresponding stratigraphic number to each feature association result, and combine it with the regional cell number of each stratigraphic grid cell to generate surface feature data corresponding to each stratigraphic grid cell.

[0025] Optionally, based on the model position coordinates of each curvature abrupt change point within the 3D model framework, and combined with preset construction type labels, construction identification data is generated, including:

[0026] Based on the first vertical coordinate range of each stratigraphic interface, the first vertical coordinate range is divided into multiple vertical coordinate intervals, and the vertical coordinate interval to which each curvature change point belongs is determined.

[0027] Calculate the distance between the position coordinates of the curvature change point in each vertical coordinate interval and the center coordinates of the corresponding stratigraphic grid cell in the 3D model frame. Take the center coordinates corresponding to the smallest distance as the model position coordinates of the corresponding curvature change point in the 3D model frame.

[0028] Calculate the difference between the vertical coordinate of each curvature abrupt change point and the center vertical coordinate of the corresponding stratigraphic grid cell, and use the model position coordinates of curvature abrupt change points with a difference less than or equal to a preset threshold as the final associated coordinates;

[0029] The final associated coordinates of each curvature abrupt change point, the corresponding stratigraphic grid cell regional cell number, the vertical coordinate interval to which it belongs, and the preset structural type label are associated to obtain structural identification data.

[0030] Optionally, based on the priority order of the comprehensive dataset, the 3D model framework, surface feature data under the normal vector constraints of each stratigraphic interface, structural identification data, and coordinate data are jointly constructed to generate a 3D geological model of the non-ferrous metal deposit, including:

[0031] The priority order of the comprehensive dataset is determined based on the influence of the normal vector, curvature value, and location coordinates of curvature change points of each stratigraphic interface on the 3D model framework.

[0032] Based on the priority order and the normal vector of each stratigraphic interface, the spatial angle of the stratigraphic mesh unit in the 3D model frame is adjusted to obtain the adjusted 3D model frame. The adjusted 3D model frame and the normal vector of each stratigraphic interface are jointly constructed to obtain the first joint model.

[0033] The texture parameters and structural identification data in the surface feature data are combined with the stratigraphic grid cells with the corresponding regional unit numbers in the first joint model to construct the second joint model.

[0034] Based on the coordinate data, the second joint model was adjusted and verified to obtain a three-dimensional geological model of the non-ferrous metal mine.

[0035] Optionally, the priority order of the integrated dataset is determined based on the influence of the normal vector, curvature value, and location coordinates of curvature abrupt change points of each stratigraphic interface on the 3D model framework, including:

[0036] Based on the degree of core role of the normal vector, curvature value, and curvature change point of each stratigraphic interface in the three-dimensional geological model, the first weight corresponding to each normal vector, the second weight corresponding to the curvature value, and the third weight corresponding to the curvature change point are determined.

[0037] Calculate the first deviation of the overall attitude of the 3D model frame when the normal vector of each stratigraphic interface changes, and combine it with the first weight to calculate the first comprehensive weight of the normal vector of each stratigraphic interface on the 3D model frame.

[0038] The second deviation of the surface detail features of the 3D model frame is calculated when the curvature value of each stratum interface changes. Combined with the second weight, the second comprehensive weight of the curvature value of each stratum interface on the 3D model frame is calculated.

[0039] The third deviation of the local structural identifier of the 3D model frame is calculated when the position coordinates of the curvature change points of each stratigraphic interface are offset, and combined with the third weight, the third comprehensive weight of the position coordinates of each curvature change point to the 3D model frame is calculated.

[0040] The priority order of the comprehensive dataset is determined based on the first comprehensive weight, the second comprehensive weight, and the third comprehensive weight.

[0041] Optionally, based on the coordinate data, the second joint model is adjusted and verified to obtain a three-dimensional geological model of the non-ferrous metal deposit, including:

[0042] Calculate the deviation between the model interface coordinates and the corresponding actual interface coordinates in the coordinate data for each stratigraphic interface in the second combined model;

[0043] The model interface coordinates of the strata interfaces in the second joint model with deviation values ​​greater than the preset deviation threshold are adjusted until all deviation values ​​are less than or equal to the preset deviation threshold, thus obtaining the adjusted second joint model.

[0044] The boundary coincidence parameters of any two adjacent stratigraphic grid cells in the adjusted second joint model are calculated, and the boundary coincidence parameters are verified according to the preset coincidence conditions to obtain the coincidence verification results.

[0045] Calculate the texture transition parameters between any two adjacent stratigraphic grid cells in the adjusted second joint model, and verify the texture transition parameters according to the preset transition conditions to obtain the transition verification results;

[0046] Based on the preset corresponding conditions, the correspondence between the structural type labels of the stratigraphic grid cells with structural identifiers and the final associated coordinates in the adjusted second joint model is determined, and the relationship analysis results are obtained.

[0047] If the overlap verification result, transition verification result, and relationship analysis result are all passed, the adjusted second joint model will be used as the three-dimensional geological model of the non-ferrous metal mine. If any result is failed, the second joint model will be readjusted until the overlap verification result, transition verification result, and relationship analysis result are all passed.

[0048] Secondly, this application provides a three-dimensional geological modeling system for non-ferrous metal deposits based on exploration line profiles, including:

[0049] The acquisition module is used to acquire the coordinate data of the stratigraphic interfaces of non-ferrous metal mines on multiple exploration profiles.

[0050] The calculation module is used to perform geometric calculations on the coordinate data to obtain the normal vector and curvature value of each stratigraphic interface. It performs numerical variation analysis on all curvature values ​​of each stratigraphic interface to identify the curvature abrupt change points of the corresponding stratigraphic interface.

[0051] The allocation module is used to assign priority identifiers to the normal vector, curvature value, and location coordinates of curvature abrupt change points of each stratigraphic interface to form a comprehensive dataset with priority identifiers.

[0052] The building module is used to construct a three-dimensional model framework based on the extreme values ​​of the three-dimensional coordinates of each exploration line profile;

[0053] The binding module is used to bind the curvature value of each stratum interface to the corresponding stratum grid cell in the 3D model framework using 3D texture mapping technology, so as to obtain the surface feature data corresponding to each stratum grid cell.

[0054] The building module is also used to generate construction identification data based on the model position coordinates of each curvature change point in the 3D model framework, combined with preset construction type labels;

[0055] The building module is also used to jointly construct a three-dimensional geological model of a non-ferrous metal mine based on the priority order of the comprehensive dataset, combining the three-dimensional model framework, surface feature data under the normal vector constraints of each stratigraphic interface, structural identification data, and coordinate data.

[0056] Thirdly, this application provides an electronic device, comprising:

[0057] Memory, used to store computer programs;

[0058] A processor is used to execute computer programs to implement the steps of a three-dimensional geological modeling method for non-ferrous metal mines based on exploration line profiles, as described in the first aspect above.

[0059] Fourthly, this application provides a computer-readable storage medium storing a computer program that, when executed by a processor, can implement the steps of the three-dimensional geological modeling method for non-ferrous metal mines based on exploration line profiles as described in the first aspect above.

[0060] This application provides a method for three-dimensional geological modeling of non-ferrous metal deposits based on exploration line profiles. By acquiring the coordinate data of stratigraphic interfaces on multiple exploration line profiles, it ensures that the model construction has a realistic geological basis; it effectively captures local morphological changes at stratigraphic interfaces, accurately locates key areas that may correspond to geological structures such as faults and joints, and avoids the loss of microstructural information due to smooth interpolation; it guarantees the overall spatial consistency and boundary integrity of the model; it enhances the interpretability and engineering guidance value of the model; and it achieves a structurally reasonable, semantically clear, and geometrically accurate overall model while preserving the details of the original profiles.

[0061] Furthermore, the priority order of the comprehensive dataset is first determined based on the influence of normal vectors, curvature values, and curvature abrupt change points on the model framework. Then, based on this priority and the normal vectors of each stratigraphic interface, the spatial angles of the stratigraphic grid units in the 3D model framework are adjusted to form the first joint model. Subsequently, the texture parameters and structural identification data in the surface feature data are fused with the stratigraphic grid units of the corresponding regional unit numbers in the first joint model to construct the second joint model. Finally, the second joint model is adjusted and verified based on the original coordinate data to output the final 3D geological model of the non-ferrous metal mine.

[0062] It solves the problem of losing curvature anomaly information due to excessive smoothing interpolation in existing schemes, and overcomes the defect that uniform texture mapping cannot distinguish the physical characteristics of different structural regions; it improves the ability to restore the micro-structural features in the exploration line profile and the accuracy of geological semantic expression. Attached Figure Description

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

[0064] Figure 1 A flowchart illustrating a three-dimensional geological modeling method for non-ferrous metal mines based on exploration line profiles, provided for an embodiment of this application;

[0065] Figure 2 A flowchart illustrating a method for generating a three-dimensional geological model, provided in an embodiment of this application;

[0066] Figure 3 This is a schematic diagram of the structure of a three-dimensional geological modeling system for non-ferrous metal mines based on exploration line profiles, provided as an embodiment of this application. Detailed Implementation

[0067] To address the problem that existing technologies for 3D geological modeling of non-ferrous metal mines often result in blurred local morphological details of strata and a lack of effective distinction between fault fracture zones and intact rock strata due to the use of global smoothing interpolation and unified texture mapping, this application proposes a 3D geological modeling method and system for non-ferrous metal mines based on exploration line profiles. This method achieves high-fidelity and high-semantic 3D geological reconstruction of minute structural features while avoiding excessive smoothing.

[0068] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0069] The core of this application is to provide a three-dimensional geological modeling method for non-ferrous metal mines based on exploration line profiles. A flowchart of one specific implementation method is shown below. Figure 1 As shown, the method includes:

[0070] Step 101: Obtain the coordinate data of the stratigraphic interfaces of the non-ferrous metal ore on multiple exploration profiles.

[0071] In this step, the exploration line profile refers to a two-dimensional cross-section corresponding to a virtual cutting line artificially set during non-ferrous metal ore exploration, used to observe the distribution of underground strata. The stratigraphic interface refers to the contact surface between strata of different ages or compositions underground.

[0072] In this embodiment of the application, the data acquisition object is first clearly defined as non-ferrous metal mine, and the data source is the stratigraphic interface on multiple exploration line profiles to obtain the coordinate data of the stratigraphic interface.

[0073] For example, in the preliminary exploration stage of three-dimensional geological modeling of a non-ferrous metal mine, the staff arranged 12 parallel exploration line profiles in the mining area, numbered L1 to L12, based on the distribution characteristics of the ore body. Each exploration line profile extends along the dip of the ore body and covers the core ore-bearing area of ​​the mining area.

[0074] Data was collected from key stratigraphic interfaces on each exploration line profile, with a focus on measuring the three-dimensional coordinates of ore-bearing stratigraphic interfaces in limestone, shale interfaces in ore-bearing strata, and sandstone interfaces in shale. During the data collection process, the unified coordinate system of the mining area was used as the benchmark to ensure that the spatial benchmark of all coordinate data was consistent. Through data processing, 156 sets of complete coordinate data for the three types of stratigraphic interfaces on the 12 exploration line profiles were obtained. Each set of data includes specific values ​​for the three dimensions of the horizontal, vertical, and vertical axes.

[0075] Step 102: Perform geometric calculations on the coordinate data to obtain the normal vector and curvature value of each stratigraphic interface. Perform numerical variation analysis on all curvature values ​​of each stratigraphic interface to identify the curvature abrupt change points of the corresponding stratigraphic interface.

[0076] In this step, the curvature abrupt change point refers to the point where a significant abrupt change in curvature value is identified after numerical variation analysis of all curvature values ​​at each stratigraphic interface.

[0077] In this embodiment of the application, firstly, geometric calculations are performed on the acquired coordinate data to obtain the normal vector and curvature value of each stratigraphic interface. The curvature value is calculated to be equal to the reciprocal of the radius of curvature of the stratigraphic interface curve. Secondly, numerical change analysis is performed on all curvature values ​​of each stratigraphic interface. By judging the magnitude of the change in curvature value, the curvature change point of the corresponding stratigraphic interface is identified.

[0078] For example, the 156 sets of coordinate data obtained in the steps are imported into professional geological data processing software. First, vector operations and curve fitting are performed on each set of coordinate data to calculate the normal vector of each stratigraphic interface. At the same time, the curvature value corresponding to each coordinate point on each stratigraphic interface is calculated. The curvature value of the ore-bearing strata interface in limestone ranges from 0.002 to 0.006, and the curvature value of the shale interface in ore-bearing strata ranges from 0.001 to 0.003.

[0079] Subsequently, all curvature values ​​of each stratigraphic interface were sorted in vertical coordinate order, and the difference between adjacent curvature values ​​was calculated. A threshold of 0.0015 was set for this difference. When the difference between the curvature value of a certain coordinate point and its adjacent points exceeded this threshold, it was identified as a curvature abrupt change point. Ultimately, four curvature abrupt change points were identified in the ore-bearing strata interfaces of the limestone and two curvature abrupt change points in the shale interfaces of the ore-bearing strata.

[0080] Step 103: Assign priority labels to the normal vector, curvature value, and location coordinates of curvature abrupt change points of each stratigraphic interface to form a comprehensive dataset with priority labels.

[0081] In this step, priority identifiers refer to the identifiers assigned to the normal vector, curvature value, and position coordinates of curvature abrupt change points of each stratigraphic interface, used to distinguish the importance of the three in the modeling process.

[0082] In this embodiment of the application, it is first determined that the objects that need to be assigned priority labels include the normal vector, curvature value and location coordinates of curvature change points of each stratigraphic interface. By reasonably assigning priority labels, these three types of data are integrated into a comprehensive dataset with priority labels. This dataset will provide a basis for determining the priority order of subsequent modeling.

[0083] For example, combining geological modeling logic, the impact of the normal vectors, curvature values, and coordinates of curvature abrupt change points at each stratigraphic interface on subsequent modeling is analyzed. Normal vectors directly determine the spatial attitude of the strata, influencing the overall model structure and playing the strongest core role; curvature values ​​determine the detailed features of the stratigraphic surface, affecting the model's surface fidelity and playing a secondary core role; the coordinates of curvature abrupt change points only affect local structural identifiers and have the weakest core role. Based on this, a priority identifier of level 1 is assigned to normal vectors, level 2 to curvature values, and level 3 to the coordinates of curvature abrupt change points. Subsequently, the normal vector data, curvature value data, and coordinates of curvature abrupt change points for each stratigraphic interface are bound to their corresponding priority identifiers, integrating them into a comprehensive dataset with priority identifiers.

[0084] Step 104: Construct a three-dimensional model framework based on the extreme values ​​of the three-dimensional coordinates of each exploration line profile.

[0085] In this step, the three-dimensional coordinate extreme values ​​refer to the maximum and minimum values ​​corresponding to the horizontal, vertical, and triangular coordinates of each exploration line profile.

[0086] Step 105: Use 3D texture mapping technology to bind the curvature value of each stratum interface to the corresponding stratum grid cell in the 3D model framework to obtain the surface feature data corresponding to each stratum grid cell.

[0087] In this step, 3D texture mapping technology refers to the technique of transforming abstract data into a visual 3D model surface texture. A stratigraphic mesh unit refers to a 3D spatial unit corresponding to a specific stratigraphic interface, divided according to a preset mesh size within the spatial coordinate boundaries of the 3D model framework.

[0088] Step 106: Based on the model position coordinates of each curvature change point in the 3D model framework, and combined with the preset construction type labels, generate construction identification data.

[0089] In this step, the preset structure type label refers to a pre-defined label used to identify the type of underground geological structure, such as faults and joints.

[0090] Step 107: Based on the priority order of the comprehensive dataset, the three-dimensional model framework, surface feature data under the normal vector constraints of each stratigraphic interface, structural identification data, and coordinate data are jointly constructed to generate a three-dimensional geological model of the non-ferrous metal mine.

[0091] In this step, the priority order refers to the order in which various types of data in the comprehensive dataset affect the 3D model framework, and is used to guide the joint construction.

[0092] This application's embodiments ensure that the model construction has a real geological basis by acquiring the coordinate data of stratigraphic interfaces on multiple exploration line profiles of non-ferrous metal mines; effectively capture local morphological changes at stratigraphic interfaces, accurately locate key areas that may correspond to geological structures such as faults and joints, and avoid the loss of microstructural information due to smooth interpolation; guarantee the overall spatial consistency and boundary integrity of the model; enhance the model's interpretability and engineering guidance value; and achieve a structurally reasonable, semantically clear, and geometrically accurate overall modeling while preserving the details of the original profiles. It solves the problems of messy data use, insufficient model detail reproduction, and ambiguous structural identification in existing technologies.

[0093] This application provides a specific embodiment. Step 104 involves constructing a three-dimensional model framework based on the extreme values ​​of the three-dimensional coordinates of each exploration line profile. This specifically includes the following steps:

[0094] Step 401: Determine the spatial coordinate boundaries of all exploration line profiles based on the extreme values ​​of the three-dimensional coordinates of each profile.

[0095] In this embodiment, the maximum and minimum values ​​of the abscissa, ordinate, and ordinate of each exploration line profile are first extracted. Then, the extreme values ​​of the abscissa, ordinate, and ordinate of all exploration line profiles are summarized. The minimum and maximum values ​​of the abscissa, ordinate, and ordinate of all profiles are taken to form the spatial coordinate boundary of all exploration line profiles, ensuring that the boundary can completely cover the spatial distribution range of all exploration line profiles.

[0096] Step 402: Based on the preset grid size and coordinate data, the three-dimensional space within the spatial coordinate boundary is divided into multiple stratigraphic grid units corresponding to each stratigraphic interface.

[0097] In this step, the preset grid size refers to the pre-defined size standard of the grid used to divide the three-dimensional space in the three dimensions of horizontal, vertical and vertical coordinates.

[0098] In this embodiment, the preset grid size is first determined, and the spatial distribution characteristics of each stratigraphic interface in the acquired coordinate data are combined with the determined spatial coordinate boundary as the delineation range. Then, the grid is evenly divided along the three dimensions of horizontal, vertical and vertical coordinates according to the preset grid size, so that each stratigraphic grid unit can accurately correspond to a specific stratigraphic interface, and finally multiple stratigraphic grid units corresponding to each stratigraphic interface are obtained.

[0099] Step 403: Based on the stratigraphic interface of each exploration line profile, the distribution range of the strata within the exploration line profile, and the exploration line number, interpolate and extend the coordinate data of adjacent exploration line profiles to obtain the attribute information of each stratigraphic grid cell. The attribute information includes the regional cell number, the regional coordinate range, and the exploration line number.

[0100] In this step, the exploration line number refers to the unique identifier assigned to each exploration line profile to distinguish between different exploration line profiles. The regional unit number refers to the unique identifier assigned to each stratigraphic grid unit to distinguish between different stratigraphic grid units. The regional coordinate range refers to the start and end intervals of the coordinates of each stratigraphic grid unit in the three dimensions of x, y, and y.

[0101] In this embodiment, the stratigraphic interface, the distribution range of the strata within the exploration line profile, and the corresponding exploration line number are first defined for each exploration line profile. Based on the coordinate data of two adjacent exploration line profiles, interpolation is used to supplement the missing coordinate data between adjacent profiles. First, the first difference between the target location and the smaller known coordinate location, the difference in the vertical coordinates of the two known coordinate values, and the difference in the horizontal coordinates of the two known coordinate locations are calculated. Then, the first difference is multiplied by the difference in the vertical coordinates and divided by the difference in the horizontal coordinates to obtain the result of the multiplication and division operation. Finally, the smaller vertical coordinate value among the two adjacent known coordinate values ​​is added to the result of the multiplication and division operation to obtain the interpolated coordinate value.

[0102] Based on the expanded complete coordinate data, a unique regional unit number is assigned to each stratigraphic grid unit, the regional coordinate range of each stratigraphic grid unit is defined, and each stratigraphic grid unit is associated with the corresponding exploration line number. Finally, the attribute information of each stratigraphic grid unit, including the regional unit number, regional coordinate range, and exploration line number, is obtained.

[0103] Step 404: Based on the attribute information of all stratigraphic grid cells, perform spatial coordinate matching and attribute consistency verification on all stratigraphic grid cells to obtain stratigraphic grid cells, and then stitch and integrate the target stratigraphic grid cells to obtain a three-dimensional model framework.

[0104] In this step, the target stratigraphic grid cell refers to the stratigraphic grid cell that meets the preset matching conditions and verification standards after spatial coordinate matching and attribute consistency verification.

[0105] In this embodiment, the attribute information of all stratigraphic grid units is first extracted. Based on the regional coordinate range in the attribute information, spatial coordinate matching is performed on all stratigraphic grid units to ensure that the regional coordinate ranges of adjacent stratigraphic grid units can be seamlessly connected. At the same time, the consistency of the regional unit number, regional coordinate range, and exploration line number in the attribute information of each stratigraphic grid unit is verified. Stratigraphic grid units that do not meet the preset matching conditions and verification standards are eliminated, and target stratigraphic grid units are selected. Finally, the target stratigraphic grid units are spliced ​​and integrated according to the regional coordinate range and attribute association relationship of each target stratigraphic grid unit to form a three-dimensional model framework with complete structure, connected coordinates, and consistent attributes.

[0106] The embodiments of this application achieve a precise correspondence between the three-dimensional model framework and the exploration line profile data, ensuring that the framework covers the spatial range of all exploration line profiles, and that the attribute information of each stratigraphic grid unit is clear and consistent. This solves the problem in the prior art where the three-dimensional model framework is not fully covered and the grid unit attributes are chaotic, resulting in an unstable modeling foundation.

[0107] This application provides a specific embodiment. Step 105 involves using three-dimensional texture mapping technology to bind the curvature value of each stratigraphic interface to the corresponding stratigraphic mesh unit in the three-dimensional model framework, in order to obtain the surface feature data corresponding to each stratigraphic mesh unit. This specifically includes the following steps:

[0108] Step 501: Based on the curvature value of each stratigraphic interface, divide the first vertical coordinate range of each stratigraphic interface to obtain multiple sub-intervals.

[0109] In this step, the first vertical coordinate range refers to the start and end interval of the coordinates of each stratigraphic interface in the vertical coordinate dimension.

[0110] In this embodiment of the application, the curvature value distribution of each stratigraphic interface is first clarified. Based on the numerical variation law of the curvature value, the range of the first vertical coordinate corresponding to each stratigraphic interface is divided according to different intervals of curvature value to ensure that the curvature values ​​in each sub-interval after division have similar characteristics, and finally multiple sub-intervals are obtained.

[0111] Step 502: Calculate the overlap between the sub-interval of each stratigraphic interface and the second vertical coordinate range of the corresponding stratigraphic grid cell in the 3D model frame. Associate the sub-intervals with an overlap greater than or equal to the preset overlap threshold with the corresponding stratigraphic grid cell to obtain multiple association relationships.

[0112] In this step, the second vertical coordinate range refers to the start and end interval of the coordinates of each stratigraphic grid cell in the vertical coordinate dimension. The preset overlap threshold is a pre-set critical value used to determine whether a sub-interval is associated with a stratigraphic grid cell.

[0113] In this embodiment, the starting and ending ranges of the vertical coordinates corresponding to the sub-intervals of each stratigraphic interface are first extracted, as well as the second vertical coordinate range of the corresponding stratigraphic grid cell in the three-dimensional model frame. The overlap between the two is calculated as follows: First, the length of the first vertical coordinate of the sub-interval is calculated, and then the length of the second vertical coordinate of the intersection of the sub-interval and the second vertical coordinate range is calculated; then, the length of the second vertical coordinate is divided by the length of the first vertical coordinate; finally, it is multiplied by 100% to obtain the overlap.

[0114] The calculated overlap is then compared with a preset overlap threshold. Sub-intervals with overlap greater than or equal to the preset overlap threshold are associated with the corresponding stratigraphic grid cells, resulting in multiple association relationships.

[0115] Step 503: According to the preset texture mapping rules, use three-dimensional texture mapping technology to determine the texture level corresponding to each curvature value, and convert each curvature value into a texture parameter according to the texture level corresponding to each curvature value.

[0116] In this step, the preset texture mapping rule refers to the pre-defined rule used to match curvature values ​​with texture levels and texture parameters.

[0117] In this embodiment of the application, the correspondence between the curvature value range and the texture level and texture parameters in the preset texture mapping rule is first clarified. Based on the rule, the curvature value range to which each curvature value belongs is determined one by one, the texture level corresponding to each curvature value is determined, and then each curvature value is converted into the corresponding texture parameter according to the preset correspondence standard between the texture level and the texture parameter, so as to ensure that each curvature value can be accurately mapped to the visualized texture parameter.

[0118] Step 504: Based on the association relationship, associate and bind the texture parameters of each stratigraphic interface with the corresponding stratigraphic grid cell to obtain multiple feature association results.

[0119] In this embodiment of the application, multiple previously obtained association relationships are first retrieved, and the stratigraphic grid unit corresponding to the sub-interval in each association relationship is identified. Then, the texture parameters of each stratigraphic interface are accurately associated and bound with the corresponding stratigraphic grid unit in the association relationship to ensure that each stratigraphic grid unit can match the corresponding texture parameters, and finally, multiple feature association results are obtained.

[0120] Step 505: Add a corresponding stratigraphic number to each feature association result, and combine it with the regional cell number of each stratigraphic grid cell to generate surface feature data corresponding to each stratigraphic grid cell.

[0121] In this step, the stratigraphic number refers to a unique identifier assigned to each stratigraphic interface to distinguish different stratigraphic interfaces.

[0122] In this embodiment of the application, firstly, a corresponding stratigraphic number is added to each feature association result to clarify the stratigraphic interface to which each feature association result belongs. Then, the feature association result with the added stratigraphic number is associated and integrated with the regional unit number of each stratigraphic grid unit so that the attribute information of each stratigraphic grid unit is completely associated with the corresponding texture parameters and stratigraphic number, and finally, the surface feature data corresponding to each stratigraphic grid unit is generated.

[0123] The embodiments of this application realize the precise spatial correlation between curvature values ​​and stratigraphic grid cells, ensuring that the surface feature data of each stratigraphic grid cell is consistent with the actual curvature features of the stratigraphic interface, providing realistic surface detail support for the three-dimensional geological model, and solving the problems of disordered binding of curvature values ​​and stratigraphic grid cells and inaccurate surface feature restoration in the prior art.

[0124] This application provides a specific embodiment. Step 106 involves generating construction identification data based on the model position coordinates of each curvature abrupt change point within the 3D model framework, combined with a preset construction type label. This specifically includes the following steps:

[0125] Step 601: Based on the first vertical coordinate range of each stratigraphic interface, divide the first vertical coordinate range into multiple vertical coordinate intervals, and determine the vertical coordinate interval to which each curvature change point belongs.

[0126] In this embodiment of the application, the first vertical coordinate range of each stratigraphic interface is first determined, and the first vertical coordinate range is evenly divided into multiple continuous vertical coordinate intervals according to the preset division standard to ensure that the range of each vertical coordinate interval is clear and does not overlap. Then, the vertical coordinate value of each curvature change point is checked one by one to determine which vertical coordinate interval it falls into, and finally the vertical coordinate interval to which each curvature change point belongs is determined.

[0127] Step 602: Calculate the distance between the position coordinates of the curvature change point in each vertical coordinate interval and the center coordinates of the corresponding stratigraphic grid cell in the 3D model frame. Take the center coordinates corresponding to the smallest distance as the model position coordinates of the corresponding curvature change point in the 3D model frame.

[0128] In this step, the center coordinates refer to the coordinates of the center point of each stratigraphic grid cell in the three dimensions of horizontal, vertical, and vertical coordinates.

[0129] In this embodiment, the position coordinates of all curvature change points within each vertical coordinate interval are first extracted, as well as the center coordinates of all stratigraphic grid cells corresponding to that vertical coordinate interval in the three-dimensional model frame. The three-dimensional distance between each curvature change point and the center coordinates of the corresponding stratigraphic grid cell is calculated. All distances corresponding to each curvature change point are compared, and the distance with the smallest value is selected. The center coordinates of the stratigraphic grid cell corresponding to the smallest value are used as the model position coordinates of the corresponding curvature change point in the three-dimensional model frame.

[0130] Step 603: Calculate the difference between the vertical coordinate of each curvature abrupt change point and the center vertical coordinate of the corresponding stratigraphic grid cell, and use the model position coordinates of curvature abrupt change points with a difference less than or equal to a preset threshold as the final associated coordinates.

[0131] In this step, the preset threshold refers to a pre-set critical value used to determine whether the difference between the vertical coordinate of the curvature abrupt change point and the center vertical coordinate of the stratigraphic grid cell meets the requirements.

[0132] In this embodiment, firstly, the vertical coordinates of each curvature abrupt change point and the center vertical coordinates of the stratigraphic grid cell to which its corresponding model position coordinates belong are extracted, and the absolute difference between the two is calculated. Secondly, the calculated absolute difference is compared with a preset threshold. If the difference is less than or equal to the preset threshold, the model position coordinates of the curvature abrupt change point are determined as the final associated coordinates. If the difference is greater than the preset threshold, other adjacent stratigraphic grid cells within the vertical coordinate interval of the curvature abrupt change point are selected again, and the absolute difference between the vertical coordinates of the abrupt change point and the center vertical coordinates of the new grid cell are calculated again. The matching and calculation operations are repeated until the difference is less than or equal to the preset threshold, and then the corresponding model position coordinates are determined as the final associated coordinates.

[0133] Step 604: Associate the final associated coordinates of each curvature change point, the regional cell number of the corresponding stratigraphic grid unit, the vertical coordinate interval to which it belongs, and the preset structural type label to obtain structural identification data.

[0134] In this embodiment of the application, the final associated coordinates of each curvature change point are first collected, as well as the regional cell number of the stratigraphic grid unit corresponding to the final associated coordinates and the vertical coordinate interval to which the curvature change point belongs. Then, these information are associated with the preset structural type labels one by one to ensure that the information of each curvature change point is completely consistent with the structural type label. Finally, all associated information is integrated to obtain structural identification data.

[0135] The embodiments of this application realize the precise spatial correlation between curvature change points and stratigraphic grid units, ensuring that the structural identification data includes complete location, unit affiliation and structural type information, providing accurate data support for the local structural identification of three-dimensional geological models, and solving the problems of mismatch between structural identification data and actual geological structural locations and incomplete information in the prior art.

[0136] This application provides a specific embodiment, such as Figure 2 As shown, step 107, based on the priority order of the comprehensive dataset, jointly constructs the 3D model framework, surface feature data under the normal vector constraints of each stratigraphic interface, structural identification data, and coordinate data to generate a 3D geological model of the non-ferrous metal deposit. This specifically includes the following steps:

[0137] Step 701: Determine the priority order of the comprehensive dataset based on the influence of the normal vector, curvature value, and position coordinates of curvature change points of each stratum interface on the 3D model framework.

[0138] In this step, the degree of influence refers to the strength of the effect of the normal vector, curvature value, and position coordinates of curvature change points of each stratum interface on the construction of the three-dimensional model framework.

[0139] Step 702: Adjust the spatial angles of the stratigraphic grid cells in the 3D model frame according to the priority order and the normal vector of each stratigraphic interface to obtain the adjusted 3D model frame. Combine the adjusted 3D model frame and the normal vector of each stratigraphic interface to construct the first joint model.

[0140] In this step, spatial angle refers to the tilt angle and azimuth angle of the stratigraphic grid cell relative to the coordinate system in three-dimensional space.

[0141] In this embodiment, the priority order of the determined integrated dataset is first clarified. Guided by this order, and combined with the dip direction and angle of the strata indicated by the normal vector of each stratum interface, the spatial angle of each stratum grid cell in the three-dimensional model frame is adjusted one by one so that the spatial orientation of the stratum grid cell is consistent with the actual attitude of the stratum interface, thus obtaining the adjusted three-dimensional model frame. Then, the adjusted three-dimensional model frame is deeply fused with the normal vector of each stratum interface to ensure that the constraint effect of the normal vector on the model frame is fully reflected, and finally the first joint model is obtained.

[0142] Step 703: Combine the texture parameters and structural identification data in the surface feature data with the stratigraphic grid cells corresponding to the regional unit numbers in the first joint model to construct the second joint model.

[0143] In this embodiment, firstly, all texture parameters included in the surface feature data are extracted, and the correlation between the stratigraphic grid cells corresponding to each texture parameter is clarified. Simultaneously, information such as the structure type label and final associated coordinates included in the structure identification data are organized. Then, using the region cell number as the matching basis, the texture parameters are mapped one by one to the stratigraphic grid cells with the same region cell number in the first joint model. At the same time, the structure identification data is also bound to the corresponding stratigraphic grid cells according to the region cell number. Through data fusion, the joint construction is completed, ultimately obtaining the second joint model.

[0144] Step 704: Adjust and verify the second joint model based on the coordinate data to obtain a three-dimensional geological model of the non-ferrous metal mine.

[0145] In this embodiment, the initially acquired coordinate data is retrieved first, and the model interface coordinates of each stratigraphic interface in the second joint model are compared one by one with the actual interface coordinates in the coordinate data to calculate the deviation between the two. Then, based on the deviation, the position and spatial angle of the stratigraphic grid unit in the second joint model are adjusted accordingly. At the same time, the binding accuracy of surface feature data and structural identification data is verified to ensure that the overall structure, surface details and local structures of the model are consistent with the actual geological conditions, and finally a three-dimensional geological model of the non-ferrous metal mine is obtained.

[0146] The embodiments of this application realize the orderly integration of various types of data into the modeling process according to their importance, ensuring that the overall occurrence, surface details and local structures of the three-dimensional geological model are accurately matched with the actual geological conditions, effectively solving the problems of disordered modeling data fusion and large deviations between the model and the actual geological conditions in the prior art.

[0147] This application provides a specific embodiment. Step 701 involves determining the priority order of the comprehensive dataset based on the influence of the normal vector, curvature value, and location coordinates of curvature abrupt change points of each stratigraphic interface on the three-dimensional model framework. This specifically includes the following steps:

[0148] Step 711: Based on the degree of core role of the normal vector, curvature value, and curvature change point of each stratigraphic interface in the three-dimensional geological model, determine the first weight corresponding to each normal vector, the second weight corresponding to the curvature value, and the third weight corresponding to the curvature change point.

[0149] In this step, the degree of core role refers to the strength of the crucial supporting role of the normal vector, curvature value, and position coordinates of curvature abrupt change points of each stratigraphic interface in the construction of the 3D geological model. The first weight refers to the quantitative value assigned to the normal vector of each stratigraphic interface, reflecting its degree of core role. The second weight refers to the quantitative value assigned to the curvature value of each stratigraphic interface, reflecting its degree of core role. The third weight refers to the quantitative value assigned to the position coordinates of the curvature abrupt change points of each stratigraphic interface, reflecting its degree of core role.

[0150] In this embodiment, the key supporting role of the normal vector of each stratigraphic interface in the overall structural formation of the three-dimensional geological model, the key supporting role of the curvature value in the restoration of surface details of the three-dimensional geological model, and the key supporting role of the position coordinates of curvature change points in the identification of local structures of the three-dimensional geological model are first analyzed to clarify the degree of core role of the three. Then, according to the strength of the degree of core role, a corresponding first weight is assigned to the normal vector, a corresponding second weight is assigned to the curvature value, and a corresponding third weight is assigned to the position coordinates of the curvature change point to ensure that the weight values ​​are positively correlated with the degree of core role.

[0151] Step 712: Calculate the first deviation of the overall attitude of the 3D model frame when the normal vector of each stratigraphic interface changes, and combine it with the first weight to calculate the first comprehensive weight of the normal vector of each stratigraphic interface on the 3D model frame.

[0152] In this step, the change in the normal vector of the stratigraphic interface refers to the change in the direction or value of the normal vector of each stratigraphic interface. The overall attitude of the 3D model framework refers to the overall orientation of the 3D model framework in 3D space, such as the tilt direction and tilt angle. The first deviation refers to the degree of deviation between the overall attitude of the 3D model framework and the original design attitude when the normal vector of the stratigraphic interface changes.

[0153] In this embodiment, the change range of the normal vector of the stratigraphic interface is first set to simulate the scenario after the change of the normal vector; then the first deviation between the overall attitude of the three-dimensional model frame after the change and the original design attitude is calculated. The calculation process of the first deviation is as follows: firstly, the square of the first difference between the changed inclination angle and the original inclination angle is calculated to obtain the first square; the square of the second difference between the changed azimuth angle and the original azimuth angle is calculated to obtain the second square; the first square and the second square are added together to obtain the sum.

[0154] Finally, the square root of the sum is taken to obtain the first deviation; then the first deviation is multiplied by the corresponding first weight to calculate the first comprehensive weight of the normal vector of each stratum interface to the three-dimensional model frame.

[0155] Step 713: Calculate the second deviation of the surface detail features of the 3D model frame when the curvature value of each stratum interface changes, and combine it with the second weight to calculate the second comprehensive weight of the curvature value of each stratum interface on the 3D model frame.

[0156] In this step, the surface detail features of the 3D model frame refer to subtle features such as roughness and texture morphology of the 3D model frame surface. The second deviation refers to the degree of deviation between the surface detail features of the 3D model frame and the original design detail features when the curvature value of the ground interface changes.

[0157] In this embodiment, the range of curvature value variation of the stratigraphic interface is first set to simulate the scenario after the curvature value changes; then, the second deviation between the surface detail features of the changed 3D model frame and the original design detail features is calculated. The calculation process for the second deviation is as follows:

[0158] First, calculate the difference between the changed surface detail feature parameters and the original surface detail feature parameters and take the absolute value to obtain the first result; then divide the first result by the original surface detail feature parameters to obtain the second result; finally, multiply the second result by 100% to obtain the second deviation; then multiply the second deviation by the corresponding second weight to finally calculate the second comprehensive weight of the curvature value of each stratigraphic interface on the three-dimensional model framework.

[0159] Step 714: Calculate the third deviation of the local structural identifier of the 3D model frame when the position coordinates of the curvature change points of each stratigraphic interface shift, and combine the third weight to calculate the third comprehensive weight of the position coordinates of each curvature change point on the 3D model frame.

[0160] In this step, the local structural identifier of the 3D model frame refers to the identification information of local geological structures such as faults and joints marked in the 3D model frame. The third deviation refers to the degree of deviation between the local structural identifier of the 3D model frame and the original design structural identifier when the position coordinates of the curvature change point of the stratum interface shift.

[0161] In this embodiment, the offset of the position coordinates of the curvature abrupt change point is first set to simulate the scene after the position coordinates are offset; then, the third deviation between the local structural identifier of the 3D model frame after the offset and the original design structural identifier is calculated. The calculation process of the third deviation is as follows:

[0162] First, calculate the first difference between the offset x-coordinate and the original x-coordinate, the second difference between the offset y-coordinate and the original y-coordinate, and the third difference between the offset y-coordinate and the original y-coordinate. Then, calculate the square of the first difference, the square of the second difference, and the square of the third difference, respectively, to obtain the first square. Finally, add the first square, the second square, and the third square and perform an open operation to obtain the third deviation. Then, multiply the third deviation by the corresponding third weight to finally calculate the third comprehensive weight of the position coordinates of each curvature change point relative to the 3D model frame.

[0163] Step 715: Determine the priority order of the comprehensive dataset based on the first comprehensive weight, the second comprehensive weight, and the third comprehensive weight.

[0164] In this embodiment of the application, the first comprehensive weight, the second comprehensive weight, and the third comprehensive weight are first collected and compared one by one, and sorted in descending order of value; then the priority order of the comprehensive dataset is determined according to the sorting result, that is, the larger the comprehensive weight value, the higher the priority of the corresponding data in the comprehensive dataset, to ensure that the priority order is consistent with the actual influence of the data on the three-dimensional model framework.

[0165] This application embodiment realizes the objective quantitative determination of the priority of the comprehensive dataset, ensuring that the priority order is consistent with the actual influence of the data on the 3D model framework, providing clear guidance for the orderly fusion of subsequent modeling data, and solving the problems of subjective judgment of the priority of the comprehensive dataset and lack of orderly constraints in data fusion in the prior art.

[0166] This application provides a specific embodiment. Step 704 involves adjusting and verifying the second joint model based on coordinate data to obtain a three-dimensional geological model of the non-ferrous metal ore. This specifically includes the following steps:

[0167] Step 721: Calculate the deviation between the model interface coordinates of each stratigraphic interface in the second combined model and the corresponding actual interface coordinates in the coordinate data.

[0168] In this step, the preset deviation threshold refers to a pre-set critical value used to determine whether the deviation between the model interface coordinates and the actual interface coordinates meets the requirements.

[0169] In this embodiment, the model interface coordinates of each stratigraphic interface in the second combined model are first extracted, and the corresponding actual interface coordinates in the coordinate data are retrieved simultaneously. The deviation between the two is calculated, and the calculation process for the deviation is as follows:

[0170] First, calculate the first difference, the second difference, and the third difference between the model interface coordinates and the actual interface coordinates in the horizontal, vertical, and vertical dimensions, respectively. Then, calculate the square of the first difference to obtain the first square, the square of the second difference to obtain the second square, and the square of the third difference to obtain the third square. Next, add the first square, the second square, and the third square and then perform a square root operation to obtain the deviation value of each stratum interface.

[0171] Step 722: Adjust the model interface coordinates of the formation interfaces in the second combined model where the deviation value is greater than the preset deviation threshold until all deviation values ​​are less than or equal to the preset deviation threshold, and obtain the adjusted second combined model.

[0172] In this embodiment, the deviation value of each stratigraphic interface is first compared with a preset deviation threshold one by one, and stratigraphic interfaces with deviation values ​​greater than the preset deviation threshold are selected. The model interface coordinates of these stratigraphic interfaces are then adjusted accordingly. During the adjustment process, the deviation value is recalculated in real time, and the optimization is repeated until the deviation values ​​of all stratigraphic interfaces are less than or equal to the preset deviation threshold. Finally, the adjusted second joint model is obtained, ensuring that the model interface coordinates and the actual interface coordinates are accurately matched.

[0173] Step 723: Calculate the boundary coincidence parameters of any two adjacent stratigraphic grid cells in the adjusted second joint model, and verify the boundary coincidence parameters according to the preset coincidence conditions to obtain the coincidence verification results.

[0174] In this step, adjacent stratigraphic mesh cells refer to stratigraphic mesh cells that are spatially adjacent and have contiguous coordinate ranges within the 3D model framework. Preset coincidence conditions refer to pre-defined standards used to determine whether the boundary coincidence parameters of adjacent stratigraphic mesh cells meet the requirements.

[0175] In this embodiment, firstly, any two adjacent stratigraphic grid cells in the adjusted second joint model are selected, and their boundary coordinate ranges are extracted. The calculation process for the boundary coincidence parameter is as follows: First, the boundary coincidence length and total boundary length of the adjacent stratigraphic grid cells are measured respectively. Then, the boundary coincidence length is divided by the total boundary length to obtain the quotient. Finally, the quotient is multiplied by 100% to obtain the boundary coincidence parameter. Then, the calculated boundary coincidence parameter is compared with the preset coincidence conditions to determine whether it meets the requirements, and the coincidence verification result is obtained to ensure that the boundaries of adjacent stratigraphic grid cells are completely connected.

[0176] Step 724: Calculate the texture transition parameters between the texture parameters of any two adjacent stratigraphic grid cells in the adjusted second joint model, and verify the texture transition parameters according to the preset transition conditions to obtain the transition verification results.

[0177] In this step, the preset transition condition refers to the pre-set standard used to determine whether the transition of texture parameters between adjacent strata grid cells is smooth.

[0178] In this embodiment, the texture parameters corresponding to any two adjacent stratigraphic grid units in the adjusted second joint model are first extracted, and the texture transition parameters between the two are calculated. The transition parameters are equal to the absolute value of the difference between the texture parameters of the two adjacent stratigraphic grid units. Then, the texture transition parameters are compared with the preset transition conditions to determine whether the transition is smooth and to obtain the transition verification result, so as to ensure that the texture connection of the model surface is natural.

[0179] Step 725: Based on the preset corresponding conditions, determine the correspondence between the structure type label of the stratigraphic grid unit with structure identifier in the adjusted second joint model and the final associated coordinate, and obtain the relationship analysis results.

[0180] In this step, the preset matching conditions refer to the pre-defined rules used to determine whether the structure type label matches the final associated coordinates. A stratigraphic grid cell with a structure identifier refers to a stratigraphic grid cell that is bound to both a structure type label and the final associated coordinates.

[0181] In this embodiment of the application, all stratigraphic grid cells with structural identifiers in the adjusted second joint model are first screened out, the structural type label and the corresponding final associated coordinates of each cell are extracted, and the correspondence between the two is analyzed one by one according to the preset corresponding conditions to determine whether the label and coordinates match, so as to obtain the relationship analysis results and ensure that the structural identifiers are accurate.

[0182] Step 726: If the overlap verification result, transition verification result, and relationship analysis result are all passed, then the adjusted second joint model is used as the three-dimensional geological model of the non-ferrous metal mine. If any result is failed, then the second joint model is readjusted until the overlap verification result, transition verification result, and relationship analysis result are all passed.

[0183] In this embodiment, the overlap verification results, transition verification results, and relationship analysis results are first summarized and comprehensively judged. If all three results pass, the adjusted second joint model is directly determined as the three-dimensional geological model of the non-ferrous metal mine. If any result fails, the process returns to the previous adjustment stage, and the coordinates, boundaries, or data binding relationships of the second joint model are re-optimized. All verification steps are executed again until all three results pass, ensuring that the final model meets all preset standards.

[0184] The embodiments of this application achieve comprehensive compliance of the three-dimensional geological model in terms of coordinate accuracy, boundary connectivity, surface smoothness, and structural identification accuracy, providing reliable model support for non-ferrous metal mineral exploration and development, and solving the problems of single-dimensional model verification and hidden defects in the prior art.

[0185] Figure 3 This application provides a schematic diagram of a specific implementation of a three-dimensional geological modeling system for non-ferrous metal mines based on exploration line profiles, as illustrated in the embodiments of this application. Figure 3 The system may include:

[0186] The acquisition module 21 is used to acquire the coordinate data of the stratigraphic interfaces of non-ferrous metal mines on multiple exploration line profiles;

[0187] The calculation module 22 is used to perform geometric calculations on the coordinate data to obtain the normal vector and curvature value of each stratigraphic interface, and to perform numerical change analysis on all curvature values ​​of each stratigraphic interface to identify the curvature abrupt change points of the corresponding stratigraphic interface.

[0188] The allocation module 23 is used to assign priority identifiers to the normal vector, curvature value, and location coordinates of curvature change points of each stratigraphic interface to form a comprehensive dataset with priority identifiers.

[0189] Module 24 is used to construct a three-dimensional model framework based on the three-dimensional coordinate extreme values ​​of each exploration line profile;

[0190] The binding module 25 is used to bind the curvature value of each stratum interface to the corresponding stratum grid unit in the 3D model frame using 3D texture mapping technology, so as to obtain the surface feature data corresponding to each stratum grid unit.

[0191] Module 24 is also used to generate construction identification data based on the model position coordinates of each curvature change point in the 3D model framework, combined with preset construction type labels.

[0192] Module 24 is also used to jointly construct a three-dimensional geological model of a non-ferrous metal mine based on the priority order of the comprehensive dataset, including the three-dimensional model framework, surface feature data under the normal vector constraints of each stratigraphic interface, structural identification data, and coordinate data.

[0193] This application provides an embodiment of a three-dimensional geological modeling system for non-ferrous metal mines based on exploration line profiles. This system is used to implement the aforementioned method for three-dimensional geological modeling of non-ferrous metal mines based on exploration line profiles. Therefore, the specific implementation of the three-dimensional geological modeling system for non-ferrous metal mines based on exploration line profiles can be found in the embodiment section of the aforementioned method for three-dimensional geological modeling of non-ferrous metal mines based on exploration line profiles. The specific implementation can be referred to the description of the corresponding embodiments, which will not be repeated here.

[0194] This application also provides an electronic device, including: a memory for storing a computer program; and a processor for executing the computer program to implement the steps of any of the above-described three-dimensional geological modeling methods for non-ferrous metal mines based on exploration line profiles.

[0195] This application also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of any of the above-described methods for three-dimensional geological modeling of non-ferrous metal mines based on exploration line profiles.

[0196] In one exemplary embodiment, the aforementioned computer-readable storage medium may include, but is not limited to, various media capable of storing computer programs, such as USB flash drives, read-only memory, random access memory, portable hard drives, magnetic disks, or optical disks.

[0197] Embodiments of the present invention also provide a computer program product, which includes a computer program that, when executed by a processor, implements the steps in any of the embodiments of the three-dimensional geological modeling method for non-ferrous metal mines based on exploration line profiles.

[0198] Those skilled in the art will further recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

[0199] The foregoing has provided a detailed description of a three-dimensional geological modeling method and system for non-ferrous metal mines based on exploration line profiles, as provided in this application. Specific examples have been used to illustrate the principles and implementation methods of this application. The descriptions of the embodiments above are merely for the purpose of helping to understand the method and its core ideas. It should be noted that those skilled in the art can make various improvements and modifications to this application without departing from its principles, and these improvements and modifications also fall within the protection scope of this application.

Claims

1. A method for three-dimensional geological modeling of non-ferrous metal deposits based on exploration line profiles, characterized in that, include: Obtain coordinate data of stratigraphic interfaces of non-ferrous metal deposits on multiple exploration profiles; Geometric calculations are performed on the coordinate data to obtain the normal vector and curvature value of each stratigraphic interface. Numerical variation analysis is performed on all curvature values ​​of each stratigraphic interface to identify the curvature abrupt change points of the corresponding stratigraphic interface. Priority labels are assigned to the normal vector, curvature value, and location coordinates of curvature abrupt change points of each stratigraphic interface to form a comprehensive dataset with priority labels; A three-dimensional model framework is constructed based on the extreme values ​​of the three-dimensional coordinates of each exploration line profile; The curvature value of each stratigraphic interface is bound to the corresponding stratigraphic mesh unit in the three-dimensional model framework using three-dimensional texture mapping technology to obtain the surface feature data corresponding to each stratigraphic mesh unit. Based on the model position coordinates of each curvature abrupt change point in the 3D model framework, and combined with the preset construction type label, construction identification data is generated; Based on the priority order of the comprehensive dataset, the three-dimensional model framework, surface feature data under the normal vector constraint of each stratigraphic interface, the structural identification data, and the coordinate data are jointly constructed to generate a three-dimensional geological model of non-ferrous metal deposits. Based on the priority order of the comprehensive dataset, the three-dimensional model framework, surface feature data under the normal vector constraints of each stratigraphic interface, the structural identification data, and the coordinate data are jointly constructed to generate a three-dimensional geological model of the non-ferrous metal deposit, including: The priority order of the comprehensive dataset is determined based on the degree of influence of the normal vector, curvature value, and position coordinates of curvature change points of each stratigraphic interface on the three-dimensional model framework. Based on the priority order and the normal vector of each stratigraphic interface, the spatial angle of the stratigraphic grid unit in the three-dimensional model frame is adjusted to obtain the adjusted three-dimensional model frame. The adjusted three-dimensional model frame and the normal vector of each stratigraphic interface are jointly constructed to obtain the first joint model. The texture parameters in the surface feature data and the structural identification data are jointly constructed with the stratigraphic grid cells corresponding to the regional unit numbers in the first joint model to obtain the second joint model. Based on the coordinate data, the second joint model is adjusted and verified to obtain a three-dimensional geological model of the non-ferrous metal mine.

2. The method according to claim 1, characterized in that, Based on the extreme values ​​of the three-dimensional coordinates of each exploration line profile, a three-dimensional model framework is constructed, including: Based on the extreme values ​​of the three-dimensional coordinates of each exploration line profile, determine the spatial coordinate boundaries of all exploration line profiles; Based on the preset grid size and the coordinate data, the three-dimensional space within the spatial coordinate boundary is divided into multiple stratigraphic grid units corresponding to each stratigraphic interface. Based on the stratigraphic interface, the distribution range of the strata within the stratigraphic profile, and the exploration line number of each exploration line profile, the coordinate data of adjacent exploration line profiles are interpolated and expanded to obtain the attribute information of each stratigraphic grid cell. The attribute information includes the regional cell number, the regional coordinate range, and the exploration line number. Based on the attribute information of all stratigraphic grid units, spatial coordinate matching is performed on all stratigraphic grid units using the regional coordinate range in the attribute information. At the same time, consistency verification is performed on the regional unit number, regional coordinate range, and exploration line number in the attribute information of each stratigraphic grid unit to obtain the target stratigraphic grid unit. The target stratigraphic grid units are then spliced ​​and integrated to obtain the three-dimensional model framework.

3. The method according to claim 1, characterized in that, The curvature value of each stratigraphic interface is bound to the corresponding stratigraphic mesh cell in the 3D model framework using 3D texture mapping technology to obtain the surface feature data corresponding to each stratigraphic mesh cell, including: Based on the curvature value of each stratigraphic interface, the range of the first vertical coordinate of each stratigraphic interface is divided to obtain multiple sub-intervals; The overlap between the sub-interval of each stratigraphic interface and the second vertical coordinate range of the corresponding stratigraphic grid cell in the three-dimensional model frame is calculated. Sub-intervals with an overlap greater than or equal to a preset overlap threshold are associated with the corresponding stratigraphic grid cells to obtain multiple association relationships. According to the preset texture mapping rules, the texture level corresponding to each curvature value is determined by the three-dimensional texture mapping technology, and each curvature value is converted into a texture parameter according to the texture level corresponding to each curvature value. Based on the aforementioned association, the texture parameters of each stratigraphic interface are associated and bound with the corresponding stratigraphic grid cell to obtain multiple feature association results; Add a corresponding stratigraphic number to each feature association result, and combine it with the regional cell number of each stratigraphic grid cell to generate surface feature data corresponding to each stratigraphic grid cell.

4. The method according to claim 1, characterized in that, Based on the model position coordinates of each curvature abrupt change point within the 3D model framework, and combined with preset construction type labels, construction identification data is generated, including: Based on the first vertical coordinate range of each stratigraphic interface, the first vertical coordinate range is divided into multiple vertical coordinate intervals, and the vertical coordinate interval to which each curvature change point belongs is determined. Calculate the distance between the position coordinates of the curvature change point in each vertical coordinate interval and the center coordinates of the corresponding stratigraphic grid cell in the 3D model frame. Take the center coordinates corresponding to the smallest distance as the model position coordinates of the corresponding curvature change point in the 3D model frame. Calculate the difference between the vertical coordinate of each curvature abrupt change point and the center vertical coordinate of the corresponding stratigraphic grid cell, and use the model position coordinates of curvature abrupt change points with a difference less than or equal to a preset threshold as the final associated coordinates; The final associated coordinates of each curvature abrupt change point, the corresponding stratigraphic grid cell regional cell number, the vertical coordinate interval to which it belongs, and the preset structural type label are associated to obtain structural identification data.

5. The method according to claim 1, characterized in that, Based on the degree of influence of the normal vector, curvature value, and location coordinates of curvature abrupt change points of each stratigraphic interface on the 3D model framework, the priority order of the comprehensive dataset is determined, including: Based on the degree of core role of the normal vector, curvature value, and curvature change point of each stratigraphic interface in the three-dimensional geological model, the first weight corresponding to each normal vector, the second weight corresponding to the curvature value, and the third weight corresponding to the curvature change point are determined. Calculate the first deviation of the overall attitude of the three-dimensional model frame when the normal vector of each stratigraphic interface changes, and combine it with the first weight to calculate the first comprehensive weight of the normal vector of each stratigraphic interface on the three-dimensional model frame. Calculate the second deviation of the surface detail features of the three-dimensional model frame when the curvature value of each stratum interface changes, and combine it with the second weight to calculate the second comprehensive weight of the curvature value of each stratum interface on the three-dimensional model frame; Calculate the third deviation of the local structural identifier of the three-dimensional model frame when the position coordinates of the curvature change points of each stratigraphic interface are offset, and combine the third weight to calculate the third comprehensive weight of the position coordinates of each curvature change point to the three-dimensional model frame. The priority order of the comprehensive dataset is determined based on the first comprehensive weight, the second comprehensive weight, and the third comprehensive weight.

6. The method according to claim 1, characterized in that, Based on the coordinate data, the second joint model is adjusted and verified to obtain a three-dimensional geological model of the non-ferrous metal mine, including: Calculate the deviation between the model interface coordinates of each stratigraphic interface in the second joint model and the corresponding actual interface coordinates in the coordinate data; The model interface coordinates of the strata interfaces in the second joint model with deviation values ​​greater than the preset deviation threshold are adjusted until all deviation values ​​are less than or equal to the preset deviation threshold, thus obtaining the adjusted second joint model. The boundary coincidence parameters of any two adjacent stratigraphic grid cells in the adjusted second joint model are calculated, and the boundary coincidence parameters are verified according to the preset coincidence conditions to obtain the coincidence verification results. Calculate the texture transition parameters between any two adjacent stratigraphic grid cells in the adjusted second joint model, and verify the texture transition parameters according to the preset transition conditions to obtain the transition verification results; Based on the preset corresponding conditions, the correspondence between the structure type labels of the stratigraphic grid cells with structure identifiers in the adjusted second joint model and the final associated coordinates is determined, and the relationship analysis results are obtained. If the overlap verification result, transition verification result, and relationship analysis result are all passed, the adjusted second joint model will be used as the three-dimensional geological model of the non-ferrous metal mine. If any result is failed, the second joint model will be readjusted until the overlap verification result, transition verification result, and relationship analysis result are all passed.

7. A three-dimensional geological modeling system for non-ferrous metal mines based on exploration line profiles, characterized in that, include: The acquisition module is used to acquire the coordinate data of the stratigraphic interfaces of non-ferrous metal mines on multiple exploration profiles. The calculation module is used to perform geometric calculations on the coordinate data to obtain the normal vector and curvature value of each stratigraphic interface, and to perform numerical change analysis on all curvature values ​​of each stratigraphic interface to identify the curvature abrupt change points of the corresponding stratigraphic interface. The allocation module is used to assign priority identifiers to the normal vector, curvature value, and location coordinates of curvature abrupt change points of each stratigraphic interface to form a comprehensive dataset with priority identifiers. The building module is used to construct a three-dimensional model framework based on the extreme values ​​of the three-dimensional coordinates of each exploration line profile; The binding module is used to bind the curvature value of each stratum interface to the corresponding stratum grid unit in the three-dimensional model framework using three-dimensional texture mapping technology, so as to obtain the surface feature data corresponding to each stratum grid unit. The construction module is also used to generate construction identification data based on the model position coordinates in the 3D model framework for each curvature abrupt change point, combined with preset construction type labels; The construction module is also used to jointly construct the three-dimensional model framework, surface feature data under the normal vector constraints of each stratigraphic interface, the structural identification data, and the coordinate data based on the priority order of the comprehensive dataset, thereby generating a three-dimensional geological model of the non-ferrous metal mine. Based on the priority order of the comprehensive dataset, the three-dimensional model framework, surface feature data under the normal vector constraints of each stratigraphic interface, the structural identification data, and the coordinate data are jointly constructed to generate a three-dimensional geological model of the non-ferrous metal deposit, including: The priority order of the comprehensive dataset is determined based on the degree of influence of the normal vector, curvature value, and position coordinates of curvature change points of each stratigraphic interface on the three-dimensional model framework. Based on the priority order and the normal vector of each stratigraphic interface, the spatial angle of the stratigraphic grid unit in the three-dimensional model frame is adjusted to obtain the adjusted three-dimensional model frame. The adjusted three-dimensional model frame and the normal vector of each stratigraphic interface are jointly constructed to obtain the first joint model. The texture parameters in the surface feature data and the structural identification data are jointly constructed with the stratigraphic grid cells corresponding to the regional unit numbers in the first joint model to obtain the second joint model. Based on the coordinate data, the second joint model is adjusted and verified to obtain a three-dimensional geological model of the non-ferrous metal mine.

8. An electronic device, characterized in that, include: Memory, used to store computer programs; A processor, configured to execute the computer program to implement the steps of a three-dimensional geological modeling method for non-ferrous metal mines based on exploration line profiles as described in any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, enables a three-dimensional geological modeling method for non-ferrous metal mines based on exploration line profiles as described in any one of claims 1 to 6.

Images