Uranium mine drilling see-through point horizontal projection batch calculation method
Technical Field
The invention belongs to a calculation method, and particularly relates to a uranium ore drilling cut-through point horizontal projection batch calculation method.
Background
The drilling ore-finding cutting-through point is the basis for drawing the reserve plane drawing, and is not only used for pushing the plane of the ore body during reserve calculation, but also used for guiding drilling construction in time. In actual reservoir plan mapping, the borehole see cut-through point is typically offset from the location envisaged by the borehole due to the inclination and azimuth offset that exists in the drilling operation. In order to re-reflect the actual mineral site on the projection map, it is often necessary to calculate the position of the drill mineral cutting point by projection. The traditional method for calculating the projection of the drilling holes is to calculate the cutting-through points of the drilling holes from the drilling holes one by one according to given drilling coordinates, drilling inclinations and mineral positions. The method consumes a large amount of manpower and material resources, if the data volume is large, a long time is required to calculate all the drilling see-through points, meanwhile, the whole process of calculation is required to be manually participated, and if the manual operation is wrong, the position of the calculated cutting-through points is inaccurate.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a uranium mine drilling see-through point horizontal projection batch calculation method.
The invention discloses a uranium mine drilling see-through point horizontal projection batch calculation method, which comprises the following steps:
step 1, preparing a drilling database, namely establishing a database structure and inputting drilling data, wherein the database comprises a drilling input information structure table and a drilling calculation information structure table;
Step 2, calling a drilling inclinometry record, calculating the coordinates of each control point, namely calling the drilling inclinometry record, calculating the related parameters of each control point by adopting a full-angle half-distance method, calculating the coordinates of each control point according to the parameters of each control point, and storing the calculation result;
Step3, calling drilling ore-finding parameters, determining upper and lower control points of an ore-finding cutting point, namely calling drilling single engineering ore-finding parameters, determining upper and lower control serial numbers of ore bodies, and storing results;
And 4, calculating coordinates of the cutting-through points of the drilling visible ores, outputting the coordinates, namely calculating the coordinates of the cutting-through points of the drilling visible ores according to the upper and lower control points of the cutting-through points of the visible ores, and outputting the projection of the cutting-through points.
The uranium mine drilling see-through point horizontal projection batch calculation method comprises the following specific steps of:
Step 1.1, designing different types of drilling information data tables based on an Access database structure, and inputting an information structure table comprising a basic information table, an inclinometry record table and a mine section information table, wherein the calculation information structure table comprises an inclinometry calculation result table and a drilling see-through point calculation table.
And 1.2, inputting drilling data into a corresponding data table in a database for storage, wherein the data table is a basic information table.
The uranium mine drilling and cutting point horizontal projection batch calculation method comprises the steps that a basic information table in the step 1.2 comprises a drilling number (holeID) and an orifice coordinate (X 0,Y0,Z0), an inclinometry record table comprises a drilling number (holeID), an inclinometry point number (N), a measuring point hole depth (H), an azimuth angle (beta) and an inclination angle (alpha), and a mine section information table comprises a drilling number (holeID), a mine layer position (H a~Hb) and a mine body number (M).
The uranium mine drilling see-through point horizontal projection batch calculation method comprises the following specific steps of:
step 2.1, calling a borehole inclinometry record, calculating the depth H i and the control length L i of each control point by adopting a full-angle half-distance method, and storing calculation results into an inclinometry calculation result table, wherein a specific calculation formula is as follows:
the control point depth calculation formula:
Control length calculation formula:
i is the serial number of the ith measuring point, H i is the depth of the ith measuring point, H i is the depth of the ith control point, L i
The ith control point controls the length.
Step 2.2, determining the coordinate system direction, wherein X is the geodetic coordinate X direction, Y is the geodetic coordinate Y direction, Z is the elevation direction, calling the drilling inclinometry record, calculating the displacement of each control section in three directions, and storing the calculation result into an inclinometry calculation result table, wherein the specific calculation formula is as follows:
The displacement of the control section in the X direction is delta X i=Licosαicosβi
The displacement of the control section in the Y direction is delta Y i=Licosαisinβi
The displacement of the control section in the Z direction is delta Z i=-Lisinαi
I is the serial number of the ith measuring point, alpha i is the dip angle of the ith measuring point, beta i is the azimuth angle of the ith measuring point, L i
The ith control point controls the length;
Delta X i is the displacement in the X direction of the ith control section, delta Y i is the displacement in the Y direction of the ith control section;
DeltaZ i is the displacement of the ith control section in the Z direction.
Step 2.3, calling the three-direction displacement calculated in step 2.2 in the inclinometer calculation result table, calculating the three-direction accumulated displacement, and storing the calculation result in the inclinometer calculation result table, wherein the specific calculation formula is as follows:
cumulative displacement ΣΔx i=△X1+△X2+△X3+…+△Xi in the X direction;
accumulated displacement ΣΔy i=△Y1+△Y2+△Y3+…+△Yi in the Y direction;
Accumulated displacement Sigma delta Z i=△Z1+△Z2+△Z3+…+△Zi in Z direction;
i is the serial number of the ith measuring point, deltaX i is the displacement of the ith control section in the X direction, deltaY i is the displacement of the ith control section in the Y direction;
DeltaZ i is the displacement of the ith control section in the Z direction.
Step 2.4, calculating the coordinates O i(Xi,Yi,Zi of each control point according to the orifice coordinates (X 0,Y0,Z0) and the three-direction accumulated displacement calculated in the step 2.3 in the inclinometry calculation result table, and storing the calculation result in the inclinometry calculation result table, wherein the calculation formula is as follows:
Xi=X0+∑ΔXi;
Yi=Y0+∑ΔYi;
Zi=Z0+∑ΔZi;
i is the serial number of the ith measuring point, X i is the X coordinate of the ith control point, Y i is the Y coordinate of the ith control point, and Z i is the Z coordinate of the ith control point;
Sigma DeltaX i th control point X is cumulatively displaced, sigma DeltaY i th control point Y is cumulatively displaced, sigma DeltaH i
The ith control point accumulates displacement in the Z-direction.
The uranium mine drilling see-through point horizontal projection batch calculation method comprises the following specific steps of:
step 3.1, calling single engineering mineral seeing parameters in a mineral segment information table, calculating the hole depth P M of mineral seeing and cutting through points of each mineral body of the drill, and storing the result in a mineral seeing and cutting through point calculation table of the drill, wherein the specific calculation formula is as follows:
p M=0.5(Ha+Hb) or H b
Wherein P M=0.5(Ha+Hb) represents that the cutting-through point of the Mth numbered ore body of the drill hole is the intersection point hole depth of the central curved surface of the ore body and the drill hole,
The cutting and passing point of the M numbered ore body of the P M=Hb drill hole is intersection point hole depth of curved surface bottom plate of ore body and drilling hole.
And 3.2, calling the hole depths and coordinates of the control points calculated in the step 2 in an inclinometer calculation result table, and determining the sequence numbers of the upper control point and the lower control point of the drilling see-through point, wherein the determination method is as follows:
let c=min| P M-Hj |;
If P M-Hj is less than 0, j is the control point number under the cut-through point of the see-through mine, and j-1 is the control point number on the cut-through point of the see-through mine;
If P M-Hj is greater than 0, j is the control point number on the cut-through point of the see-through mine, and j+1 is the control point number under the cut-through point of the see-through mine.
In the judgment, P M is the depth of a cutting point of a drilling M-shaped ore body, H j is the depth of a j-th control point, and j is the sequence number of an upper control point or a lower control point of the cutting point.
Step 3.3, according to the sequence numbers of the upper control point and the lower control point of the ore cutting point calculated in the step 3.2, the horizontal coordinates and the hole depth of the upper control point and the lower control point in an inclinometer calculation result table are called, the result is stored in a drilling ore cutting point calculation table, and Record = { holeID, M and P M,X Upper part ,Y Upper part ,H Upper part ,X Lower part(s) ,Y Lower part(s) ,H Lower part(s) };
In the storage record, holeID is a drilling number, M is a mineral body number, P M is the depth of a cutting-through point of the mineral body with M number, X Upper part 、Y Upper part is the horizontal coordinate of a control point on the cutting-through point, X Lower part(s) ,Y Lower part(s) is the horizontal coordinate of a control point below the cutting-through point, H Upper part is the depth of a control point on the cutting-through point, and H Lower part(s) is the depth of a control point below the cutting-through point.
The uranium mine drilling see-through point horizontal projection batch calculation method comprises the following specific steps of:
And 4.1, calling a storage record in the drilling see-through point calculation table in the step 3.3, calculating horizontal projection coordinates (X M,YM) of the drilling see-through point, and storing the result into the drilling see-through point calculation table.
The calculation formula of horizontal projection coordinates of the drilling see-through point is as follows:
XM=(X Lower part(s) -X Upper part )*(PM-H Upper part )/(H Lower part(s) -H Upper part )+X Upper part
YM=(Y Lower part(s) -Y Upper part )*(PM-H Upper part )/(H Lower part(s) -H Upper part )+Y Upper part
X M、YM is a horizontal projection coordinate of a cutting-through point, X Upper part 、Y Upper part is a horizontal coordinate of an upper control point of the cutting-through point, X Lower part(s) 、Y Lower part(s) is a horizontal seat of a lower control point of the cutting-through point, and P M is the depth of a cutting-through point of a drilling M-shaped ore body;
H Upper part is the depth of the control point hole above the cut-through point, and H Lower part(s) is the depth of the control point hole below the cut-through point.
And 4.2, outputting the horizontal projection coordinates of the drilling see-through points calculated in the step 4.2 in batches, and outputting output= { holeID, M, X M,YM }.
According to the uranium deposit drilling see-through point horizontal projection batch calculation method, the database in the step 1 is a uranium deposit drilling data table and a data structure which are established based on the Access database.
The invention has the remarkable effects that (1) the invention provides the uranium mine drilling see-through point horizontal projection batch calculation method, which can automatically calculate the coordinates of drilling cut-through points, reduce errors in the process of manually participating in data processing, increase calculation accuracy and greatly improve the drawing efficiency of a reserve plan.
(2) The invention provides a horizontal projection batch calculation method for uranium mine drilling see-through points, which breaks through single-hole one-by-one calculation in the traditional calculation process, realizes batch calculation of drilling see-through point coordinates and output, and saves a large amount of calculation time.
(3) The invention provides a uranium mine drilling and cutting point horizontal projection batch calculation method, wherein intermediate calculation data and result data are stored in a database, and the database is stored in a form of a table, so that the data can be conveniently and rapidly inquired and exported.
Drawings
FIG. 1 horizontal projection calculation flow chart of drilling see-through point
FIG. 2Access database structure
FIG. 3 is an exemplary diagram of a projected spatial relationship of borehole survey points
Detailed Description
The invention is described in further detail below with reference to the drawings and the single-hole examples.
The technical scheme of the invention is that the uranium mine drilling see-through point horizontal projection batch calculation method specifically comprises the following steps (figure 1):
step 1, preparing a drilling database, namely establishing a database structure and inputting drilling data, wherein the database comprises a drilling input information structure table and a drilling calculation information structure table;
The specific steps of the step 1 are as follows:
Step 1.1, designing a drilling different-category information data table (figure 2) based on an Access database structure, wherein the input information structure table comprises a basic information table (DBJBXX), an inclinometry record table (DBCXJL) and an ore section information table (DBKDXX), and the calculation information structure table comprises an inclinometry calculation result table (DBCXJG) and a drilling see-through point calculation table (DBQCD).
Step 1.2, inputting drilling data into a corresponding data table in a database for storage, wherein a basic information table (DBJBXX) comprises a drilling number (holeID) and an orifice coordinate (X 0,Y0,Z0), an inclinometry record table (DBCXJL) comprises a drilling number (holeID), an inclinometry point number (N), a measuring point hole depth (H), an azimuth angle (beta) and an inclination angle (alpha), and a mine section information table (DBKDXX) comprises a drilling number (holeID), a mine layer position (H a~Hb), a mine body number (M) and the like.
Step 2, calling a drilling inclinometry record, calculating the coordinates of each control point, namely calling the drilling inclinometry record, calculating the related parameters of each control point by adopting a full-angle half-distance method, calculating the coordinates of each control point according to the parameters of each control point, and storing the calculation result;
The specific steps of the step 2 are as follows:
Step 2.1, calling a drilling inclinometry record table (DBCXJL), calculating the depth H i and the control length L i of each control point by adopting a full-angle half-distance method, and storing calculation results into an inclinometry calculation result table (DBCXJG), wherein the specific calculation formula is as follows:
the control point depth calculation formula:
Control length calculation formula:
i is the serial number of the ith measuring point, H i is the depth of the ith measuring point, and H i is the depth of the ith control point;
l i ith control point controls length.
Step 2.2, determining the coordinate system direction, wherein X is the geodetic coordinate X direction, Y is the geodetic coordinate Y direction, Z is the elevation direction (figure 3), calling a drilling inclinometry record table (DBCXJL), calculating the displacement of each control section in three directions, and storing the calculation result into an inclinometry calculation result table (DBCXJG), wherein the specific calculation formula is as follows:
The displacement of the control section in the X direction is delta X i=Licosαicosβi
The displacement of the control section in the Y direction is delta Y i=Licosαisinβi
The displacement of the control section in the Z direction is delta Z i=-Lisinαi
I is the serial number of the ith measuring point, alpha i is the dip angle of the ith measuring point, beta i is the azimuth angle of the ith measuring point, L i
The ith control point controls the length;
Delta X i is the displacement in the X direction of the ith control section, delta Y i is the displacement in the Y direction of the ith control section;
DeltaZ i is the displacement of the ith control section in the Z direction.
Step 2.3, calling the three-direction displacement calculated in step 2.2 in the inclinometer calculation result table (DBCXJG), calculating the three-direction accumulated displacement, and storing the calculation result in the inclinometer calculation result table (DBCXJG), wherein the specific calculation formula is as follows:
cumulative displacement ΣΔx i=△X1+△X2+△X3+…+△Xi in the X direction;
accumulated displacement ΣΔy i=△Y1+△Y2+△Y3+…+△Yi in the Y direction;
Accumulated displacement Sigma delta Z i=△Z1+△Z2+△Z3+…+△Zi in Z direction;
i is the serial number of the ith measuring point, deltaX i is the displacement of the ith control section in the X direction, deltaY i is the ith
Controlling the Y-direction displacement of the section;
DeltaZ i is the displacement of the ith control section in the Z direction.
Step 2.4, calculating the coordinates O i(Xi,Yi,Zi of each control point according to the three-direction accumulated displacement calculated in step 2.3 in the orifice coordinates (X 0,Y0,Z0) and the inclinometry calculation result table (DBCXJG), and storing the calculation result in the inclinometry calculation result table (DBCXJG), wherein the calculation formula is as follows:
Xi=X0+∑ΔXi;
Yi=Y0+∑ΔYi;
Zi=Z0+∑ΔZi;
i is the serial number of the ith measuring point, X i is the X coordinate of the ith control point, Y i is the Y coordinate of the ith control point, and Z i is the Z coordinate of the ith control point;
Sigma delta X i ith control point X cumulative displacement sigma delta Y i ith control point Y cumulative displacement;
Sigma delta H i ith control point Z-direction accumulated displacement.
Step 3, calling drilling ore-finding parameters, determining ore-finding cutting-through point control points, namely calling drilling single-project ore-finding parameters, determining ore body up-down control serial numbers, and storing results;
the specific steps of the step 3 are as follows:
Step 3.1, calling single-engineering mineral seeing parameters in a mineral segment information table (DBKDXX), calculating the hole depth P M of mineral seeing and cutting through points of each mineral body of the drill hole, and storing the result in a mineral seeing and cutting through point calculation table (DBQCD), wherein the specific calculation formula is as follows:
p M=0.5(Ha+Hb) or H b
Wherein P M=0.5(Ha+Hb) represents that the cutting-through point of the Mth numbered ore body of the drill hole is the intersection point hole depth of the central curved surface of the ore body and the drill hole,
The cutting and passing point of the M numbered ore body of the P M=Hb drill hole is intersection point hole depth of curved surface bottom plate of ore body and drilling hole.
And 3.2, calling the hole depths and coordinates of the control points calculated in the step 2 in a table of an inclinometer calculation result table (DBCXJG), and determining the sequence numbers of the upper control point and the lower control point of the drilling ore-finding cutting-through point, wherein the judging method is as follows:
let c=min| P M-Hj |;
If P M-Hj is less than 0, j is the control point number under the cut-through point of the see-through mine, and j-1 is the control point number on the cut-through point of the see-through mine;
If P M-Hj is greater than 0, j is the control point number on the cut-through point of the see-through mine, and j+1 is the control point number under the cut-through point of the see-through mine.
In the judgment, P M is the depth of a cutting point of a drilling M-shaped ore body, H j is the depth of a j-th control point, and j is the sequence number of an upper control point or a lower control point of the cutting point.
Step 3.3, according to the sequence numbers of the upper control point and the lower control point of the ore cutting point calculated in the step 3.2, calling the horizontal coordinates and the hole depth of the upper control point and the lower control point in an inclinometry calculation result table (DBCXJG), storing the result in a drilling ore cutting point calculation table (DBQCD), and storing records record= { holeID, M and P M,X Upper part ,Y Upper part ,H Upper part ,X Lower part(s) ,Y Lower part(s) ,H Lower part(s) };
In the storage record, holeID is a drilling number, M is a mineral body number, P M is the depth of a cutting-through point of the mineral body with M number, X Upper part 、Y Upper part is the horizontal coordinate of a control point on the cutting-through point, X Lower part(s) ,Y Lower part(s) is the horizontal coordinate of a control point below the cutting-through point, H Upper part is the depth of a control point on the cutting-through point, and H Lower part(s) is the depth of a control point below the cutting-through point.
And 4, calculating coordinates of the cutting-through points of the drilling visible ores, outputting the coordinates, namely calculating the coordinates of the cutting-through points of the drilling visible ores according to the upper and lower control points of the cutting-through points of the visible ores, and outputting the projection of the cutting-through points.
The specific steps of the step 4 are as follows:
and 4.1, calling a storage record in a drilling see-through point calculation table (DBQCD) in the step 3.3, calculating horizontal projection coordinates (X M,YM) of the drilling see-through point, and storing the result into the drilling see-through point calculation table (DBQCD).
The calculation formula of horizontal projection coordinates of the drilling see-through point is as follows:
XM=(X Lower part(s) -X Upper part )*(PM-H Upper part )/(H Lower part(s) -H Upper part )+X Upper part
YM=(Y Lower part(s) -Y Upper part )*(PM-H Upper part )/(H Lower part(s) -H Upper part )+Y Upper part
X M、YM is the horizontal projection coordinate of the cutting point, X Upper part 、Y Upper part is the horizontal coordinate of the upper control point of the cutting point, X Lower part(s) 、Y Lower part(s) is the horizontal seat of the lower control point of the cutting point, P M is the hole depth of the cutting point of the M-shaped ore body for drilling, H Upper part is the hole depth of the upper control point of the cutting point, and H Lower part(s) is the hole depth of the lower control point of the cutting point.
And 4.2, outputting the horizontal projection coordinates of the drilling see-through points calculated in the step 4.2 in batches, and outputting output= { holeID, M, X M,YM }.
The present invention has been described in detail with reference to the drawings and the embodiments, but the present invention is not limited to the embodiments described above, and various changes can be made within the knowledge of those skilled in the art without departing from the spirit of the present invention. The invention may be practiced otherwise than as specifically described.