A batch calculation method for horizontal projection of uranium drilling penetration points

Abstract

The present invention belongs to a calculation method, and specifically relates to a batch calculation method for the horizontal projection of uranium borehole penetration points. It includes: step 1, preparing a borehole database: establishing a database structure and entering borehole data; step 2, calling the borehole inclination measurement record to calculate the coordinates of each control point: calling the borehole inclination measurement record and storing the calculation result; step 3, calling the borehole penetration parameters to determine the upper and lower control points of the penetration point: calling the single borehole project penetration parameters to determine the upper and lower control numbers of the ore body, and storing the result; step 4, calculating the coordinates of the borehole penetration point and outputting the coordinates: calculating the coordinates of the borehole penetration point based on the upper and lower control points of the penetration point, and outputting the projection of the penetration point. The significant effects of the present invention are: improving the drawing efficiency of the reserve plan; saving a lot of calculation time; facilitating and quickly querying and exporting data.

Description

Technical Field

[0001] The invention belongs to a calculation method, and in particular relates to a batch calculation method for horizontal projections of uranium ore drilling holes and ore penetration points. Background Art

[0002] The drill hole penetration point is the basis for drawing the reserve plan map. It is not only used for extrapolation of the ore body plane during reserve calculation, but also can be used to guide drilling construction in a timely manner. In the actual reserve plan map drawing, due to the deviation of the inclination and azimuth during drilling construction, the drill hole penetration point usually deviates from the assumed position of the drill hole. In order to reflect the actual penetration point on the projection map, it is usually necessary to project and calculate the position of the drill hole penetration point. The previous method of borehole projection calculation is to calculate the drill hole penetration point one by one according to the given borehole coordinates, borehole inclination measurement and ore penetration position. This method consumes a lot of manpower and material resources. If the amount of data is large, it takes a long time to calculate all the drill hole penetration points. At the same time, manual participation is required throughout the calculation. If manual operation errors occur, the calculated penetration point position will be inaccurate. Summary of the invention

[0003] The present invention aims at the defects of the prior art and provides a method for batch calculating the horizontal projection of the ore-penetrating points of uranium ore drilling holes.

[0004] The present invention is implemented as follows: a method for batch calculating horizontal projections of uranium ore drilling penetration points, the method comprising the following steps:

[0005] Step 1, prepare drilling database: establish database structure and input drilling data, the database includes drilling input information structure table and drilling calculation information structure table;

[0006] Step 2, calling the borehole inclinometer record and calculating the coordinates of each control point: calling the borehole inclinometer record, using the full-angle half-distance method to calculate the relevant parameters of each control point, and calculating the coordinates of each control point according to the parameters of each control point, and storing the calculation results;

[0007] Step 3, call the drilling ore-seeing parameters to determine the upper and lower control points of the ore-seeing and cutting point: call the drilling single project ore-seeing parameters to determine the upper and lower control numbers of the ore body, and store the results;

[0008] Step 4, calculate the coordinates of the drill hole's ore-cutting point and output the coordinates: calculate the coordinates of the drill hole's ore-cutting point based on the upper and lower control points of the ore-cutting point, and output the projection of the cutting point.

[0009] The method for batch calculating horizontal projections of uranium ore drilling penetration points as described above, wherein the specific steps of step 1 are as follows:

[0010] Step 1.1, based on the Access database structure, design different types of drilling information data tables, and input information structure tables including basic information table, inclinometer record table, and ore section information table; calculation information structure tables include inclinometer calculation result table and drilling hole penetration point calculation table.

[0011] Step 1.2, the drilling data is entered into the corresponding data table in the database for storage, and the data table is a basic information table.

[0012] The method for batch calculating horizontal projections of uranium ore drill hole penetration points as described above, wherein the basic information table in step 1.2 includes the borehole number (holeID) and the hole mouth coordinates (X0, Y0, Z0); the inclination measurement record table includes the borehole number (holeID), the inclination measurement point number (N), the hole depth (h), the azimuth (β), and the inclination (α) of the measurement point; the ore section information table includes the borehole number (holeID), the ore layer position (H a ~H b ), ore body number (M).

[0013] The method for batch calculating horizontal projections of uranium ore drilling penetration points as described above, wherein the specific steps of step 2 are as follows:

[0014] Step 2.1: Call the borehole inclinometer record and use the full-angle half-distance method to calculate the depth H of each control point. i and control length L i , and store the calculation results in the inclinometer calculation result table. The specific calculation formula is as follows:

[0015] Control point depth calculation formula:

[0016] Control length calculation formula:

[0017] i is the serial number of the i-th measuring point; h i is the depth of the i-th measuring point; H i is the depth of the i-th control point; L i

[0018] The i-th control point controls the length.

[0019] Step 2.2, determine the direction of the coordinate system, X is the geodetic coordinate X direction, Y is the geodetic coordinate Y direction, Z is the elevation direction, call the borehole inclination measurement record, calculate the displacement in the three directions of each control section, and store the calculation results in the inclination measurement calculation result table. The specific calculation formula is as follows:

[0020] Control segment X-direction displacement: △X i =L i cosα i cosβ i

[0021] Control segment displacement in Y direction: △Y i =L i cosα i sinβ i

[0022] Control section Z direction displacement: △Z i =-L i sinα i

[0023] i is the serial number of the i-th measuring point; α i is the inclination angle of the i-th measuring point; β i is the azimuth of the i-th measuring point; L i

[0024] The i-th control point controls the length;

[0025] △X i is the displacement in the X direction of the i-th control segment; △Y i is the displacement in the Y direction of the i-th control segment;

[0026] △Z i is the displacement in the Z direction of the ith control segment.

[0027] Step 2.3, call the three-direction displacements calculated in step 2.2 in the inclinometer calculation result table, calculate the three-direction cumulative displacements, and store the calculation results in the inclinometer calculation result table. The specific calculation formula is as follows:

[0028] Cumulative displacement ∑ΔX in X direction i =△X1+△X2+△X3+…+△X i ;

[0029] Cumulative displacement in Y direction ∑ΔY i =△Y1+△Y2+△Y3+…+△Y i ;

[0030] Cumulative displacement ∑ΔZ in Z direction i =△Z1+△Z2+△Z3+…+△Z i ;

[0031] i is the number of the i-th measuring point; △X i is the displacement in the X direction of the i-th control segment; △Y i is the displacement in the Y direction of the i-th control segment;

[0032] △Z i is the displacement in the Z direction of the ith control segment.

[0033] Step 2.4, calculate the coordinates of each control point based on the hole coordinates (X0, Y0, Z0) and the cumulative displacement in three directions calculated in step 2.3 in the inclinometric calculation result table. i (X i , Y i , Z i ), and store the calculation results in the inclinometer calculation result table. The calculation formula is as follows:

[0034] X i =X0+∑ΔX i ;

[0035] Y i =Y0+∑ΔY i ;

[0036] Z i =Z0+∑ΔZ i ;

[0037] i is the serial number of the i-th measuring point; X i is the X coordinate of the i-th control point; Y i is the Y coordinate of the i-th control point; Z i is the Z coordinate of the i-th control point;

[0038] ∑△X i Cumulative displacement of the i-th control point in the X direction; ∑△Y i Cumulative displacement of the i-th control point in Y direction; ∑△H i

[0039] Cumulative displacement of the i-th control point in the Z direction.

[0040] The method for batch calculating horizontal projections of uranium ore drilling penetration points as described above, wherein the specific steps of step 3 are as follows:

[0041] Step 3.1: call the drilling single project ore-seeing parameters in the ore section information table, and calculate the hole depth P of each ore body in the drilling hole. M The results are stored in the drilling penetration point calculation table. The specific calculation formula is as follows:

[0042] P M =0.5(H a +H b ) or H b

[0043] Where: P M =0.5(H a +H b ) indicates that the penetration point of the borehole No. M is the intersection point of the central curved surface of the ore body and the depth of the borehole.

[0044] P M =H bThe point at which the borehole No. M cuts through the ore body is the depth of the intersection of the curved bottom plate of the ore body and the borehole.

[0045] Step 3.2, call the hole depth and coordinates of each control point calculated in step 2 in the inclinometric calculation result table, and determine the upper and lower control point numbers of the drill hole penetration point. The determination method is as follows:

[0046] Let C = min | P M -H j |;

[0047] If P M -H j <0, j is the number of the control point below the ore-seeking point, and j-1 is the number of the control point above the ore-seeking point;

[0048] If P M -H j >0, then j is the sequence number of the upper control point of the ore-seeking and cutting point, and j+1 is the sequence number of the lower control point of the ore-seeking and cutting point.

[0049] In the judgment formula: P M H is the hole depth of the drilling hole M ore body; j is the depth of the jth control point; j is the sequence number of the upper or lower control point of the cutting point.

[0050] Step 3.3, according to the upper and lower control point serial numbers of the mine cutting point calculated in step 3.2, call the horizontal coordinates and hole depths of the upper and lower control points in the inclinometric calculation result table, and store the results in the drilling hole mine cutting point calculation table, storing the record Record = {holeID, M, P M , X 上 , Y 上 , H 上 , X 下 , Y 下 , H 下};

[0051] In the storage record, holeID is the drilling hole number, M is the ore body number, P is M The depth of the hole at the cutting point of the Mth ore body, X 上 , Y 上 is the horizontal coordinate of the control point on the cutting point, X 下 , Y 下 H is the horizontal coordinate of the control point under the cutting point, 上 H is the hole depth of the control point on the cutting point; 下 The hole depth of the control point below the cutting point.

[0052] The method for batch calculating horizontal projections of uranium ore drilling penetration points as described above, wherein the specific steps of step 4 are as follows:

[0053] Step 4.1, call the stored record in the drilling hole penetration point calculation table in step 3.3, and calculate the horizontal projection coordinates (X M , Y M ) and store the results in the drill hole penetration point calculation table.

[0054] The calculation formula for the horizontal projection coordinates of the drilling hole penetration point is as follows:

[0055] X M =(X 下 -X 上 )*(P M -H 上 ) / (H 下 -H 上 )+X 上

[0056] Y M =(Y 下 -Y 上 )*(P M -H 上 ) / (H 下 -H 上 )+Y 上

[0057] X M , Y M Horizontal projection coordinates of the cutting point; X 上 , Y 上 X is the horizontal coordinate of the control point on the cutting point; 下 , Y 下 P is the horizontal position of the control point below the cutting point; M The hole depth of the drilling hole M ore body;

[0058] H 上 H is the hole depth of the control point on the cutting point; 下 The hole depth of the control point below the cutting point.

[0059] Step 4.2, batch output the horizontal projection coordinates of the drill hole penetration points calculated in step 4.2, output Output = {holeID, M, X M , Y M}.

[0060] As described above, a method for batch calculating horizontal projections of uranium ore drilling penetration points is provided, wherein the database in step 1 is a uranium ore drilling data table and data structure established based on an Access database.

[0061] The remarkable effects of the present invention are as follows: (1) The present invention provides a method for batch calculating the horizontal projections of uranium ore drill hole penetration points, which can automatically calculate the coordinates of the drill hole penetration points, reduce the errors in the process of manual participation in data processing, increase the calculation accuracy, and greatly improve the drawing efficiency of the reserve plan.

[0062] (2) The present invention provides a method for batch calculation of horizontal projections of uranium ore drilling penetration points, which breaks through the traditional calculation process of single hole calculation one by one, realizes batch calculation of drilling hole penetration point coordinates and output, and saves a lot of calculation time.

[0063] (3) The present invention provides a method for batch calculation of horizontal projections of uranium ore drilling penetration points, wherein intermediate calculation data and result data are stored in a database, and the database is stored in a table form, which is convenient for rapid query and export of data. BRIEF DESCRIPTION OF THE DRAWINGS

[0064] Figure 1 Flow chart of horizontal projection calculation of drilling hole penetration point

[0065] Figure 2 Access Database Structure

[0066] Figure 3 Example diagram of the spatial relationship of drilling point projection DETAILED DESCRIPTION

[0067] The present invention is further described in detail below with reference to the accompanying drawings and single hole examples.

[0068] The technical solution of the present invention is as follows: A method for batch calculating the horizontal projection of the uranium ore drilling point, the method specifically comprises the following steps: Figure 1 ):

[0069] Step 1, prepare drilling database: establish database structure and input drilling data, the database includes drilling input information structure table and drilling calculation information structure table;

[0070] The specific steps of step 1 are as follows:

[0071] Step 1.1, based on the Access database structure, design the data tables for different types of drilling information ( Figure 2 ), the input information structure table includes the basic information table (DBJBXX), the inclinometer record table (DBCXJL), and the mining section information table (DBKDXX); the calculation information structure table includes the inclinometer calculation result table (DBCXJG) and the drilling hole penetration point calculation table (DBQCD).

[0072] Step 1.2, the drilling data is entered into the corresponding data table in the database for storage. The basic information table (DBJBXX) includes the drilling number (holeID) and the coordinates of the hole mouth (X0, Y0, Z0); the inclination record table (DBCXJL) includes the drilling number (holeID), the inclination point number (N), the hole depth (h), the azimuth (β), and the inclination (α); the mining section information table (DBKDXX) includes the drilling number (holeID), the ore layer position (H a ~H b ), ore body number (M), etc.

[0073] Step 2, calling the borehole inclinometer record and calculating the coordinates of each control point: calling the borehole inclinometer record, using the full-angle half-distance method to calculate the relevant parameters of each control point, and calculating the coordinates of each control point according to the parameters of each control point, and storing the calculation results;

[0074] The specific steps of step 2 are as follows:

[0075] Step 2.1, call the borehole inclinometer record table (DBCXJL), use the full angle half distance method to calculate the depth H of each control point i and control length L i The calculation results are stored in the inclination calculation result table (DBCXJG). The specific calculation formula is as follows:

[0076] Control point depth calculation formula:

[0077] Control length calculation formula:

[0078] i is the serial number of the i-th measuring point; h i is the depth of the i-th measuring point; H i is the depth of the i-th control point;

[0079] L i The i-th control point controls the length.

[0080] Step 2.2, determine the direction of the coordinate system, X is the geodetic X direction, Y is the geodetic Y direction, and Z is the elevation direction ( Figure 3 ), call the drilling inclination record table (DBCXJL), calculate the displacement of each control section in three directions, and store the calculation results in the inclination calculation result table (DBCXJG). The specific calculation formula is as follows:

[0081] Control segment X-direction displacement: △X i =L i cosα i cosβ i

[0082] Control segment displacement in Y direction: △Y i=L i cosα i sinβ i

[0083] Control section Z direction displacement: △Z i =-L i sinα i

[0084] i is the serial number of the i-th measuring point; α i is the inclination angle of the i-th measuring point; β i is the azimuth of the i-th measuring point; L i

[0085] The i-th control point controls the length;

[0086] △X i is the displacement in the X direction of the i-th control segment; △Y i is the displacement in the Y direction of the i-th control segment;

[0087] △Z i is the displacement in the Z direction of the ith control segment.

[0088] Step 2.3, call the three-direction displacements calculated in step 2.2 in the inclinometer calculation result table (DBCXJG), calculate the three-direction cumulative displacements, and store the calculation results in the inclinometer calculation result table (DBCXJG). The specific calculation formula is as follows:

[0089] Cumulative displacement ∑ΔX in X direction i =△X1+△X2+△X3+…+△X i ;

[0090] Cumulative displacement in Y direction ∑ΔY i =△Y1+△Y2+△Y3+…+△Y i ;

[0091] Cumulative displacement ∑ΔZ in Z direction i =△Z1+△Z2+△Z3+…+△Z i ;

[0092] i is the number of the i-th measuring point; △X i is the displacement in the X direction of the i-th control segment; △Y i For the i

[0093] Control segment displacement in Y direction;

[0094] △Z i is the displacement in the Z direction of the ith control segment.

[0095] Step 2.4, calculate the coordinates of each control point based on the hole coordinates (X0, Y0, Z0) and the three-direction cumulative displacement calculated in step 2.3 in the inclinometric calculation result table (DBCXJG). i (X i , Y i , Z i ), and store the calculation results in the inclinometer calculation result table (DBCXJG). The calculation formula is as follows:

[0096] X i =X0+∑ΔX i ;

[0097] Y i =Y0+∑ΔY i ;

[0098] Z i =Z0+∑ΔZ i ;

[0099] i is the serial number of the i-th measuring point; X i is the X coordinate of the i-th control point; Y i is the Y coordinate of the i-th control point; Z i is the Z coordinate of the i-th control point;

[0100] ∑△X i Cumulative displacement of the i-th control point in the X direction; ∑△Y i Cumulative displacement of the i-th control point in Y direction;

[0101] ∑△H i Cumulative displacement of the i-th control point in the Z direction.

[0102] Step 3, call the drilling ore-seeing parameters and determine the ore-seeing and cutting point control points: call the drilling single project ore-seeing parameters, determine the upper and lower control numbers of the ore body, and store the results;

[0103] The specific steps of step 3 are as follows:

[0104] Step 3.1, call the drilling single project ore-seeing parameters in the ore section information table (DBKDXX), and calculate the hole depth P of each ore body in the drilling hole. M The results are stored in the Drilling Qualification and Cutting Point Calculation Table (DBQCD). The specific calculation formula is as follows:

[0105] P M =0.5(H a +H b ) or H b

[0106] Where: P M =0.5(H a +H b) indicates that the penetration point of the borehole No. M is the intersection point of the central curved surface of the ore body and the depth of the borehole.

[0107] P M =H b The point at which the borehole No. M cuts through the ore body is the depth of the intersection of the curved bottom plate of the ore body and the borehole.

[0108] Step 3.2, call the hole depth and coordinates of each control point calculated in step 2 in the inclinometer calculation result table (DBCXJG), determine the upper and lower control point numbers of the drill hole penetration point, and the determination method is as follows:

[0109] Let C = min | P M -H j |;

[0110] If P M -H j <0, j is the number of the control point below the ore-seeking point, and j-1 is the number of the control point above the ore-seeking point;

[0111] If P M -H j >0, then j is the sequence number of the upper control point of the ore-seeking and cutting point, and j+1 is the sequence number of the lower control point of the ore-seeking and cutting point.

[0112] In the judgment formula: P M H is the hole depth of the drilling hole M ore body; j is the depth of the jth control point; j is the sequence number of the upper or lower control point of the cutting point.

[0113] Step 3.3, according to the upper and lower control point serial numbers of the mine cutting point calculated in step 3.2, call the horizontal coordinates and hole depths of the upper and lower control points in the inclinometric calculation result table (DBCXJG), and store the results in the drilling hole mine cutting point calculation table (DBQCD), and store the record Record = {holeID, M, P M , X 上 , Y 上 , H 上 , X 下 , Y 下 , H 下};

[0114] In the storage record, holeID is the drilling hole number, M is the ore body number, P is M The depth of the hole at the cutting point of the Mth ore body, X 上 , Y 上 is the horizontal coordinate of the control point on the cutting point, X 下 , Y 下 H is the horizontal coordinate of the control point under the cutting point, 上 H is the hole depth of the control point on the cutting point; 下 The hole depth of the control point below the cutting point.

[0115] Step 4, calculate the coordinates of the drill hole's ore-cutting point and output the coordinates: calculate the coordinates of the drill hole's ore-cutting point based on the upper and lower control points of the ore-cutting point, and output the projection of the cutting point.

[0116] The specific steps of step 4 are as follows:

[0117] Step 4.1, call the stored record in the drill hole penetration point calculation table (DBQCD) in step 3.3, and calculate the horizontal projection coordinates (X M , Y M ) and store the results in the Drill Hole Qualification and Cutting Point Calculation Table (DBQCD).

[0118] The calculation formula for the horizontal projection coordinates of the drilling hole penetration point is as follows:

[0119] X M =(X 下 -X 上 )*(P M -H 上 ) / (H 下 -H 上 )+X 上

[0120] Y M =(Y 下 -Y 上 )*(P M -H 上 ) / (H 下 -H 上 )+Y 上

[0121] X M , Y M Horizontal projection coordinates of the cutting point; X 上 , Y 上 X is the horizontal coordinate of the control point on the cutting point; 下 , Y 下 P is the horizontal position of the control point below the cutting point; M H is the hole depth of the drilling hole M ore body; 上 H is the hole depth of the control point on the cutting point; 下 The hole depth of the control point below the cutting point.

[0122] Step 4.2, batch output the horizontal projection coordinates of the drill hole penetration points calculated in step 4.2, output Output = {holeID, M, X M , Y M}.

[0123] The present invention is described in detail above with reference to the accompanying drawings and embodiments, but the present invention is not limited to the above embodiments, and various changes can be made within the knowledge of ordinary technicians in the field without departing from the purpose of the present invention. The contents not described in detail in the present invention can adopt the existing technology.

Claims

1. A method for batch calculation of horizontal projections of uranium ore drilling penetration points, characterized in that: The method comprises the following steps: Step 1, prepare drilling database: establish database structure and input drilling data, the database includes drilling input information structure table and drilling calculation information structure table; Step 2, calling the borehole inclinometer record and calculating the coordinates of each control point: calling the borehole inclinometer record, using the full-angle half-distance method to calculate the relevant parameters of each control point, and calculating the coordinates of each control point according to the parameters of each control point, and storing the calculation results; Step 3, call the drilling ore-seeing parameters to determine the upper and lower control points of the ore-seeing and cutting point: call the drilling single project ore-seeing parameters to determine the upper and lower control numbers of the ore body, and store the results; Step 4, calculate the coordinates of the drill hole's ore-penetrating point and output the coordinates: calculate the coordinates of the drill hole's ore-penetrating point based on the upper and lower control points of the ore-penetrating point, and output the projection of the penetration point. The specific steps of step 1 are as follows: Step 1.1, based on the Access database structure, design different types of drilling information data tables, and enter the information structure table including basic information table, inclination record table, and mining section information table; The calculation information structure table includes the inclination measurement calculation result table and the drilling hole penetration point calculation table. Step 1.2, enter the drilling data into the corresponding data table in the database for storage, which is the basic information table. The basic information table in step 1.2 includes the borehole number holeID and the hole mouth coordinates (X0, Y0, Z0); the inclination record table includes the borehole number holeID, the inclination point number N, the hole depth h, the azimuth β, and the inclination α; the ore section information table includes the borehole number holeID, the ore layer position H a and H b , ore body number M, The specific steps of step 2 are as follows: Step 2.1: Call the borehole inclinometer record and use the full-angle half-distance method to calculate the depth H of each control point. i and control length L i , and store the calculation results in the inclinometer calculation result table. The specific calculation formula is as follows: Control point depth calculation formula: Control length calculation formula: i is the serial number of the i-th measuring point; h i is the depth of the i-th measuring point; H i is the depth of the i-th control point; L i is the control length of the i control point, Step 2.2, determine the direction of the coordinate system, X is the geodetic coordinate X direction, Y is the geodetic coordinate Y direction, Z is the elevation direction, call the borehole inclination measurement record, calculate the displacement in the three directions of each control section, and store the calculation results in the inclination measurement calculation result table. The specific calculation formula is as follows: Control segment X-direction displacement: △X i =L i cosα i cosβ i Control segment displacement in Y direction: △Y i =L i cosα i sinβ i Control section Z direction displacement: △Z i =-L i sinα i α i is the inclination angle of the i-th measuring point; β i is the azimuth of the i-th measuring point; △X i is the displacement in the X direction of the i-th control segment; △Y i is the displacement in the Y direction of the i-th control segment; △Z i is the displacement in the Z direction of the ith control segment, Step 2.3, call the three-direction displacements calculated in step 2.2 in the inclinometer calculation result table, calculate the three-direction cumulative displacements, and store the calculation results in the inclinometer calculation result table. The specific calculation formula is as follows: Cumulative displacement ∑ΔX in X direction i =△X1+△X2+△X3+…+△X i ; Cumulative displacement in Y direction ∑ΔY i =△Y1+△Y2+△Y3+…+△Y i ; Cumulative displacement ∑ΔZ in Z direction i =△Z1+△Z2+△Z3+…+△Z i ; Step 2.4, calculate the coordinates of each control point based on the hole coordinates (X0, Y0, Z0) and the cumulative displacement in three directions calculated in step 2.3 in the inclinometric calculation result table. i (X i , Y i , Z i ), and store the calculation results in the inclinometer calculation result table. The calculation formula is as follows: X i =X0+∑ΔX i ; Y i =Y0+∑ΔY i ; WITH i =Z0+∑ΔZ i ; X i is the X coordinate of the i-th control point; Y i is the Y coordinate of the i-th control point; Z i is the Z coordinate of the i-th control point; ∑△X i is the cumulative displacement of the i-th control point in the X direction; ∑△Y i is the cumulative displacement of the i-th control point in the Y direction; ∑ΔZ i is the cumulative displacement of the i-th control point in the Z direction, The specific steps of step 3 are as follows: Step 3.1: call the drilling single project ore-seeing parameters in the ore section information table, and calculate the hole depth P of each ore body in the drilling hole. M The results are stored in the drilling penetration point calculation table. The specific calculation formula is as follows: P M =0.5(H a +H b ) or H b Where: P M =0.5(H a +H b ) indicates that the penetration point of the borehole No. M is the intersection point of the central curved surface of the ore body and the depth of the borehole. P M =H b Indicates that the penetration point of the borehole No. M is the intersection of the curved bottom plate of the ore body and the depth of the borehole. Step 3.2, call the hole depth and coordinates of each control point calculated in step 2 in the inclinometric calculation result table, and determine the upper and lower control point numbers of the drill hole penetration point. The determination method is as follows: Let C = min | P M -H j |; If P M -H j <0, j is the number of the control point below the ore-seeking point, and j-1 is the number of the control point above the ore-seeking point; If P M -H j >0, then j is the number of the upper control point of the ore-seeking and cutting point, and j+1 is the number of the lower control point of the ore-seeking and cutting point. In the judgment formula: H j is the depth of the jth control point; j is the number of the upper or lower control point of the cutting point, Step 3.3, according to the upper and lower control point serial numbers of the mine cutting point calculated in step 3.2, call the horizontal coordinates and hole depths of the upper and lower control points in the inclinometric calculation result table, and store the results in the drilling hole mine cutting point calculation table, storing the record Record = {holeID, M, P M , X 上 , Y 上 , H 上 , X 下 , Y 下 , H 下 }; Storage record, X 上 , Y 上 is the horizontal coordinate of the control point on the cutting point, X 下 , Y 下 H is the horizontal coordinate of the control point under the cutting point, 上 H is the hole depth of the control point on the cutting point; 下 The hole depth of the control point under the cutting point, The specific steps of step 4 are as follows: Step 4.1, call the stored record in the drilling hole penetration point calculation table in step 3.3, and calculate the horizontal projection coordinates (X M , Y M ), and store the results in the drill hole penetration point calculation table. The calculation formula for the horizontal projection coordinates of the drilling hole penetration point is as follows: X M =(X 下 -X 上 )*(P M -H 上 ) / (H 下 -H 上 )+X 上 Y M =(Y 下 -Y 上 )*(P M -H 上 ) / (H 下 -H 上 )+Y 上 Step 4.2, batch output the horizontal projection coordinates of the drill hole penetration points calculated in step 4.1, output Output = {holeID, M, X M , Y M }.

Images