CN109885936B - Tracking and identifying method and system for geological boundary in two-dimensional geological section diagram - Google Patents
Tracking and identifying method and system for geological boundary in two-dimensional geological section diagram Download PDFInfo
- Publication number
- CN109885936B CN109885936B CN201910132520.0A CN201910132520A CN109885936B CN 109885936 B CN109885936 B CN 109885936B CN 201910132520 A CN201910132520 A CN 201910132520A CN 109885936 B CN109885936 B CN 109885936B
- Authority
- CN
- China
- Prior art keywords
- geological
- geological boundary
- boundary
- coordinates
- identification code
- 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
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
The invention discloses a method and a system for tracking and identifying geological boundary lines in a two-dimensional geological section chart, which belong to the technical field of civil engineering, wherein the method for tracking and identifying geological boundary lines in the two-dimensional geological section chart comprises the following steps: s1, preparing basic data; s2, standardizing geological information; s3, identifying a drilling datum point; s4, extracting geological boundary lines; s5, geological boundary tracking; a) Comparing the geological boundary with the relative positions of the drill holes; b) Correcting the identification code of the starting point of the geological boundary; c) Correcting the identification code of the geological boundary end point; d) Tracking geological boundary lines; s6, coordinate conversion; s7, identifying attribute information; s8, outputting standard data; by adopting the technical scheme, the method for tracking and integrating the geological boundary in the two-dimensional geological section map and further identifying the geological boundary has the advantages of strong controllability, higher efficiency and capability of meeting the requirement of large-scale application.
Description
Technical Field
The invention belongs to the technical field of civil engineering, and particularly relates to a method and a system for tracking and identifying geological boundary lines in a two-dimensional geological section.
Background
Three-dimensional geologic modeling is a technique for geologic research by combining geologic information with visualization tools in a three-dimensional environment using computer technology.
Methods and techniques for three-dimensional geologic modeling are not yet mature due to the complexity of engineering geologic conditions, limitations in the degree of cognition. Most of the existing three-dimensional geological modeling methods are based on two-dimensional engineering geological data, such as two-dimensional geological sections, and the like, geological boundary lines and attribute information thereof are extracted, further deduction calculation is performed to form a geological interface, and a three-dimensional geological model is finally constructed. The geological boundary line in the two-dimensional geological section is comprehensively judged by a geological engineer in combination with drilling information, regional geological conditions and the like, the traditional database has no geological boundary line data with accurate geological attribute information, the geological boundary lines of the same stratum are not organically connected, and the geological boundary line and the attribute information thereof need to be manually extracted during three-dimensional geological modeling, so that the workload is high and the efficiency is low. The invention provides a method capable of accurately tracking and integrating geological boundary lines in a two-dimensional geological section chart and further identifying the geological boundary lines, which has the advantages of strong controllability and higher efficiency of man-machine interaction operation, and meets the requirement of large-scale application.
Disclosure of Invention
Aiming at the problems existing in the prior art, the invention provides a method and a system for tracking and identifying geological boundary lines in a two-dimensional geological section chart, wherein the method for tracking and integrating the geological boundary lines in the two-dimensional geological section chart and further identifying the geological boundary lines has the advantages of strong controllability, higher efficiency and capability of meeting the requirement of large-scale application.
One of the purposes of the present invention is to provide a method for tracking and identifying geological boundary in a two-dimensional geological section, which at least comprises the following steps:
s1, preparing basic data;
the basic data specifically comprise engineering geological data, including an engineering geological data backup database, an engineering geological section chart in a CAD format, an engineering geological histogram and an engineering geological report;
s2, standardizing geological information;
arranging exploration data in an engineering geological data backup database according to the following sequence: drilling numbers, drilling hole coordinates (x, y, z), stratum layering node coordinates (x, y, z), stratum layer numbers and attribute information;
s3, identifying a drilling datum point;
extracting the drilling number and the orifice elevation in the engineering geological section map in CAD format, establishing drilling reference points, recording the drilling number and the orifice elevation coordinates (alpha) 1 ，β 1 )；
S4, extracting geological boundary lines;
extracting geological boundary lines in the engineering geological section map in CAD format, recording the numbers of the geological boundary lines, and the number of nodes and the coordinates (alpha) of the nodes on the geological boundary lines 2 ，β 2 ) Node coordinates α by node 2 Arranging records from small to large;
s5, geological boundary tracking;
attribute information defining geological boundaries: the stratum information on one side above the geological boundary is attribute information of the geological boundary;
a) Comparing the geological boundary with the relative positions of the drill holes;
for the kth drill hole in the drill hole reference point of step S3, the orifice elevation coordinate (α 1 ，β 1 ) Coordinates (α) of the start point and the end point in each geological boundary node extracted in step S4 2 ，β 2 ) Comparing;
if alpha is 1 And alpha is 2 Identical and beta 1 And beta 2 Different, the start point or the end point of the geological boundary line is indicated to be positioned at the geological layering position below the drilling hole, and the identification code of the geological layering position is marked as 1;
in other cases, indicating that the start point or the end point of the geological boundary is not at the geological stratification position below the drilling hole orifice, and temporarily marking the identification code of the geological boundary;
circulating all the drilling holes, judging the position relation between the drilling holes and the geological boundary, and marking the identification codes of the starting points or the ending points of the drilling holes;
b) Correcting the identification code of the starting point of the geological boundary;
A start point for the ith geological boundary;
if the starting point has no identification code, marking the identification code of the starting point as 0;
if the starting point identification code is 1, searching all other geological boundary starting points, finding out geological boundary overlapped with the starting point coordinates of the ith geological boundary, counting the number m of the geological boundary, and correcting the starting point identification code according to the following rule: if m=1, the start point identification code of the ith geological boundary is not corrected and is still 1; if m >1, assuming that the starting points of m geological boundaries are overlapped at the point O, comparing the slope ratio of a section of straight line (OP, OS) adjacent to the starting point O of the geological boundary (OP, OSQ), wherein the starting point identification code of the geological boundary (OP) where the straight line (OP) with the maximum slope ratio is positioned is not corrected, is still 1, and the starting point identification codes of other geological boundaries (OSQ) are corrected to be 2;
circularly judging the starting points of all geological boundaries, and marking and correcting the identification codes (one of 0,1 and 2) of the geological boundaries;
c) Correcting the identification code of the geological boundary end point;
endpoint for the ith geological boundary:
if the terminal does not have the identification code, marking the terminal identification code as 0;
if the terminal identification code is 1, searching all other geological boundary terminals, finding out geological boundary overlapped with the terminal coordinates of the ith geological boundary, counting the number n of the geological boundary, and correcting the terminal identification code according to the following rule; if n=1, the endpoint identification code of the ith geological boundary is not corrected and is still 1; if n >1, assuming that the end points of n geological boundaries are overlapped at a point C, comparing the slope ratio of a section of straight line (KC, BC) of the geological boundaries (AKC, BC) close to the end point C, wherein the end point identification code of the geological boundary (AKC) where the straight line (KC) with the largest slope ratio is positioned is not corrected, and is still 1, and the end point identification codes of other geological Boundaries (BC) are corrected to be 2;
Circularly judging the end points of all geological boundaries, and marking and correcting the identification codes (one of 0,1 and 2) of the geological boundaries;
d) Tracking geological boundary lines;
the start point and the end point of each geological boundary have identification codes through the steps a), b) and c);
tracking integration rules: according to the origin coordinates alpha of all geological boundary lines 2 Arranging geological boundaries from small to large, and arranging the starting point coordinates alpha of the geological boundaries 2 Minimal initial tracking integration; if the identification code of the end point of one geological boundary line and the start point of the other geological boundary line are 1 and the coordinates are coincident, connecting the two geological boundary lines end to end, and integrating the two geological boundary lines into a new geological boundary line; tracking the end point of the new geological boundary, and comparing and tracking the new geological boundary with the start points of all other non-integrated geological boundaries again according to the rules;
tracking termination rules: the end point identification code of the original existing or newly integrated geological boundary is not 1, or the start point coordinates of all the rest geological boundaries are not coincident;
tracking and integrating all geological boundaries according to the rules to form a plurality of new geological boundaries;
s6, coordinate conversion;
s61, comparing the drilling numbers and establishing stepsThe one-to-one correspondence of the borehole aperture coordinates in the borehole reference point of step S3 to the aperture coordinates in the geological information normalized in step S2, i.e., α 1 And (x, y), beta 1 Correspondence with z;
s62, tracking node coordinates on the integrated geological boundary line in the step S5, and converting the node coordinates into coordinates in standardized geological information according to the corresponding relation;
s7, identifying attribute information;
for the geological boundary line with the transformed coordinates in the step S6, a node with an identification code of 1 is found, the coordinates of the node are compared with the coordinates of the stratum layering nodes in the step S2, and stratum layer numbers and attribute information with coincident coordinates are given to the geological boundary line;
s8, outputting standard data;
outputting geological boundary with attribute information according to a certain format, and using the geological boundary as three-dimensional geological modeling basic data, for example, according to the following sequence: geological boundary serial number, node coordinates (x, y, z), stratum layer number, attribute information.
Another object of the present invention is to provide a system for tracking and identifying geological boundary in two-dimensional geological section, which at least comprises:
a basic data preparation module;
the basic data specifically comprise engineering geological data, including an engineering geological data backup database, an engineering geological section chart in a CAD format, an engineering geological histogram and an engineering geological report;
a geological information standardization module;
arranging exploration data in an engineering geological data backup database according to the following sequence: drilling numbers, drilling hole coordinates (x, y, z), stratum layering node coordinates (x, y, z), stratum layer numbers and attribute information;
A drilling reference point identification module;
extracting the drilling number and the orifice elevation in the engineering geological section map in CAD format, establishing drilling reference points, recording the drilling number and the orifice elevation coordinates (alpha) 1 ，β 1 )；
A geological boundary extraction module;
extracting geological boundary lines in the engineering geological section map in CAD format, recording the numbers of the geological boundary lines, and the number of nodes and the coordinates (alpha) of the nodes on the geological boundary lines 2 ，β 2 ) Node coordinates α by node 2 Arranging records from small to large;
a geological boundary tracking module;
attribute information defining geological boundaries: the stratum information on one side above the geological boundary is attribute information of the geological boundary;
a) Comparing the geological boundary with the relative positions of the drill holes;
for the kth drill hole in the drill hole reference point of step S3, the orifice elevation coordinate (α 1 ，β 1 ) Coordinates (α) of the start point and the end point in each geological boundary node extracted in step S4 2 ，β 2 ) Comparing;
if alpha is 1 And alpha is 2 Identical and beta 1 And beta 2 Different, the start point or the end point of the geological boundary line is indicated to be positioned at the geological layering position below the drilling hole, and the identification code of the geological layering position is marked as 1;
in other cases, indicating that the start point or the end point of the geological boundary is not at the geological stratification position below the drilling hole orifice, and temporarily marking the identification code of the geological boundary;
Circulating all the drilling holes, judging the position relation between the drilling holes and the geological boundary, and marking the identification codes of the starting points or the ending points of the drilling holes;
b) Correcting the identification code of the starting point of the geological boundary;
a start point for the ith geological boundary;
if the starting point has no identification code, marking the identification code of the starting point as 0;
if the starting point identification code is 1, searching all other geological boundary starting points, finding out geological boundary overlapped with the starting point coordinates of the ith geological boundary, counting the number m of the geological boundary, and correcting the starting point identification code according to the following rule: if m=1, the start point identification code of the ith geological boundary is not corrected and is still 1; if m >1, assuming that the starting points of m geological boundaries are overlapped at the point O, comparing the slope ratio of a section of straight line (OP, OS) adjacent to the starting point O of the geological boundary (OP, OSQ), wherein the starting point identification code of the geological boundary (OP) where the straight line (OP) with the maximum slope ratio is positioned is not corrected, is still 1, and the starting point identification codes of other geological boundaries (OSQ) are corrected to be 2;
circularly judging the starting points of all geological boundaries, marking and correcting the identification codes of all geological boundaries, wherein the corrected identification codes are one of 0,1 and 2;
c) Correcting the identification code of the geological boundary end point;
endpoint for the ith geological boundary:
if the terminal does not have the identification code, marking the terminal identification code as 0;
If the terminal identification code is 1, searching all other geological boundary terminals, finding out geological boundary overlapped with the terminal coordinates of the ith geological boundary, counting the number n of the geological boundary, and correcting the terminal identification code according to the following rule; if n=1, the endpoint identification code of the ith geological boundary is not corrected and is still 1; if n >1, assuming that the end points of n geological boundaries are overlapped at a point C, comparing the slope ratio of a section of straight line (KC, BC) of the geological boundaries (AKC, BC) close to the end point C, wherein the end point identification code of the geological boundary (AKC) where the straight line (KC) with the largest slope ratio is positioned is not corrected, and is still 1, and the end point identification codes of other geological Boundaries (BC) are corrected to be 2;
circularly judging the end points of all geological boundaries, marking and correcting the identification codes of all geological boundaries, wherein the corrected identification codes are one of 0,1 and 2;
d) Tracking geological boundary lines;
the start point and the end point of each geological boundary have identification codes through the steps a), b) and c);
tracking integration rules: according to the origin coordinates alpha of all geological boundary lines 2 Arranging geological boundaries from small to large, and arranging the starting point coordinates alpha of the geological boundaries 2 Minimal initial tracking integration; if the identification code of the end point of one geological boundary line and the start point of the other geological boundary line are 1 and the coordinates are coincident, connecting the two geological boundary lines end to end, and integrating the two geological boundary lines into a new geological boundary line; tracking the end point of the new geological boundary, and comparing and tracking the new geological boundary with the start points of all other non-integrated geological boundaries again according to the rules;
Tracking termination rules: the end point identification code of the original existing or newly integrated geological boundary is not 1, or the start point coordinates of all the rest geological boundaries are not coincident;
tracking and integrating all geological boundaries according to the rules to form a plurality of new geological boundaries;
a coordinate conversion module;
s61, comparing the drilling numbers to establish a one-to-one correspondence between the drilling hole coordinates in the drilling reference point in the step S3 and the hole coordinates in the geological information standardized in the step S2, namely alpha 1 And (x, y), beta 1 Correspondence with z;
s62, tracking node coordinates on the integrated geological boundary line in the step S5, and converting the node coordinates into coordinates in standardized geological information according to the corresponding relation;
an attribute information identification module;
for the geological boundary line with the transformed coordinates in the step S6, a node with an identification code of 1 is found, the coordinates of the node are compared with the coordinates of the stratum layering nodes in the step S2, and stratum layer numbers and attribute information with coincident coordinates are given to the geological boundary line;
a standard data output module;
outputting geological boundary with attribute information according to a certain format, and using the geological boundary as three-dimensional geological modeling basic data, for example, according to the following sequence: geological boundary serial number, node coordinates (x, y, z), stratum layer number, attribute information.
It is a further object of the present invention to provide a computer program for implementing a method for tracking and identifying geological boundary in a two-dimensional geological section.
The invention aims at providing an information data processing terminal for realizing a tracking identification method of geological boundary lines in a two-dimensional geological section.
It is a fifth object of the present invention to provide a computer readable storage medium comprising instructions which, when run on a computer, cause the computer to perform a method of tracking and identifying geological boundaries in a two-dimensional geological profile.
In summary, the invention has the advantages and positive effects that:
by adopting the technical scheme, the invention provides a method capable of accurately tracking and integrating the geological boundary lines in the two-dimensional geological section map and further identifying the geological boundary lines, aiming at the problem that the geological boundary lines and the geological attribute information thereof in the two-dimensional engineering geological data and the geological boundary lines of the same stratum are not organically integrated during three-dimensional geological modeling. The method combines the layer-connecting rule of the geological boundary, defines and judges the identification code of the geological boundary, further tracks and integrates the geological boundary, and finally identifies the attribute information of the geological boundary by comparing the identification code with standardized geological information; the number of geological boundaries is effectively reduced, the identification efficiency is greatly improved, the man-machine interaction operation is high in controllability, the large-scale application requirements are met, and a solid foundation is laid for building three-dimensional geological modeling and production and application thereof.
Drawings
FIG. 1 is a flow chart of a preferred embodiment of the present invention;
FIG. 2 is a diagram illustrating the correction of the start code of a geological boundary according to a preferred embodiment of the present invention;
FIG. 3 is a schematic diagram of endpoint identification code correction of geological boundaries in accordance with a preferred embodiment of the present invention;
the figure indicates:
A. b, C, O, P, Q: stratum layering node
AKC, BC, OP, OSQ: geological boundary
Detailed Description
The present invention will be described in further detail with reference to the following examples in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
In the three-dimensional geological modeling process, the existing two-dimensional engineering geological data are fully utilized, but the traditional database has no geological boundary data with accurate geological attribute information, the geological boundary and the geological attribute information thereof in the two-dimensional engineering geological data are manually extracted, and the geological boundary of the same stratum is organically integrated and connected, so that the workload is high and the efficiency is low.
Referring to fig. 1 to 3, a method for tracking and identifying geological boundary in a two-dimensional geological section includes:
(1) Basic data preparation:
engineering geological data: the method mainly comprises an engineering geological data backup library, an engineering geological section chart (CAD format), an engineering geological histogram, an engineering geological report and the like.
(2) Geological information standardization:
arranging exploration data in an engineering geological data backup database according to the following sequence: borehole number, borehole orifice coordinates (x, y, z), stratigraphic layering node coordinates (x, y, z), stratigraphic layering number (i.e., stratigraphic layering node), and attribute information (geological information such as geotechnical name, geologic age, geologic origin, weathering degree, bearing capacity, etc.).
(3) Drilling reference point identification:
extracting the drilling number and the orifice elevation in the engineering geological section map in CAD format, establishing drilling reference points, recording the drilling number and the orifice elevation coordinates (alpha) 1 ，β 1 )；
(4) Geological boundary extraction:
extracting geological boundary lines in the engineering geological section map in CAD format, recording the numbers of the geological boundary lines, and the number of nodes and the coordinates (alpha) of the nodes on the geological boundary lines 2 ，β 2 ) Node coordinates α by node 2 Arranging records from small to large;
(5) Geological boundary tracking:
Attribute information defining geological boundaries: the stratum information on the upper side of the geological boundary is attribute information of the geological boundary. a) Comparing the geological boundary with the relative position of the drilling hole;
for the kth drill hole in the drill hole reference point of step S3, the orifice elevation coordinate (α 1 ，β 1 ) Coordinates (α) of the start point and the end point in each geological boundary node extracted in step S4 2 ，β 2 ) Comparing;
if alpha is 1 And alpha is 2 Identical and beta 1 And beta 2 Different, the start point or the end point of the geological boundary line is indicated to be positioned at the geological layering position below the drilling hole, and the identification code of the geological layering position is marked as 1;
in other cases, the start point or the end point of the geological boundary is indicated not to be at the geological layering position below the drilling hole, and the identification code of the geological boundary is temporarily not marked;
circulating all the drilling holes, judging the position relation between the drilling holes and the geological boundary, and marking the identification codes of the starting point or the end point of the geological boundary;
b) Identification code correction of geological boundary origin
Start point for the ith geological boundary:
if the starting point has no identification code, marking the identification code of the starting point as 0;
if the starting point identification code is 1, searching all other geological boundary starting points, finding out geological boundary overlapped with the starting point coordinates of the ith geological boundary, counting the number m of the geological boundary (including the ith geological boundary), and correcting the starting point identification code according to the following rule: if m=1, the start point identification code of the ith geological boundary is not corrected and is still 1; if m >1, assuming that the starting points of m (m=2) geological boundary lines are overlapped at the point O, comparing the slope ratio of a section of straight line (OP, OS) adjacent to the starting point O of the geological boundary lines (OP, OSQ), wherein the starting point identification code of the geological boundary line (OP) where the straight line (OP) with the maximum slope ratio is positioned is not corrected, and is still 1, and the starting point identification codes of other geological boundary lines (OSQ) are corrected to be 2;
And circularly judging the starting points of all geological boundaries, and marking and correcting the identification codes (one of three types 0,1 and 2).
Illustrating:
it is known that: the O-point coordinates are (. Alpha 4 ，β 4 ) P point coordinates (alpha 5 ，β 5 ) S point coordinates (alpha 6 ，β 6 )
Then: OP slope ratio K OP ＝(β 5 -β 4 )/(|α 5 -α 4 |)
OS slope ratio K OS ＝(β 6 -β 4 )/(|α 6 -α 4 |)
K OP >K OS Therefore, the start point identification code of the geological boundary OP is not corrected, but is still 1, and the start point identification code of the geological boundary OSQ is corrected to 2.
Note that: the ratio is mathematically large and the ratio is negative, the smaller the absolute value, the larger the ratio.
c) Identification code correction of geological boundary end point
Endpoint for the ith geological boundary:
if the terminal does not have the identification code, marking the terminal identification code as 0;
if the terminal identification code is 1, searching all other geological boundary terminals, finding out geological boundary lines overlapped with the terminal coordinates of the ith geological boundary line, counting the number n of the geological boundary lines (including the ith geological boundary line), and correcting the terminal identification code according to the following rule; if n=1, the endpoint identification code of the ith geological boundary is not corrected and is still 1; if n >1, assuming that the end points of n (n=2) geological boundaries overlap at the point C, comparing the slope ratio of a section of straight line (KC, BC) of the geological boundary (AKC, BC) adjacent to the end point C, wherein the end point identification code of the geological boundary (AKC) where the straight line (KC) with the maximum slope ratio is positioned is not corrected, and is still 1, and the end point identification codes of other geological Boundary (BC) are corrected to be 2;
And circularly judging the end points of all geological boundaries, and marking and correcting the identification codes (one of 0,1 and 2) of the geological boundaries.
Illustrating:
it is known that: the C point coordinates are (alpha) 7 ，β 7 ) K Point coordinates (alpha) 8 ，β 8 ) Point B coordinates (alpha 9 ，β 9 )
Then: KC slope ratio K KC ＝(β 8 -β 7 )/(|α 8 -α 7 |)
BC slope ratio K BC ＝(β 9 -β 7 )/(|α 9 -α 7 |)
K KC >K BC Therefore, the end point identification code of the geological boundary AKC is not corrected, but is still 1, and the end point identification code of the geological boundary BC is corrected to 2.
Note that: the ratio is mathematically large and the ratio is negative, the smaller the absolute value, the larger the ratio.
d) Geological boundary tracing
The start point and the end point of each geological boundary have identification codes through the steps a), b) and c).
Tracking integration rules: according to the origin coordinates alpha of all geological boundary lines 2 Arranging geological boundaries from small to large, and arranging the starting point coordinates alpha of the geological boundaries 2 Minimal initial tracking integration; if the identification code of the end point of one geological boundary line and the start point of the other geological boundary line are 1 and the coordinates are coincident, connecting the two geological boundary lines end to end, and integrating the two geological boundary lines into a new geological boundary line; tracking the end point of the new geological boundary, and comparing the end point with the start points of all other non-integrated geological boundaries again according to the rules to track and integrate.
Tracking termination rules: the end point identification code of the original existing or newly integrated geological boundary is not 1, or the start point coordinates of all the rest geological boundaries are not coincident.
And tracking and integrating all geological boundaries according to the rules to form a plurality of new geological boundaries.
Illustrating:
this step has been completed: according to the origin coordinates alpha of all geological boundary lines 2 The geological boundaries are arranged from small to large. Then from the geological boundary origin coordinate alpha 2 Minimal start tracking integration:
the new ID of the ith geological boundary is given as num_1,
if the end point identification code is 0 or 2, the tracking integration is terminated, and the geological boundary is ended;
if the end point identification code is 1, tracking and searching all other geological boundaries, if a certain geological boundary has a starting point of 1 and the coordinates are the same as the starting point coordinates of the ith geological boundary, merging the geological boundary into the ith geological boundary, and tracking the end point identification code of the newly merged geological boundary, if the end point identification code is 1, continuing tracking, and if the end point identification code is 0 or 2, ending tracking;
the remaining geological boundaries are integrated by cyclic tracking.
(6) Coordinate conversion:
a) Comparing the drilling numbers to establish a one-to-one correspondence between the drilling hole coordinates in the drilling reference point in the step S3 and the hole coordinates in the geological information standardized in the step S2, namely alpha 1 And (x, y), beta 1 Correspondence with z;
b) The node coordinates on the geological boundary line which is tracked and integrated in the step (5) are converted into coordinates in standardized geological information according to the corresponding relation;
(7) And (3) identifying attribute information:
and (3) for the geological boundary after the coordinate conversion in the step (6), finding out a node with an identification code of 1, comparing the node coordinate with the stratum layering node coordinate in the step (2), and giving the geological boundary the stratum number and attribute information with the overlapped coordinates.
(8) Standard data output:
outputting geological boundary with attribute information according to a certain format, and using the geological boundary as three-dimensional geological modeling basic data, for example, according to the following sequence: geological boundary number, node coordinates (x, y, z), stratum number (i.e., stratum layered nodes), attribute information (geological information such as geotechnical name, geological age, geological origin, weathering degree, bearing capacity, etc.).
A second preferred embodiment is a system for implementing a method for tracking and identifying geological boundary in a two-dimensional geological section, comprising:
a basic data preparation module:
engineering geological data: the method mainly comprises an engineering geological data backup library, an engineering geological section chart (CAD format), an engineering geological histogram, an engineering geological report and the like.
Geological information standardization module:
Arranging exploration data in an engineering geological data backup database according to the following sequence: borehole number, borehole orifice coordinates (x, y, z), stratigraphic layering node coordinates (x, y, z), stratigraphic layering number (i.e., stratigraphic layering node), and attribute information (geological information such as geotechnical name, geologic age, geologic origin, weathering degree, bearing capacity, etc.).
A drilling datum point identification module:
engineering land of CAD formatThe drilling number and the orifice elevation in the mass section are extracted, a drilling reference point is established, the drilling number is recorded, and the orifice elevation coordinate (alpha 1 ，β 1 )；
A geological boundary extraction module:
extracting geological boundary lines in the engineering geological section map in CAD format, recording the numbers of the geological boundary lines, and the number of nodes and the coordinates (alpha) of the nodes on the geological boundary lines 2 ，β 2 ) Node coordinates α by node 2 Arranging records from small to large;
geological boundary tracking module:
attribute information defining geological boundaries: the stratum information on the upper side of the geological boundary is attribute information of the geological boundary. a) Comparison of geological boundary and relative position of borehole
For the kth drill hole in the drill hole reference point of step S3, the orifice elevation coordinate (α 1 ，β 1 ) Coordinates (α) of the start point and the end point in each geological boundary node extracted in step S4 2 ，β 2 ) Comparing;
if alpha is 1 And alpha is 2 Identical and beta 1 And beta 2 Different, the start point or the end point of the geological boundary line is indicated to be positioned at the geological layering position below the drilling hole, and the identification code of the geological layering position is marked as 1;
in other cases, the start point or the end point of the geological boundary is indicated not to be at the geological layering position below the drilling hole, and the identification code of the geological boundary is temporarily not marked;
circulating all the drilling holes, judging the position relation between the drilling holes and the geological boundary, and marking the identification codes of the starting point or the end point of the geological boundary;
b) Identification code correction of geological boundary origin
Start point for the ith geological boundary:
if the starting point has no identification code, marking the identification code of the starting point as 0;
if the starting point identification code is 1, searching all other geological boundary starting points, finding out geological boundary overlapped with the starting point coordinates of the ith geological boundary, counting the number m of the geological boundary (including the ith geological boundary), and correcting the starting point identification code according to the following rule: if m=1, the start point identification code of the ith geological boundary is not corrected and is still 1; if m >1, assuming that the starting points of m (m=2) geological boundary lines are overlapped at the point O, comparing the slope ratio of a section of straight line (OP, OS) adjacent to the starting point O of the geological boundary lines (OP, OSQ), wherein the starting point identification code of the geological boundary line (OP) where the straight line (OP) with the maximum slope ratio is positioned is not corrected, and is still 1, and the starting point identification codes of other geological boundary lines (OSQ) are corrected to be 2;
And circularly judging the starting points of all geological boundaries, and marking and correcting the identification codes (one of three types 0,1 and 2).
Illustrating:
it is known that: the O-point coordinates are (. Alpha 4 ，β 4 ) P point coordinates (alpha 5 ，β 5 ) S point coordinates (alpha 6 ，β 6 )
Then: OP slope ratio K OP ＝(β 5 -β 4 )/(|α 5 -α 4 |)
OS slope ratio K OS ＝(β 6 -β 4 )/(|α 6 -α 4 |)
K OP >K OS Therefore, the start point identification code of the geological boundary OP is not corrected, but is still 1, and the start point identification code of the geological boundary OSQ is corrected to 2.
Note that: the ratio is mathematically large and the ratio is negative, the smaller the absolute value, the larger the ratio.
c) Identification code correction of geological boundary end point
Endpoint for the ith geological boundary:
if the terminal does not have the identification code, marking the terminal identification code as 0;
if the terminal identification code is 1, searching all other geological boundary terminals, finding out geological boundary lines overlapped with the terminal coordinates of the ith geological boundary line, counting the number n of the geological boundary lines (including the ith geological boundary line), and correcting the terminal identification code according to the following rule; if n=1, the endpoint identification code of the ith geological boundary is not corrected and is still 1; if n >1, assuming that the end points of n (n=2) geological boundaries overlap at the point C, comparing the slope ratio of a section of straight line (KC, BC) of the geological boundary (AKC, BC) adjacent to the end point C, wherein the end point identification code of the geological boundary (AKC) where the straight line (KC) with the maximum slope ratio is positioned is not corrected, and is still 1, and the end point identification codes of other geological Boundary (BC) are corrected to be 2;
And circularly judging the end points of all geological boundaries, and marking and correcting the identification codes (one of 0,1 and 2) of the geological boundaries.
Illustrating:
it is known that: the C point coordinates are (alpha) 7 ，β 7 ) K Point coordinates (alpha) 8 ，β 8 ) Point B coordinates (alpha 9 ，β 9 )
Then: KC slope ratio K KC ＝(β 8 -β 7 )/(|α 8 -α 7 |)
BC slope ratio K BC ＝(β 9 -β 7 )/(|α 9 -α 7 |)
K KC >K BC Therefore, the end point identification code of the geological boundary AKC is not corrected, but is still 1, and the end point identification code of the geological boundary BC is corrected to 2.
Note that: the ratio is mathematically large and the ratio is negative, the smaller the absolute value, the larger the ratio.
d) Geological boundary tracing
The start point and the end point of each geological boundary have identification codes through the steps a), b) and c).
Tracking integration rules: according to the origin coordinates alpha of all geological boundary lines 2 Arranging geological boundaries from small to large, and arranging the starting point coordinates alpha of the geological boundaries 2 Minimal initial tracking integration; if the identification code of the end point of one geological boundary line and the start point of the other geological boundary line are 1 and the coordinates are coincident, connecting the two geological boundary lines end to end, and integrating the two geological boundary lines into a new geological boundary line; tracking the end point of the new geological boundary, and comparing the end point with the start points of all other non-integrated geological boundaries again according to the rules to track and integrate.
Tracking termination rules: the end point identification code of the original existing or newly integrated geological boundary is not 1, or the start point coordinates of all the rest geological boundaries are not coincident.
And tracking and integrating all geological boundaries according to the rules to form a plurality of new geological boundaries.
Illustrating:
this step has been completed: according to the origin coordinates alpha of all geological boundary lines 2 The geological boundaries are arranged from small to large. Then from the geological boundary origin coordinate alpha 2 Minimal start tracking integration:
the new ID of the ith geological boundary is given as num_1,
if the end point identification code is 0 or 2, the tracking integration is terminated, and the geological boundary is ended;
if the end point identification code is 1, tracking and searching all other geological boundaries, if a certain geological boundary has a starting point of 1 and the coordinates are the same as the starting point coordinates of the ith geological boundary, merging the geological boundary into the ith geological boundary, and tracking the end point identification code of the newly merged geological boundary, if the end point identification code is 1, continuing tracking, and if the end point identification code is 0 or 2, ending tracking;
the remaining geological boundaries are integrated by cyclic tracking.
And a coordinate conversion module:
a) Comparing the drilling numbers to establish a one-to-one correspondence between the drilling hole coordinates in the drilling reference point in the step S3 and the hole coordinates in the geological information standardized in the step S2, namely alpha 1 And (x, y), beta 1 Correspondence with z;
b) The node coordinates on the geological boundary line which is tracked and integrated in the step (5) are converted into coordinates in standardized geological information according to the corresponding relation;
attribute information identification module:
and (3) for the geological boundary after the coordinate conversion in the step (6), finding out a node with an identification code of 1, comparing the node coordinate with the stratum layering node coordinate in the step (2), and giving the geological boundary the stratum number and attribute information with the overlapped coordinates.
Standard data output module:
outputting geological boundary with attribute information according to a certain format, and using the geological boundary as three-dimensional geological modeling basic data, for example, according to the following sequence: geological boundary number, node coordinates (x, y, z), stratum number (i.e., stratum layered nodes), attribute information (geological information such as geotechnical name, geological age, geological origin, weathering degree, bearing capacity, etc.).
A third preferred embodiment is a computer program for implementing a method for tracking and identifying a geological boundary in a two-dimensional geological section, the method for tracking and identifying a geological boundary in a two-dimensional geological section comprising the steps of:
(1) Basic data preparation:
engineering geological data: the method mainly comprises an engineering geological data backup library, an engineering geological section chart (CAD format), an engineering geological histogram, an engineering geological report and the like.
(2) Geological information standardization:
arranging exploration data in an engineering geological data backup database according to the following sequence: borehole number, borehole orifice coordinates (x, y, z), stratigraphic layering node coordinates (x, y, z), stratigraphic layering number (i.e., stratigraphic layering node), and attribute information (geological information such as geotechnical name, geologic age, geologic origin, weathering degree, bearing capacity, etc.).
(3) Drilling reference point identification:
extracting the drilling number and the orifice elevation in the engineering geological section map in CAD format, establishing drilling reference points, recording the drilling number and the orifice elevation coordinates (alpha) 1 ，β 1 )；
(4) Geological boundary extraction:
extracting geological boundary lines in the engineering geological section map in CAD format, recording the numbers of the geological boundary lines, and the number of nodes and the coordinates (alpha) of the nodes on the geological boundary lines 2 ，β 2 ) Node coordinates α by node 2 Arranging records from small to large;
(5) Geological boundary tracking:
attribute information defining geological boundaries: the stratum information on the upper side of the geological boundary is attribute information of the geological boundary. a) Comparison of geological boundary and relative position of borehole
For the kth drill hole in the drill hole reference point of step S3, the orifice elevation coordinate (α 1 ，β 1 ) Starting point sum in each geological boundary node extracted in step S4 Coordinates of the endpoint (alpha) 2 ，β 2 ) Comparing;
if alpha is 1 And alpha is 2 Identical and beta 1 And beta 2 Different, the start point or the end point of the geological boundary line is indicated to be positioned at the geological layering position below the drilling hole, and the identification code of the geological layering position is marked as 1;
in other cases, the start point or the end point of the geological boundary is indicated not to be at the geological layering position below the drilling hole, and the identification code of the geological boundary is temporarily not marked;
circulating all the drilling holes, judging the position relation between the drilling holes and the geological boundary, and marking the identification codes of the starting point or the end point of the geological boundary; b) Identification code correction of geological boundary origin
Start point for the ith geological boundary:
if the starting point has no identification code, marking the identification code of the starting point as 0;
if the starting point identification code is 1, searching all other geological boundary starting points, finding out geological boundary overlapped with the starting point coordinates of the ith geological boundary, counting the number m of the geological boundary (including the ith geological boundary), and correcting the starting point identification code according to the following rule: if m=1, the start point identification code of the ith geological boundary is not corrected and is still 1; if m >1, assuming that the starting points of m (m=2) geological boundary lines are overlapped at the point O, comparing the slope ratio of a section of straight line (OP, OS) adjacent to the starting point O of the geological boundary lines (OP, OSQ), wherein the starting point identification code of the geological boundary line (OP) where the straight line (OP) with the maximum slope ratio is positioned is not corrected, and is still 1, and the starting point identification codes of other geological boundary lines (OSQ) are corrected to be 2;
And circularly judging the starting points of all geological boundaries, and marking and correcting the identification codes (one of three types 0,1 and 2).
Illustrating:
it is known that: the O-point coordinates are (. Alpha 4 ，β 4 ) P point coordinates (alpha 5 ，β 5 ) S point coordinates (alpha 6 ，β 6 )
Then: OP slope ratio K OP ＝(β 5 -β 4 )/(|α 5 -α 4 |)
OS slope ratio K OS ＝(β 6 -β 4 )/(|α 6 -α 4 |)
K OP >K OS Therefore, the start point identification code of the geological boundary OP is not corrected, but is still 1, and the start point identification code of the geological boundary OSQ is corrected to 2.
Note that: the ratio is mathematically large and the ratio is negative, the smaller the absolute value, the larger the ratio.
c) Identification code correction of geological boundary end point
Endpoint for the ith geological boundary:
if the terminal does not have the identification code, marking the terminal identification code as 0;
if the terminal identification code is 1, searching all other geological boundary terminals, finding out geological boundary lines overlapped with the terminal coordinates of the ith geological boundary line, counting the number n of the geological boundary lines (including the ith geological boundary line), and correcting the terminal identification code according to the following rule; if n=1, the endpoint identification code of the ith geological boundary is not corrected and is still 1; if n >1, assuming that the end points of n (n=2) geological boundaries overlap at the point C, comparing the slope ratio of a section of straight line (KC, BC) of the geological boundary (AKC, BC) adjacent to the end point C, wherein the end point identification code of the geological boundary (AKC) where the straight line (KC) with the maximum slope ratio is positioned is not corrected, and is still 1, and the end point identification codes of other geological Boundary (BC) are corrected to be 2;
And circularly judging the end points of all geological boundaries, and marking and correcting the identification codes (one of 0,1 and 2) of the geological boundaries.
Illustrating:
it is known that: the C point coordinates are (alpha) 7 ，β 7 ) K Point coordinates (alpha) 8 ，β 8 ) Point B coordinates (alpha 9 ，β 9 )
Then: KC slope ratio K KC ＝(β 8 -β 7 )/(|α 8 -α 7 |)
BC slope ratio K BC ＝(β 9 -β 7 )/(|α 9 -α 7 |)
K KC >K BC Therefore, the end point identification code of the geological boundary AKC is not corrected, but is still 1, and the end point identification code of the geological boundary BC is corrected to 2.
Note that: the ratio is mathematically large and the ratio is negative, the smaller the absolute value, the larger the ratio.
d) Geological boundary tracing
The start point and the end point of each geological boundary have identification codes through the steps a), b) and c).
Tracking integration rules: according to the origin coordinates alpha of all geological boundary lines 2 Arranging geological boundaries from small to large, and arranging the starting point coordinates alpha of the geological boundaries 2 Minimal initial tracking integration; if the identification code of the end point of one geological boundary line and the start point of the other geological boundary line are 1 and the coordinates are coincident, connecting the two geological boundary lines end to end, and integrating the two geological boundary lines into a new geological boundary line; tracking the end point of the new geological boundary, and comparing the end point with the start points of all other non-integrated geological boundaries again according to the rules to track and integrate.
Tracking termination rules: the end point identification code of the original existing or newly integrated geological boundary is not 1, or the start point coordinates of all the rest geological boundaries are not coincident.
And tracking and integrating all geological boundaries according to the rules to form a plurality of new geological boundaries.
Illustrating:
this step has been completed: according to the origin coordinates alpha of all geological boundary lines 2 The geological boundaries are arranged from small to large. Then from the geological boundary origin coordinate alpha 2 Minimal start tracking integration:
the new ID of the ith geological boundary is given as num_1,
if the end point identification code is 0 or 2, the tracking integration is terminated, and the geological boundary is ended;
if the end point identification code is 1, tracking and searching all other geological boundaries, if a certain geological boundary has a starting point of 1 and the coordinates are the same as the starting point coordinates of the ith geological boundary, merging the geological boundary into the ith geological boundary, and tracking the end point identification code of the newly merged geological boundary, if the end point identification code is 1, continuing tracking, and if the end point identification code is 0 or 2, ending tracking;
the remaining geological boundaries are integrated by cyclic tracking.
(6) Coordinate conversion:
a) Comparing the drilling numbers to establish a one-to-one correspondence between the drilling hole coordinates in the drilling reference point in the step S3 and the hole coordinates in the geological information standardized in the step S2, namely alpha 1 And (x, y), beta 1 Correspondence with z;
b) The node coordinates on the geological boundary line which is tracked and integrated in the step (5) are converted into coordinates in standardized geological information according to the corresponding relation;
(7) And (3) identifying attribute information:
and (3) for the geological boundary after the coordinate conversion in the step (6), finding out a node with an identification code of 1, comparing the node coordinate with the stratum layering node coordinate in the step (2), and giving the geological boundary the stratum number and attribute information with the overlapped coordinates.
(8) Standard data output:
outputting geological boundary with attribute information according to a certain format, and using the geological boundary as three-dimensional geological modeling basic data, for example, according to the following sequence: geological boundary number, node coordinates (x, y, z), stratum number (i.e., stratum layered nodes), attribute information (geological information such as geotechnical name, geological age, geological origin, weathering degree, bearing capacity, etc.).
In a fourth preferred embodiment, an information data processing terminal for implementing a method for tracking and identifying geological boundary in a two-dimensional geological section is provided. The method for tracking and identifying the geological boundary in the two-dimensional geological section comprises the following steps:
(1) Basic data preparation:
engineering geological data: the method mainly comprises an engineering geological data backup library, an engineering geological section chart (CAD format), an engineering geological histogram, an engineering geological report and the like.
(2) Geological information standardization:
arranging exploration data in an engineering geological data backup database according to the following sequence: borehole number, borehole orifice coordinates (x, y, z), stratigraphic layering node coordinates (x, y, z), stratigraphic layering number (i.e., stratigraphic layering node), and attribute information (geological information such as geotechnical name, geologic age, geologic origin, weathering degree, bearing capacity, etc.).
(3) Drilling reference point identification:
extracting the drilling number and the orifice elevation in the engineering geological section map in CAD format, establishing drilling reference points, recording the drilling number and the orifice elevation coordinates (alpha) 1 ，β 1 )；
(4) Geological boundary extraction:
extracting geological boundary lines in the engineering geological section map in CAD format, recording the numbers of the geological boundary lines, and the number of nodes and the coordinates (alpha) of the nodes on the geological boundary lines 2 ，β 2 ) Node coordinates α by node 2 Arranging records from small to large;
(5) Geological boundary tracking:
attribute information defining geological boundaries: the stratum information on the upper side of the geological boundary is attribute information of the geological boundary. a) Comparison of geological boundary and relative position of borehole
For the kth drill hole in the drill hole reference point of step S3, the orifice elevation coordinate (α 1 ，β 1 ) Coordinates (α) of the start point and the end point in each geological boundary node extracted in step S4 2 ，β 2 ) Comparing;
if alpha is 1 And alpha is 2 Identical and beta 1 And beta 2 Different, the start point or the end point of the geological boundary line is indicated to be positioned at the geological layering position below the drilling hole, and the identification code of the geological layering position is marked as 1;
in other cases, the start point or the end point of the geological boundary is indicated not to be at the geological layering position below the drilling hole, and the identification code of the geological boundary is temporarily not marked;
circulating all the drilling holes, judging the position relation between the drilling holes and the geological boundary, and marking the identification codes of the starting point or the end point of the geological boundary;
b) Identification code correction of geological boundary origin
Start point for the ith geological boundary:
if the starting point has no identification code, marking the identification code of the starting point as 0;
if the starting point identification code is 1, searching all other geological boundary starting points, finding out geological boundary overlapped with the starting point coordinates of the ith geological boundary, counting the number m of the geological boundary (including the ith geological boundary), and correcting the starting point identification code according to the following rule: if m=1, the start point identification code of the ith geological boundary is not corrected and is still 1; if m >1, assuming that the starting points of m (m=2) geological boundary lines are overlapped at the point O, comparing the slope ratio of a section of straight line (OP, OS) adjacent to the starting point O of the geological boundary lines (OP, OSQ), wherein the starting point identification code of the geological boundary line (OP) where the straight line (OP) with the maximum slope ratio is positioned is not corrected, and is still 1, and the starting point identification codes of other geological boundary lines (OSQ) are corrected to be 2;
And circularly judging the starting points of all geological boundaries, and marking and correcting the identification codes (one of three types 0,1 and 2).
Illustrating:
it is known that: the O-point coordinates are (. Alpha 4 ，β 4 ) P point coordinates (alpha 5 ，β 5 ) S point coordinates (alpha 6 ，β 6 )
Then: OP slope ratio K OP ＝(β 5 -β 4 )/(|α 5 -α 4 |)
OS slope ratio K OS ＝(β 6 -β 4 )/(|α 6 -α 4 |)
K OP >K OS Therefore, the start point identification code of the geological boundary OP is not corrected, but is still 1, and the start point identification code of the geological boundary OSQ is corrected to 2.
Note that: the ratio is mathematically large and the ratio is negative, the smaller the absolute value, the larger the ratio.
c) Identification code correction of geological boundary end point
Endpoint for the ith geological boundary:
if the terminal does not have the identification code, marking the terminal identification code as 0;
if the terminal identification code is 1, searching all other geological boundary terminals, finding out geological boundary lines overlapped with the terminal coordinates of the ith geological boundary line, counting the number n of the geological boundary lines (including the ith geological boundary line), and correcting the terminal identification code according to the following rule; if n=1, the endpoint identification code of the ith geological boundary is not corrected and is still 1; if n >1, assuming that the end points of n (n=2) geological boundaries overlap at the point C, comparing the slope ratio of a section of straight line (KC, BC) of the geological boundary (AKC, BC) adjacent to the end point C, wherein the end point identification code of the geological boundary (AKC) where the straight line (KC) with the maximum slope ratio is positioned is not corrected, and is still 1, and the end point identification codes of other geological Boundary (BC) are corrected to be 2;
And circularly judging the end points of all geological boundaries, and marking and correcting the identification codes (one of 0,1 and 2) of the geological boundaries.
Illustrating:
it is known that: the C point coordinates are (alpha) 7 ，β 7 ) K Point coordinates (alpha) 8 ，β 8 ) Point B coordinates (alpha 9 ，β 9 )
Then: KC slope ratio K KC ＝(β 8 -β 7 )/(|α 8 -α 7 |)
BC slope ratio K BC ＝(β 9 -β 7 )/(|α 9 -α 7 |)
K KC >K BC Therefore, the end point identification code of the geological boundary AKC is not corrected, but is still 1, and the end point identification code of the geological boundary BC is corrected to 2.
Note that: the ratio is mathematically large and the ratio is negative, the smaller the absolute value, the larger the ratio.
d) Geological boundary tracing
The start point and the end point of each geological boundary have identification codes through the steps a), b) and c).
Tracking integration rules: according to the origin coordinates alpha of all geological boundary lines 2 Arranging geological boundaries from small to large, and arranging the starting point coordinates alpha of the geological boundaries 2 Minimal initial tracking integration; if the identification code of the end point of one geological boundary line and the start point of the other geological boundary line are 1 and the coordinates are coincident, connecting the two geological boundary lines end to end, and integrating the two geological boundary lines into a new geological boundary line; tracking the end point of the new geological boundary, and comparing the end point with the start points of all other non-integrated geological boundaries again according to the rules to track and integrate.
Tracking termination rules: the end point identification code of the original existing or newly integrated geological boundary is not 1, or the start point coordinates of all the rest geological boundaries are not coincident.
And tracking and integrating all geological boundaries according to the rules to form a plurality of new geological boundaries.
Illustrating:
this step has been completed: according to the origin coordinates alpha of all geological boundary lines 2 The geological boundaries are arranged from small to large. Then from the geological boundary origin coordinate alpha 2 Minimal start tracking integration:
the new ID of the ith geological boundary is given as num_1,
if the end point identification code is 0 or 2, the tracking integration is terminated, and the geological boundary is ended;
if the end point identification code is 1, tracking and searching all other geological boundaries, if a certain geological boundary has a starting point of 1 and the coordinates are the same as the starting point coordinates of the ith geological boundary, merging the geological boundary into the ith geological boundary, and tracking the end point identification code of the newly merged geological boundary, if the end point identification code is 1, continuing tracking, and if the end point identification code is 0 or 2, ending tracking;
the remaining geological boundaries are integrated by cyclic tracking.
(6) Coordinate conversion:
a) Comparing the drilling numbers to establish a one-to-one correspondence between the drilling hole coordinates in the drilling reference point in the step S3 and the hole coordinates in the geological information standardized in the step S2, namely alpha 1 And (x, y), beta 1 Correspondence with z;
b) The node coordinates on the geological boundary line which is tracked and integrated in the step (5) are converted into coordinates in standardized geological information according to the corresponding relation;
(7) And (3) identifying attribute information:
and (3) for the geological boundary after the coordinate conversion in the step (6), finding out a node with an identification code of 1, comparing the node coordinate with the stratum layering node coordinate in the step (2), and giving the geological boundary the stratum number and attribute information with the overlapped coordinates.
(8) Standard data output:
outputting geological boundary with attribute information according to a certain format, and using the geological boundary as three-dimensional geological modeling basic data, for example, according to the following sequence: geological boundary number, node coordinates (x, y, z), stratum number (i.e., stratum layered nodes), attribute information (geological information such as geotechnical name, geological age, geological origin, weathering degree, bearing capacity, etc.).
A fifth preferred embodiment is a computer-readable storage medium comprising instructions that, when executed on a computer, cause the computer to perform a method of tracking and identifying a geological boundary in a two-dimensional geological profile, the method of tracking and identifying a geological boundary in a two-dimensional geological profile comprising the steps of: the method for tracking and identifying the geological boundary in the two-dimensional geological section comprises the following steps:
(1) Basic data preparation:
engineering geological data: the method mainly comprises an engineering geological data backup library, an engineering geological section chart (CAD format), an engineering geological histogram, an engineering geological report and the like.
(2) Geological information standardization:
arranging exploration data in an engineering geological data backup database according to the following sequence: borehole number, borehole orifice coordinates (x, y, z), stratigraphic layering node coordinates (x, y, z), stratigraphic layering number (i.e., stratigraphic layering node), and attribute information (geological information such as geotechnical name, geologic age, geologic origin, weathering degree, bearing capacity, etc.).
(3) Drilling reference point identification:
extracting the drilling number and the orifice elevation in the engineering geological section map in CAD format, establishing drilling reference points, recording the drilling number and the orifice elevation coordinates (alpha) 1 ，β 1 )；
(4) Geological boundary extraction:
extracting geological boundary lines in the engineering geological section map in CAD format, recording the numbers of the geological boundary lines, and the number of nodes and the coordinates (alpha) of the nodes on the geological boundary lines 2 ，β 2 ) Node coordinates α by node 2 Arranging records from small to large;
(5) Geological boundary tracking:
attribute information defining geological boundaries: the stratum information on the upper side of the geological boundary is attribute information of the geological boundary. a) Comparison of geological boundary and relative position of borehole
For the kth drill hole in the drill hole reference point of step S3, the orifice elevation coordinate (α 1 ，β 1 ) Coordinates (α) of the start point and the end point in each geological boundary node extracted in step S4 2 ，β 2 ) Comparing;
if alpha is 1 And alpha is 2 Identical and beta 1 And beta 2 Different, the start point or the end point of the geological boundary line is indicated to be positioned at the geological layering position below the drilling hole, and the identification code of the geological layering position is marked as 1;
in other cases, the start point or the end point of the geological boundary is indicated not to be at the geological layering position below the drilling hole, and the identification code of the geological boundary is temporarily not marked;
circulating all the drilling holes, judging the position relation between the drilling holes and the geological boundary, and marking the identification codes of the starting point or the end point of the geological boundary;
b) Identification code correction of geological boundary origin
Start point for the ith geological boundary:
if the starting point has no identification code, marking the identification code of the starting point as 0;
if the starting point identification code is 1, searching all other geological boundary starting points, finding out geological boundary overlapped with the starting point coordinates of the ith geological boundary, counting the number m of the geological boundary (including the ith geological boundary), and correcting the starting point identification code according to the following rule: if m=1, the start point identification code of the ith geological boundary is not corrected and is still 1; if m >1, assuming that the starting points of m (m=2) geological boundary lines are overlapped at the point O, comparing the slope ratio of a section of straight line (OP, OS) adjacent to the starting point O of the geological boundary lines (OP, OSQ), wherein the starting point identification code of the geological boundary line (OP) where the straight line (OP) with the maximum slope ratio is positioned is not corrected, and is still 1, and the starting point identification codes of other geological boundary lines (OSQ) are corrected to be 2;
And circularly judging the starting points of all geological boundaries, and marking and correcting the identification codes (one of three types 0,1 and 2).
Illustrating:
it is known that: the O-point coordinates are (. Alpha 4 ，β 4 ) P point coordinates (alpha 5 ，β 5 ) S point coordinates (alpha 6 ，β 6 )
Then: OP slope ratio K OP ＝(β 5 -β 4 )/(|α 5 -α 4 |)
OS slope ratio K OS ＝(β 6 -β 4 )/(|α 6 -α 4 |)
K OP >K OS Therefore, the start point identification code of the geological boundary OP is not corrected, but is still 1, and the start point identification code of the geological boundary OSQ is corrected to 2.
Note that: the ratio is mathematically large and the ratio is negative, the smaller the absolute value, the larger the ratio.
c) Identification code correction of geological boundary end point
Endpoint for the ith geological boundary:
if the terminal does not have the identification code, marking the terminal identification code as 0;
if the terminal identification code is 1, searching all other geological boundary terminals, finding out geological boundary lines overlapped with the terminal coordinates of the ith geological boundary line, counting the number n of the geological boundary lines (including the ith geological boundary line), and correcting the terminal identification code according to the following rule; if n=1, the endpoint identification code of the ith geological boundary is not corrected and is still 1; if n >1, assuming that the end points of n (n=2) geological boundaries overlap at the point C, comparing the slope ratio of a section of straight line (KC, BC) of the geological boundary (AKC, BC) adjacent to the end point C, wherein the end point identification code of the geological boundary (AKC) where the straight line (KC) with the maximum slope ratio is positioned is not corrected, and is still 1, and the end point identification codes of other geological Boundary (BC) are corrected to be 2;
And circularly judging the end points of all geological boundaries, and marking and correcting the identification codes (one of 0,1 and 2) of the geological boundaries.
Illustrating:
it is known that: the C point coordinates are (alpha) 7 ，β 7 ) K Point coordinates (alpha) 8 ，β 8 ) Point B coordinates (alpha 9 ，β 9 )
Then: KC slope ratio K KC ＝(β 8 -β 7 )/(|α 8 -α 7 |)
BC slope ratio K BC ＝(β 9 -β 7 )/(|α 9 -α 7 |)
K KC >K BC Therefore, the end point identification code of the geological boundary AKC is not corrected, but is still 1, and the end point identification code of the geological boundary BC is corrected to 2.
Note that: the ratio is mathematically large and the ratio is negative, the smaller the absolute value, the larger the ratio.
d) Geological boundary tracing
The start point and the end point of each geological boundary have identification codes through the steps a), b) and c).
Tracking integration rules: according to the origin coordinates alpha of all geological boundary lines 2 Arranging geological boundaries from small to large, and arranging the starting point coordinates alpha of the geological boundaries 2 Minimal initial tracking integration; if the identification code of the end point of one geological boundary line and the start point of the other geological boundary line are 1 and the coordinates are coincident, connecting the two geological boundary lines end to end, and integrating the two geological boundary lines into a new geological boundary line; tracking the end point of the new geological boundary, and comparing the end point with the start points of all other non-integrated geological boundaries again according to the rules to track and integrate.
Tracking termination rules: the end point identification code of the original existing or newly integrated geological boundary is not 1, or the start point coordinates of all the rest geological boundaries are not coincident.
And tracking and integrating all geological boundaries according to the rules to form a plurality of new geological boundaries.
Illustrating:
this step has been completed: according to the origin coordinates alpha of all geological boundary lines 2 The geological boundaries are arranged from small to large. Then from the geological boundary origin coordinate alpha 2 Minimal start tracking integration:
the new ID of the ith geological boundary is given as num_1,
if the end point identification code is 0 or 2, the tracking integration is terminated, and the geological boundary is ended;
if the end point identification code is 1, tracking and searching all other geological boundaries, if a certain geological boundary has a starting point of 1 and the coordinates are the same as the starting point coordinates of the ith geological boundary, merging the geological boundary into the ith geological boundary, and tracking the end point identification code of the newly merged geological boundary, if the end point identification code is 1, continuing tracking, and if the end point identification code is 0 or 2, ending tracking;
the remaining geological boundaries are integrated by cyclic tracking.
(6) Coordinate conversion:
a) Comparing the drilling numbers to establish a one-to-one correspondence between the drilling hole coordinates in the drilling reference point in the step S3 and the hole coordinates in the geological information standardized in the step S2, namely alpha 1 And (x, y), beta 1 Correspondence with z;
b) The node coordinates on the geological boundary line which is tracked and integrated in the step (5) are converted into coordinates in standardized geological information according to the corresponding relation;
(7) And (3) identifying attribute information:
and (3) for the geological boundary after the coordinate conversion in the step (6), finding out a node with an identification code of 1, comparing the node coordinate with the stratum layering node coordinate in the step (2), and giving the geological boundary the stratum number and attribute information with the overlapped coordinates.
(8) Standard data output:
outputting geological boundary with attribute information according to a certain format, and using the geological boundary as three-dimensional geological modeling basic data, for example, according to the following sequence: geological boundary number, node coordinates (x, y, z), stratum number (i.e., stratum layered nodes), attribute information (geological information such as geotechnical name, geological age, geological origin, weathering degree, bearing capacity, etc.).
In the above embodiments, it may be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When used in whole or in part, is implemented in the form of a computer program product comprising one or more computer instructions. When loaded or executed on a computer, produces a flow or function in accordance with embodiments of the present invention, in whole or in part. The computer may be a general purpose computer, a special purpose computer, a computer network, or other programmable apparatus. The computer instructions may be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another computer-readable storage medium, for example, the computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center by a wired (e.g., coaxial cable, fiber optic, digital Subscriber Line (DSL), or wireless (e.g., infrared, wireless, microwave, etc.) means. The computer readable storage medium may be any available medium that can be accessed by a computer or a data storage device such as a server, data center, etc. that contains an integration of one or more available media. The usable medium may be a magnetic medium (e.g., floppy Disk, hard Disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium (e.g., solid State Disk (SSD)), etc.
The foregoing description of the preferred embodiments of the invention is not intended to be limiting, but rather is intended to cover all modifications, equivalents, and alternatives falling within the spirit and principles of the invention.
Claims (4)
1. A tracking and identifying method of geological boundary in a two-dimensional geological section diagram is characterized by comprising the following steps of: at least comprises the following steps
S1, preparing basic data;
the basic data specifically comprise engineering geological data, including an engineering geological data backup database, an engineering geological section chart in a CAD format, an engineering geological histogram and an engineering geological report;
s2, standardizing geological information;
arranging exploration data in an engineering geological data backup database according to the following sequence: drilling numbers, drilling hole coordinates (x, y, z), stratum layering node coordinates (x, y, z), stratum layer numbers and attribute information;
s3, identifying a drilling datum point;
extracting the drilling number and the orifice elevation in the engineering geological section map in CAD format, establishing drilling reference points, recording the drilling number and the orifice elevation coordinates (alpha) 1 ，β 1 )；
S4, extracting geological boundary lines;
extracting geological boundary lines in the engineering geological section map in CAD format, recording the numbers of the geological boundary lines, and the number of nodes and the coordinates (alpha) of the nodes on the geological boundary lines 2 ，β 2 ) Node coordinates α by node 2 Arranging records from small to large;
s5, geological boundary tracking;
attribute information defining geological boundaries: the stratum information on one side above the geological boundary is attribute information of the geological boundary;
a) Comparing the geological boundary with the relative positions of the drill holes;
for the kth drill hole in the drill hole reference point of step S3, the orifice elevation coordinate (α 1 ，β 1 ) Coordinates (α) of the start point and the end point in each geological boundary node extracted in step S4 2 ，β 2 ) Comparing;
if alpha is 1 And alpha is 2 Identical and beta 1 And beta 2 Different, the start point or the end point of the geological boundary line is indicated to be positioned at the geological layering position below the drilling hole, and the identification code of the geological layering position is marked as 1;
in other cases, the start point or the end point of the geological boundary is indicated not to be at the geological layering position below the drilling hole, and the identification code of the geological boundary is temporarily not marked;
circulating all the drilling holes, judging the position relation between the drilling holes and the geological boundary, and marking the identification codes of the starting point or the end point of the geological boundary;
b) Correcting the identification code of the starting point of the geological boundary;
a start point for the ith geological boundary;
if the starting point has no identification code, marking the identification code of the starting point as 0;
if the starting point identification code is 1, searching all other geological boundary starting points, finding out geological boundary overlapped with the starting point coordinates of the ith geological boundary, counting the number m of the geological boundary, and correcting the starting point identification code according to the following rule: if m=1, the start point identification code of the ith geological boundary is not corrected and is still 1; if m >1, assuming that the starting points of m geological boundaries are overlapped at the point O, comparing the slope ratio of a section of straight line (OP, OS) adjacent to the starting point O of the geological boundary (OP, OSQ), wherein the starting point identification code of the geological boundary (OP) where the straight line (OP) with the maximum slope ratio is positioned is not corrected, is still 1, and the starting point identification codes of other geological boundaries (OSQ) are corrected to be 2;
Circularly judging the starting points of all geological boundaries, and marking a correction starting point identification code which is one of 0,1 and 2;
c) Correcting the identification code of the geological boundary end point;
endpoint for the ith geological boundary:
if the terminal does not have the identification code, marking the terminal identification code as 0;
if the terminal identification code is 1, searching all other geological boundary terminals, finding out geological boundary overlapped with the terminal coordinates of the ith geological boundary, counting the number n of the geological boundary, and correcting the terminal identification code according to the following rule; if n=1, the endpoint identification code of the ith geological boundary is not corrected and is still 1; if n >1, assuming that the end points of n geological boundaries are overlapped at a point C, comparing the slope ratio of a section of straight line (KC, BC) of the geological boundaries (AKC, BC) close to the end point C, wherein the end point identification code of the geological boundary (AKC) where the straight line (KC) with the largest slope ratio is positioned is not corrected, and is still 1, and the end point identification codes of other geological Boundaries (BC) are corrected to be 2;
circularly judging the end points of all geological boundaries, and marking a corrected end point identification code which is one of 0,1 and 2;
d) Tracking geological boundary lines;
the start point and the end point of each geological boundary have identification codes through the steps a), b) and c);
Tracking integration rules: according to the origin coordinates alpha of all geological boundary lines 2 Arranging geological boundaries from small to large, and arranging the starting point coordinates alpha of the geological boundaries 2 Minimal initial tracking integration; if the identification code of the end point of one geological boundary line and the start point of the other geological boundary line are 1 and the coordinates are coincident, connecting the two geological boundary lines end to end, and integrating the two geological boundary lines into a new geological boundary line; tracking the end point of the new geological boundary, and comparing and tracking the new geological boundary with the start points of all other non-integrated geological boundaries again according to the rules;
tracking termination rules: the end point identification code of the original existing or newly integrated geological boundary is not 1, or the start point coordinates of all the rest geological boundaries are not coincident;
tracking and integrating all geological boundaries according to the rules to form a plurality of new geological boundaries;
s6, coordinate conversion;
s61, comparing the drilling numbers to establish a one-to-one correspondence between the drilling hole coordinates in the drilling reference point in the step S3 and the hole coordinates in the geological information standardized in the step S2, namely alpha 1 And (x, y), beta 1 Correspondence with z;
s62, tracking node coordinates on the integrated geological boundary line in the step S5, and converting the node coordinates into coordinates in standardized geological information according to the corresponding relation;
S7, identifying attribute information;
for the geological boundary line with the transformed coordinates in the step S6, a node with an identification code of 1 is found, the coordinates of the node are compared with the coordinates of the stratum layering nodes in the step S2, and stratum layer numbers and attribute information with coincident coordinates are given to the geological boundary line;
s8, outputting standard data;
outputting geological boundary with attribute information according to a certain format, and using the geological boundary as three-dimensional geological modeling basic data, for example, according to the following sequence: geological boundary serial number, node coordinates (x, y, z), stratum layer number, attribute information.
2. A system for tracking and identifying geological boundary in a two-dimensional geological section, which is characterized in that: at least comprises:
a basic data preparation module;
the basic data specifically comprise engineering geological data, including an engineering geological data backup database, an engineering geological section chart in a CAD format, an engineering geological histogram and an engineering geological report;
a geological information standardization module;
arranging exploration data in an engineering geological data backup database according to the following sequence: drilling numbers, drilling hole coordinates (x, y, z), stratum layering node coordinates (x, y, z), stratum layer numbers and attribute information;
a drilling reference point identification module;
Extracting the drilling number and the orifice elevation in the engineering geological section map in CAD format, establishing drilling reference points, recording the drilling number and the orifice elevation coordinates (alpha) 1 ，β 1 )；
A geological boundary extraction module;
extracting geological boundary lines in the engineering geological section map in CAD format, recording the numbers of the geological boundary lines, and the number of nodes and the coordinates (alpha) of the nodes on the geological boundary lines 2 ，β 2 ) Node coordinates α by node 2 Arranging records from small to large;
a geological boundary tracking module;
attribute information defining geological boundaries: the stratum information on one side above the geological boundary is attribute information of the geological boundary;
a) Comparing the geological boundary with the relative positions of the drill holes;
for the kth drill hole in the drill hole reference point of step S3, the orifice elevation coordinate (α 1 ，β 1 ) Coordinates (α) of the start point and the end point in each geological boundary node extracted in step S4 2 ，β 2 ) Comparing;
if alpha is 1 And alpha is 2 Identical and beta 1 And beta 2 Different, the start point or the end point of the geological boundary line is indicated to be positioned at the geological layering position below the drilling hole, and the identification code of the geological layering position is marked as 1;
in other cases, the start point or the end point of the geological boundary is indicated not to be at the geological layering position below the drilling hole, and the identification code of the geological boundary is temporarily not marked;
circulating all the drilling holes, judging the position relation between the drilling holes and the geological boundary, and marking the identification codes of the starting point or the end point of the geological boundary;
b) Correcting the identification code of the starting point of the geological boundary;
a start point for the ith geological boundary;
if the starting point has no identification code, marking the identification code of the starting point as 0;
if the starting point identification code is 1, searching all other geological boundary starting points, finding out geological boundary overlapped with the starting point coordinates of the ith geological boundary, counting the number m of the geological boundary, and correcting the starting point identification code according to the following rule: if m=1, the start point identification code of the ith geological boundary is not corrected and is still 1; if m >1, assuming that the starting points of m geological boundaries are overlapped at the point O, comparing the slope ratio of a section of straight line (OP, OS) adjacent to the starting point O of the geological boundary (OP, OSQ), wherein the starting point identification code of the geological boundary (OP) where the straight line (OP) with the maximum slope ratio is positioned is not corrected, is still 1, and the starting point identification codes of other geological boundaries (OSQ) are corrected to be 2;
circularly judging the starting points of all geological boundaries, and marking a correction starting point identification code which is one of 0,1 and 2;
c) Correcting the identification code of the geological boundary end point;
endpoint for the ith geological boundary:
if the terminal does not have the identification code, marking the terminal identification code as 0;
if the terminal identification code is 1, searching all other geological boundary terminals, finding out geological boundary overlapped with the terminal coordinates of the ith geological boundary, counting the number n of the geological boundary, and correcting the terminal identification code according to the following rule; if n=1, the endpoint identification code of the ith geological boundary is not corrected and is still 1; if n >1, assuming that the end points of n geological boundaries are overlapped at a point C, comparing the slope ratio of a section of straight line (KC, BC) of the geological boundaries (AKC, BC) close to the end point C, wherein the end point identification code of the geological boundary (AKC) where the straight line (KC) with the largest slope ratio is positioned is not corrected, and is still 1, and the end point identification codes of other geological Boundaries (BC) are corrected to be 2;
Circularly judging the end points of all geological boundaries, and marking a corrected end point identification code which is one of 0,1 and 2;
d) Tracking geological boundary lines;
the start point and the end point of each geological boundary have identification codes through the steps a), b) and c);
tracking integration rules: according to the origin coordinates alpha of all geological boundary lines 2 Arranging geological boundaries from small to large, and arranging the starting point coordinates alpha of the geological boundaries 2 Minimal start trackingIntegrating; if the identification code of the end point of one geological boundary line and the start point of the other geological boundary line are 1 and the coordinates are coincident, connecting the two geological boundary lines end to end, and integrating the two geological boundary lines into a new geological boundary line; tracking the end point of the new geological boundary, and comparing and tracking the new geological boundary with the start points of all other non-integrated geological boundaries again according to the rules;
tracking termination rules: the end point identification code of the original existing or newly integrated geological boundary is not 1, or the start point coordinates of all the rest geological boundaries are not coincident;
tracking and integrating all geological boundaries according to the rules to form a plurality of new geological boundaries;
a coordinate conversion module;
s61, comparing the drilling numbers to establish a one-to-one correspondence between the drilling hole coordinates in the drilling reference point in the step S3 and the hole coordinates in the geological information standardized in the step S2, namely alpha 1 And (x, y), beta 1 Correspondence with z;
s62, tracking node coordinates on the integrated geological boundary line in the step S5, and converting the node coordinates into coordinates in standardized geological information according to the corresponding relation;
an attribute information identification module;
for the geological boundary line with the transformed coordinates in the step S6, a node with an identification code of 1 is found, the coordinates of the node are compared with the coordinates of the stratum layering nodes in the step S2, and stratum layer numbers and attribute information with coincident coordinates are given to the geological boundary line;
a standard data output module;
outputting geological boundary with attribute information according to a certain format, and using the geological boundary as three-dimensional geological modeling basic data, for example, according to the following sequence: geological boundary serial number, node coordinates (x, y, z), stratum layer number, attribute information.
3. An information data processing terminal for realizing the tracking and identifying method of geological boundary in the two-dimensional geological section of claim 1.
4. A computer readable storage medium comprising instructions which, when run on a computer, cause the computer to perform the tracking identification method of geological boundaries in a two-dimensional geological section as claimed in claim 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910132520.0A CN109885936B (en) | 2019-02-22 | 2019-02-22 | Tracking and identifying method and system for geological boundary in two-dimensional geological section diagram |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910132520.0A CN109885936B (en) | 2019-02-22 | 2019-02-22 | Tracking and identifying method and system for geological boundary in two-dimensional geological section diagram |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109885936A CN109885936A (en) | 2019-06-14 |
| CN109885936B true CN109885936B (en) | 2023-07-14 |
Family
ID=66929000
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910132520.0A Active CN109885936B (en) | 2019-02-22 | 2019-02-22 | Tracking and identifying method and system for geological boundary in two-dimensional geological section diagram |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109885936B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1773554A (en) * | 2004-11-09 | 2006-05-17 | 中国石油天然气集团公司 | A two-dimensional closed geological structure model building method |
| JP2006154046A (en) * | 2004-11-26 | 2006-06-15 | Katsuhiko Mima | Land data management system and network data collation system |
| CN105957146A (en) * | 2016-04-29 | 2016-09-21 | 铁道第三勘察设计院集团有限公司 | Linear engineering three-dimensional geological modeling method |
| CN106777391A (en) * | 2017-02-21 | 2017-05-31 | 河海大学 | Geologic section modeling method based on drill hole information and knowledge reasoning technology |
| CN107194966A (en) * | 2017-05-17 | 2017-09-22 | 中国铁路设计集团有限公司 | A kind of discrimination method of two-dimentional geological section geological boundary and its attribute information |
| CN107909627A (en) * | 2017-11-17 | 2018-04-13 | 中国铁路设计集团有限公司 | Geology two dimension based on common layer concept connects layer method |
| CN108109203A (en) * | 2017-11-14 | 2018-06-01 | 中国市政工程中南设计研究总院有限公司 | Based on the geological boundary cloud data extracting method and system under BIM environment |
-
2019
- 2019-02-22 CN CN201910132520.0A patent/CN109885936B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1773554A (en) * | 2004-11-09 | 2006-05-17 | 中国石油天然气集团公司 | A two-dimensional closed geological structure model building method |
| JP2006154046A (en) * | 2004-11-26 | 2006-06-15 | Katsuhiko Mima | Land data management system and network data collation system |
| CN105957146A (en) * | 2016-04-29 | 2016-09-21 | 铁道第三勘察设计院集团有限公司 | Linear engineering three-dimensional geological modeling method |
| CN106777391A (en) * | 2017-02-21 | 2017-05-31 | 河海大学 | Geologic section modeling method based on drill hole information and knowledge reasoning technology |
| CN107194966A (en) * | 2017-05-17 | 2017-09-22 | 中国铁路设计集团有限公司 | A kind of discrimination method of two-dimentional geological section geological boundary and its attribute information |
| CN108109203A (en) * | 2017-11-14 | 2018-06-01 | 中国市政工程中南设计研究总院有限公司 | Based on the geological boundary cloud data extracting method and system under BIM environment |
| CN107909627A (en) * | 2017-11-17 | 2018-04-13 | 中国铁路设计集团有限公司 | Geology two dimension based on common layer concept connects layer method |
Non-Patent Citations (1)
| Title |
|---|
| 基于交叉折剖面的高精度三维地质模型快速构建方法研究;屈红刚等;《北京大学学报(自然科学版)》;20081120(第06期);正文第1-2章 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109885936A (en) | 2019-06-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104504047B (en) | A kind of SOLID MINERAL RESOURCES reserve estimate system based on dual-layer data storehouse | |
| CN111563103B (en) | Method and system for detecting data blood relationship | |
| CN103969682A (en) | Drilling and logging data depth matching method and system | |
| CN109271697B (en) | System and method for processing interrupt data in AutoCAD based on VBA, and computer program | |
| CN105183811A (en) | Informatization data collection system for exploration drilling catalog | |
| CN107992426A (en) | A kind of software error localization method excavated based on Frequent tree mining and processing unit | |
| CN113190910B (en) | Slice analysis system and method based on BIM three-dimensional geological tunnel model | |
| CN111429471A (en) | Geological disaster information management system and method | |
| CN109885936B (en) | Tracking and identifying method and system for geological boundary in two-dimensional geological section diagram | |
| CN107194966B (en) | Identification method for geological boundary and attribute information of two-dimensional geological section map | |
| CN104462288B (en) | A kind of similarity of paths analysis method and system | |
| CN107013207B (en) | Method for acquiring underground well position information of complex fault block oil reservoir sand body | |
| CN113987639A (en) | Intelligent three-dimensional cable addressing analysis system and method based on subway BIM model | |
| CN110956693B (en) | Three-dimensional geological modeling method | |
| CN103049583B (en) | The conversion method of earth monitor data and translation interface | |
| CN114123190A (en) | Method and device for determining target region to which ammeter belongs, electronic equipment and storage medium | |
| CN108694744B (en) | Geological structure surface form comprehensive interpretation method based on multiple outcrops | |
| CN107992590B (en) | Big data system beneficial to information comparison | |
| CN111783373A (en) | Topology analysis method of PSCAD simulation model | |
| CN114971255B (en) | Automatic parameter performance detection system and method | |
| CN107832348B (en) | Method for processing network data stream from intelligent acquisition terminal to cloud | |
| CN111752940B (en) | Sand body communication relation data acquisition method and device | |
| CN109272042B (en) | Small-scale crack modeling method, device, equipment and medium based on PR model | |
| CN117632773B (en) | Data interaction anomaly detection method and detection port based on computer terminal | |
| CN113418511B (en) | Digital mapping system and mapping method for water conservancy and hydropower engineering geology |
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 |