[go: up one dir, main page]

CN114329301B - A batch calculation method for horizontal projection of uranium drilling penetration points - Google Patents

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

Info

Publication number
CN114329301B
CN114329301B CN202011047965.8A CN202011047965A CN114329301B CN 114329301 B CN114329301 B CN 114329301B CN 202011047965 A CN202011047965 A CN 202011047965A CN 114329301 B CN114329301 B CN 114329301B
Authority
CN
China
Prior art keywords
point
drilling
calculation
ore
control point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011047965.8A
Other languages
Chinese (zh)
Other versions
CN114329301A (en
Inventor
刘洋
张明林
祝洪涛
王良玉
王青
李继木
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Nuclear Inner Mongolia Mining Investment Co ltd
Original Assignee
China Nuclear Inner Mongolia Mining Investment Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Nuclear Inner Mongolia Mining Investment Co ltd filed Critical China Nuclear Inner Mongolia Mining Investment Co ltd
Priority to CN202011047965.8A priority Critical patent/CN114329301B/en
Publication of CN114329301A publication Critical patent/CN114329301A/en
Application granted granted Critical
Publication of CN114329301B publication Critical patent/CN114329301B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Earth Drilling (AREA)

Abstract

本发明属于计算方法,具体涉及一种铀矿钻孔见矿切穿点水平投影批量计算方法。它包括:步骤1,准备钻孔数据库:建立数据库结构,并录入钻孔数据;步骤2,调用钻孔测斜记录,计算各控制点坐标:调用钻孔测斜记录,并将计算结果存储;步骤3,调用钻孔见矿参数,确定见矿切穿点上下控制点:调用钻孔单工程见矿参数,确定矿体上下控制序号,并将结果存储;步骤4,计算钻孔见矿切穿点坐标,并输出坐标:依据见矿切穿点上下控制点,计算钻孔见矿切穿点坐标,并将切穿点投影输出。本发明的显著效果是:提高储量平面图的绘图效率;节省大量计算时间;方便快速查询导出数据。

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

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.

Claims (1)

1.一种铀矿钻孔见矿切穿点水平投影批量计算方法,其特征在于,该方法所包括的如下步骤:1. A method for batch calculation of horizontal projections of uranium ore drilling penetration points, characterized in that the method comprises the following steps: 步骤1,准备钻孔数据库:建立数据库结构,并录入钻孔数据,该数据库包含钻孔录入信息结构表和钻孔计算信息结构表;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; 步骤2,调用钻孔测斜记录,计算各控制点坐标:调用钻孔测斜记录,采用全角半距法计算各控制点相关参数,并根据各控制点参数计算各控制点坐标,并将计算结果存储;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; 步骤3,调用钻孔见矿参数,确定见矿切穿点上下控制点:调用钻孔单工程见矿参数,确定矿体上下控制序号,并将结果存储;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; 步骤4,计算钻孔见矿切穿点坐标,并输出坐标:依据见矿切穿点上下控制点,计算钻孔见矿切穿点坐标,并将切穿点投影输出,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. 所述步骤1的具体步骤如下:The specific steps of step 1 are as follows: 步骤1.1,以Access数据库结构为基础,设计钻孔不同类别信息数据表,录入信息结构表包括基本信息表、测斜记录表、矿段信息表;计算信息结构表包含测斜计算结果表、钻孔见矿切穿点计算表,Step 1.1, based on the Access database structure, design different types of drilling information data tables, the input information structure table includes the basic information table, the inclination record table, and the mining section information table; the calculation information structure table includes the inclination calculation result table and the drilling penetration point calculation table. 步骤1.2,将钻孔数据录入到数据库中对应的数据表中存储,该数据表为基本信息表,Step 1.2, enter the drilling data into the corresponding data table in the database for storage, which is the basic information table. 所述的步骤1.2中的基本信息表包含钻孔号holeID、孔口坐标(X0,Y0,Z0);测斜记录表包含钻孔号holeID、测斜点号N、测点孔深h、方位角β、倾角α;矿段信息表包含钻孔号holeID、矿层位置Ha和Hb、矿体编号M,The basic information table in step 1.2 includes the borehole number holeID and the hole mouth coordinates (X 0 , Y 0 , Z 0 ); 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 positions H a and H b , and the ore body number M. 所述步骤2的具体步骤如下:The specific steps of step 2 are as follows: 步骤2.1,调用钻孔测斜记录,采用全角半距法,计算各控制点深度Hi和控制长度Li,并将计算结果存储到测斜计算结果表中,具体计算公式如下:Step 2.1, call the borehole inclinometer record, use the full-angle half-distance method to calculate the depth Hi and control length Li of each control point, 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为第i测点序号;hi为第i测点深度;Hi为第i控制点深度;Li为i控制点控制长度,i is the serial number of the i-th measuring point; hi is the depth of the i-th measuring point; Hi is the depth of the i-th control point; Li is the control length of the i-th control point, 步骤2.2,确定坐标系方向,X为大地坐标X方向,Y为大地坐标Y方向,Z为高程方向,调用钻孔测斜记录,计算各控制段三方向的位移量,并将计算结果存储到测斜计算结果表中,具体计算公式如下: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: 控制段X方向位移量:△Xi=Licosαicosβi Control segment displacement in X direction: △X i =L i cos α i cos β i 控制段Y方向位移量:△Yi=Licosαisinβi Displacement of the control segment in the Y direction: △Y i =L i cosα i sinβ i 控制段Z方向位移量:△Zi=-Lisinαi Displacement of the control section in the Z direction: △Z i = -L i sinα i αi为第i测点倾角;βi为第i测点方位角;α i is the inclination angle of the i-th measuring point; β i is the azimuth angle of the i-th measuring point; △Xi为第i控制段X方向位移量;△Yi为第i控制段Y方向位移量;△Zi为第i控制段Z方向位移量,Xi is the displacement of the ith control segment in the X direction; △ Yi is the displacement of the ith control segment in the Y direction; △ Zi is the displacement of the ith control segment in the Z direction. 步骤2.3,调用测斜计算结果表中步骤2.2计算的三方向位移量,计算三方向累计位移量,并将计算结果存储到测斜计算结果表中,具体计算公式如下: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: X方向累计位移∑ΔXi=△X1+△X2+△X3+…+△XiCumulative displacement in X direction ∑ΔX i = ΔX 1 + ΔX 2 + ΔX 3 + … + ΔX i ; Y方向累计位移∑ΔYi=△Y1+△Y2+△Y3+…+△YiCumulative displacement in Y direction ∑ΔY i = ΔY 1 + ΔY 2 + ΔY 3 + … + ΔY i ; Z方向累计位移∑ΔZi=△Z1+△Z2+△Z3+…+△ZiThe cumulative displacement in the Z direction ∑ΔZ i = ΔZ 1 + ΔZ 2 + ΔZ 3 + … + ΔZ i ; 步骤2.4,依据孔口坐标(X0,Y0,Z0)和测斜计算结果表中步骤2.3计算的三方向累计位移量,计算各控制点坐标Oi(Xi,Yi,Zi),并将计算结果存储到测斜计算结果表中,计算公式如下:Step 2.4, calculate the coordinates of each control point O i (X i , Yi i , Zi ) based on the orifice coordinates (X 0 , Y 0 , Z 0 ) and the cumulative displacements in three directions calculated in step 2.3 in the inclinometer calculation result table, and store the calculation results in the inclinometer calculation result table. The calculation formula is as follows: Xi=X0+∑ΔXi Xi = X0 + ∑ΔXi ; Yi=Y0+∑ΔYi Yi = Y0 + ∑ΔYi ; Zi=Z0+∑ΔZi Zi = Z0 + ∑ΔZi ; Xi为第i控制点X坐标;Yi为第i控制点Y坐标;Zi为第i控制点Z坐标; Xi is the X coordinate of the ith control point; Yi is the Y coordinate of the ith control point; Zi is the Z coordinate of the ith control point; ∑△Xi为第i控制点X向累计位移;∑△Yi为第i控制点Y向累计位移;∑ΔZi为第i控制点Z向累计位移,∑△X i is the cumulative displacement of the ith control point in the X direction; ∑△Y i is the cumulative displacement of the ith control point in the Y direction; ∑ΔZ i is the cumulative displacement of the ith control point in the Z direction, 所述步骤3的具体步骤如下:The specific steps of step 3 are as follows: 步骤3.1,调用矿段信息表中钻孔单工程见矿参数,计算钻孔各矿体见矿切穿点孔深PM,并将结果存储在钻孔见矿切穿点计算表中,具体计算公式如下:Step 3.1, call the drilling single project ore-seeing parameters in the ore section information table, calculate the hole depth PM of each ore body of the drilling hole, and store the result in the drilling hole ore-seeing and cutting point calculation table. The specific calculation formula is as follows: PM=0.5(Ha+Hb)或Hb PM = 0.5 (H a + H b ) or H b 式中:PM=0.5(Ha+Hb)表示钻孔第M编号矿体切穿点为矿体中心曲面与钻孔的交点孔深,Where: PM = 0.5 ( Ha + Hb ) means that the penetration point of the Mth ore body is the intersection of the central curved surface of the ore body and the depth of the drill hole. PM=Hb表示钻孔第M编号矿体切穿点为矿体曲面底板与钻孔的交点孔深, PM = Hb means that the penetration point of the Mth ore body is the intersection of the bottom plate of the ore body and the depth of the drill hole. 步骤3.2,调用测斜计算结果表中步骤2计算的各控制点孔深和坐标,确定钻孔见矿切穿点上下控制点序号,判定方法如下: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: 令C=min|PM-Hj|;Let C = min |P M - H j |; 若PM-Hj<0,则j为见矿切穿点下控制点序号,j-1为见矿切穿点上控制点序号;If P M -H j <0, then j is the number of the lower control point of the ore-seeking and cutting point, and j-1 is the number of the upper control point of the ore-seeking and cutting point; 若PM-Hj>0,则j为见矿切穿点上控制点序号,j+1为见矿切穿点下控制点序号,If PM - Hj >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. 判定式中:Hj为第j控制点深度;j为切穿点上或下控制点序号,In the judgment formula: Hj is the depth of the jth control point; j is the sequence number of the upper or lower control point of the cutting point, 步骤3.3,根据步骤3.2计算的矿切穿点上、下控制点序号,调用测斜计算结果表中的上、下控制点的水平坐标和孔深,并将结果存储在钻孔见矿切穿点计算表,存储记录Record={holeID,M,PM,X,Y,H,X,Y,H};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 mine cutting point calculation table, storing the record Record = {holeID, M, PM , Xup , Yup , Hup , Xdown , Ydown , Hdown }; 存储记录中,X、Y为切穿点上控制点水平坐标,X,Y为切穿点下控制点水平坐标,H为切穿点上控制点孔深;H为切穿点下控制点孔深,In the storage record, X -up and Y -up are the horizontal coordinates of the upper control point of the cutting point, X -down and Y -down are the horizontal coordinates of the lower control point of the cutting point, H -up is the hole depth of the upper control point of the cutting point, and H -down is the hole depth of the lower control point of the cutting point. 所述步骤4的具体步骤如下:The specific steps of step 4 are as follows: 步骤4.1,调用步骤3.3在钻孔见矿切穿点计算表中的存储记录,计算钻孔见矿切穿点水平投影坐标(XM,YM),并将结果存储到钻孔见矿切穿点计算表中,Step 4.1, call the stored record in the drill hole penetration point calculation table in step 3.3, calculate the horizontal projection coordinates (X M , Y M ) of the drill hole penetration point, and store the result 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: XM=(X-X)*(PM-H)/(H-H)+X X M = (X down - X up ) * (P M - H up ) / (H down - H up ) + X up YM=(Y-Y)*(PM-H)/(H-H)+Y Y M = (Y down - Y up ) * (P M - H up ) / (H down - H up ) + Y up 步骤4.2,批量输出步骤4.1算的钻孔见矿切穿点水平投影坐标,输出Output={holeID,M,XM,YM}。Step 4.2, batch output the horizontal projection coordinates of the drill hole penetration points calculated in step 4.1, and output Output = {holeID, M, X M , Y M }.
CN202011047965.8A 2020-09-29 2020-09-29 A batch calculation method for horizontal projection of uranium drilling penetration points Active CN114329301B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011047965.8A CN114329301B (en) 2020-09-29 2020-09-29 A batch calculation method for horizontal projection of uranium drilling penetration points

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011047965.8A CN114329301B (en) 2020-09-29 2020-09-29 A batch calculation method for horizontal projection of uranium drilling penetration points

Publications (2)

Publication Number Publication Date
CN114329301A CN114329301A (en) 2022-04-12
CN114329301B true CN114329301B (en) 2025-02-14

Family

ID=81011136

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011047965.8A Active CN114329301B (en) 2020-09-29 2020-09-29 A batch calculation method for horizontal projection of uranium drilling penetration points

Country Status (1)

Country Link
CN (1) CN114329301B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119124065A (en) * 2024-07-01 2024-12-13 大连理工大学 A slope displacement monitoring and early warning method and inclinometer

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110648378A (en) * 2019-09-05 2020-01-03 核工业二四三大队 Electronic drawing method for horizontal projection drawing of sandstone-type uranium ore body

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7089266B2 (en) * 2003-06-02 2006-08-08 The Board Of Trustees Of The Leland Stanford Jr. University Computer systems and methods for the query and visualization of multidimensional databases
CN111488328B (en) * 2019-12-17 2023-06-13 核工业二四三大队 Automatic drawing method for uranium ore geological drilling core geophysical prospecting and recording curve

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110648378A (en) * 2019-09-05 2020-01-03 核工业二四三大队 Electronic drawing method for horizontal projection drawing of sandstone-type uranium ore body

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
地质勘探工程数据模型的可视化研究;任佳 等;软件导刊;20101030(第10期);第123-125 *

Also Published As

Publication number Publication date
CN114329301A (en) 2022-04-12

Similar Documents

Publication Publication Date Title
CN105631753B (en) Horizontal well oil reservoir profile modeling method based on stratum attitude
JP6550358B2 (en) Construction time prediction system for construction machinery
CN107264724A (en) It is used for the baseline design method that fitting-out piping is installed in the anti-state segmentation of cabin curved surface
CN112859935A (en) Leveling control method for multi-point supporting platform
CN112576257A (en) Tunnel blast hole lofting positioning method
CN101078627A (en) On-line calibration method for shield machine automatic guiding system based on optical fiber gyro and PSD laser target
KR102604505B1 (en) Drilling device and drilling management system including the same
CN109614743A (en) Excavator and bucket positioning method, electronic equipment, storage medium
CN114329301B (en) A batch calculation method for horizontal projection of uranium drilling penetration points
CN108571287B (en) The hole trajectory control system of Kernel-based methods control
CN110555559B (en) A method for predicting land subsidence under the condition of insufficient mining of thick loose layer
CN106638318B (en) Locating method of anchor point at end of cable-stayed bridge and installation method of cable conduit
CN110530358A (en) Car body navigation positional device and navigation system and method
CN114411840B (en) Land leveling control method and device and excavator
CN109196168A (en) Equipment control device and Work machine
CN103925854B (en) A kind of guided missile level measurement method based on three-dimensional coordinate transfer principle
CN216052188U (en) Guiding device of shaft heading machine
CN118065924A (en) An intelligent tunneling construction method based on tunnel boring machine
CN104931027B (en) Projection paying-off method for decorative pattern
CN104635528B (en) Diclinic cutting section accuracy control method
CN118601536A (en) Open-pit drilling rig and control method thereof
CN216049921U (en) High-precision large-water-depth underwater engineering positioning system
Shen et al. Development of virtual laser target board for tunnel boring machine guidance control
WO2024037120A1 (en) Precise positioning system for hydraulic impact hammer
CN214502334U (en) Blast hole arrangement device for tunnel surrounding rock three-arm drill jumbo

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant