CN116228998A - Automatic generation method of cut geological profile based on geological space database - Google Patents

Abstract

The application relates to an automatic generation method of a cut geological profile based on a geological space database, which comprises the following steps: establishing a geological map space database; modeling is carried out on the geological map space database, and the geological map space database comprises attribute models of boundary points, association models of the boundary points and corresponding boundary lines, projection point corresponding relation models of the boundary points and map cutting section views, left and right geological body attribute models corresponding to the positions of the boundary points and graphic parameter models; acquiring section lines drawn by a user in a graphic work area and based on a geological map; based on the section lines, performing cutting section calculation to obtain parameters corresponding to each model; and obtaining drawing parameters of the cut geological section according to the parameters corresponding to the models, and automatically drawing by utilizing the drawing parameters to obtain the cut geological section. The method and the device solve the problem that in the prior art, for areas without drilling data, a corresponding three-dimensional geological model cannot be established, and then a corresponding map-cut geological section map cannot be obtained.

Description

Automatic generation method of cut geological profile based on geological space database

Technical Field

The application relates to the technical field of geological exploration, in particular to an automatic generation method of a cut geological section map based on a geological space database.

Background

The map cut geological section is a geological section in a specified direction which is formed by a projection method according to a certain scale on a topographic geological map and according to geographic and geological elements, and the map cut geological section can intuitively reflect the spatial distribution characteristics of stratum, lithology and structure in a map area. The map cut geological section has important significance in mineral exploration and geology research.

With the continuous advancement of geological mineral work informatization in recent years, the application of three-dimensional geological modeling technology has become more and more popular, and the mode of making a graph cut section based on a three-dimensional geological model is also receiving more and more attention. However, actually measured drilling data is required for building the three-dimensional geologic model, and for areas without drilling data, such as villages or mountains, the corresponding three-dimensional geologic model cannot be built, and then the corresponding map cut geologic profile cannot be obtained. Improvements are therefore urgently needed.

Disclosure of Invention

In order to solve the problem that in the prior art, for areas without drilling data, a corresponding three-dimensional geological model cannot be established and then a corresponding cut geological section map cannot be obtained, the application provides an automatic generation method of the cut geological section map based on a geological space database.

The automatic generation method of the cut geological section map based on the geological space database adopts the following technical scheme:

an automatic generation method of a cut geological profile based on a geological space database, comprising the following steps:

establishing a geological map space database;

modeling is carried out on the geological map space database, and the geological map space database comprises attribute models of boundary points, association models of the boundary points and corresponding boundary lines, projection point corresponding relation models of the boundary points and map cutting section views, left and right geological body attribute models corresponding to the positions of the boundary points and graphic parameter models;

acquiring section lines drawn by a user in a graphic work area and based on a geological map;

based on the section lines, performing cutting section calculation to obtain parameters corresponding to each model;

obtaining drawing parameters of the cut geological section according to the parameters corresponding to the models, wherein the drawing parameters comprise: the times of stratum, layer sequence, contact relation, distribution and extension conditions of new and old stratum, anticline and syncline conditions, axial plane shape, pivot state, turning end shape, extending direction of folds, length-width ratio of folds and formation times of folds;

and automatically drawing by using the drawing parameters to obtain a cut geological section.

Preferably, the geological map space database comprises element types, namely, a topographic contour line, a geological body, a geological boundary line, a fault and a producing state;

the method further comprises the steps of:

setting a data layer and corresponding attribute data items based on a geological map, and collecting corresponding data; the data layer and the corresponding attribute data items comprise: a contour map layer, wherein the corresponding attribute data item is an elevation value; the geologic body map layer corresponds to the attribute data items such as geologic age, geologic symbols and lithology description; the geological boundary map layer is provided with corresponding attribute data items which are trends and dip angles; the fault layer is provided with a tendency and an inclination angle corresponding to the attribute data item; a yield layer, wherein the corresponding attribute data items are trends and dip angles;

and acquiring parameters corresponding to the model based on the acquired data of the data layer and the corresponding attribute data item.

Preferably, the calculating the cutting profile based on the section line, to obtain parameters corresponding to each model includes:

performing intersection analysis of the cutting section line and the contour line of the terrain to obtain the elevation value of the intersection point, simultaneously obtaining the elevation values of the cutting line end point and the turning point, and writing the elevation value into the attribute model of the boundary point; wherein, the attribute model of the boundary point comprises: category code number, coordinate X, coordinate Y, elevation, primitive ID, accumulated length, label, stratum sequence, and association ID;

and/or

Performing intersection analysis of the cutting section line and the geologic body to obtain intersection point position information of a boundary line of the geologic body, through which the cutting line passes; writing an attribute model of the boundary point of the geological body through calculation; the attribute model of the boundary point of the geologic body comprises: category code number, coordinate X, coordinate Y, elevation, primitive ID, cumulative length, label, type, remark, formation sequence, association ID;

and/or

Performing intersection analysis of the cutting section line and the fault, acquiring intersection point position information of the fault boundary line penetrated by the cutting line, and writing an attribute model of the intersection point through calculation; the attribute model of the intersection point comprises: category code, serial number, coordinate X, coordinate Y, elevation, primitive ID, dip, true dip, cumulative length, label, type of fault, remark, association ID.

Preferably, when the intersection analysis of the cutting section line and the geologic body is performed, the method further comprises: judging whether an intersection point of a boundary line of the geological body penetrated by the cutting line is a fault boundary point or not; if the fault boundary point is the fault boundary point, acquiring a fault name, putting the fault name into a type data item, and then acquiring inclination and true dip angle information of the fault boundary according to the boundary where the fault is located; if the inclination and true dip angle information exists, writing the inclination and true dip angle information into an inclination and true dip angle data item; similarly, if the geological boundary point is the geological boundary point, the contact relation of the geological boundary, the geological body geological code and the stratum sequence information are written into the data item, and if the inclination and true dip angle information exists, the contact relation, the geological body geological code and the stratum sequence information are also written into the inclination and true dip angle data item; if no trend and dip angle information exists, obtaining the trend of the occurrence meeting the conditions, writing the true dip angle value into the trend and true dip angle data item; and searching a yield closest to the intersection point, and if the yield point falls in the stratum corresponding to the intersection point, obtaining the yield which meets the condition.

Preferably, the anticline and syncline conditions of the stratum are judged by the following methods:

assigning sequence numbers as new and old stratum relation sequence codes; wherein, the older stratum of the age is, the larger the sequence number value is;

analyzing stratum crossing the cutting section line to obtain sequence numbers of each stratum;

judging the nucleus and the wing through the mutation condition of the stratum sequence number, and further judging the anticline and syncline conditions: if the stratum sequence number is changed from small to large and is changed into small, judging that the stratum is anticline according to the new principle of the old wing of the nucleus; if the stratum sequence number is changed from large to small and is larger, the stratum sequence number is judged to be syncline according to the principle that the new wing of the nucleus is old.

Preferably, by drawing a graph cut section line and adopting a space analysis intersecting technology, calculating the intersection point of the graph cut section line and a geological body arc section and a geological boundary line; and obtaining the geological unit attributes of the geological body around the intersection point, thereby obtaining the geological times of the boundary around the boundary point, the stratum sequence and the contact relation of the boundary where the boundary point is located.

Preferably, the distribution and extension conditions of the new stratum and the old stratum are automatically judged by the following method:

sequentially segmenting the intersection points of the graph cut section line and the geological boundary line at turning points of the graph cut section line, and sequencing the intersection points in sequence from left to right and from top to bottom, so that the spatial distribution condition of new and old strata is obtained through the attribute of the intersection points; the stratum sequence of the new stratum and the old stratum is obtained by the geologist according to the stratum table of the China regional year or the international geological year representation and sequencing the geological conditions of the working area; assigning a sequence number as a new and old stratum relation sequence code, wherein the older stratum is the older, and the numerical value of the assigned sequence number is larger;

calculating the relative occurrence of the intersection point by using the set buffer zone, and obtaining occurrence representing the geological boundary where the intersection point is located by a nearest distance method, thereby obtaining occurrence attributes and obtaining trend and dip angle values of new and old strata; or according to the correlation model of the boundary points and the corresponding boundary lines, if the corresponding boundary lines have the trend and inclination angle attributes, uniformly converting the trend and inclination angle of the boundary lines into numerical values, and then calculating the apparent inclination angle to obtain the stratum extension condition, wherein if the apparent inclination angle is 90 degrees, the stratum output state is vertical; if the apparent dip angle is 0 degrees, the output state of the stratum is horizontal; if the apparent dip angle epsilon (0,90), the formation output state is inclined towards the right end of the sectional view; if the apparent dip angle e (-90,0), the formation is tilted toward the left end of the profile.

Preferably, the axial plane occurrence information is obtained by:

judging the relation between the inclination of the two wings and the inclination angle, and if the inclination of the two wings is opposite and the inclination angles are nearly equal, indicating that the axial surfaces are vertical; if the inclination angles of the two wings are not equal, the inclination of the axial plane is indicated.

Preferably, the pivot state is obtained by:

judging the trend state of the two wings, and if the trend directions of the two wings are parallel, the hinge is horizontal; if the two wings meet or are curved in an arc shape, the hinge is inclined; wherein, the tip of the intersection curve of the two wings and the same stratum boundary line is the tilting direction of the anticline hinge and faces the lifting direction of the anticline hinge.

Preferably, the form of the turning end is judged by the tight degree of the folds, namely, the turning end is tightly closed as the stratum distance of the two wings is closer, and conversely, the turning end is wider and more gentle.

The distance between two points is calculated through the adjacent intersection points of the section line and the crossed geological boundary line, and the conditions of the width, the smoothness and the like can be known through a series of calculation, the judgment of the dip angle and stratum age, the calculation and the judgment of the form and the comprehensive result, so that the graph drawing of the section view is convenient.

Preferably, the point, line and plane space database of the cross section is established while generating the graph cut cross section; the automatic drawing by using the drawing parameters to obtain the cut geological section comprises the following steps:

calculating and processing the section layering boundary line, namely firstly calculating the visual inclination angle and the true inclination angle, then drawing a graphic sketch, firstly projecting each intersection point of a section line and a topographic line on a geological map onto the section map according to a given scale, and forming a topographic section curve; then, each intersection point of the section line and the geological boundary line on the geological map is projected onto a section curve of the topographic map; layering boundaries are drawn according to formation trends and dip angles near each point.

Preferably, the automatically mapping by using the mapping parameters to obtain the map-cut geological section further comprises:

obtaining complete boundary boundaries forming the sectional view according to the calculated sectional view sketch, wherein the boundary boundaries comprise vertical lines of projection points of starting and ending endpoints of the sectional view, stratum dividing lines, topographic curves and lower base lines of the sectional view; through space analysis, the automatic line approaching, automatic line cutting, redundant line deleting, automatic calculation and treatment of back oblique line, oblique line and water bottom boundary line are adopted, and then the line is converted into an arc segment; automatically checking and processing topological relations among arc segments, and automatically forming cross-section plane primitives of each geological unit according to the enclosed closed arc segments; then, corresponding graphic elements are automatically assigned according to the geological map related attributes, and the establishment of a profile map attribute library is completed; the color of the original corresponding stratum is given to the primitive as the display parameter of the primitive, thereby completing the coloring work of the section; finally, according to the relevant attribute of the projection point, judging and selecting a proper position, and automatically marking the position including stratum code, fracture name, water system name and azimuth angle finishing of the section view.

In summary, the present application includes at least one of the following beneficial technical effects:

1. according to the method, on the premise that no drilling data exists, the existing geological data is fully utilized, modeling is conducted based on a geological map space database, and related parameters of an attribute model of a boundary point, a correlation model of the boundary point and a corresponding boundary, a projection point corresponding relation model of the boundary point and a map-cut section, a left geological attribute model and a right geological attribute model corresponding to the position of the boundary point and a graphic parameter model are obtained, so that drawing parameters are obtained, and the method comprises the following steps: the formation times, layer sequence, contact relation, distribution and extension conditions of new and old formations, anticline and syncline conditions, axial surface appearance, pivot state, turning end shape, extending direction of folds, length ratio of folds, formation times of folds, and finally, automatically drawing by using the drawing parameters to obtain a cut geological section. According to the method and the device, the map cut geological section map is directly obtained based on the geological map of the plane, and the problem that in the prior art, for areas without drilling data, a corresponding three-dimensional geological model cannot be built, and then the corresponding map cut geological section map cannot be obtained is solved. In addition, the technical scheme of this application convenient operation, the practicality is strong, can improve drawing efficiency and drawing accuracy moreover.

2. In the application, the sequence number is used as a new stratum relation sequence code and an old stratum relation sequence code; then analyzing the stratum intersected with the cutting section line to obtain the sequence number of each stratum; then judging the nucleus and the wing through the mutation condition of the stratum sequence number, and further judging the anticline and syncline conditions: if the stratum sequence number is changed from small to large and is changed into small, judging that the stratum is anticline according to the new principle of the old wing of the nucleus; if the stratum sequence number is changed from large to small and is larger, judging to be syncline according to the principle that the new wing of the core is old; compared with the prior art, the method for judging the core and the wings in a symmetrical mode, and further judging the anticline and the syncline, the accuracy of anticline and syncline judgment is greatly improved, and the accuracy of automatic drawing is further improved.

Drawings

FIG. 1 is a method flow diagram of one embodiment of the present application.

FIG. 2 is a comparison of a cut geological profile (bottom) obtained using the method of the present application with a cut geological profile (top) obtained using the prior art;

FIG. 3 is a schematic view of the attribute parameters of various points of a map cut geological section and map cut geological section projection obtained using the method of the present application.

FIG. 4 is a schematic drawing of mapping parameters of a map cut geological profile of the present application.

Detailed Description

The present application is described in further detail below in conjunction with fig. 1-4.

The embodiment of the application discloses an automatic generation method of a cut geological section based on a geological space database, as shown in fig. 1, comprising the following steps:

s1, establishing a geological map space database; specifically, on the premise of no drilling data, the existing geological data can be fully utilized, and a geological map space database containing elements such as topographic contour lines, geologic bodies, geological boundaries, faults, occurrence and the like is established by adopting GIS software;

the method further comprises the steps of:

s11, setting a data layer and corresponding attribute data items based on a geological map, and collecting corresponding data; the data layer and the corresponding attribute data items comprise: a contour map layer, wherein the corresponding attribute data item is an elevation value; the geologic body map layer corresponds to the attribute data items such as geologic age, geologic symbols and lithology description; the geological boundary map layer is characterized in that corresponding attribute data items are trend and dip angle (the content can be optional or not, and trend and dip angle data in the occurrence map layer can be adopted); a fault layer, wherein the corresponding attribute data items are trend and dip angle (the content can be optional or not, and trend and dip angle data in a zone layer can be adopted); a yield layer, wherein the corresponding attribute data items are trends and dip angles;

s12, acquiring parameters corresponding to the model based on the acquired data of the data layer and the corresponding attribute data item;

s2, modeling is carried out on the geological map space database, wherein the geological map space database comprises attribute models of boundary points, association models of the boundary points and corresponding boundary lines, projection point corresponding relation models of the boundary points and map cutting section views, left and right geological body attribute models corresponding to positions of the boundary points and graphic parameter models; the modeling described herein is not modeling in a conventional sense, but rather is a data item to be acquired to create a cut geological profile after drawing a section line on the geological map

S3, acquiring section lines drawn by a user in the graphic work area and based on a geological map, wherein the section lines can be multi-point broken lines;

s4, based on the section lines, performing cutting section calculation to obtain parameters corresponding to the models; comprising the following steps:

performing intersection analysis of the cutting section line and the contour line of the terrain to obtain the elevation value of the intersection point, simultaneously (according to adjacent analysis and extension inference) obtaining the elevation values of the cutting line end point and the turning point, and writing into the attribute model of the boundary point; wherein, the attribute model of the boundary point comprises: category code number, coordinate X, coordinate Y, elevation, primitive ID, accumulated length, label (geological symbol or geological name required to be labeled in section), type, remark (descriptive information such as lithology of the geologic body is recorded), stratum sequence (new and old sequence of stratum age is recorded), and association ID;

and/or

Performing intersection analysis of the cutting section line and the geologic body to obtain intersection point position information of a boundary line of the geologic body, through which the cutting line passes; writing an attribute model of the boundary point of the geological body through calculation; the attribute model of the boundary point of the geologic body comprises: category code number, coordinate X, coordinate Y, elevation, primitive ID (tendency, true dip angle, visual dip angle, section line azimuth angle), accumulated length, label (geological symbol or geological name required to be labeled by section diagram), type, remark (record descriptive information of lithology of geologic body, etc.), stratum sequence (record stratum age new and old sequence number), and association ID;

and/or

Performing intersection analysis of the cutting section line and the fault, acquiring intersection point position information of the fault boundary line penetrated by the cutting line, and writing an attribute model of the intersection point through calculation; the attribute model of the intersection point comprises: class code (the class code can be set as follows: 0-section cutting line control point, 1-elevation control point, 2-nearby ground object marking point, 3-section cutting line and fault intersection point, 4-geologic body demarcation point; the class code is 3 here), serial number, coordinate X, coordinate Y, elevation, primitive ID, tendency, true dip angle, accumulated length, label (fault name or code number required to record the section drawing), fault type (fault type name recorded), remark and associated ID;

when the intersecting analysis of the cutting section line and the geologic body is carried out, the method further comprises the following steps: judging whether an intersection point of a boundary line of the geological body penetrated by the cutting line is a fault boundary point or not; if the fault boundary point is the fault boundary point, acquiring a fault name, putting the fault name into a type data item, and then acquiring inclination and true dip angle information of the fault boundary according to the boundary where the fault is located; if the inclination and true dip angle information exists, writing the inclination and true dip angle information into an inclination and true dip angle data item; similarly, if the geological boundary point is the geological boundary point, the contact relation of the geological boundary, the geological body geological code and the stratum sequence information are written into the data item, and if the inclination and true dip angle information exists, the contact relation, the geological body geological code and the stratum sequence information are also written into the inclination and true dip angle data item; if no trend and dip angle information exists, obtaining the trend of the occurrence meeting the conditions, writing the true dip angle value into the trend and true dip angle data item; wherein, the closest occurrence of the intersection is searched from occurrence (such as formation occurrence, formation inversion occurrence, horizontal occurrence, volcanic formation occurrence, invasion contact occurrence, etc.), and if the occurrence falls in the formation corresponding to the intersection, the occurrence is a satisfactory occurrence.

It should be noted that because the geologic database filling tendency is often done in two ways, one is a specific data value and the other is an interval code. In order to be able to perform the calculation, the trend data items are here collectively numerical. The conversion method comprises the following steps: "S-E" -135, "E" -180, "N" -0.0, "W" -270, "NEE" -70.0, "NNW" -340.0, "NW" -315.0, "NWW" -290.0, "SE" -135.0, "SEE" -115.0, "SSW" -200.0, "SW" -135.0, "SWW" -250.0, "SSE" -160.0, "NNE" -20.0.

S5, drawing parameters of a cut geological section are obtained according to the parameters corresponding to the models, as shown in FIG. 4, and the method comprises the following steps: the times of stratum, layer sequence, contact relation, distribution and extension conditions of new and old stratum, anticline and syncline conditions, axial plane shape, pivot state, turning end shape, extending direction of folds, length-width ratio of folds and formation times of folds; the drawing parameters are directly obtained or calculated through corresponding data items in the limit point attribute model; the age of forming the folds can be obtained through new and old stratum relations, the length-width ratio of the folds can be set according to actual conditions, and the extending direction of the folds can be obtained according to inclination angles.

Specifically, the anticline and syncline conditions of the stratum are judged by the following methods:

s51, assigning sequence numbers as new and old stratum relation sequence codes; wherein, the older stratum of the age is, the larger the sequence number value is; wherein, the new and old stratum sequence can be confirmed by geologist in advance according to stratum (geologic age) table of China regional ages or international geologic age representation;

s52, analyzing stratum crossing the cutting section line to obtain sequence numbers of each stratum;

s53, judging the nucleus and the wing through the mutation condition of the stratum sequence number, and further judging the anticline and syncline conditions: if the stratum sequence number is changed from small to large and is changed into small, judging that the stratum is anticline according to the new principle of the old wing of the nucleus; if the stratum sequence number is changed from large to small and is larger, the stratum sequence number is judged to be syncline according to the principle that the new wing of the nucleus is old.

Optionally, in order to accurately obtain the intersection point of the graph cut section line and the geological boundary line, and rapidly obtain geological units and related attribute contents on the left and right sides of the boundary line point through the position of the intersection point, so as to draw an accurate section view, and calculate the intersection point of the graph cut section line, the geological body arc section and the geological boundary line by drawing the graph cut section line and adopting a space analysis intersection technology; and obtaining the geological unit attributes of the geological body around the intersection point, thereby obtaining the geological times of the boundary around the boundary point, the stratum sequence and the contact relation of the boundary where the boundary point is located.

Optionally, the distribution and extension conditions of the new stratum and the old stratum are automatically judged by the following method:

sequentially segmenting the intersection points of the graph cut section line and the geological boundary line at turning points of the graph cut section line, and sequencing the intersection points in sequence from left to right and from top to bottom, so that the spatial distribution condition of new and old strata is obtained through the attribute of the intersection points; the stratum sequence of the new stratum and the old stratum is obtained by the geologist according to the stratum table of the China regional year or the international geological year representation and sequencing the geological conditions of the working area; assigning a sequence number as a new and old stratum relation sequence code, wherein the older stratum is the older, and the numerical value of the assigned sequence number is larger;

calculating the relative occurrence of the intersection point by using the set buffer zone (without limiting the distance), and obtaining occurrence of geological boundary lines where the intersection point is located (such as formation occurrence) by a nearest distance method, thereby obtaining occurrence attributes and obtaining trend and dip angle values of new and old formations; or according to the correlation model of the boundary points and the corresponding boundary lines, if the corresponding boundary lines have the trend and inclination angle attributes, uniformly converting the trend and inclination angle of the boundary lines into numerical values (corresponding to the numerical value B in the table below), and calculating the magnitude of the apparent inclination angle to obtain the stratum extension condition, wherein if the apparent inclination angle is 90 degrees, the output state of the stratum is vertical; if the apparent dip angle is 0 degrees, the output state of the stratum is horizontal; if the apparent dip angle epsilon (0,90), the formation output state is inclined towards the right end of the sectional view; if the apparent dip angle e (-90,0), the formation is tilted toward the left end of the profile.

Specifically, the apparent inclination angle calculation formula is: tan (a/180 pi) =tan (B/180 pi) ×cos (V/180 pi); wherein, the angle between the A-view inclination angle, the B-true inclination angle, the V-section direction and the stratum (or fault) trend, V = section azimuth angle-stratum trend; (the apparent dip angle A epsilon (-90, 90) calculated from the above formula, the absolute value represents the magnitude of the apparent dip angle;

the tilting mode is mainly determined by V:

if V epsilon (-90, 90), cos (V/180. Times. PI) >0, tilt toward the right end of the cross-section;

if V epsilon (-180, -90) U (90, 180), cos (V/180. Pi) <0, tilt to the left end of the cross section;

if V e { -90,90}, cos (V/180 pi) =0, horizontal, the direction of the section cut line is consistent with the strike of the stratum or fault; but when b=90, tan (B/180 pi) = infinity, where there is always a=90, the formation or fault stands up, independent of V.

Therefore, the above 4 cases can be separated from the value of A, and the stratum extension situation can be obtained.

A takes on the value	tan(A/180*PI)	Output status	Determining factor
				90	∞	Vertical stand	B＝90
0	0	Horizontal level	V∈{-90,90}
				∈(0,90)	>0	Inclined to the right end of the sectional view	V∈(-90,90),
∈(-90,0)	<0	Inclined to the left end of the sectional view	V∈(-180,-90)U(90,180)

Optionally, the axial plane occurrence information is obtained by: judging the relation between the inclination of the two wings and the inclination angle, and if the inclination of the two wings is opposite and the inclination angles are nearly equal, indicating that the axial surfaces are vertical; if the inclination angles of the two wings are not equal, the inclination of the axial plane is indicated.

According to the above calculation and discrimination results, the form of the wrinkles can be obtained, and the wrinkles can be distinguished according to the axial plane shape and the two wing shape:

erect wrinkles: the axial surface is nearly vertical, the inclination of the two wings is opposite, and the inclination angles are nearly equal;

oblique fold: the axial surface is inclined, the two wings are inclined towards the same direction, and the stratum of one wing is inverted;

lying fold: the axial plane is nearly horizontal, one wing stratum is normal, and the other wing stratum is inverted;

inverted fold: the axial plane is inclined, two wings are inclined, the rock strata of the two wings tend to be the same, the inclination angles are equal or unequal, the sequence of one wing rock stratum is normal, and the sequence of the other wing rock stratum is inverted. The inclined shaft surfaces of the inclined and inverted wrinkles are consistent with the tendency of the slow wings (the slow wings are relatively gentle wings of the two wings, and the slow wings are judged by viewing the inclination angle), so that the drawing can be performed according to the rule.

Optionally, the pivot state is obtained by: judging the trend state of the two wings, and if the trend directions of the two wings are parallel, the hinge is horizontal; if the two wings meet or are curved in an arc shape, the hinge is inclined; wherein, the tip of the intersection curve of the two wings and the same stratum boundary line is the tilting direction of the anticline hinge and faces the lifting direction of the anticline hinge.

The specific description is as follows-horizontal fold: the hinge extends approximately horizontally, and the two wing rock strata run approximately in parallel and are symmetrically distributed; and (5) tilting wrinkles: the hinge is inclined towards one end, and the trend of the rock stratum with two wings is arc-shaped. For anticlockwise upward, the surrounding tip points in the direction of tipping; for an upward incline, the surrounding opening points in the direction of the hinge's incline.

In addition, the form of the turning end can be judged by the tight degree of the folds, namely, the turning end is tightly closed when the stratum distance of the two wings is closer, and conversely, the turning end is wider and is gentle. The distance between two points is calculated through the adjacent intersection points of the section line and the crossed geological boundary line, and the conditions of the width and the smoothness are known through a series of calculation, the inclination angle, the discrimination of stratum age, the calculation and discrimination of the form (such as anticline, syncline, axes, axial planes, two wings and the like) and the comprehensive result, so that the graph drawing of the sectional view is convenient.

And S6, automatically drawing by using the drawing parameters to obtain a cut geological section.

Creating a graph cut section and simultaneously establishing a point, line and plane space database of the section; the automatic drawing by using the drawing parameters to obtain the cut geological section comprises the following steps:

s61, creating a graph cut cross section and simultaneously establishing a point, line and plane space database of the cross section; the automatic drawing by using the drawing parameters to obtain the cut geological section comprises the following steps:

calculating and processing the section layering boundary line, namely firstly calculating the visual inclination angle and the true inclination angle, then drawing a graphic sketch (a sketch without filling colors), firstly projecting each intersection point of a section line on a geological map and a topographic line (i.e. a contour line) onto the section map according to a given scale to form a topographic section curve; then, each intersection point of the section line on the geological map and the geological boundary line (the boundary line of the stratum, the non-integrated line, the fault line and the like) is projected onto the section curve of the topographic map; the layering boundary is drawn according to the stratum tendency and the dip angle near each point (when the profile and the trend are obliquely crossed, the layering boundary is drawn according to the apparent dip angle of the profile direction).

If the quality problem exists in the graph generated by calculation, such as calculation errors caused by inaccuracy of original data or calculation methods of computer software, such as a buffer method, the method can obtain all the shapes in the buffer, then judge whether the obtained shapes are suitable and utilized according to the requirements, a nearest distance method and the like, and the obtained data are inaccurate due to incomplete comprehensive consideration of certain characteristics in calculation, the data are corrected by an editing tool provided by a geological map space database, and after corrected drawing parameters are obtained (the parameters are stored in intersection point attribute data items calculated before), the section layering boundary line is recalculated and processed.

S62, obtaining complete boundary boundaries which form the sectional view according to the calculated sectional view sketch, wherein the boundary boundaries comprise vertical lines (elevation lines) of projection points of starting and ending endpoints of the sectional view, stratum dividing lines, terrain curves and lower bottom lines (horizontal bottom lines) of the sectional view; through space analysis, the automatic line approaching, automatic line cutting, redundant line deleting, automatic calculation and treatment of back oblique line, oblique line and water bottom boundary line are adopted, and then the line is converted into an arc segment; automatically checking and processing topological relations among arc segments, and automatically forming cross-section plane primitives of each geological unit according to the enclosed closed arc segments (reducing human intervention); then, according to the geological map related attributes (including geological unit attributes such as geological code, geological body name, geological body related description and the like, and also including parameter characteristic attributes of geological body primitives, geological body surface colors, filling patterns and the like, the contents are obtained and assigned to attribute data items of intersection points of graph section lines and boundary lines in the process of calculating the graph section map) are automatically assigned to corresponding primitives, and the establishment of a section map attribute library is rapidly completed; the color of the original corresponding stratum is given to the primitive as the display parameter of the primitive, and the coloring work of the section view is completed rapidly; (the generated sectional view automatically establishes the topological relation of the sectional view according to a space topological processing method, ensures the consistency of lines and arc segments and the consistency of primitive topologies), and finally judges and selects proper positions according to the related attribute of the projection points, and automatically marks the positions including stratum codes (or names), broken names, (river and the like) water system names and azimuth finishing of the sectional view.

The inventor designs and realizes a tool set integrating automatic generation, interactive editing, data management and output printing of geological profiles based on the technology contents of the application by adopting C++ development language and applying a geological map space database based on the technologies such as a MapGIS platform, geographical information and the like. FIG. 2 is a comparison of a cut geological profile (bottom) obtained using the method of the present application with a cut geological profile (top) obtained using the prior art; FIG. 3 is a schematic view of the attribute parameters of various points of a map cut geological section and map cut geological section projection obtained using the method of the present application. Through application inspection, the tool set is efficient and practical, the result drawing is standard and attractive, the actual application requirements of the geological structure complex region can be met, and the expected target is achieved.

Claims (10)

1. The automatic generation method of the cut geological section based on the geological space database is characterized by comprising the following steps:

establishing a geological map space database;

modeling is carried out on the geological map space database, and the geological map space database comprises attribute models of boundary points, association models of the boundary points and corresponding boundary lines, projection point corresponding relation models of the boundary points and map cutting section views, left and right geological body attribute models corresponding to the positions of the boundary points and graphic parameter models;

acquiring section lines drawn by a user in a graphic work area and based on a geological map;

based on the section lines, performing cutting section calculation to obtain parameters corresponding to each model;

obtaining drawing parameters of the cut geological section according to the parameters corresponding to the models, wherein the drawing parameters comprise: the times of stratum, layer sequence, contact relation, distribution and extension conditions of new and old stratum, anticline and syncline conditions, axial plane shape, pivot state, turning end shape, extending direction of folds, length-width ratio of folds and formation times of folds;

and automatically drawing by using the drawing parameters to obtain a cut geological section.

2. The automatic generation method of a cut geological profile based on a geological space database according to claim 1, wherein: the geological map space database comprises element types, namely a topographic contour line, a geologic body, a geological boundary line, a fault and a yield;

the method further comprises the steps of:

setting a data layer and corresponding attribute data items based on a geological map, and collecting corresponding data; the data layer and the corresponding attribute data items comprise: a contour map layer, wherein the corresponding attribute data item is an elevation value; the geologic body map layer corresponds to the attribute data items such as geologic age, geologic symbols and lithology description; the geological boundary map layer is provided with corresponding attribute data items which are trends and dip angles; the fault layer is provided with a tendency and an inclination angle corresponding to the attribute data item; a yield layer, wherein the corresponding attribute data items are trends and dip angles;

and acquiring parameters corresponding to the model based on the acquired data of the data layer and the corresponding attribute data item.

3. The automatic generation method of a cut geological profile based on a geological space database according to claim 1, wherein: and performing cutting profile calculation based on the section lines to obtain parameters corresponding to each model, wherein the parameters comprise:

performing intersection analysis of the cutting section line and the contour line of the terrain to obtain the elevation value of the intersection point, simultaneously obtaining the elevation values of the cutting line end point and the turning point, and writing the elevation value into the attribute model of the boundary point; wherein, the attribute model of the boundary point comprises: category code number, coordinate X, coordinate Y, elevation, primitive ID, cumulative length, label, type, remark, formation sequence, association ID;

and/or

Performing intersection analysis of the cutting section line and the geologic body to obtain intersection point position information of a boundary line of the geologic body, through which the cutting line passes; writing an attribute model of the boundary point of the geological body through calculation; the attribute model of the boundary point of the geologic body comprises: category code number, coordinate X, coordinate Y, elevation, primitive ID, cumulative length, label, type, remark, formation sequence, association ID;

and/or

Performing intersection analysis of the cutting section line and the fault, acquiring intersection point position information of the fault boundary line penetrated by the cutting line, and writing an attribute model of the intersection point through calculation; the attribute model of the intersection point comprises: category code, serial number, coordinate X, coordinate Y, elevation, primitive ID, dip, true dip, cumulative length, label, type of fault, remark, association ID.

4. A method for automatically generating a cut-to-map geologic profile based on a geospatial database as defined in claim 3, wherein: when the intersecting analysis of the cutting section line and the geologic body is carried out, the method further comprises the following steps: judging whether an intersection point of a boundary line of the geological body penetrated by the cutting line is a fault boundary point or not; if the fault boundary point is the fault boundary point, acquiring a fault name, putting the fault name into a type data item, and then acquiring inclination and true dip angle information of the fault boundary according to the boundary where the fault is located; if the inclination and true dip angle information exists, writing the inclination and true dip angle information into an inclination and true dip angle data item; similarly, if the geological boundary point is the geological boundary point, the contact relation of the geological boundary, the geological body geological code and the stratum sequence information are written into the data item, and if the inclination and true dip angle information exists, the contact relation, the geological body geological code and the stratum sequence information are also written into the inclination and true dip angle data item; if no trend and dip angle information exists, obtaining the trend of the occurrence meeting the conditions, writing the true dip angle value into the trend and true dip angle data item; and searching a yield closest to the intersection point, and if the yield point falls in the stratum corresponding to the intersection point, obtaining the yield which meets the condition.

5. The automatic generation method of the cut geological section based on the geological space database according to claim 1, wherein the anticline and syncline conditions of the stratum are judged by the following methods:

assigning sequence numbers as new and old stratum relation sequence codes; wherein, the older stratum of the age is, the larger the sequence number value is;

analyzing stratum crossing the cutting section line to obtain sequence numbers of each stratum;

judging the nucleus and the wing through the mutation condition of the stratum sequence number, and further judging the anticline and syncline conditions: if the stratum sequence number is changed from small to large and is changed into small, judging that the stratum is anticline according to the new principle of the old wing of the nucleus; if the stratum sequence number is changed from large to small and is larger, the stratum sequence number is judged to be syncline according to the principle that the new wing of the nucleus is old.

6. The method for automatically generating a cut geological profile based on a geological space database according to claim 1 or 5, wherein: calculating the intersection point of the graph cut section line and the geological body arc section and the geological boundary line by drawing the graph cut section line and adopting a space analysis intersection technology; and obtaining the geological unit attributes of the geological body around the intersection point, thereby obtaining the geological times of the boundary around the boundary point, the stratum sequence and the contact relation of the boundary where the boundary point is located.

7. The automatic generation method of a cut geological profile based on a geological space database according to claim 1, wherein: the distribution and extension conditions of new and old strata are automatically judged by the following method:

sequentially segmenting the intersection points of the graph cut section line and the geological boundary line at turning points of the graph cut section line, and sequencing the intersection points in sequence from left to right and from top to bottom, so that the spatial distribution condition of new and old strata is obtained through the attribute of the intersection points; the stratum sequence of the new stratum and the old stratum is obtained by the geologist according to the stratum table of the China regional year or the international geological year representation and sequencing the geological conditions of the working area; assigning a sequence number as a new and old stratum relation sequence code, wherein the older stratum is the older, and the numerical value of the assigned sequence number is larger;

calculating the relative occurrence of the intersection point by using the set buffer zone, and obtaining occurrence representing the geological boundary where the intersection point is located by a nearest distance method, thereby obtaining occurrence attributes and obtaining trend and dip angle values of new and old strata; or according to the correlation model of the boundary points and the corresponding boundary lines, if the corresponding boundary lines have the trend and inclination angle attributes, uniformly converting the trend and inclination angle of the boundary lines into numerical values, and then calculating the apparent inclination angle to obtain the stratum extension condition, wherein if the apparent inclination angle is 90 degrees, the stratum output state is vertical; if the apparent dip angle is 0 degrees, the output state of the stratum is horizontal; if the apparent dip angle epsilon (0,90), the formation output state is inclined towards the right end of the sectional view; if the apparent dip angle E (-90,0), the stratum is inclined towards the left end of the section;

and/or

The axial face attitude information is obtained by the following method:

judging the relation between the inclination of the two wings and the inclination angle, and if the inclination of the two wings is opposite and the inclination angles are nearly equal, indicating that the axial surfaces are vertical; if the inclination angles of the two wings are unequal, the inclination of the axial plane is indicated;

and/or

The hub state is obtained by:

judging the trend state of the two wings, and if the trend directions of the two wings are parallel, the hinge is horizontal; if the two wings meet or are curved in an arc shape, the hinge is inclined; wherein, the tip of the intersection curve of the two wings and the same stratum boundary line is the tilting direction of the anticline hinge and faces the lifting direction of the anticline hinge.

8. The automatic generation method of a graph cut section based on a geospatial database according to claim 1, characterized in that: the form of the turning end is judged by the tight degree of the fold, namely, the turning end is tightly closed as the stratum distance of the two wings is closer, and conversely, the turning end is wider and is gentle.

9. The automatic generation method of a graph cut section based on a geospatial database according to claim 1, characterized in that: creating a graph cut section and simultaneously establishing a point, line and plane space database of the section; the automatic drawing by using the drawing parameters to obtain the cut geological section comprises the following steps:

calculating and processing the section layering boundary line, namely firstly calculating the visual inclination angle and the true inclination angle, then drawing a graphic sketch, firstly projecting each intersection point of a section line and a topographic line on a geological map onto the section map according to a given scale, and forming a topographic section curve; then, each intersection point of the section line and the geological boundary line on the geological map is projected onto a section curve of the topographic map; layering boundaries are drawn according to formation trends and dip angles near each point.

10. The automatic generation method of a graph cut profile based on a geospatial database according to claim 9, characterized in that: the automatic drawing by using the drawing parameters to obtain the cut geological section further comprises the following steps:

obtaining complete boundary boundaries forming the sectional view according to the calculated sectional view sketch, wherein the boundary boundaries comprise vertical lines of projection points of starting and ending endpoints of the sectional view, stratum dividing lines, topographic curves and lower base lines of the sectional view; through space analysis, the automatic line approaching, automatic line cutting, redundant line deleting, automatic calculation and treatment of back oblique line, oblique line and water bottom boundary line are adopted, and then the line is converted into an arc segment; automatically checking and processing topological relations among arc segments, and automatically forming cross-section plane primitives of each geological unit according to the enclosed closed arc segments; then, corresponding graphic elements are automatically assigned according to the geological map related attributes, and the establishment of a profile map attribute library is completed; the color of the original corresponding stratum is given to the primitive as the display parameter of the primitive, thereby completing the coloring work of the section; finally, according to the relevant attribute of the projection point, judging and selecting a proper position, and automatically marking the position including stratum code, fracture name, water system name and azimuth angle finishing of the section view.

Images