CN117910114A - Construction method of complex geologic model and earth and stone digital analysis and calculation method - Google Patents

Abstract

The invention belongs to the field of reconnaissance design, and discloses a construction method of a complex geologic model and a digital analysis and calculation method of earth and stone, wherein the construction method of the complex geologic model comprises the following steps: calibrating key information of a geological survey histogram, constructing a survey hole soil layer database, creating digital survey column holes, generating triangular mesh surfaces, constructing survey column hole groups, constructing geological solid models of all the survey column hole groups, assembling an integral geological solid model and grouping according to soil layer categories. The method automatically generates any complicated digital geological entity model based on the reading of the geological survey histogram, overcomes the defects of large manual participation, low efficiency and non-universal model format of the existing modeling method, breaks through the limitation of rough construction of the geological entity model and lack of corresponding geological attribute information, and simultaneously achieves the purposes of automatically organizing geological data and constructing the model.

Description

Construction method of complex geologic model and earth and stone digital analysis and calculation method

Technical Field

The invention belongs to the field of reconnaissance design, and particularly relates to a construction method of a complex geological model and a digital analysis and calculation method of earth and stone.

Background

Geological investigation is an important and indispensable front-end work in construction engineering, and is a main basis and a necessary condition for the engineering investigation design work. The method is characterized by finding out, analyzing and evaluating the geology, geographic environmental characteristics and geotechnical engineering conditions of a construction site according to the requirements of the construction engineering, providing reasonable basic suggestions and compiling construction engineering investigation files (short for a geological investigation report).

At present, a geological survey report is usually about 20 meters and an exploration point, and the result is mainly represented by a text report and a CAD plane drawing (an exploration point layout diagram, a section diagram and a histogram), so that the soil layer is limited in reflection condition, and the geological soil layer condition at a position of a non-exploration point cannot be intuitively reflected as engineering design data. When the situation is met, a simple three-dimensional geological model is established according to the section view and the histogram of the existing exploration points in the past possibly according to the continuous assumption of the soil level or is estimated approximately by the experience of a mature engineer, but the problem is that the method is effective for the simple geological situation, but serious misjudgment is inevitably caused when the complex geological situation is discontinuously distributed such as cracks, interlayers, karst caves and the like, so that design errors are caused; in addition, the exploration point can be directly supplemented for direct exploration, and the method is most reliable, but has the problem that a large amount of manpower, material resources and time cost can be consumed; and finally, special three-dimensional geological modeling software can be used, which is mainly used in the fields of geological exploration, petroleum and mineral development, has strong modeling capability and higher precision, but has the problems that the software has high use cost (the software is expensive and is only used by special persons), the requirements on data format are strict, and the generated three-dimensional geological model is not convenient to flow downwards, so that the software has very few applications in general construction engineering, particularly civil construction field, and is not as easy to directly detect by adopting direct supplementary exploration points in many cases.

Most construction projects have problems of earth and stone excavation and backfill, and especially for projects with underground buildings or structures, earth and stone analysis is also a very important work in the process of investigation and design. In the past, the analysis is generally performed by CAD plug-in units by adopting approximation algorithms such as irregular triangular network method, square grid method, contour line method, average section method, average elevation method and the like. These methods have advantages and disadvantages, but both suffer from the following two problems:

1. the ratio of earthwork to stone (simply called soil-stone ratio) cannot be obtained, and the accuracy of engineering cost is affected;

2. The loose coefficients of the soil and the stone cannot be effectively considered, and large statistical deviation of engineering quantity is easily caused.

For the first problem, in the past, the soil-to-stone ratio is generally calculated by utilizing on-site soil sample screening and according to the proportion of the mass of the soil and the stone, and the method has large workload and is easy to cause the calculation distortion of the soil-to-stone ratio due to sample problems.

For the second problem, the average value is estimated by experienced engineers, and the method has large human factors and doubtful reliability, and is easy to dispute.

Based on the method, the invention provides a construction method for realizing any complex geological entity model by only reading the geological survey histogram, and simultaneously provides a method for realizing earth-rock digital analysis and calculation by simulating earth-rock entity excavation so as to solve the problems.

Disclosure of Invention

The invention aims to solve the problems of large labor participation, low adaptability, low model construction efficiency, high use cost of professional three-dimensional geological modeling software and low model multiplexing rate of the existing complex geological model construction method and the problems that soil Dan Bi cannot be obtained and soil and stone loose coefficients cannot be effectively assigned by the existing soil and stone analysis method. Therefore, the invention provides a construction method of a complex geological model and a digitized analysis and calculation method of the earth and stone side, which can construct any complex geological entity model by only reading a geological survey histogram and realize digitized analysis and calculation of the earth and stone side by simulating earth and stone side entity excavation.

The technical scheme adopted for solving the technical problems is as follows:

in order to solve the problems of large labor participation, low adaptability, low model construction efficiency, high use cost of professional three-dimensional geological modeling software in construction engineering and low model multiplexing rate of the traditional complex geological model construction method, the invention provides a complex geological model construction method, which comprises the following steps:

S1, calibrating key information of a geodesic histogram: importing a geological survey histogram in modeling software, and calibrating key information in any geological survey histogram;

S2, constructing a exploratory hole soil layer database: analyzing all key information in all the geological survey column diagrams through the key information calibrated in the S1, and constructing a survey pore soil layer database;

S3, creating a digital exploration column hole: calling a exploration hole soil layer database in the S2, creating digital exploration column holes, filling corresponding soil layer information in the digital exploration column holes, and enabling the digital exploration column holes to correspond to the actual exploration holes one by one;

S4, generating a triangular mesh surface and constructing an exploration column hole group: extracting axes of the digital exploration column holes in the step S3, wherein the axes are composed of line segments, each line segment represents a soil layer at a corresponding position, extracting orifice points (namely upper end points of the axes) of all exploration holes, adopting Delaunay algorithm to form the exploration orifice point positions into triangular mesh surfaces, dividing each digital exploration column hole in the S3 into a plurality of groups through vertex codes of each triangular mesh, wherein each triangular mesh corresponds to one digital exploration column hole group;

s5, constructing a geological solid model of each exploration column hole group: according to the mode in S4, each digital exploration column hole group is essentially composed of 3 groups of line segments, according to line segment attribute matching, a surface is composed of points or lines, and then a closed body is composed of the surfaces, so that a digital geological entity model of each exploration column hole group can be constructed;

S6, assembling the whole geological entity model and grouping according to soil layer categories: and traversing each triangular grid in the S4 by the method in the S5 to complete the construction of all the digital geological solid models, and finally storing the digital geological solid models in groups according to soil layer categories.

In a preferred embodiment, in S1, for the critical information calibration, a closed frame of the arbitrary survey bar chart is selected, and the closed frame is required to contain all the information in the survey bar chart and cannot contain the information of other survey bar charts, and the critical information content includes the survey hole number, the plane coordinates of the survey hole, the elevation of the survey hole, the numbers of each soil layer and the layering thickness of the soil layers, and other information can be calculated from the critical information.

In a preferred embodiment, in S2, the exploratory soil layer database is constructed as follows:

S21, selecting English letters (such as 'X' or 'X', 'Y' or 'Y') in plane coordinates of the exploration holes as specific marks in the calibrated key information, and recording the positions of the English letters in CAD as positioning reference points of all calibration key information closed frames;

s22, searching and matching the data information in all the geological survey bar charts through the specific marks in S21, finding the positions of the corresponding specific marks, and mapping all the calibration key information closing frames to the positions of the specific marks through the positioning reference points;

s23, identifying data information in each closed line frame, using the exploratory hole numbers as indexes, storing soil layer key information of each exploratory hole in a tree structure mode, and constructing an exploratory hole soil layer database;

In the exploration hole soil layer database, each exploration hole stores exploration hole numbers and model units, exploration hole points formed by exploration hole coordinates and exploration hole port elevations, each soil layer number (formed by ' identification ' +main layer number + ' "+sub-layer number, if the number is ① ₂, the number is stored as TuC1.2, the soil layer numbers and the soil layer types are in one-to-one correspondence), soil layer layering thickness, and each soil layer top elevation obtained by calculation of the exploration hole elevations and the soil layer layering thickness. This way, the identification of stored data is accurate, efficient, and unaffected by the lot, number, and location of the survey histogram.

In a preferred embodiment, in S3, the digital survey column hole is created by: using the exploration hole numbers as indexes, and creating digital exploration column holes through exploration hole opening points, soil layering thicknesses and top mark heights of all soil layers stored in an exploration hole soil layer database;

The height of each section of digital exploration column hole is consistent with the thickness of the soil layer, the attribute of the column hole is consistent with the soil layer number (soil layer type) and is expressed by a specific color corresponding to the soil layer type, and the diameter of the column hole can be the actual exploration hole diameter (for clarity of display, a larger diameter, such as 1200mm, can be adopted). This approach is equivalent to mapping real exploration holes into digital space, and also facilitates later clearer viewing.

In a preferred embodiment, in S4, a Delaunay algorithm is adopted to store a complete and continuous grid surface (without gaps) formed by triangular grids, wherein the grid surface is formed by a plurality of exploration orifice points (if the exploration orifice points are n, the number of the triangular grids is not less than n-2, and n is more than or equal to 3) as an original ground surface of an exploration area; the exploration hole numbers correspond to exploration hole points and are stored in groups of 3, so that corresponding exploration hole data can be extracted according to the exploration hole number index in S5.

In a preferred embodiment, in S6, a cyclic traversal is adopted, which is specifically as follows:

firstly, selecting the group with the largest line number in the 3 groups of line segments as an initial line segment group, and starting traversing all line segments from top to bottom in sequence; each cycle starts from the line segment traversed by the initial line segment group, then sequentially goes to the 2 nd and 3 rd groups to search whether the line segments with the same attribute exist, if yes, the entity element group is added, if no, the traversed line segment end points are added into the entity element group, and then the unmatched line segments in the 2 nd and 3 rd groups are traversed in a similar way; in the traversal process, the following rules need to be guaranteed:

all line segments of the 3 groups are added into the entity element group;

ensuring that no repeated line segments exist in the entity element group;

The endpoints joining the set of physical elements cannot cross the line segment on a sequential basis.

Thus, in the entity element group, all are arrays composed of line segments or points and containing 3 elements, there are the following 4 cases: (line, wire); (line, dot); (line, dot, line); (line, dot). The most basic pentahedron or tetrahedron geological solid model can be formed by the 4 element groups, and the attribute of the geological solid model is determined by the sideline section. By the method, any complicated digital geologic solid model can be constructed, and the situation that gaps and overlapping exists in the geologic solid model formed by the method can be verified by a 'negative evidence method' on the premise of guaranteeing the 3 rules. The three-dimensional model formats circulated in the market all support the digital geological entity model constructed by the most basic topological relation, which is convenient for subsequent diversion application.

Noteworthy are: the construction mode of the exploration hole soil layer database can be replaced by directly reading the key information data obtained by the exploration hole soil layer database and storing the key information data as the soil layer database, namely S1 and S2 are replaced according to the construction mode, and the rest steps are unchanged.

In order to solve the problems that the existing earth-rock analysis method cannot obtain the earth-rock ratio and cannot effectively assign the earth-rock loose coefficient, the invention provides a digital analysis and calculation method for the earth-rock, which comprises the following steps:

s1', constructing a geological entity model: constructing a geologic solid model by adopting the construction method of the complex geologic model of any one of claims 1 to 6;

s2', calibrating the attribute of the geological entity model: grouping and calibrating the attribute of the geological entity model according to each soil layer type;

S3', generating building object blocks: building object blocks are generated by building ranges and elevation of the base surface;

S4', constructing a residual geological solid model and an excavated geological solid model: traversing each geological entity model in groups, and judging the relation between each geological entity and a building object block in the following manner:

When disjoint, incorporating the geologic entity into a set of remaining geologic entity models;

When intersecting, there are two cases: the geological entity is completely inside the building block, and the geological entity is incorporated into the excavated geological entity model group; and the geological entity part is overlapped with the building object block, the geological entity obtained by the Boolean intersection operation is included in the excavated geological entity model group, and the geological entity obtained by the Boolean difference operation is included in the rest geological entity model group, and the geological entity obtained by the Boolean operation is endowed with the attribute before the Boolean operation.

Grouping and storing the residual geological entity model and the excavated geological entity model according to soil layer types;

S5', calculating a statistical filling amount: extracting the center points of all triangular grids by adopting the triangular grid surface generated in the S4', judging the relation with the bottom surface of the building object block, and calculating statistical filling when the center point is below the bottom surface of the building object block and indicating filling engineering;

s6', statistically analyzing the excavation amount of each soil layer and corresponding parameters: and according to the results of S2 'and S4', the excavation amount, the soil-to-stone ratio and the comprehensive loosening coefficient of each soil layer are obtained through statistical analysis.

In a preferred embodiment, S2' is specifically: assigning each soil layer property parameter x _i according to the geological survey report, wherein i in the parameter is a soil layer number for distinguishing soil layer types, and when x _i =1, the parameter is rock; when x _i =0, it is earth;

And according to each soil layer type, assigning a corresponding loosening coefficient s _i, wherein i in the parameter is a soil layer number for distinguishing the soil layer type.

In a preferred embodiment, in S5', when it is determined that a filling project exists, dividing the bottom surface of the building block into a plurality of grids, obtaining the center point of each grid, projecting downward one by one, and if the building block falls on the triangular grid surface generated in S4, recording the area of the grid as a _i, where the calculation formula of the filling TF is as follows:

TF＝∑A_i×cos(α_i)×l_i；

Wherein i is the filling grid number: 1.2, 3, …, n; alpha _i is the included angle between the normal line of the grid surface and the Z axis; and l _i is the distance between the center point and the falling point.

In a preferred embodiment, S6' is calculated as follows:

setting n groups of the excavated geological entity model groups, wherein each group comprises m geological entities, and each geological entity volume is The excavation amount WF ^k of each soil layer is:

The soil-to-stone ratio TSB is:

The integrated loose coefficients SSXS are:

Where k is the number of the group of the excavated geological entity model: 1.2, 3, …, n; j is the number of each geological entity in the geological entity model group: 1.2, 3, …, m.

Noteworthy are: in S3', for the case that the substrate is a slope, the method can be solved by firstly creating the slope according to the initial elevation of the slope and then stretching the slope vertically (Z-axis) to generate building blocks, and other steps are unchanged.

Compared with the prior art, the invention has the beneficial effects that:

1. On the premise of not increasing extra cost, the invention can automatically generate any complicated digital geological entity model based on the read geological survey histogram, overcomes the defects of large manual participation, low efficiency and non-universal model format of the existing modeling method, breaks through the limitation of rough construction of the geological entity model and lack of corresponding geological attribute information, and simultaneously realizes the purposes of automatically organizing geological data and constructing the model.

2. The invention also realizes detailed acquisition of the detailed excavation quantity, the soil-to-stone ratio, the comprehensive loosening coefficient and the like of each soil layer by simulating the soil-stone entity excavation, and greatly improves the efficiency and the accuracy of the soil-stone engineering quantity statistics compared with the prior art.

Drawings

For a clearer description of the technical solutions of the embodiments of the present invention, the drawings that are needed in the description of the embodiments will be briefly introduced below, it being obvious that the drawings in the description below are only some embodiments of the present invention, and that other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art, wherein:

FIG. 1 is a flow frame diagram of a method for constructing a complex geologic model and a method for digitally analyzing and calculating earthwork;

FIG. 2 is a schematic diagram of key information of a marker survey histogram;

FIG. 3 is a diagram of analytical survey histogram stored data;

FIG. 4 is a schematic diagram of creating a digital survey string hole;

FIG. 5 is a schematic diagram of a digitized geologic solid model generated from reading data packets;

FIG. 6 is a schematic diagram of a single digital survey column well bore set construction digitized geological entity;

FIG. 7 is a schematic diagram of the construction of an integral digitized geologic solid model;

FIG. 8 is a schematic diagram of base data;

FIG. 9 is a schematic diagram of a digitized analysis of earth and rock excavation performed by reading substrate data;

FIG. 10 is a schematic diagram of a simulation of earth excavation;

fig. 11 is a schematic diagram of analysis and statistics results of excavation of earth and stone.

Wherein: fig. 2 (1) is a diagram frame of arbitrary survey histogram; (2) numbering the exploration holes; (3) exploratory hole coordinates; (4) the exploration orifice elevation; (5) numbering soil layers; and (6) the layering thickness of the soil layer.

Detailed Description

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 order to realize the technical scheme of the invention, a set grasshopper of functional implementation batteries (called gh batteries for short) are formed through a rhinoceros CAD platform and adopting Python script language after secondary development. The technical scheme of the invention will be clearly and completely described below with reference to the gh battery and the attached drawings. It will be apparent that the described embodiments are only some, but not all, embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.

Example 1: the invention provides a construction method of a complex geological model, which automatically generates a digital geological entity model based on reading a geological survey histogram, and is carried out according to the following steps when being implemented:

First, the key information of the histogram is surveyed by the mark, as shown in fig. 2: and adopting a closed wire frame to frame and select key information such as exploration hole numbers, exploration hole plane coordinates (XY plane coordinates), exploration hole orifice elevations, each soil layer number, soil layer layering thickness and the like in the geological exploration histogram.

And secondly, mapping key information selected by a frame in the closed line frame to all the geological survey data by adopting a columnar geological survey gh battery, analyzing the data information in the geological survey data, and storing the geological survey data in a tree structure according to rules in the scheme in the invention, as shown in figure 3. To facilitate the review of the geosurvey data, the cylindrical geosurvey gh battery creates a three-dimensional digital survey cylinder hole in rhinoceros CAD space, as shown in fig. 4.

Thirdly, three-dimensional geological gh batteries are adopted to read tree-shaped geological survey data, the tree-shaped geological survey data are grouped according to the vertex numbers of triangular meshes, and geological entity models are generated by sequentially grouping, as shown in the figure 5, and the three-dimensional geological gh batteries are as follows: {8} in the geologic entity data storage structure represents a 9 th set of geologic entity models, where there are 629 TuC geologic entities, and so on. When the geologic solid models are generated in groups, the points are connected into lines according to rules in the scheme in the content of the invention, the lines form planes in a surrounding mode, and the planes are sealed into a body to generate a single group of geologic solid models, as shown in figure 6. Then traversing all groups in turn, the entire geologic solid model is generated, as shown in FIG. 7.

Example 2: the invention provides a method for digitally analyzing and calculating an earth and stone side, which realizes the digital analysis and calculation of the earth and stone side by simulating the excavation of the earth and stone side entity, and is concretely implemented according to the following steps:

In a first step, the extent and elevation of the building base is defined on the original ground surface (or the current ground) of the exploration area, as shown in fig. 8. Building object blocks are generated by building base ranges and elevations according to rules in the scheme described in the present disclosure. The building base range is represented by a closed line frame, and can support any plurality of closed line frames, namely when the base range has a plurality of elevations, the base range is divided into a plurality of closed line frames, and then the closed line frames are assigned to be realized by different elevations.

And secondly, assigning corresponding property parameters and loosening coefficients according to each soil layer type, as shown in figure 9. The assignment format is [ soil layer category, loose coefficient, property parameter ]. The property parameters and the loose coefficients can be obtained by inquiring in related technical documents according to soil layer types.

Thirdly, adopting a soil and stone excavation analysis gh battery, and accessing the geological entity model generated in the embodiment 1 to perform soil and stone excavation analysis, as shown in fig. 9: the earth-rock excavation analysis of the underside of the gh battery box in fig. 9 is time consuming, wherein the three-dimensional geological generation takes 10.3s and the earth-rock entity excavation simulation analysis takes 12.4s.

The first and second steps in this embodiment are actually implemented as an interface between "base data" and "soil parameters" of the earth-rock excavation analysis gh battery. This gh cell simulates solid earth excavation to obtain each "excavated geologic entity" and "remaining geologic entities" in a boolean operation, as shown in fig. 10. Statistical parameters are then calculated to guide the engineering design through rules described in the schemes described in this summary, as shown in fig. 11.

The foregoing description is only illustrative of the present invention and is not intended to limit the scope of the invention, and all equivalent structures or equivalent processes or direct or indirect application in other related arts are included in the scope of the present invention.

Claims (10)

1. A construction method of a complex geological model is characterized by comprising the following steps: the method comprises the following steps:

S1, calibrating key information of a geodesic histogram: importing a geological survey histogram in modeling software, and calibrating key information in any geological survey histogram;

S2, constructing a exploratory hole soil layer database: analyzing all key information in all the geological survey column diagrams through the key information calibrated in the S1, and constructing a survey pore soil layer database;

S3, creating a digital exploration column hole: calling a exploration hole soil layer database in the S2, creating digital exploration column holes, filling corresponding soil layer information in the digital exploration column holes, and enabling the digital exploration column holes to correspond to the actual exploration holes one by one;

S4, generating a triangular mesh surface and constructing an exploration column hole group: extracting axes of the digital exploration column holes in the step S3, wherein the axes are formed by line segments, each line segment represents a soil layer at a corresponding position, extracting orifice points of all exploration holes, adopting a Delaunay algorithm to form the exploration orifice points into triangular mesh surfaces, dividing each digital exploration column hole in the step S3 into a plurality of groups through vertex coding of each triangular mesh, and each triangular mesh corresponds to one digital exploration column hole group;

s5, constructing a geological solid model of each exploration column hole group: according to the mode in S4, each digital exploration column hole group is essentially composed of 3 groups of line segments, according to line segment attribute matching, a surface is composed of points or lines, and then a closed body is composed of the surfaces, so that a digital geological entity model of each exploration column hole group can be constructed;

S6, assembling the whole geological entity model and grouping according to soil layer categories: and traversing each triangular grid in the S4 by the method in the S5 to complete the construction of all the digital geological solid models, and finally storing the digital geological solid models in groups according to soil layer categories.

2. The method for constructing a complex geologic model according to claim 1, wherein: in S1, when the key information is marked, a closed frame of any survey bar chart is selected, and the closed frame is required to contain all information in the survey bar chart and cannot contain information of other survey bar charts, and the key information content includes the survey hole number, the plane coordinates of the survey hole, the elevation of the survey hole, the numbers of each soil layer and the layering thickness of the soil layers.

3. The method for constructing a complex geologic model according to claim 2, wherein: s2, constructing a exploratory hole soil layer database according to the following method:

S21, selecting English letters in plane coordinates of the exploration holes in the calibrated key information as specific marks, and recording the positions of the English letters in CAD as positioning reference points of all calibration key information closed frames;

s22, searching and matching the data information in all the geological survey bar charts through the specific marks in S21, finding the positions of the corresponding specific marks, and mapping all the calibration key information closing frames to the positions of the specific marks through the positioning reference points;

s23, identifying data information in each closed line frame, using the exploratory hole numbers as indexes, storing soil layer key information of each exploratory hole in a tree structure mode, and constructing an exploratory hole soil layer database;

In the exploration hole soil layer database, each exploration hole stores an exploration hole number and a model unit, an exploration hole point formed by exploration hole coordinates and exploration hole orifice elevation, each soil layer number, soil layer layering thickness and each soil layer top elevation obtained by calculation of the exploration hole orifice elevation and the soil layer layering thickness.

4. A method of constructing a complex geologic model as claimed in claim 3, wherein: in S3, the creation mode of the digital exploration column hole is as follows: using the exploration hole numbers as indexes, and creating digital exploration column holes through exploration hole opening points, soil layering thicknesses and top mark heights of all soil layers stored in an exploration hole soil layer database;

The height of each section of digital exploration column hole is consistent with the thickness of the soil layer, the attribute of the column hole is consistent with the number of the soil layer and expressed by a specific color corresponding to the category of the soil layer, and the diameter of the column hole can be the actual exploration hole diameter.

5. The method for constructing a complex geologic model according to claim 1, wherein: s6, performing in a cyclic traversal mode, wherein the method specifically comprises the following steps:

firstly, selecting the group with the largest line number in the 3 groups of line segments as an initial line segment group, and starting traversing all line segments from top to bottom in sequence; each cycle starts from the line segment traversed by the initial line segment group, then sequentially goes to the 2 nd and 3 rd groups to search whether the line segments with the same attribute exist, if yes, the entity element group is added, if no, the traversed line segment end points are added into the entity element group, and then the unmatched line segments in the 2 nd and 3 rd groups are traversed in a similar way; in the traversal process, the following rules need to be guaranteed:

all line segments of the 3 groups are added into the entity element group;

ensuring that no repeated line segments exist in the entity element group;

The endpoints joining the set of physical elements cannot cross the line segment on a sequential basis.

6. A digitalized analysis and calculation method for earthwork is characterized in that: the method comprises the following steps:

s1', constructing a geological entity model: constructing a geologic solid model by adopting the construction method of the complex geologic model of any one of claims 1 to 6;

s2', calibrating the attribute of the geological entity model: grouping and calibrating the attribute of the geological entity model according to each soil layer type;

S3', generating building object blocks: building object blocks are generated by building ranges and elevation of the base surface;

S4', constructing a residual geological solid model and an excavated geological solid model: traversing each geological entity model in groups, and judging the relation between each geological entity and a building object block in the following manner:

When disjoint, incorporating the geologic entity into a set of remaining geologic entity models;

When the geological entity is completely inside the building block, the geological entity is incorporated into the excavated geological entity model group;

When the geological entity part is overlapped with the building object block, firstly adopting the geological entity obtained by Boolean intersection operation to incorporate the excavated geological entity model group, and then adopting the geological entity obtained by Boolean difference operation to incorporate the residual geological entity model group;

grouping and storing the residual geological entity model and the excavated geological entity model according to soil layer types;

S5', calculating a statistical filling amount: extracting the center points of all triangular grids by adopting the triangular grid surface generated in the S4', judging the relation with the bottom surface of the building object block, and calculating statistical filling when the center point is below the bottom surface of the building object block and indicating filling engineering;

s6', statistically analyzing the excavation amount of each soil layer and corresponding parameters: and according to the results of S2 'and S4', the excavation amount, the soil-to-stone ratio and the comprehensive loosening coefficient of each soil layer are obtained through statistical analysis.

7. The method for digitally analyzing and calculating the earthwork according to claim 6, wherein the method comprises the following steps: the S2' is specifically as follows: assigning each soil layer property parameter x _i according to the geological survey report, wherein i in the parameter is a soil layer number for distinguishing soil layer types, and when x _i =1, the parameter is rock; when x _i =0, it is earth;

And according to each soil layer type, assigning a corresponding loosening coefficient s _i, wherein i in the parameter is a soil layer number for distinguishing the soil layer type.

8. The method for digitally analyzing and calculating the earthwork according to claim 6, wherein the method comprises the following steps: in S3', the building object blocks are space entities, when the building is an underground building, the building object blocks take actual heights, otherwise, the building object blocks take relatively larger heights, and soil layer excavation insufficiency is avoided.

9. The method for digitally analyzing and calculating the earthwork according to claim 7, wherein the method comprises the steps of: in S5', when it is judged that filling engineering exists, dividing the bottom surface of a building block into a plurality of grids, obtaining the center point of each grid, projecting downwards one by one, and if the building block falls on the triangular grid surface generated in S4, recording the grid area as A _i, wherein the calculation formula of filling TF is as follows:

TF＝∑A_i×cos(α_i)×l_i；

Wherein i is the filling grid number: 1.2, 3, …, n; alpha _i is the included angle between the normal line of the grid surface and the Z axis; and l _i is the distance between the center point and the falling point.

10. The method for digitally analyzing and calculating the earthwork according to claim 9, wherein the method comprises the steps of: in S6', the following is calculated:

setting n groups of the excavated geological entity model groups, wherein each group comprises m geological entities, and each geological entity volume is The excavation amount WF ^k of each soil layer is:

The soil-to-stone ratio TSB is:

The integrated loose coefficients SSXS are:

Where k is the number of the group of the excavated geological entity model: 1.2, 3, …, n; j is the number of each geological entity in the geological entity model group: 1.2, 3, …, m.