Ore body virtual sectioning and drilling data generation method and ore body contour line intelligent generation method

Abstract

The invention discloses an ore body virtual sectioning and drilling data generation method and an ore body contour line intelligent generation method. The method comprises the steps that a three-dimensional model file of an ore body model and corresponding sectioning / drilling parameters are obtained; extracting all the three-dimensional entities and converting the three-dimensional entities into a polyhedral data structure, calculating a minimum bounding rectangle of the polyhedral data structure and performing model alignment; performing equal-interval sectioning along the long side direction of the model and extracting to obtain a corresponding contour line; performing contour line rendering and contour boundary calculation on each section, and performing calculation according to drilling parameters to obtain a drilling position so as to realize simulation of a drilling path; a sectioning result file and a corresponding image file are generated, and ore body virtual sectioning and drilling data generation are completed. According to the method, the generation of the ore body virtual sectioning and drilling data and the training of the corresponding ore body contour line intelligent generation model can be realized, and the scheme is higher in reliability and better in accuracy.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of mines, and particularly relates to a method for virtually cutting a mineral body and generating drilling data and a method for intelligently generating a mineral body contour line. BACKGROUND

[0002] In the field of mineral exploration technology, the delineation of a mineral body contour line has always been a key work in the field of mineral exploration. Therefore, the accurate and reliable delineation of a mineral body contour line is of great significance.

[0003] At present, the traditional delineation scheme of a mineral body contour line still adopts a manual delineation scheme; that is, a geological engineer delineates and draws a mineral body contour line according to the obtained cutting and drilling data of the mineral body, according to the geological rules and working experience. However, such a manual delineation scheme is not only time-consuming and laborious, but also has relatively poor reliability and accuracy.

[0004] In recent years, artificial intelligence schemes have developed rapidly, and therefore some researchers have applied artificial intelligence schemes to the field of mineral body contour line generation and have achieved certain results. However, the performance of an artificial intelligence scheme is directly related to the training data of the corresponding mineral body contour line; the more data and the higher the data quality, the better the performance of the corresponding artificial intelligence scheme obtained by training. However, the real data of a mineral body contour line (including drilling data and corresponding contour line data) is extremely scarce at present, and the data quality is poor, which seriously affects the application of artificial intelligence schemes in the field of mineral body contour line generation. SUMMARY

[0005] One of the purposes of the present application is to provide a method for virtually cutting a mineral body and generating drilling data, which has high reliability and good accuracy.

[0006] The second purpose of the present application is to provide a method for intelligently generating a mineral body contour line, which includes the method for virtually cutting a mineral body and generating drilling data.

[0007] The method for virtually cutting a mineral body and generating drilling data provided by the present application includes the following steps:

[0008] S1. Obtain a three-dimensional model file of a mineral body model and corresponding cutting / drilling parameters;

[0009] S2. Based on the three-dimensional model file obtained in step S1, extract all three-dimensional entities and convert them into a polyhedral data structure, calculate the minimum enclosing rectangle of the polyhedral data structure, and align the model;

[0010] S3. Equally interval cut along the long side direction of the model, and extract the corresponding contour line;

[0011] S4. For each profile, rendering of the contour line and calculation of the contour boundary are performed, and the drilling position is calculated according to the drilling parameters to realize simulation of the drilling path;

[0012] S5. According to the obtained data information, a profile cutting result file and a corresponding image file are generated, and the generation of the virtual profile cutting and drilling data of the ore body is completed.

[0013] The step S1 specifically comprises the following steps:

[0014] A file in DMF format of a three-dimensional model of the ore body is obtained;

[0015] The obtained profile cutting / drilling parameters include a profile interval , a drilling interval , and a pixel accuracy in the width direction .

[0016] The step S2 comprises the following steps:

[0017] The three-dimensional model file obtained in the step S1 is loaded, the three-dimensional model file is regarded as a layered database, and all layers are traversed to extract all three-dimensional entities;

[0018] Each three-dimensional entity is converted into a polyhedral data structure, corresponding boundary information is calculated, and is projected onto an XY plane to obtain a two-dimensional point set;

[0019] A minimum enclosing rectangle of the two-dimensional point set and four corresponding corner points are calculated, and a unit vector corresponding to each of the long side and the short side of the minimum enclosing rectangle is set;

[0020] A rotation matrix is constructed to map the original coordinate system of the three-dimensional model file to a target coordinate system taking the minimum enclosing rectangle as a reference to realize alignment of the model.

[0021] The step S2 specifically comprises the following steps:

[0022] The three-dimensional model file obtained in the step S1 is loaded; the three-dimensional model file is regarded as a layered database D; all layers are traversed to extract a three-dimensional entity list , wherein is the total number of three-dimensional entities;

[0023] Each three-dimensional entity is converted into a corresponding polyhedral data structure ; boundary of the polyhedral data structure is calculated, wherein , wherein is the minimum value of the x coordinates of all points in the polyhedral data structure ,​ for polyhedral data structure the maximum value of x coordinate of all points in for polyhedral data structure the minimum value of y coordinate of all points in for polyhedral data structure the maximum value of y coordinate of all points in for polyhedral data structure the minimum value of z coordinate of all points in for polyhedral data structure the maximum value of z coordinate of all points in

[0024] projecting to XY plane to get two-dimensional point set is wherein is the i-th point in two-dimensional point set is the x-axis coordinate of is the i-axis coordinate of N is the total number of points in two-dimensional point set

[0025] calculating the minimum enclosing rectangle of two-dimensional point set the four corner points of the minimum enclosing rectangle are calculated by minimum area rectangle algorithm

[0026] defining the first vector , the second vector wherein represents the coordinate of the point; at the same time, it is stipulated that the modulus of is the long side of the minimum enclosing rectangle, the modulus of is the short side of the minimum enclosing rectangle; normalizing the first vector and the second vector to get the first unit vector and the second unit vector ;

[0027] setting the origin of the minimum enclosing rectangle wherein is the set vector and ;

[0028] in order to align the model, setting the rotation matrix ​​​​​​​This maps the original coordinate system of the 3D model file to the target coordinate system based on the minimum bounding rectangle; in the original coordinate system, the origin... The coordinates are x-axis unit vector for y-axis unit vector for z-axis unit vector for In the target coordinate system, the origin The coordinates are x-axis Represented as y-axis Represented as z-axis Represented as ;

[0029] The rotation matrix is ​​calculated using the following formula. :

[0030] In the formula The coordinates of the point in the target coordinate system; These are the coordinates of the point in the original coordinate system;

[0031] The resulting rotation matrix Multiply by a two-dimensional point set The points in the two-dimensional point set Transform the points in the target coordinate system to obtain a two-dimensional point set in the target coordinate system. ;

[0032] Based on the two-dimensional point set of the target coordinate system Calculations yielded The minimum bounding rectangle as The primary bounding box of the target coordinate system; the primary bounding box of the target coordinate system is expanded outward by a set margin to obtain the target coordinate system bounding box. .

[0033] Step S3 includes the following steps:

[0034] According to the set profile interval The number of cutting planes is calculated along the long side of the model.

[0035] The model cutting function is used to cut the polyhedron data structure at equal intervals, and the corresponding polyline set is extracted.

[0036] Insert a rectangular boundary into the polyline set and extract the corresponding contour line.

[0037] The step S3 specifically comprises the following steps:

[0038] cutting along the long side direction of the model at equal intervals; the number of cutting planes is expressed as wherein is the length of the long side of the model, is the floor symbol; meanwhile, the number of cutting planes is set as The minimum value is 1;

[0039] the initial distance is expressed as The initial distance is defined as the distance from the first cutting plane to the starting point of the long side of the model;

[0040] For the jth cutting plane, the value of j is The origin of the cutting plane is expressed as The normal vector is expressed as wherein is the origin of the minimum bounding rectangle, is the length of the short side of the minimum bounding rectangle;

[0041] The working plane is set as is expressed as wherein represents any point in space, and the points satisfying the equation constitute the working plane ;

[0042] The model cutting function is used for cutting to obtain the output polyhedral data structure When cutting, the forward cutting is set, and the tolerance is 0;

[0043] The multi-segment line set is extracted from the polyhedral data structure , The element in is a point sequence;

[0044] The orthogonal projection of the center point c on the cutting plane is calculated by the following formula:

[0045] wherein is the projection center point; is the working plane normal vector; represents the modulus of , that is, the length of the vector;

[0046] Finally, the multi-segment of the rectangular boundary is obtained, and the corresponding contour line is extracted.

[0047] The step S4 specifically comprises the following steps.

[0048] For each profile, the resolution and the pixel of the image are calculated.

[0049] For each profile, the contour image RGB channel matrix, the drilling image RGB channel matrix and the combined image RGB channel matrix are constructed.

[0050] For each element in the multi-segment line set , the line segment endpoints are traversed to map to the pixel coordinates, so as to realize the rendering of the contour line. The contour boundary is calculated, and the number of drillings is calculated according to the contour boundary and the set drilling interval.

[0051] For each drilling, the start point and the end point of the vertical line segment are set, and the intersection set of the vertical line segment and all contour lines is solved.

[0052] According to the obtained intersection set, the intersection points are connected to complete the generation and rotation of the drilling segment.

[0053] Finally, the drilling line is obtained and saved to complete the simulation of the drilling path.

[0054] The step S4 specifically comprises the following steps.

[0055] For each profile, the image resolution is calculated: the image width is

[0056] , the image height is , the width direction pixel size is , and the height direction pixel size is ; wherein is the maximum value of the x-axis coordinate of the target coordinate system frame , the minimum value of the x-axis coordinate of the target coordinate system frame , the maximum value of the y-axis coordinate of the target coordinate system frame , the minimum value of the y-axis coordinate of the target coordinate system frame , the image width , and the image height ; the symbol is the upward rounding symbol. ​​​​​

[0057] For each profile, construct the RGB channel matrix of the contour image. RGB channel matrix of borehole image and combined image RGB channel matrix ;

[0058] Rendering outlines: For polyline sets Each element in traversal endpoints of all line segments And mapped to the pixel x-coordinate and pixel ordinate : , ,in For the endpoints of the line segment Quantity, For the endpoints of the line segment Quantity, intermediate variable and Then, the Bresenham algorithm is used to draw lines of the specified color, completing the rendering of the outline.

[0059] Calculate the contour boundary for ,in For the outline boundary exist Minimum value in the direction, For the outline boundary exist Maximum value in direction, For the outline boundary exist Minimum value in the direction, For the outline boundary exist Maximum value in the direction; calculate the number of boreholes. for And set the number of holes The minimum value is 1;

[0060] Starting drilling distance Represented as Starting drilling distance Define a two-dimensional coordinate system along the X direction of the cutting plane, starting from the minimum coordinates. Initial offset distance to the first borehole, starting borehole distance Used to determine the starting point of the drilling path;

[0061] For the m-th borehole, the starting point s of the vertical line segment is defined as... The endpoint e is ,in It is the second intermediate variable and ;

[0062] Calculate the set of intersection points of the vertical line segment with all contour lines. ;Will Along the Y direction, starting from the starting point s, sort the points to obtain the set of intersection points. ;

[0063] Make a judgment:

[0064] like If the number is even, connect the intersection points in pairs to draw drill segments of the set color, and add gray line segments between adjacent pairs and at the boundaries; then, add drills at the left and right boundaries with a set margin to ensure the integrity of the covered area; where, Indicates the number of intersections between the vertical line segment and the outline;

[0065] like If the number is odd, it is considered an anomaly, and the tolerance is increased along the X direction of the cutting plane. - ) / 100, and recalculate;

[0066] Finally, the drill line is obtained and saved to complete the simulation of the drill path.

[0067] The present invention also provides an intelligent method for generating ore body contour lines, which includes the aforementioned method for virtual ore body sectioning and borehole data generation, comprising the following steps:

[0068] S1. Using the aforementioned method for virtual ore body sectioning and borehole data generation, a training dataset is generated;

[0069] S2. Based on the diffusion model and Unet network, construct an initial model for intelligent generation of ore body contour lines;

[0070] S3. Using the training dataset obtained in step S1, train the initial intelligent generation model of ore body contour lines constructed in step S2 to obtain the intelligent generation model of ore body contour lines.

[0071] S4. Using the intelligent generation model of the ore body contour obtained in step S3, intelligently generate the contour of the actual target ore body.

[0072] The method for generating virtual ore body profile and borehole data, as well as the intelligent generation method for ore body contour lines provided by this invention, through the conversion of three-dimensional models, profile simulation, borehole simulation, rendering, and image generation, can not only realize the generation of virtual ore body profile and borehole data and the training of the corresponding intelligent generation model for ore body contour lines, but also the solution of this invention has higher reliability and better accuracy. Attached Figure Description

[0073] Figure 1 The method flow diagram of the data generation method of the present application.

[0074] Figure 2 The virtual sectioning diagram of the ore body model of the embodiment of the data generation method of the present application; wherein Figure 2 (a) is a diagram of the three-dimensional model of the ore body, Figure 2 (b) is a diagram of the virtual sectioning of the ore body model.

[0075] Figure 3 The generation effect diagram of the No. 3 sectioning surface of the embodiment of the data generation method of the present application; wherein, Figure 3 (a) is a diagram of the sectioning contour line of the ore body, Figure 3 (b) is a diagram of the virtual drilling generation process, Figure 3 (c) is a diagram of the virtual drilling generation result.

[0076] Figure 4 The generation effect diagram of the No. 6 sectioning surface of the embodiment of the data generation method of the present application; wherein, Figure 4 (a) is a diagram of the sectioning contour line of the ore body, Figure 4 (b) is a diagram of the virtual drilling generation process, Figure 4 (c) is a diagram of the virtual drilling generation result.

[0077] Figure 5 The method flow diagram of the contour line intelligent generation method of the present application.

[0078] Figure 6 The comparison diagram of the contour line intelligent generation method of the present application and the artificial interpretation contour line; wherein, Figure 6 (a) is a diagram of the drilling of drilling 1, Figure 6 (b) is a diagram of the contour line of drilling 1 generated by the present application, Figure 6 (c) is a diagram of the contour line of drilling 1 interpreted by a professional engineer, Figure 6 (d) is a diagram of the drilling of drilling 2, Figure 6 (e) is a diagram of the contour line of drilling 2 generated by the present application, Figure 6 (f) is a diagram of the contour line of drilling 2 interpreted by a professional engineer. DETAILED DESCRIPTION

[0079] As Figure 1 shown is the method flow diagram of the data generation method of the present application: the ore body virtual sectioning and drilling data generation method disclosed by the present application comprises the following steps:

[0080] S1. Obtain the three-dimensional model file of the ore body model and the corresponding sectioning / drilling parameters; specifically comprising the following steps:

[0081] obtaining a three-dimensional model file of the ore body model in DMF format;

[0082] obtaining sectioning / drilling parameters, including section interval , drilling interval , and width direction pixel precision .

[0083] S2. Based on the three-dimensional model file obtained in step S1, extract all three-dimensional entities and convert them into a polyhedral data structure, calculate the minimum bounding rectangle of the polyhedral data structure, and align the model; including the following steps:

[0084] Load the three-dimensional model file obtained in step S1, treat the three-dimensional model file as a layered database, and traverse all layers to extract all three-dimensional entities;

[0085] Convert each three-dimensional entity into a polyhedral data structure, calculate the corresponding boundary information, and project it onto the XY plane to obtain a two-dimensional point set;

[0086] Calculate the minimum bounding rectangle and the corresponding four corner points of the two-dimensional point set, and set the unit vectors corresponding to the long and short sides of the minimum bounding rectangle, respectively;

[0087] Construct a rotation matrix to map the original coordinate system of the three-dimensional model file to the target coordinate system based on the minimum bounding rectangle to achieve model alignment.

[0088] In specific implementation, the following steps can be used:

[0089] Load the three-dimensional model file obtained in step S1; treat the three-dimensional model file as a layered database D; traverse all layers to extract a list of all three-dimensional entities for i = 1 to N , where N is the total number of three-dimensional entities;

[0090] Convert each three-dimensional entity into a corresponding polyhedral data structure ; calculate the boundary of the polyhedral data structure as , where is the minimum x-coordinate value of all points in the polyhedral data structure , is the maximum x-coordinate value of all points in the polyhedral data structure , is the minimum y-coordinate value of all points in the polyhedral data structure , is the maximum y-coordinate value of all points in the polyhedral data structure The maximum y-coordinate of all points. For polyhedral data structures The minimum z-coordinate of all points. For polyhedral data structures The maximum z-coordinate of all points;

[0091] Will Projecting onto the XY plane (i.e., setting the Z-axis coordinates of the left and right points to 0) yields a two-dimensional point set. for ,in Two-dimensional point set The i-th point in for x-axis coordinates for The i-axis coordinates, where N is a two-dimensional point set. The total number of points in;

[0092] Calculate a two-dimensional point set The minimum bounding rectangle; the four corner points of the minimum bounding rectangle are calculated using the minimum area rectangle algorithm. The coordinates;

[0093] Define the first vector for The second vector for ,in Represents the coordinates of a point; and also stipulates that, model Let the longest side of the smallest enclosing rectangle be . model The short side of the smallest enclosing rectangle; the first vector Second vector Normalization yields the first unit vector. Second unit vector ;

[0094] Set the origin of the minimum bounding rectangle for ,in For the given vector and ;

[0095] To align the model, set the rotation matrix. This maps the original coordinate system of the 3D model file to the target coordinate system based on the minimum bounding rectangle; in the original coordinate system, the origin... The coordinates are x-axis unit vector for y-axis unit vector for , z-axis unit vector is ; in the target coordinate system, the coordinates of the origin are , the x-axis is represented as , the y-axis is represented as , and the z-axis is represented as ;

[0096] The rotation matrix is calculated using the following formula:

[0097] In the formula, is the coordinates of the point in the target coordinate system; is the coordinates of the point in the original coordinate system; in the implementation, the direct axis alignment method can be used;

[0098] The obtained rotation matrix is multiplied by the points in the two-dimensional point set , and the points in the two-dimensional point set are converted to the target coordinate system to obtain the target coordinate system two-dimensional point set ;

[0099] According to the target coordinate system two-dimensional point set , the minimum bounding rectangle of is calculated as the target coordinate system primary bounding box of ; the target coordinate system primary bounding box is expanded outward by a set margin to obtain the target coordinate system bounding box ; for example, if the set margin is 10%, the target coordinate system primary bounding box is expanded outward by to ensure that potential cutting errors are accommodated, , is the maximum value of the x-axis coordinates of the points in the target coordinate system primary bounding box, is the minimum value of the x-axis coordinates of the points in the target coordinate system primary bounding box, is the maximum value of the y-axis coordinates of the points in the target coordinate system primary bounding box, is the minimum value of the y-axis coordinates of the points in the target coordinate system primary bounding box.

[0100] S3. Equidistant cutting is performed along the long edge direction of the model, and the corresponding contour lines are extracted; including the following steps:

[0101] According to the set section interval and the direction of the long edge of the model, the number of cutting planes is calculated; here, the direction of the long edge of the model is defined as the bounding box The direction of the longer side;

[0102] The model cutting function is used to cut the polyhedron data structure at equal intervals, and the corresponding polyline set is extracted.

[0103] Insert a rectangular boundary into the polyline set and extract the corresponding outline; the rectangular boundary is defined as an extended form (boundary box) based on the minimum bounding rectangle (MER). The generated closed polyline is defined on each cutting plane, and its purpose is to provide a reference boundary for the cutting results.

[0104] In practice, the following steps can be taken:

[0105] Perform equally spaced cuts along the long side of the model. Here, the long side of the model is defined as the bounding box obtained in step S2. The direction of the long side; the number of cutting planes Represented as ,in The length of the longer side of the model. Sets the floor sign; also sets the number of cutting planes. The minimum value is 1;

[0106] Starting distance Represented as The starting distance Defined as the distance from the first cutting plane to the starting point of the long side of the model;

[0107] For the j-th cutting plane, the value of j is... The origin of the cutting plane Represented as normal vector Represented as ;in The origin of the smallest enclosing rectangle. The length of the shorter side of the smallest enclosing rectangle;

[0108] Set the work plane , represented as ,in Represent any point in space that satisfies the equation The points constitute the working plane ;

[0109] The model is cut using a cutting function to obtain the output polyhedral data structure. The cutting is set to a positive cutting direction with a tolerance of 0; the model cutting function is defined as the thModelFunc.CutModelData function in DmPyBindInterface.

[0110] extracting a polygonal data structure from a polyhedral data structure , from a polyhedral data structure as a point sequence;

[0111] The orthogonal projection of the center point c on the cutting plane is calculated using the following formula: :

[0112] wherein is the projection center point; is the working plane normal vector; denotes the modulus, i.e. the length of the vector;

[0113] Finally, the polygonal boundary is obtained by acquiring the multi-segment line and extracting the corresponding contour line;

[0114] In this step, the working plane is defined to determine the model cutting position and generate the contour line, and the two-dimensional contour multi-segment line is extracted by calculating the focal point of the working plane and the model cutting function, and the reference frame is set by defining the working plane to provide the plane coordinate system for subsequent drilling simulation and image rendering;

[0115] In this step, the purpose of calculating is to set the positioning reference point to provide a geometric center reference point for the contour line of each cutting plane, to ensure that the spatial positioning of the contour line is consistent with the overall coordinate system of the model, and to standardize the output: the projection center point is used to assist in generating consistent DMF files and image outputs, to ensure the relative position accuracy of the contour line and the boundary multi-segment line on the plane.

[0116] S4. For each section, the rendering of the contour line and the calculation of the contour boundary are performed, and the drilling position is calculated according to the drilling parameters to realize the simulation of the drilling path; including the following steps:

[0117] For each section, the resolution and pixels of the image are calculated;

[0118] For each section, the contour image RGB channel matrix, the drilling image RGB channel matrix, and the combined image RGB channel matrix are constructed;

[0119] For each element in the multi-segment line set , the line segment endpoints are traversed to map to the pixel coordinates to realize the rendering of the contour line;

[0120] Calculate the contour boundary, and according to the contour boundary and the set drilling interval , the number of drillings is calculated;

[0121] For each drilling, set the starting point and ending point of the vertical line segment, and solve the intersection set of the vertical line segment with all contour lines;

[0122] According to the obtained intersection set, the connection of the intersection points is carried out to complete the generation and rotation of the drilling segment;

[0123] Finally, the drilling line is obtained and saved to complete the simulation of the drilling path.

[0124] In specific implementation, the following steps can be adopted:

[0125] For each profile, calculate the image resolution: image width is , image height is , width direction pixel size is , height direction pixel size is ; wherein is the maximum value of the x-axis coordinate of the target coordinate system frame , is the minimum value of the x-axis coordinate of the target coordinate system frame , is the maximum value of the y-axis coordinate of the target coordinate system frame , is the minimum value of the y-axis coordinate of the target coordinate system frame , is the image width, is the image height; is the upward rounding symbol;

[0126] For each profile, construct the contour image RGB channel matrix , drilling image RGB channel matrix and combined image RGB channel matrix ;

[0127] Render the contour line: for each element in the multi-segment line set , traverse all end points of the line segment in , and map to pixel horizontal coordinate and pixel vertical coordinate : , , wherein endpoints of the line segment component, endpoints of the line segment component, intermediate variable and ; then, a line of a set color (preferably a red line) is drawn by using the Bresenham algorithm (i.e. the line function of OpenCV); the rendering of the contour line is completed;

[0128] the contour boundary is calculated as , wherein is the contour boundary in the direction of , the minimum value, is the contour boundary in the direction of , the maximum value, is the contour boundary in the direction of , the minimum value, is the contour boundary in the direction of , the maximum value; the number of drill holes is calculated as , and the number of drill holes is set to be at least 1;

[0129] the starting drill hole distance is represented as , the starting drill hole distance defines a two-dimensional X direction along the section plane, from the minimum coordinate to the initial offset distance of the first drill hole, the starting drill hole distance is used to determine the starting point of the drill hole path;

[0130] for the mth drill hole, the starting point s of the vertical line segment is defined as , and the ending point e is , wherein is a second intermediate variable and ;

[0131] the intersection set of the vertical line segment and all contour lines is calculated as , in the specific implementation, the Intersect function can be used to find the intersection of the parameterized line segments, and the tolerance is set to be ; the is sorted along the Y direction (the Y direction refers to the vertical axis of the three-dimensional ore body model on the section plane , which corresponds to the Y coordinate axis of the model (in the aligned coordinate system)), starting from the starting point s of ;

[0132] Judgment is made:

[0133] If is even, then the paired connection intersection points are drawn to draw the drilling section of the set color (preferably red), and the adjacent pairs and the boundary part are supplemented with a gray line segment; then, the drilling of the left and right boundaries with a set margin (preferably 5%) is supplemented to ensure the integrity of the coverage area; wherein, represents the number of intersection points of the vertical line segment and the contour line;

[0134] If is odd, it is considered abnormal, the tolerance ( - ) / 100 is increased along the X direction of the section plane, and recalculation is performed;

[0135] Finally, the drilling line is obtained and saved to complete the simulation of the drilling path.

[0136] S5. According to the obtained data information, a section result file and a corresponding image file are generated, and the generation of the virtual section of the ore body and the drilling data is completed; in specific implementation, the generated files include a section result file (DMF format) and a JSON state file, and three pictures are generated for each section: a contour line image, a drilling image, and a contour line-drilling combination image.

[0137] The method of the application is further described below in combination with an embodiment:

[0138] Because the drilling data and the contour line data are relatively sparse in the production exploration of the mine, they cannot meet the training requirements of the artificial intelligence model; therefore, the scheme of the application can generate data by using the constructed three-dimensional ore body model, perform virtual sectioning on the ore body model, generate virtual contour line data, and generate virtual drilling according to the formulated rules on the basis of the generated contour line, and the two kinds of data are used in cooperation to meet the training requirements.

[0139] As Figure 2 shown, Figure 2 (a) is a three-dimensional ore body model that has been constructed, Figure 2 (b) is a schematic diagram of the virtual sectioning of the ore body model according to the set parameters; a total of 7 section planes are formed in the diagram, and the distance between the section planes is 50.

[0140] Figure 3 is the result of the virtual sectioning of the model according to the No. 3 section plane, Figure 3 (a) is the generated ore body contour line, Figure 3 (c) is a drilling diagram generated according to the set spacing, in which red represents the ore body part and gray represents the non-ore body part. Figure 4 is the result of the virtual sectioning of the model according to the No. 6 section plane.

[0141] By Figures 2-4 It can be seen that the data generation method of the present application has the characteristics of high precision and data consistency, and by using accurate plane sectioning and pixel coordinate mapping, the generated two-dimensional contour line and drilling path data have high resolution and geometric consistency, which is suitable for machine learning training. Moreover, the method of the present application has the characteristics of flexibility and standardization, supports adjustable section interval and drilling number, and outputs standard format (DMF and JPEG), which is convenient for integration into existing geological software and AI workflow.

[0142] As Figure 5 The method flowchart of the contour line intelligent generation method of the present application is shown in the figure: the contour line intelligent generation method of the present application including the virtual sectioning of the ore body and the drilling data generation method comprises the following steps:

[0143] S1. Using the virtual sectioning of the ore body and the drilling data generation method to generate a training data set;

[0144] S2. Based on the diffusion model and the Unet network, an initial model for intelligent generation of ore body contour lines is constructed. In specific implementation, the model proposed by Lvmin Zhang in the paper "Adding Conditional Control to Text-to-Image Diffusion Models" in 2023 can be used;

[0145] S3. Using the training data set obtained in step S1, the initial model for intelligent generation of ore body contour lines constructed in step S2 is trained to obtain a model for intelligent generation of ore body contour lines;

[0146] S4. Using the model for intelligent generation of ore body contour lines obtained in step S3, the actual contour line intelligent generation of the target ore body is carried out.

[0147] The following will combine an embodiment to explain the effect of the data generation method of the present application:

[0148] The contour line intelligent generation model formed by this method can efficiently and high-quality interpret the contour line based on the drilling data; Figure 6 The comparison between the contour line intelligently generated based on the drilling data and the contour line interpreted by human is shown. It can be seen that the contour line generated by the model not only conforms to the shape of the drilling, but also is basically the same as the contour line shape interpreted by professional engineers, fully illustrating the effectiveness and accuracy of the data generation method and the contour line intelligent generation method of the present application.

[0149] In addition, the ore body virtual sectioning and drilling data generation method of the present application can also be used for geological structure analysis and underground engineering design.

Claims

1. A method for virtual ore body sectioning and borehole data generation, characterized in that... The method comprises the following steps: S1. Obtain a three-dimensional model file of the ore body model and corresponding sectioning / drilling parameters; S2. Based on the three-dimensional model file obtained in step S1, extract all three-dimensional entities and convert them into a polyhedral data structure, calculate the minimum enclosing rectangle of the polyhedral data structure, and align the model; S3. Perform equal-interval sectioning along the long edge direction of the model, and extract the corresponding profile lines; S4. For each section, perform rendering of the profile lines and calculation of the profile boundary, and calculate the drilling position according to the drilling parameters to realize simulation of the drilling path; S5. According to the obtained data information, generate a sectioning result file and a corresponding image file, and complete generation of the ore body virtual sectioning and drilling data.

2. The method of claim 1, wherein The step S1 specifically comprises the following steps: Obtain a three-dimensional model DMF format file of the ore body model; Acquired sectioning / drilling parameters, including section interval , drilling interval , and width direction pixel accuracy .

3. The method of claim 2, wherein The step S2 comprises the following steps: Load the three-dimensional model file obtained in step S1, treat the three-dimensional model file as a layered database, traverse all layers to extract all three-dimensional entities; Convert each three-dimensional entity into a polyhedral data structure, calculate the corresponding boundary information, and project it onto the XY plane to obtain a two-dimensional point set; Calculate the minimum enclosing rectangle and the corresponding four corner points of the two-dimensional point set, and set the unit vectors corresponding to the long side and the short side of the minimum enclosing rectangle respectively; Construct a rotation matrix to map the original coordinate system of the three-dimensional model file to the target coordinate system with the minimum enclosing rectangle as the reference to realize alignment of the model.

4. The method of claim 3, wherein The step S2 specifically comprises the following steps: Loading the three-dimensional model file obtained in step S1; regarding the three-dimensional model file as a layered database D; traversing all layers to extract a list of all three-dimensional entities For wherein is the total number of three-dimensional entities; Each three-dimensional entity Transformed into the corresponding polyhedral data structure ; Computation of polyhedral data structures boundary for ,in For polyhedral data structures The minimum x-coordinate of all points. For polyhedral data structures The maximum x-coordinate of all points. For polyhedral data structures The minimum y-coordinate of all points. For polyhedral data structures The maximum y-coordinate of all points. For polyhedral data structures The minimum z-coordinate of all points. For polyhedral data structures The maximum z-coordinate of all points; Will Projecting onto the XY plane yields a two-dimensional point set. for ,in Two-dimensional point set The i-th point in for x-axis coordinates for The i-axis coordinates, where N is a two-dimensional point set. The total number of points in; Computing the minimum enclosing rectangle of a two-dimensional point set ; computing the coordinates of the four corner points of the minimum enclosing rectangle by a minimum area rectangle algorithm ​ Definition of the first vector is , the second vector is where denotes the coordinates of the points; it is also provided that the modulus of is the long side of the smallest enclosing rectangle, the modulus of is the short side of the smallest enclosing rectangle; the first vector and the second vector are normalized, obtaining the first unit vector and the second unit vector ; Setting the origin of the minimum bounding rectangle For where is a set vector and ; To align the model, set the rotation matrix , which maps the original coordinate system of the three-dimensional model file to the target coordinate system with the minimum bounding rectangle as the reference. In the original coordinate system, the coordinates of the origin are , the unit vector of the x-axis is , the unit vector of the y-axis is , and the unit vector of the z-axis is ; in the target coordinate system, the coordinates of the origin are , the x-axis is represented as , the y-axis is represented as , and the z-axis is represented as ; The rotation matrix is calculated using the following equation : wherein is the coordinate of the point in the target coordinate system; is the coordinate of the point in the original coordinate system; The resulting rotation matrix is multiplied by the points in the two-dimensional point set , converting the points in the two-dimensional point set into the target coordinate system, resulting in the target coordinate system two-dimensional point set ; According to the target coordinate system two-dimensional point set , the minimum enclosing rectangle of is calculated as the target coordinate system primary frame of ; the target coordinate system primary frame is expanded outward by a set margin to obtain the target coordinate system frame .

5. The method of claim 4, wherein The step S3 comprises the following steps: According to the set profile interval The number of profile planes is calculated in the longitudinal direction of the model. Perform equal-interval cutting using a model sectioning function to obtain the output polyhedral data structure, and extract the corresponding multi-segment line set; Insert the rectangular boundary in the multi-segment line set and extract the corresponding profile line.

6. The method of claim 5, wherein The step S3 specifically comprises the following steps: cutting along the length direction of the model at equal intervals; the number of cutting planes is expressed as wherein is the length of the long side of the model, is a floor function; meanwhile, the number of cutting planes is set as the minimum value is 1; starting distance is represented as , the starting distance is defined as the distance of the starting point of the long side of the first cut plane distance model; For the j-th section plane, j takes values from 1 to n , the origin of the section plane is denoted as , the normal vector is denoted as ; wherein is the origin of the minimum bounding rectangle, is the length of the short side of the minimum bounding rectangle; Setting a work plane , is represented as where represents any point in space that satisfies the equation The points that satisfy the equation form the work plane Adopting the model section function to cut, get the output polyhedral data structure , when cutting, set to positive cutting, and the tolerance is 0; Extracting a set of polylines from a polyhedral data structure , , , as a sequence of points; The orthogonal projection of the center point c on the cutting plane is calculated using the following formula : wherein is the projection center point; is the work plane normal vector; denotes the norm of Finally, obtain the multi-segment line of the rectangular boundary and extract the corresponding profile line.

7. The method of claim 6, wherein The step S4 comprises the following steps: For each section, calculate the resolution and pixels of the image; For each section, construct a profile image RGB channel matrix, a drilling image RGB channel matrix, and a combined image RGB channel matrix; For polyline sets Each element in The endpoints of the line segments are traversed to map to pixel coordinates in order to render the outline. calculating a contour boundary, and calculating a number of drill holes based on the contour boundary and a set drill hole spacing , calculating a number of drill holes; For each drilling, set the starting point and ending point of the vertical line segment, and solve the intersection set of the vertical line segment with all profile lines; According to the obtained intersection set, connect the intersections to complete generation and rotation of the drilling segment; Finally, obtain the drilling line and save it to complete simulation of the drilling path.

8. The method of claim 7, wherein The step S4 specifically comprises the following steps: For each profile, calculate the image resolution: image width for Image height for Pixel size in width direction for Pixel size in the height direction for ;in Border of the target coordinate system The maximum value of the x-axis coordinate. Border of the target coordinate system The minimum value of the x-axis coordinate. Border of the target coordinate system The maximum value of the y-axis coordinate. Border of the target coordinate system The minimum value of the y-axis coordinate. Image width, Image height; The rounding up symbol; For each profile, a contour image RGB channel matrix is constructed , a bore image RGB channel matrix , and a combined image RGB channel matrix ; Rendering outlines: For polyline sets Each element in traversal endpoints of all line segments And mapped to the pixel x-coordinate and pixel ordinate : , ,in For the endpoints of the line segment Quantity, For the endpoints of the line segment Quantity, intermediate variable and Then, the Bresenham algorithm is used to draw lines of the specified color, completing the rendering of the outline. calculating the contour boundary for wherein is the contour boundary is the minimum value in the direction, is the contour boundary is the maximum value in the direction, is the contour boundary is the minimum value in the direction, is the contour boundary is the maximum value in the direction; calculating the number of drill holes for and setting the number of drill holes to at least 1; Starting hole distance is represented as , starting hole distance defines a two-dimensional X direction along the plane of the cut from the minimum coordinate to the initial offset distance of the first hole, starting hole distance for determining the starting point of the hole path; For the mth borehole, define the start point s of the vertical line segment as and the end point e as where is a second intermediate variable and ; Compute the intersection set of the vertical line segment and all contour lines ; if Sort the intersection set in Y direction from the start point s of the vertical line segment ; Judge: If If the number of intersection points is even, pair the intersection points to draw the drilling section of the set color, and supplement the gray line section between adjacent pairs and the boundary part; then, supplement the drilling of the left and right boundaries with the set margin to ensure the integrity of the coverage area; wherein, represents the number of intersection points of the vertical line segment and the contour line; If If the number is odd, it is considered abnormal, increase the tolerance in the X direction along the section plane ( - ) / 100, and recalculate; Finally, obtain the drilling line and save it to complete simulation of the drilling path.

9. A method for automatically generating the contour lines of a mineral body, comprising the method for virtually sectioning a mineral body and generating borehole data according to any one of claims 1 to 8, characterized in that The method comprises the following steps: S1. Use the ore body virtual sectioning and drilling data generation method in any one of claims 1-8 to generate a training data set; S2. Based on a diffusion model and a Unet network, construct an initial model for intelligent generation of ore body profile lines; S3. Use the training data set obtained in step S1 to train the initial model for intelligent generation of ore body profile lines constructed in step S2 to obtain an intelligent generation model for ore body profile lines; S4. Using the intelligent generation model of the ore body contour obtained in step S3, intelligently generate the contour of the actual target ore body.

Images