CN117152301B - Geological boundary drawing system based on geological point occurrence and coordinates - Google Patents
Geological boundary drawing system based on geological point occurrence and coordinates Download PDFInfo
- Publication number
- CN117152301B CN117152301B CN202311424841.0A CN202311424841A CN117152301B CN 117152301 B CN117152301 B CN 117152301B CN 202311424841 A CN202311424841 A CN 202311424841A CN 117152301 B CN117152301 B CN 117152301B
- Authority
- CN
- China
- Prior art keywords
- point
- geological
- geological boundary
- boundary
- lines
- 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
Links
- 239000002245 particle Substances 0.000 claims abstract description 47
- 238000000034 method Methods 0.000 claims description 49
- 238000005553 drilling Methods 0.000 claims description 29
- 238000013507 mapping Methods 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 13
- 238000010276 construction Methods 0.000 claims description 4
- 238000011161 development Methods 0.000 claims description 2
- 238000011835 investigation Methods 0.000 abstract 1
- 239000011229 interlayer Substances 0.000 description 18
- 239000010410 layer Substances 0.000 description 10
- 239000011435 rock Substances 0.000 description 9
- 238000012876 topography Methods 0.000 description 8
- 238000001514 detection method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000005452 bending Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 3
- 239000000523 sample Substances 0.000 description 2
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000003306 harvesting Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a geological boundary drawing system based on geological point occurrence and coordinates. The system belongs to the technical field of engineering geological investigation, and comprises a geological point information input module, an auxiliary line generation module and a geological boundary generation module; the geological point information input module is used for inputting geological point information; the auxiliary line generation module is used for generating auxiliary lines according to the geological point information and preset contour lines and assigning values to the elevations of the auxiliary lines; and the geological boundary generation module is used for acquiring the dew point of the geological boundary formed by intersecting the contour lines with the same elevation in the control range of the ground particles and the auxiliary lines, and connecting the dew points of the geological boundary to generate the geological boundary. The invention greatly improves the working efficiency and the accuracy of drawing the drawing.
Description
Technical Field
The invention relates to the technical field of engineering geological exploration, in particular to a geological boundary drawing system based on geological point occurrence and coordinates.
Background
Geological boundaries refer to the intersection of various geological points such as strata, rock mass, structural layers and the like in a geological map with the ground surface. The accurate drawing of geological boundary is the premise and basis for evaluating the geological conditions of site engineering and reasonably determining the exploration work.
In engineering geological mapping, the geological boundary can be accurately drawn by directly tracking the coordinates of the geological points in the areas with better exposure conditions of the geological points, but in the areas with wide coverage layer distribution and poor exposure conditions, the exposure is usually manually disclosed, after the occurrence of the geological points of the observation points is determined, the geological boundary is extended to two sides along the stratum trend by 1/2 line distance according to a V-shaped method. Because the V-shaped rule is only expressed by the geological phenomenon, the exposed boundary line of the geological point is influenced by the occurrence and the topography, and the bending degree of the exposed curve can be drawn only by depending on the experience of geological personnel. Therefore, the geological boundary between the observation points is less accurate.
Disclosure of Invention
The purpose of the invention is that: a geological boundary drawing system based on geological point occurrence and coordinates is provided. The invention greatly improves the working efficiency and the accuracy of drawing the drawing.
The technical scheme of the invention is as follows: a geological boundary drawing system based on geological point occurrence and coordinates comprises a geological point information input module, an auxiliary line generation module and a geological boundary generation module; the geological point information input module is used for inputting geological point information; the auxiliary line generation module is used for generating auxiliary lines according to the geological point information and preset contour lines and assigning values to the elevations of the auxiliary lines; the geological boundary generation module is used for acquiring geological boundary dew points formed by intersecting the contour lines with the same elevation in the control range of the ground particles and the auxiliary lines, and connecting the dew points of the geological boundary to generate the geological boundary;
the geological point information comprises coordinates, elevation, inclination and inclination angle of the geological points;
the auxiliary line is generated as follows:
according to the formulaL=Δh·cot(α) Calculating pay-off distanceLIn the followingΔhAt the same time, the distance between the contour lines is equal to the distance between the contour lines,αis an inclination angle; along the trend of the ground mass point, generating a coordinate passing through the geological pointP 0 (X,Y) Auxiliary line of (2)lThe method comprises the steps of carrying out a first treatment on the surface of the Based on auxiliary lineslThen according to the pay-off distanceLGenerating a series parallel tolWherein the auxiliary lines of inclination along the geological point are numbered in sequence:l-1,l-2,…,l-nthe method comprises the steps of carrying out a first treatment on the surface of the Auxiliary lines opposite to the tendency of ground particles are numbered in sequence:l+1,l+2,…,l+n；
the assignment process of the auxiliary line elevation is as follows:
elevation of ground pointHAssigning values to auxiliary linesl；
Will beH-Δh,H-2·Δh,…,H-n·ΔhAssignment to auxiliary linel-1,l-2,…,l-n；
Will beH+Δh,H+2·Δh,…,H+n·ΔhAssignment to auxiliary linel+1,l+2,…,l+n；
The process of connecting the geological boundary to the dew point is as follows: acquiring a starting point of a geological boundary; based on the starting point, connecting dew points of all geological boundary lines according to a shortest path method;
the starting point is obtained as follows: make the following stepsConstructing a point set by using all obtained geological boundary exposure pointsP n Point-to-point setP n The dew points of each geological boundary line in the map are ordered according to the elevation, and the found out dew point construction point set of each geological boundary line with the largest elevation is usedP k The method comprises the steps of carrying out a first treatment on the surface of the Separately calculating the geological pointsP 0 And point setP k Distance between dew points of each geological boundaryd i ，d i The dew point of the geological boundary corresponding to the maximum value of (2) is taken as a starting point;
the shortest path method is specifically as follows: starting from a starting point, starting from a set of pointsP n Finding out another geological boundary nearest to the starting point to find out the dew point as a 2 nd point; from the point setP n Finding out the dew point of another geological boundary closest to the 2 nd point and excluding the starting point as the 3 rd point; from the point setP n Finding out the nearest geological boundary line from the 3 rd point, excluding the starting point, the 2 nd point and the 3 rd point, and taking the dew point as the 4 th point; analogize to traverse the set of pointsP n The dew points of all the geological boundary lines are found out, the dew points of all the geological boundary lines forming the geological boundary lines are found out, and the obtained geological boundary lines are sequentially connected.
In the geological boundary drawing system based on the geological occurrence and coordinates, when the local mass point is an underground mass point, the geological point is revealed by drilling, at this time, the coordinates of the ground mass point are the hole position coordinates of the drilling hole, and the elevation of the ground mass point is equal to the difference between the hole elevation of the drilling hole and the depth of the ground mass point revealed by the drilling hole.
In the geological boundary drawing system based on the geological point occurrence and the coordinates, the system is integrated into an AutoCAD platform in a secondary development plug-in mode in an AutoLISP environment of AutoCAD.
The invention has the advantages that: according to the invention, quick and accurate drawing of the geological boundary based on the occurrence and coordinates of the geological points is realized, the line-placing distance is not required to be calculated manually in the drawing process, the complex auxiliary line is not required to be drawn manually to determine the intersection point of the rock stratum boundary and the contour line, the line of the dew point of the geological boundary is not required to be connected manually, the geological boundary can be automatically generated only by knowing the information of the occurrence of the coordinates and Gao Chengji of the ground points, and the drawing efficiency is greatly improved.
The connecting line between the projection points of the intersection points of the different elevation trend lines of the rock stratum and the corresponding elevation topography contour lines on the horizontal plane is the geological boundary line. The geological boundary is plotted according to the pay-off distance, and is obtained by plotting based on the relation between the contour line and the projection of the geological point on the plane. The contour lines in the plan view only reflect the topography of the ground surface, and the contour lines below the ground surface are absent. Therefore, conventional pay-off distance mapping methods do not allow for the mapping of geological boundaries of subsurface earth particles.
According to the coordinates of the holes and the heights of the holes at the positions of the ground points, the coordinates and the heights of the ground points exposed in the holes are calculated by combining the depths of the ground points exposed in the holes: the coordinates of the ground particles are the hole site coordinates of the drill holes, and the elevation of the ground particles is equal to the difference between the hole elevation of the drill holes and the depth of the ground particles disclosed by the drill holes. Using the ground particles revealed by the drilling holes as reference points, and generating an auxiliary line passing through the coordinates of the ground points along the trend of the ground particleslThen a series of parallel lines are made with the paying-off distance as the intervallAnd is high Cheng Fuzhi. The Gao Chengda auxiliary line above the hole elevation belongs to the auxiliary line below the ground surface between the hole elevation and the elevation of the ground particles, but belongs to the trend line of the ground particles disclosed by the drilling holes at different elevations both above and below the ground surface. At the moment, the intersection point of the auxiliary line and the contour line, which is obtained by CAD automatic calculation and has the same elevation as the drilling hole, is the dew point of the geological point revealed by the drilling on the geological boundary line of the earth surface; then, according to the method of determining the geological boundary by using the dew point of the geological boundary by using the earth surface particles, the geological boundary of the earth surface particles is obtained, so that the aim of exposing the earth particles based on drilling and programming the geological boundary is fulfilled. Based on the technical principle, the invention can not only realize the tracking of geological boundary lines of known geological points exposed on the earth surface, but also realize the drawing of geological boundary lines of geological points such as faults, interlayers and the like exposed in drilling, namely, the invention can realize the drawing of geological boundary lines of earth surface and underground earth particles.
The method can rapidly and accurately infer the specific exposed boundary of the ground particles such as stratum and rock stratum based on the observed geologic point occurrence and coordinates, and improves the efficiency and accuracy of geologic mapping. Meanwhile, for monoclinic rock stratum with poor outcrop and stable production, a more complete geological map can be rapidly drawn according to geological boundary information revealed by drilling.
The geological boundary line drawn by the invention is based on the pay-off distance principle, the outcrop boundary line drawn is based on accurate calculation to obtain the intersection point of the geological boundary line and each contour line, and then the intersection point is connected. Therefore, the connecting line between the ground particles automatically bends along with the bending change of the contour lines at the positions of the gullies and the ridges, and the quantification of the V-shaped rule when the geological boundary line passes through the gullies or the ridges is effectively realized.
According to the invention, geological boundary drawing is integrated into an AUTOCAD through a plug-in form, so that the process of geological boundary drawing is realized, traditional manual calculation, assignment and connection are handed to a computer, automation is realized, the problem that geological boundary defects of ground particles disclosed by a borehole cannot be compiled and drawn by a traditional method is solved, automatic drawing of geological boundary of ground surface mapping points and ground particles disclosed by the borehole is realized, and the working efficiency and drawing accuracy are greatly improved.
In traditional engineering geological mapping, in areas with better exposed geological points, geological boundaries can be accurately drawn by directly tracking coordinates of particles in the mapping areas. In the area with widely distributed covering layers and poor outcrop conditions, manual exposure is usually needed, after the occurrence of the geological point of the observation point is determined, the geological boundary is extended to two sides along the stratum trend by 1/2 line distance according to the V-shaped method. The V-shaped rule is only expressed by the geological phenomenon, the exposed line of the ground particles is affected by the occurrence and the topography, and the bending degree of the exposed line can only be drawn by relying on the experience of geology staff, so that the accuracy of geological boundary between observation points is poor. According to the water and electricity engineering geological mapping regulations (NB-T10074-2018), the distance between mapping geological points is 2 cm-3 cm on the scale chart of the corresponding chart. Taking 1:1000 geological map mapping as an example, the distance between mapping points in an actual site is 200-300 m, the geological boundary observed at the mapping points is accurate, but the geological boundary between two mapping points (about 200-300 m distance) is determined according to a V-shaped rule, and the accuracy is poor due to the fact that the geological boundary is drawn depending on experience of geological personnel. The invention relates to a method for quantitatively drawing geological boundary between two known ground particles, which is characterized in that the geological boundary of an inclined rock stratum is drawn according to a line-setting distance, specifically, the method is based on the relation between contour lines and projection of geological points on a plane, deduces based on a strict mathematical formula, accurately determines the dew point of the geological boundary on each contour line between the two mapping points through a drawing method, and is equivalent to adding a series of accurate geological mapping points at different elevations between the two mapping points, thereby improving the precision and accuracy of drawing the geological boundary.
The rock mass at the fault and weak interlayer part is broken and is easy to be weathered, and is mostly covered by the covering layer. There is usually no good outcrop in the building site, and in order to find out the concrete position, character and other basic geological information of faults or weak interlayers, a great number of tank probes or vertical shafts need to be excavated in the whole site area to provide guidance advice for engineering construction or treatment. Because the thickness of the fault layer and the weak interlayer is generally not large and is covered by the covering layer in the field range, a large number of groove probes are required to be longitudinally and transversely arranged for determination, and the workload is large. According to the invention, the outcrop point of the known fault or weak interlayer can be quickly and accurately determined, the spreading direction of the fault or weak interlayer on a topographic map is drawn, and a small amount of groove detection can be arranged on site for verification after the extending direction is determined, so that unnecessary workload such as groove detection and the like is effectively reduced. The method comprises the following steps: according to the invention, fault and weak interlayer points with good exposure conditions can be found around the site area, the boundary line of the fault or the weak interlayer of the site area can be rapidly obtained on the indoor topographic map, and then the groove detection or the vertical well is arranged near the fault or the weak interlayer, so that the groove detection and the vertical well are arranged more pertinently, and the engineering quantity of the groove detection or the vertical well is reduced. Specifically, the process of reducing unnecessary workload such as groove detection by using the method comprises the following steps:
(1) Searching a specific position or section of a fault or a weak interlayer near or around a site, and acquiring the coordinate, elevation and attitude basic information of the dew point through a measuring instrument;
(2) The dew point of the fault or the weak interlayer is on the topographic map, the exposed boundary line of the fault or the weak interlayer is rapidly drawn on the topographic map by using the method, and coordinates are obtained;
(3) Introducing boundary coordinates of faults or interlayers into a measuring instrument, and lofting by using the measuring instrument on site;
(4) The position of the fault or weak interlayer is verified by manually excavating a trench or shaft at or near the presumed fault or interlayer boundary.
Traditionally, auxiliary lines are drawn in CAD by array commands of CAD, a series of equidistant parallel auxiliary lines are generated according to line spacing, but the auxiliary lines are free of elevation attributes, each auxiliary line needs to be manually high Cheng Fuzhi through attribute buttons, and the auxiliary lines are quite large in workload and easy to make mistakes. After the invention is integrated into CAD through the plug-in form, through simple interface operation: selecting the position of the ground particles, inputting the elevation and the occurrence of the ground particles, inputting the interval between the contour lines, and clicking the button for generating the auxiliary line (see figure 3), thereby realizing the automation of the assignment of the auxiliary line and the elevation, and greatly improving the working efficiency and the accuracy. Taking 50 auxiliary lines and a height Cheng Fuzhi as an example, the traditional method needs 2-3 minutes, and the operation of the invention only needs 3-5 seconds.
In general, it is relatively easy to find the intersection point of the same elevation of the contour line and the auxiliary line in the CAD drawing, and in the case of a sparse contour line with a small scale and a small number of auxiliary lines, manually select the intersection point of the auxiliary line with the same elevation and the contour line as the geological boundary line to generate the dew point. However, for large scales, the contour lines are very dense, and the dew point of the geological boundary is very difficult to obtain by the traditional method. If in a certain geological boundary line sketching, the number of intersection points of 61 auxiliary lines and contour lines exceeds 4000, but only 24 intersection points are geological boundary line dew points, and among 4000 intersection points of CAD, 24 geological boundary lines dew points are found by a manual comparison elevation attribute method, so that the method is time-consuming and labor-consuming and is very prone to error. In addition, the traditional geological boundary line drawing is that the dew point of the geological boundary line found in the steps is manually found by manpower, and discrete points are sequentially connected by using a multi-section line. The method is convenient for connecting lines with fewer geological exposure points, and if the exposure points are more, the method also has the problem of time and labor waste. The auxiliary line is automatically generated by the auxiliary line generating module, the elevation of the auxiliary line is automatically assigned, and the geological boundary dew point in the control range of the geological point is automatically obtained by the geological boundary generating module, so that the efficiency and the accuracy are greatly improved.
Generally, it is relatively easy to manually connect a series of point sets, because a person can determine the sequence of connecting points, but it is relatively difficult for a computer to automatically connect geological boundaries. In order to overcome the technical problem, the invention realizes the automatic connection of dew points of all geological boundaries based on a shortest path method, and simultaneously considers that the automatic connection efficiency of a point set by a shortest path is lower, and the invention optimizes and adjusts the algorithm: the start point of the geological boundary is searched before connecting the lines. The calculation efficiency is greatly improved through the optimization adjustment.
According to the conventional method, 5 to 10 minutes are required to complete the geological boundary drawing of 1 known earth particle (earth surface earth particle), and 30 to 40 minutes are generally required to complete the geological boundary drawing of 1 borehole revealing earth particle (underground earth particle). The invention only needs 0.5 to 1 minute to finish the geological boundary drawing time of the earth surface and underground earth particles, thereby greatly improving the working efficiency and the drawing accuracy of the drawing.
Drawings
FIG. 1 is a step diagram of the invention for mapping geological boundaries in CAD;
FIG. 2 is a diagram of a plug-in interface of the present invention after integration into a CAD;
FIG. 3 is an interface operational step for implementing auxiliary line generation and elevation assignment automation;
FIG. 4 is a schematic diagram of determining a shortest path computation starting point;
FIG. 5 is a diagram of a conventional method for determining geological boundaries based on borehole reveal geological information;
FIG. 6 illustrates determining geological boundaries according to the present process based on borehole reveal geological information in accordance with the present invention;
FIG. 7 shows the dew point of the geological boundary in the measured topography;
FIG. 8 is a diagram of a load plug-in acquiring earth particle coordinates and inputting an elevation interface;
FIG. 9 is an auxiliary line automatically drawn based on elevation and pay-off distance;
FIG. 10 is a sketch of geological boundaries generated from occurrence and coordinates;
FIG. 11 is a final geological boundary generated by clicking the clear auxiliary line button;
FIG. 12 is an exposed pattern of geological boundaries in a plane;
FIG. 13 is a plot of the spatial spread projection of an inclined formation versus contour.
Detailed Description
The invention is further illustrated by the following figures and examples, which are not intended to be limiting.
Example 1. The invention forms a geological boundary line to show dew point by intersecting the contour lines with the same elevation in the control range of the ground particles and the auxiliary lines, thereby generating the geological boundary line, and the technical principle is as follows:
1) As shown in fig. 12, the principle of geological map programming shows that the straight line connecting the intersection point of the geological boundary exposed by a certain ground particle and the same elevation contour is the trend line of the ground particle. That is, a series of lines with different heights are made on a certain ground point according to the height of the contour line and the line-releasing distance, as shown in fig. 12. The intersection point of the trend line and the corresponding terrain contour line is the dew point of the geological boundary line of the geological point on the geological boundary line, gao Chengji of the dew point of the geological boundary line represents the ground elevation and also represents the elevation of the geological point layer, and the geological boundary line dew points are sequentially connected by smooth curves to form the geological boundary line of the stratum.
2) As shown in FIG. 13, the plane ABCD is inclinedαAs can be seen from FIG. 13, for the dip aspect of the stable attitude, the dip angle isαIs higher by one contour line along the inclination directionDistance of%Δh) The horizontal distance of the projection distance between two adjacent trend lines on the same horizontal plane (i.e. the paying-off distanceL) Has the following relationship:L=Δh·cot(α)。
3) According to the projection correlations between the trend lines of different heights and the topographic lines on the topographic map, for the known geological points, a series of parallel lines parallel to the trend of the rock stratum are made on the topographic map along the trend direction, and the parallel line distance is thatLThus, the trend lines of different elevations of the rock stratum can be obtained, and the intersection point of the trend lines and the topography lines of the corresponding elevations is the dew point of the geological boundary line of the layer surface on the elevation. The geological boundary of the ground particles is obtained by connecting the points from low to high or from high to low by using a smooth curve.
Based on the technical principle, the invention designs a geological boundary drawing system based on the occurrence and the coordinates of ground particles, which is shown in figures 1-6 and comprises a geological point information input module, an auxiliary line generation module and a geological boundary generation module; the geological point information input module is used for inputting geological point information; the auxiliary line generation module is used for generating auxiliary lines according to the geological point information and preset contour lines and assigning values to the elevations of the auxiliary lines; and the geological boundary generation module is used for acquiring the dew point of the geological boundary formed by intersecting the contour lines with the same elevation in the control range of the ground particles and the auxiliary lines, and connecting the dew points of the geological boundary to generate the geological boundary.
The geologic point information includes coordinates, elevation, inclination, and dip of the earth's points.
When the local mass point is an underground mass point, the geological point is revealed through the drilling, the coordinate of the local mass point is the hole position coordinate of the drilling, and the elevation of the local mass point is equal to the difference between the hole elevation of the drilling and the depth of the local mass point revealed by the drilling.
The generation process of the auxiliary line (trend line) is as follows:
according to the formulaL=Δh·cot(α) Calculating pay-off distanceLIn the followingΔhAt the same time, the distance between the contour lines is equal to the distance between the contour lines,αis an inclination angle; along the trend of the ground mass point, generating a coordinate passing through the geological pointP 0 (X,Y) Auxiliary line of (2)lThe method comprises the steps of carrying out a first treatment on the surface of the Based on auxiliary lineslThen according to the pay-off distanceLGenerating a series parallel tolWherein the auxiliary lines of inclination along the geological point are numbered in sequence:l-1,l-2,…,l-nthe method comprises the steps of carrying out a first treatment on the surface of the Auxiliary lines opposite to the tendency of ground particles are numbered in sequence:l+1,l+2,…,l+n. Azimuth of the strike of a ground particle = trend of a ground particleβ+90°. The wire-releasing distance is as follows: the horizontal distance of the projection interval of two adjacent trend lines on the same horizontal plane; the trend line is a straight line connecting the geological boundary line and the intersection point of the same elevation contour line.
The assignment process of the auxiliary line elevation is as follows:
elevation of ground pointHAssigning values to auxiliary linesl；
Will beH-Δh,H-2·Δh,…,H-n·ΔhAssignment to auxiliary linel-1,l-2,…,l-n；
Will beH+Δh,H+2·Δh,…,H+n·ΔhAssignment to auxiliary linel+1,l+2,…,l+n。
The aforementioned process of connecting the geological boundary to the dew point is as follows: acquiring a starting point of a geological boundary; and connecting the dew points of the geological boundary lines according to the shortest path method based on the starting point.
The aforementioned acquisition process of the start point is as follows (see fig. 4): constructing a set of points using all geological boundary exposure points obtainedP n Point-to-point setP n The dew points of each geological boundary line in the map are ordered according to the elevation, and the found out dew point construction point set of each geological boundary line with the largest elevation is usedP k (in figure 4,P 1 、P 2 belongs to the point setP k ) The method comprises the steps of carrying out a first treatment on the surface of the Separately calculating the geological pointsP 0 And point setP k Distance between dew points of each geological boundaryd i ，d i Maximum value of (2)The dew point of the corresponding geological boundary is taken as a starting point;
the shortest path method is specifically as follows:
starting from a starting point, starting from a set of pointsP n Finding out another geological boundary nearest to the starting point to find out the dew point as a 2 nd point; from the point setP n Finding out the dew point of another geological boundary closest to the 2 nd point and excluding the starting point as the 3 rd point; from the point setP n Finding out the nearest geological boundary line from the 3 rd point, excluding the starting point, the 2 nd point and the 3 rd point, and taking the dew point as the 4 th point; analogize to traverse the set of pointsP n The dew points of all the geological boundary lines are found out, the dew points of all the geological boundary lines forming the geological boundary lines are found out, and the obtained geological boundary lines are sequentially connected.
Conventionally, determining the surface exposure boundary of a drilling exposure structural surface according to drilling information, and specifically, drawing a series of sections through drilling points, and drawing an engineering geological section according to the depth of the drilling exposure geological boundary and the formation yield, wherein the intersection elevation of the drilling exposure geological boundary and the surface is the surface exposure point of the geological boundary. As shown in fig. 5: the ZK1 drilling elevation is 935m, the drilling depth is 25m, lithology layering is disclosed at 20m, and the stratum yield is 80 DEG & lt 40 deg. According to the geological boundary line disclosed by drilling, the surface exposure boundary line is drawn, firstly, a series of sections S1, S2, S3 and S4 passing through the drilling are drawn, and the surface exposure elevation point of the interlayer is determined according to engineering geological sections. According to the invention, besides the geological boundary can be generated according to the known points of the earth surface, the exposed boundary of the earth surface of the interlayer can be generated according to the interlayer information disclosed by the drilling. Taking fig. 6 as an example, the heights of the geological boundary lines at S1, S2, S3 and S4 according to the present invention are calculated to be completely consistent with the conventional method (the method of fig. 5), which illustrates that the geological boundary lines generated according to the present invention based on the geological information revealed by the borehole are completely correct.
The invention operates steps in AutoCAD as follows:
1) The aforementioned system is integrated into an AutoCAD platform in the form of a secondarily developed plug-in the AutoCAD's autocalap environment, the plug-in interface being seen in fig. 2.
2) And opening a topographic map of the geological boundary to be drawn and ground particles of the field surveying and mapping collection point in the AutoCAD.
3) Loading the plug-in the model space; the specific loading mode is as follows: the "load" is entered in the AutoCAD command line and returned, after which the path in which the "geological boundary" drawing plug-in is located is selected in the pop-up dialog box, and the main program file in the path folder is selected.
4) Activating a geological boundary drawing plug-in a model space after loading; the activation method is to input "geological boundary" to draw the plug-in name (i.e. main program file name), such as "czv", input "czv" in the AutoCAD command line, and open the program interface as carriage returns, as shown in fig. 2 below.
5) Clicking the 'select place particle position' button to select the position of the known place particle in the topographic map, the program automatically acquires the coordinates (X, Y) of the point, and the coordinates can be modified according to the actual engineering situation.
6) Inputting the elevation of the geological point according to the field harvest conditionH) And birth information (trend, tilt angle α).
7) Contour line spacing for inputting contour line of known topographyΔh。
8) Clicking the "generate auxiliary line" button, making a series of parallel lines parallel to the formation strike in the dip direction on the topography as auxiliary lines, and raising Cheng Fuzhi the auxiliary strike lines. Wherein the distance between the auxiliary lines isL=Δh·cot(α) The elevation of the auxiliary line passing through the geological point is consistent with the elevation of the same-place particle of the geological point to be H, and the elevation of the auxiliary line along the direction of inclination of the geological point is sequentially reducedΔhThe auxiliary line heights with opposite directions along the direction of the tendency of the geological point are sequentially increasedΔhHereby the height Cheng Fuzhi of the auxiliary line is completed.
9) Clicking the 'acquire topographic map layer' button causes the program to acquire the layer in which the contour line is located.
10 Clicking a button of 'select a yield auxiliary line', selecting the generated auxiliary line, clicking a button of 'select a contour line', and selecting a contour line which needs to be determined at the intersection point of the geological boundary line in the engineering range.
11 The radius of the drawing circle is input, and the size of the drawing intersection point of the auxiliary line and the contour line of the terrain after drawing is adjusted.
12 The intersection point of the trend auxiliary line and the contour line with the same elevation is the stratum dew point, the 'imaging' button is clicked, the program automatically generates the intersection point of the trend auxiliary line and the contour line of the terrain, the generated points are represented by circles, and a series of points are connected to represent the exposed boundary line of the ground particles on the ground surface.
13 Finally clicking the 'clear auxiliary line' button, automatically deleting auxiliary contour lines generated in the middle process of drawing and crossing points of the exposed boundary lines of the stratum, only retaining the final geological boundary line and displaying the geological boundary line with different colors.
Example 2. In the following, a dew point of a known measured geological point is taken as an example, and the exposed boundary of the geological point is drawn as shown in fig. 7. The known geological point gives a dew point number D12, the yield is: 80 angle 40, coordinates: x=524574.940, y=3079383.214, z= 864.000. The scale of the topographic map is 1:2000, and the interval between the contour lines is 1m.
The card is loaded in the CAD and the coordinates of D12 are obtained, filling in its corresponding shape and elevation, as shown in fig. 8.
Clicking the generate auxiliary line button automatically calculates the line spacing to draw the auxiliary line as shown in fig. 9.
Selecting auxiliary lines and contour lines, clicking the map forming button, and generating a geological boundary sketch is shown in fig. 10.
Clicking on the clear auxiliary line button clears the auxiliary line, creating the final geological boundary as shown in FIG. 11.
Claims (3)
1. The geological boundary drawing system based on the geological point occurrence and coordinates is characterized by comprising a ground point information input module, an auxiliary line generation module and a geological boundary generation module; the geological point information input module is used for inputting geological point information; the auxiliary line generation module is used for generating auxiliary lines according to the geological point information and preset contour lines and assigning values to the elevations of the auxiliary lines; the geological boundary generation module is used for acquiring geological boundary dew points formed by intersecting the contour lines with the same elevation in the control range of the ground particles and the auxiliary lines, and connecting the dew points of the geological boundary to generate the geological boundary;
the geological point information comprises coordinates, elevation, inclination and inclination angle of the geological points;
the auxiliary line is generated as follows:
according to the formulaL=Δh·cot(α) Calculating pay-off distanceLIn the followingΔhAt the same time, the distance between the contour lines is equal to the distance between the contour lines,αis an inclination angle; along the trend of the ground mass point, generating a coordinate passing through the geological pointP 0 (X,Y) Auxiliary line of (2)lThe method comprises the steps of carrying out a first treatment on the surface of the Based on auxiliary lineslThen according to the pay-off distanceLGenerating a series parallel tolWherein the auxiliary lines along the trends of the ground particles are numbered sequentially:l-1,l-2,…,l-nthe method comprises the steps of carrying out a first treatment on the surface of the Auxiliary lines opposite to the tendency of ground particles are numbered in sequence:l+1,l+2,…,l+n；
the assignment process of the auxiliary line elevation is as follows:
elevation of ground pointHAssigning values to auxiliary linesl；
Will beH-Δh,H-2·Δh,…,H-n·ΔhAssignment to auxiliary linel-1,l-2,…,l-n；
Will beH+Δh,H+2·Δh,…,H+n·ΔhAssignment to auxiliary linel+1,l+2,…,l+n；
The process of connecting the geological boundary to the dew point is as follows: acquiring a starting point of a geological boundary; based on the starting point, connecting dew points of all geological boundary lines according to a shortest path method;
the starting point is obtained as follows: constructing a set of points using all geological boundary exposure points obtainedP n Point-to-point setP n The dew points of each geological boundary line in the map are ordered according to the elevation, and the found out dew point construction point set of each geological boundary line with the largest elevation is usedP k The method comprises the steps of carrying out a first treatment on the surface of the Separately calculating the geological pointsP 0 Point of attachmentCollection setP k Distance between dew points of each geological boundaryd i ，d i The dew point of the geological boundary corresponding to the maximum value of (2) is taken as a starting point;
the shortest path method is specifically as follows: starting from a starting point, starting from a set of pointsP n Finding out another geological boundary nearest to the starting point to find out the dew point as a 2 nd point; from the point setP n Finding out the dew point of another geological boundary closest to the 2 nd point and excluding the starting point as the 3 rd point; from the point setP n Finding out the nearest geological boundary line from the 3 rd point, excluding the starting point, the 2 nd point and the 3 rd point, and taking the dew point as the 4 th point; analogize to traverse the set of pointsP n The dew points of all the geological boundary lines are found out, the dew points of all the geological boundary lines forming the geological boundary lines are found out, and the obtained geological boundary lines are sequentially connected.
2. The geological boundary drawing system according to claim 1, wherein when the local mass point is an underground mass point, the geological point is revealed by drilling, the coordinate of the local mass point is the hole position coordinate of the drilling hole, and the elevation of the local mass point is equal to the difference between the hole elevation of the drilling hole and the depth of the local mass point revealed by the drilling hole.
3. The geological boundary mapping system based on geologic point yield and coordinates of claim 1, wherein the system is integrated into an AutoCAD platform in the form of a plug-in for secondary development in an AutoCAD's AutoCAD environment.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311424841.0A CN117152301B (en) | 2023-10-31 | 2023-10-31 | Geological boundary drawing system based on geological point occurrence and coordinates |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311424841.0A CN117152301B (en) | 2023-10-31 | 2023-10-31 | Geological boundary drawing system based on geological point occurrence and coordinates |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117152301A CN117152301A (en) | 2023-12-01 |
| CN117152301B true CN117152301B (en) | 2024-02-06 |
Family
ID=88903139
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311424841.0A Active CN117152301B (en) | 2023-10-31 | 2023-10-31 | Geological boundary drawing system based on geological point occurrence and coordinates |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117152301B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104408755A (en) * | 2014-10-30 | 2015-03-11 | 中国电建集团成都勘测设计研究院有限公司 | Geologic plane graph drawing method based on structural plane attitude and geologic observation spot coordinate |
| CN106600661A (en) * | 2016-12-20 | 2017-04-26 | 黄河勘测规划设计有限公司 | Method of accurately producing arc geological section |
| CN107065009A (en) * | 2017-04-18 | 2017-08-18 | 中国石油大学(华东) | A kind of buried hill oil gas field underground geological map preparation method based on seismic data |
| JP2018004494A (en) * | 2016-07-04 | 2018-01-11 | 株式会社奥村組 | Method for predicting geological boundary surface or fault surface |
| CN113539051A (en) * | 2021-06-24 | 2021-10-22 | 南京师范大学 | Geological map-based stratum boundary point-by-point rock stratum attitude obtaining method and device |
| CN113819894A (en) * | 2021-09-29 | 2021-12-21 | 中水北方勘测设计研究有限责任公司 | Point-line lofting method for field geological survey by using high-precision oblique photography model |
| CN115272619A (en) * | 2022-09-28 | 2022-11-01 | 东华理工大学南昌校区 | Geological map connecting method |
| CN116228998A (en) * | 2023-01-06 | 2023-06-06 | 广东佛山地质工程勘察院 | Automatic generation method of cut geological profile based on geological space database |
-
2023
- 2023-10-31 CN CN202311424841.0A patent/CN117152301B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104408755A (en) * | 2014-10-30 | 2015-03-11 | 中国电建集团成都勘测设计研究院有限公司 | Geologic plane graph drawing method based on structural plane attitude and geologic observation spot coordinate |
| JP2018004494A (en) * | 2016-07-04 | 2018-01-11 | 株式会社奥村組 | Method for predicting geological boundary surface or fault surface |
| CN106600661A (en) * | 2016-12-20 | 2017-04-26 | 黄河勘测规划设计有限公司 | Method of accurately producing arc geological section |
| CN107065009A (en) * | 2017-04-18 | 2017-08-18 | 中国石油大学(华东) | A kind of buried hill oil gas field underground geological map preparation method based on seismic data |
| CN113539051A (en) * | 2021-06-24 | 2021-10-22 | 南京师范大学 | Geological map-based stratum boundary point-by-point rock stratum attitude obtaining method and device |
| CN113819894A (en) * | 2021-09-29 | 2021-12-21 | 中水北方勘测设计研究有限责任公司 | Point-line lofting method for field geological survey by using high-precision oblique photography model |
| CN115272619A (en) * | 2022-09-28 | 2022-11-01 | 东华理工大学南昌校区 | Geological map connecting method |
| CN116228998A (en) * | 2023-01-06 | 2023-06-06 | 广东佛山地质工程勘察院 | Automatic generation method of cut geological profile based on geological space database |
Non-Patent Citations (8)
| Title |
|---|
| "V"字型法则定量化研究及其地质意义;郭淑芳 等;黄金;第23卷(第11期);第10-12页 * |
| "V"字型法则约束下基于GIS的地层界线计算机辅助连接;李伟忠 等;中国区域地质(第03期);第300-303页 * |
| "V"字形法则的量化指标及其应用;崔小东;水力发电;第34卷(第11期);第39-40, 44页 * |
| Stochastic Distance Based Geological Boundary Modeling with Realistic Curvilinear Features;Jeff B. Boisvert;Mathematical Geosciences;第45卷;第651-665页 * |
| 剖面法编联地质界线在地质填图中的应用;王辉 等;西北大学学报(自然科学版);第38卷(第06期);第1001-1005页 * |
| 解析几何法在倾斜结构面地质图编制中的应用;刘悦忠;贵州地质;第17卷(第03期);第196-199页 * |
| 走向线在地质工作中的应用;张金英;有色矿冶;第19卷(第01期);第1-4, 16页 * |
| 马鞍桥金矿控矿构造规律及成矿预测;魏宽义 等;西安工程学院学报;第23卷(第01期);第20-24页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117152301A (en) | 2023-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2944757A1 (en) | Interactive well pad plan | |
| WO2023185734A1 (en) | Mine digital model construction method based on three-dimensional point cloud database | |
| CN109933867A (en) | A kind of geotechnical engineering investigation scheme dynamic optimization method | |
| CN109143361B (en) | Method for compiling carbonate rock stratum paleogeology map based on sequence stratigraphy | |
| CN105184864A (en) | Site stratum three-dimensional geological structure model generation method for natural foundation replacement quantities calculation | |
| CN105899749A (en) | Stratigraphic and structural interpretation of deviated and horizontal wellbores | |
| CN103806478A (en) | Method for measuring deformation of underground diaphragm wall | |
| NO322922B1 (en) | Method for determining position uncertainty when drilling wells | |
| CN104809266A (en) | Spline based accurate predicating method for face coal seam occurrence condition | |
| CN101625422B (en) | Method for monitoring time-lapse microgravity oil and gas reservoirs by fixed field method | |
| CN109472866B (en) | A kind of modeling method of big geology Building Information Model | |
| CN112434360B (en) | Three-dimensional modeling and profile mapping method for complex rock-soil body | |
| JP2018004494A (en) | Method for predicting geological boundary surface or fault surface | |
| NO20121440A1 (en) | Rotations from gradient directions | |
| CN112597572A (en) | BIM technology-based pile foundation soil penetration depth determination method | |
| CN110073246B (en) | Improved method relating to quality control | |
| CN107633556B (en) | Method for quantitatively obtaining uncertainty of three-dimensional ore deposit geological model | |
| CN110764142A (en) | Method for assisting seismic data interpretation | |
| CN108256259B (en) | The method and device of seat earth slip casting target area is determined based on Fuzzy Level Analytic Approach | |
| CN106886584B (en) | Underground space development and utilization current situation estimation method based on various urban geographic data | |
| CN117152301B (en) | Geological boundary drawing system based on geological point occurrence and coordinates | |
| CN113376695A (en) | Full waveform inversion method suitable for complex collapse column of coal seam floor | |
| CN114820969B (en) | Three-dimensional geological model construction method | |
| CN113420348B (en) | Method for rapidly drawing uranium mine exploration line profile map | |
| CN113887046B (en) | Coal mine tunnel modeling method based on three-dimensional geologic body |
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 |