BIM (Building Information Modeling)-based method for measuring and calculating earthwork amount of water

Abstract

The invention relates to the technical field of computer-aided building design, in particular to a BIM (Building Information Modeling)-based water conservancy project earthwork engineering quantity measuring and calculating method. The method comprises the following steps: constructing a three-dimensional model, and determining an operation area and a zero line; dividing the three-dimensional model based on different preset side length values to obtain an initial grid; according to the model height difference, screening to obtain a mutation point; determining a construction strength correction value according to the area distribution and height design change of the operation region where the mutation point is located; determining a reference value coefficient of the initial grid according to the number of sudden change points on each side of the initial grid, the construction strength correction values of the sudden change points and the number of zero points, and further determining the optimal side length; and realizing square grid division based on the optimal side length, and carrying out engineering quantity measurement and calculation. According to the method, the complexity of the measurement and calculation process can be avoided, the self-adaption degree of grid division is improved, and the measurement and calculation precision is enhanced.

Description

Technical Field

[0001] The present invention relates to the technical field of computer-aided architectural design, and particularly to a method for calculating earthwork volume of water conservancy projects based on BIM. Background Art

[0002] The technological advancement of water conservancy projects enables people to manage water resources more effectively. Site leveling is one of the basic preparatory works before the construction of some water conservancy projects (such as drainage systems, transportation systems, etc.), providing necessary preparations for subsequent various construction activities.

[0003] Existing methods use the grid method to pre-calculate the earthwork volume before site leveling. However, when pre-calculating a site with slightly larger undulations (higher irregularity), since there may be both areas that need to be filled and areas that need to be excavated in the same grid, and the grid method only uses the data of the grid corners for pre-calculation, the accuracy of the calculated results obtained in this way is often not satisfactory. For example, if a certain grid area is sunken in the middle and slightly higher around, it actually needs to be filled, but calculating only based on the data of the grid corners will result in a result that it needs to be excavated. In order to obtain a more accurate result, it is necessary to continuously reduce the side length of the grid to improve the calculation accuracy. In this way, the calculation process becomes extremely complex and the calculation efficiency is poor. Summary of the Invention

[0004] In order to solve the technical problems in the related art that only the data of the grid corners are used for pre-calculation, resulting in insufficient accuracy of the calculated results and poor calculation efficiency, the present invention provides a method for calculating earthwork volume of water conservancy projects based on BIM. The specific technical solutions adopted are as follows: The present invention proposes a method for calculating earthwork volume of water conservancy projects based on BIM. The method includes: Construct a three-dimensional model of the site to be measured, and combine it with the design elevation surface to determine the operation area. The operation area includes the excavation area and the filling area. The plane line where the design elevation surface overlaps with the three-dimensional model surface is the zero line; divide the three-dimensional model based on the grid method and different preset side length values to obtain the initial grid; Periodically obtain different analysis points on any side of any initial grid; screen out the mutation points during the construction process according to the model height difference of the analysis points; determine the construction intensity correction value according to the area distribution and height design change of the operation area where the mutation points are located; Take the intersection point of the initial grid side and the zero line as the zero point; determine the reference value coefficient of the initial grid according to the number of mutation points on each side of the initial grid, the construction intensity correction value of the mutation points, and the number of zero points. Combine all the reference value coefficients with the same preset side length value to determine the optimal side length; Implement grid division based on the optimal side length and conduct engineering quantity measurement.

[0005] Furthermore, the three-dimensional model area that is not on the same horizontal plane as the design elevation plane is used as the operation area.

[0006] Furthermore, according to the model height difference of the analysis points, mutation points during the construction process are screened, including: On the initial grid side, confirm the other two analysis points closest to any analysis point as adjacent points; Calculate the absolute value of the model height difference between the analysis point and the adjacent points respectively, average and normalize to obtain the mutation degree; Take the analysis points with the mutation degree greater than the preset mutation threshold as mutation points.

[0007] Furthermore, according to the area distribution and height design change of the operation area where the mutation points are located, determine the construction intensity correction value, including: According to the comparison of the operation area where the mutation points are located, the initial grid, and the area of the operation area within the initial grid, determine the area influence index; Analyze the height to be constructed in the operation area where the mutation points are located within the initial grid to determine the height influence index; Combine the height influence index and the area influence index to determine the construction intensity correction value.

[0008] Furthermore, according to the comparison of the operation area where the mutation points are located, the initial grid, and the area of the operation area within the initial grid, determine the area influence index, including: Take the area of the operation area where the mutation points are located within the initial grid as the first area, the area of the operation area where the mutation points are located as the second area, and the area of the initial grid as the third area; Calculate the ratio of the first area to the second area as the first ratio and the ratio of the first area to the third area as the second ratio respectively; Take the sum value of the first ratio and the second ratio as the area influence index.

[0009] Furthermore, analyze the height to be constructed in the operation area where the mutation points are located within the initial grid to determine the height influence index, including: Average and determine a preset number of height points from the operation area where the mutation points are located within the initial grid; Take the absolute value of the difference between the model height and the design elevation plane height at the same point as the construction height; Take the average value of the construction heights of all height points as the height influence index.

[0010] Furthermore, combine the height influence index and the area influence index to determine the construction intensity correction value, including: Calculate the product value of the area influence index and the height influence index, and normalize it as the construction intensity correction value of the mutation point.

[0011] Furthermore, according to the number of mutation points on each side of the initial grid, the construction intensity correction value of the mutation point, and the number of zero points, determine the reference value coefficient of the initial grid, including: Take the average value of the construction intensity correction values of all mutation points on the same side as the side intensity coefficient; Take the ratio of the number of zero points on the same side to the total number of zero points on the four sides of the initial grid as the zero point proportion; Calculate the product value of the side intensity coefficient and the zero point proportion as the side reference index on the side of the initial grid; Take the opposite number of the sum of the side reference indexes of all sides of the initial grid, and normalize it as the reference value coefficient.

[0012] Furthermore, combine all the reference value coefficients with the same preset side length value to determine the optimal side length, including: Calculate the average value of the reference value coefficients of all initial grids with the same preset side length value, and take the product value of the average value and the preset side length value as the reference coefficient; Take the preset side length value with the largest reference coefficient as the optimal side length.

[0013] Furthermore, the preset side length values include 1 meter, 3 meters, and 5 meters.

[0014] The present invention has the following beneficial effects: The present invention determines the operation area by constructing a three-dimensional model of the site to be measured and combining it with the designed elevation plane. The three-dimensional model is divided by different preset side length values to obtain an initial grid. By analyzing different preset side length values, it is convenient to select a more appropriate grid division situation, so that when calculating the engineering quantity, the reliability is stronger. Then, in any side of any initial grid, different analysis points are periodically obtained; according to the model height difference of the analysis points, mutation points during the construction process are screened; according to the area distribution and height design change of the operation area where the mutation points are located, the construction intensity correction value is determined. Through two dimensions of area distribution and height distribution, the construction intensity correction analysis at the position of the mutation points is realized. The construction intensity correction value can represent the terrain complexity value at the corresponding mutation point position; then, according to the number of mutation points on each side of the initial grid, the construction intensity correction value of the mutation points, and the number of zero points, the reference value coefficient of the initial grid is determined. Combining all the reference value coefficients with the same preset side length value, the optimal side length is determined, that is, analyzing from the mutation point to the edge of the grid where the mutation point is located, and further analyzing the complexity of the initial grid itself, obtaining the reference value coefficient of the initial grid, and analyzing and obtaining the optimal side length to realize the grid division and calculate the engineering quantity. Compared with directly calculating through the corner points of a fixed grid, in the embodiment of the present invention, multiple points on the grid edges are analyzed to determine a more reliable grid division, which can avoid the complexity of the calculation process. At the same time, the adaptability of the grid division is improved, and the calculation accuracy is enhanced. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] In order to more clearly illustrate the technical solutions and advantages in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required to be used in the description of the embodiments or the prior art. Obviously, the following drawings are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0016] Figure 1 It is a flowchart of a method for calculating the earthwork quantity of a water conservancy project based on BIM provided by an embodiment of the present invention; Figure 2 It is a schematic diagram of a three-dimensional model provided by an embodiment of the present invention; Figure 3 It is a schematic diagram of the zero point line in the three-dimensional model provided by an embodiment of the present invention; Figure 4 It is a three-dimensional analysis diagram of the construction intensity correction value provided by an embodiment of the present invention; Figure 5 It is a distribution diagram of grid edge point positions provided by an embodiment of the present invention; Figure 6The optimal side length decision analysis diagram provided by an embodiment of the present invention; Figure 7 The schematic diagram of grid division provided by an embodiment of the present invention. Detailed implementation manners

[0017] In order to further elaborate on the technical means and effects adopted by the present invention to achieve the predetermined invention purpose, the following combines the accompanying drawings and preferred embodiments to elaborate in detail on a method for calculating the earthwork volume of a water conservancy project based on BIM proposed according to the present invention, its specific implementation manners, structures, features and effects. In the following description, different "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. In addition, the specific features, structures or characteristics in one or more embodiments can be combined in any suitable form.

[0018] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by those skilled in the technical field to which the present invention belongs.

[0019] The following specifically describes the specific solution of a method for calculating the earthwork volume of a water conservancy project based on BIM provided by the present invention with reference to the accompanying drawings.

[0020] Please refer to Figure 1 , which shows a flowchart of a method for calculating the earthwork volume of a water conservancy project based on BIM provided by an embodiment of the present invention. The method includes: S101: Construct a three-dimensional model of the site to be measured, and combine it with the design elevation surface to determine the operation area, which includes the excavation area and the filling area. The plane line where the design elevation surface overlaps with the three-dimensional model surface is the zero line; divide the three-dimensional model based on the grid method and different preset side length values to obtain the initial grid.

[0021] The existing method uses the grid method to predict the earthwork volume before site leveling. However, when predicting and calculating a site with a slightly larger undulation (higher irregularity), since there may be both areas that need to be filled and areas that need to be excavated in the same grid, and the grid method only uses the data at the grid corners to predict and calculate, the accuracy of the calculated result is often not satisfactory. For example, if a certain grid area is sunken in the middle and slightly higher around, it actually needs to be filled, but calculating only based on the data at the grid corners will result in a result that it needs to be excavated. In order to obtain a more accurate result, it is necessary to continuously reduce the side length of the grid to improve the calculation accuracy. In this way, the calculation process becomes extremely complicated.

[0022] Therefore, in the embodiment of the present invention, the side length of the grid is mainly analyzed according to the terrain to determine the adaptive side length of the grid, which improves the adaptability of grid division while reducing the calculation amount, so as to enhance the calculation accuracy.

[0023] In the embodiments of the present invention, a Building Information Modeling (BIM) tool can be used to establish a three-dimensional model of a virtual construction project based on the site information for which the earthwork volume needs to be measured by using digital technology. Refer to Figure 2 , Figure 2 which is a schematic diagram of the three-dimensional model provided by an embodiment of the present invention. It should be noted that during the process of constructing the three-dimensional model, preprocessing operations can also be performed on the model, including but not limited to cleaning the model data, removing duplicate data, noise, and outliers; appropriately reducing the model complexity, optimizing the number of model polygons, etc. The processed three-dimensional model is used for subsequent analysis and calculation.

[0024] After determining the three-dimensional model, the plane of the site terrain in the three-dimensional model is the original plane, and the designed elevation plane is the plane corresponding to after the project is completed. Directly overlap the original plane at the same position and the same height with the designed elevation plane, and take the three-dimensional model area that is not on the same horizontal plane as the designed elevation plane as the operation area. The part where the original plane is lower is the filling area, and the part where the original plane is higher is the excavation area; and the plane line where the designed elevation plane overlaps with the original plane of the three-dimensional model surface is the zero-point line. Refer to Figure 3 , Figure 3 which is a schematic diagram of the zero-point line in the three-dimensional model provided by an embodiment of the present invention.

[0025] In the embodiments of the present invention, the preset side length values include 1 meter, 5 meters, and 10 meters, and the grid division can be realized according to different preset side lengths to obtain the initial grid corresponding to the preset side length.

[0026] S102: Periodically obtain different analysis points on any side of any initial grid; screen out the mutation points during the construction process according to the model height difference of the analysis points; determine the construction intensity correction value according to the area distribution and height design change of the operation area where the mutation points are located.

[0027] In the grid method, the positive and negative values of the construction height of the grid corner points are used to represent whether excavation or filling is required here. A positive value indicates that excavation is required, and a negative value indicates that filling is required. No matter what grid side length is selected, there may still be a phenomenon that the values of the four corner points cannot correctly represent the construction situation required in the internal area. In the embodiments of the present invention, by adding other key points on the grid side, such as the zero point of positive and negative conversion, the mutation point of construction height mutation, etc., the construction information required in the internal area can be more accurately expressed, avoiding information loss.

[0028] In the embodiments of the present invention, in order to analyze the state of the initial grid itself, different analysis points can be periodically and additionally selected on the edges of the initial grid. For example, an analysis point is selected every 10 cm on the edge of the initial grid, and there is no limitation in this regard. It should be noted that the analysis points are the point positions on the surface of the three-dimensional model, and their specific positions are the grid point positions.

[0029] Among them, the mutation point is an analysis point where the surrounding properties change violently, mainly reflected in the violent change of height; since the model height of the analysis point determines the horizontal height of the ground, and the designed elevation plane is the height to be constructed, the screening of mutation points can be realized according to the height changes of different analysis points.

[0030] Furthermore, in some embodiments of the present invention, mutation points during the construction process are screened according to the model height differences of the analysis points, including: identifying the other two analysis points closest to any analysis point on the edge of the initial grid as adjacent points; respectively calculating the absolute values of the model height differences between the analysis point and the adjacent points, averaging and normalizing them to obtain the mutation degree; and taking the analysis points with the mutation degree greater than the preset mutation threshold as mutation points.

[0031] That is, when the model height changes more greatly, it indicates that a relatively large change occurs in the corresponding terrain, and the mutation degree of the position of this analysis point is higher. In the embodiments of the present invention, the mutation degree is obtained through normalization processing, which is convenient for subsequent threshold analysis. This normalization processing can specifically be, for example, linear normalization, such as maximum-minimum linear normalization, or normalization processing based on the Sigmod activation function, and there is no limitation in this regard.

[0032] In the embodiments of the present invention, the preset mutation threshold is the threshold value of the mutation degree. The preset mutation threshold can be set to 0.6. That is to say, when the mutation degree is greater than 0.6, the corresponding analysis point is taken as a mutation point. See Table 1, which is the mutation degree analysis table: Table 1 The mutation point indicates that a relatively large terrain height change occurs correspondingly. Therefore, it is necessary to analyze the operation area where the mutation point is located. In the embodiments of the present invention, it is mainly through two dimensions: height analysis and area analysis.

[0033] Furthermore, in some embodiments of the present invention, the construction intensity correction value is determined according to the area distribution and height design change of the operation area where the mutation point is located, including: determining the area influence index according to the comparison of the operation area where the mutation point is located, the initial grid, and the area of the operation area within the initial grid; analyzing the height to be constructed in the operation area where the mutation point is located within the initial grid to determine the height influence index; and combining the height influence index and the area influence index to determine the construction intensity correction value.

[0034] First, for the analysis of area, it can be understood that the construction height of each mutation point is calculated only based on the data at that point, and it cannot well represent the situation of the site around the mutation point. For a certain mutation point, the construction situation required in a small neighborhood range can be represented by the data of its own construction height; while the construction situation required in a slightly larger range around it needs to be comprehensively represented by the data within the operation area where the mutation point is located, that is, the analysis is achieved through the area comparison of the operation area.

[0035] Further, in some embodiments of the present invention, according to the operation area where the mutation point is located, the initial grid, and the area comparison of the operation area within the initial grid, an area influence index is determined, including: taking the area of the operation area where the mutation point is located within the initial grid as the first area, the area of the operation area where the mutation point is located as the second area, and the area of the initial grid as the third area; calculating the ratio of the first area to the second area as the first ratio, and the ratio of the first area to the third area as the second ratio; taking the sum value of the first ratio and the second ratio as the area influence index, and the corresponding calculation formula is: In the formula, K represents the area influence index, represents the area of the operation area where the th mutation point is located, that is, the first area; L represents the side length of the initial grid, represents the area of the initial grid, that is, the third area; represents the area of the operation area where the mutation point is located, that is, the second area. represents the second ratio, represents the first ratio.

[0036] That is to say, calculate the ratio of the area of the operation area where the mutation point is located within the initial grid to the total area of the operation area where the mutation point is located as the first ratio. The larger the value of the first ratio, it means that the corresponding operation area belongs to the initial grid, and the influence of the operation area where the mutation point is located on the initial grid is greater.

[0037] Among them, calculate the ratio of the area of the operation area where the mutation point is located within the initial grid to the area of the initial grid as the second ratio. The larger the second ratio, it also means that the operation area accounts for a larger proportion in the initial grid, and the influence of the operation area where the mutation point is located on the initial grid is greater.

[0038] Therefore, directly calculate the sum value of the first ratio and the second ratio as the area influence index. In the embodiments of the present invention, the area influence index characterizes the influence index of the operation area where the mutation point is located on the initial grid, and the larger its value, the greater the influence of the corresponding mutation point on the initial grid.

[0039] Further, in some embodiments of the present invention, according to the height to be constructed in the operation area where the mutation point is located in the initial grid, a height influence index is determined, including: averaging and determining a preset number of height points in the operation area where the mutation point is located in the initial grid; taking the absolute value of the difference between the model height and the design elevation height at the same point as the construction height; taking the average value of the construction heights of all height points as the height influence index.

[0040] Among them, averaging and determining means uniformly determining height points within the same interval range in the operation area. The preset number can be, for example, 10, mainly analyzing the overall level of the operation area itself.

[0041] In the embodiments of the present invention, taking the average value of the construction heights of all height points as the height influence index means that the greater the construction height, the greater the absolute value of the difference between the model height and the design elevation height at the corresponding height point, and the greater the construction intensity to be carried out, that is, the greater the influence of the mutation point on the initial grid.

[0042] In summary, in the embodiments of the present invention, by combining the height influence index and the area influence index, a construction intensity correction value is determined, including: calculating the product value of the area influence index and the height influence index, and normalizing it as the construction intensity correction value of the mutation point.

[0043] See Figure 4 , Figure 4 which is a three-dimensional analysis diagram of the construction intensity correction value provided by an embodiment of the present invention. Among them, the area influence index, the height influence index, and the correction value are respectively used as three dimensions to more intuitively represent the parameter changes of the construction intensity correction value for excavation and filling.

[0044] By combining the area influence and the height influence for analysis, the construction intensity correction value of the mutation point to the initial grid can be accurately judged. The larger this value is, the greater the construction intensity of the mutation point to the initial grid, and it can also represent the construction complexity within the initial grid.

[0045] S103: Taking the points where the initial grid edges intersect the zero line as zeros; according to the number of mutation points on each edge of the initial grid, the construction intensity correction value of the mutation points, and the number of zeros, determining the reference value coefficient of the initial grid, and combining all the reference value coefficients with the same preset side length value to determine the optimal side length.

[0046] The edge of the initial grid has intersections with the zero line, and this intersection can be called a zero. The zero on the edge of the initial grid represents the boundary of the corresponding operation area. The number of zeros on the edge of the initial grid can, to a certain extent, represent the complexity of the shape within the initial grid itself. See Figure 5 , Figure 5The grid edge point distribution diagram provided by an embodiment of the present invention. The grid edge includes corner points, mutation points, and zero points. The more mutation points and zero points are distributed on the edge, the higher the corresponding complexity level.

[0047] Determine the reference value coefficient of the initial grid according to the number of mutation points on each side of the initial grid, the construction intensity correction value of the mutation points, and the number of zero points, including: taking the average value of the construction intensity correction values of all mutation points on the same side as the edge intensity coefficient; taking the ratio of the number of zero points on the same side to the total number of zero points on the four sides of the initial grid as the zero point ratio; calculating the product value of the edge intensity coefficient and the zero point ratio as the edge reference index on the edge of the initial grid; taking the opposite number of the sum value of the edge reference indexes of all edges of the initial grid and performing normalization processing as the reference value coefficient.

[0048] For the analysis of the same initial grid, there are 4 edges in the same initial grid, and each edge may have mutation points. Therefore, taking the average value of the construction intensity correction values of all mutation points on the same side as the edge intensity coefficient, and when there are no mutation points on the same side, the corresponding edge intensity coefficient is 0.

[0049] Among them, taking the ratio of the number of zero points on the same side to the total number of zero points on the four sides of the initial grid as the zero point ratio. The larger the zero point ratio, it also indicates that the situation of the corresponding edge is more complex and the weight of the edge intensity coefficient is higher. Therefore, the zero point ratio can be used as the weight to calculate the product value of the edge intensity coefficient and the zero point ratio as the edge reference index on the edge of the initial grid.

[0050] Since the larger the value of the edge reference index, it indicates that the mutation point situation of the corresponding initial grid edge is more complex, and further indicates that the terrain in the initial grid fluctuates more during construction. Therefore, the reference significance of the initial grid division is lower. In the embodiment of the present invention, the reference value coefficient is used to characterize the reference significance during the initial grid division. Taking the opposite number of the sum value of the edge reference indexes of all edges of the initial grid and performing normalization processing as the reference value coefficient.

[0051] After determining the reference value coefficient of each initial grid itself, analysis can be carried out for different preset side length values. Combining the reference value coefficients of all initial grids with the same preset side length value, determine the optimal side length, including: calculating the average value of the reference value coefficients of all initial grids with the same preset side length value, taking the product value of the average value and the preset side length value as the reference coefficient; taking the preset side length value with the largest reference coefficient as the optimal side length.

[0052] For the analysis of the reference value coefficient for different preset side lengths, as shown in Table 2: Table 2 See Figure 6 ,Figure 6 The optimal side length decision analysis diagram provided by an embodiment of the present invention constructs a two-dimensional rectangular coordinate system with the preset side length value as the abscissa and the reference coefficient as the ordinate. The optimal side length is determined by the numerical values of the reference coefficients under different preset side length values.

[0053] Sort all the initial grids with the same preset side length value to obtain the average value of the reference value coefficients. This average value can initially characterize the reliability of the division. However, since the smaller the side length, the more complex the overall calculation and the smaller the reference significance, the preset side length value is also used as a parameter and supplemented to the overall analysis link. The product value of the average value and the preset side length value is used as the reference coefficient, and the reference coefficient is used as the preferred analysis coefficient for the final preset side length value. The preset side length value with the largest reference coefficient is used as the optimal side length.

[0054] S104: Implement the grid division based on the optimal side length and perform the engineering quantity measurement.

[0055] In the embodiment of the present invention, after determining the optimal side length, the grid method division is performed based on the optimal side length. Refer to Figure 7 , Figure 7 which is the grid division schematic diagram provided by an embodiment of the present invention; the engineering quantity measurement and allocation of different grids are realized based on the final division result.

[0056] The present invention determines the operation area by constructing a three-dimensional model of the site to be measured and combining it with the designed elevation plane. The three-dimensional model is divided by different preset side length values to obtain an initial grid. By analyzing different preset side length values, a more suitable grid division can be selected, so that when calculating the engineering quantity, the reliability is stronger. Then, at any side of any initial grid, different analysis points are periodically obtained; according to the model height difference of the analysis points, the mutation points during the construction process are screened; according to the area distribution and height design change of the operation area where the mutation points are located, the construction intensity correction value is determined. Through two dimensions of area distribution and height distribution, the construction intensity correction analysis at the position of the mutation point is realized. The construction intensity correction value can represent the terrain complexity value at the corresponding mutation point position. After that, according to the number of mutation points on each side of the initial grid, the construction intensity correction value of the mutation points, and the number of zero points, the reference value coefficient of the initial grid is determined. Combining all the reference value coefficients with the same preset side length value, the optimal side length is determined, that is, the analysis is carried out from the mutation point to the edge of the grid where the mutation point is located, and further to the complexity of the initial grid itself, the reference value coefficient of the initial grid is obtained, and the optimal side length is analyzed to realize the grid division and calculate the engineering quantity. Compared with directly calculating through the corner points of a fixed grid, in the embodiment of the present invention, multiple points on the grid sides are analyzed to determine a more reliable grid division, which can avoid the complexity of the calculation process. At the same time, the adaptability of the grid division is improved, and the calculation accuracy is enhanced.

[0057] It should be noted that the above sequence of embodiments of the present invention is only for description and does not represent the advantages or disadvantages of the embodiments. The processes depicted in the drawings do not necessarily require the specific order or continuous order shown to achieve the desired result. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0058] Each embodiment in this specification is described in a progressive manner. The same or similar parts between each embodiment can be referred to each other. Each embodiment focuses on the differences from other embodiments.

Claims

1. A method for calculating the earthwork volume of water conservancy projects based on BIM, characterized in that, The method includes: Construct a three-dimensional model of the site to be measured and calculated, and combine it with the designed elevation plane to determine the operation area. The operation area includes an excavation area and a filling area. The plane line where the designed elevation plane overlaps with the three-dimensional model plane is the zero-point line; divide the three-dimensional model based on the grid method and different preset side lengths to obtain the initial grid; Periodically obtain different analysis points on any side of any initial grid; screen out the mutation points during the construction process according to the model height difference of the analysis points; determine the construction intensity correction value according to the area distribution and height design change of the operation area where the mutation points are located; Take the intersection point of the initial grid side and the zero-point line as the zero point; determine the reference value coefficient of the initial grid according to the number of mutation points on each side of the initial grid, the construction intensity correction value of the mutation points, and the number of zero points, and combine all the reference value coefficients with the same preset side length value to determine the optimal side length; Implement grid division based on the optimal side length and conduct engineering quantity measurement.

2. The method for calculating the earthwork volume of a water conservancy project based on BIM according to claim 1, characterized in that Regard the three-dimensional model area that is not on the same horizontal plane as the designed elevation plane as the operation area.

3. A method for calculating the earthwork volume of a water conservancy project based on BIM according to claim 1, characterized in that, Screen out the mutation points during the construction process according to the model height difference of the analysis points, including: Confirm the other two analysis points closest to any analysis point on the initial grid side as adjacent points; Calculate the absolute value of the model height difference between the analysis point and the adjacent points respectively, average and normalize it to obtain the mutation degree; Regard the analysis points with the mutation degree greater than the preset mutation threshold as mutation points.

4. The method for calculating the earthwork volume of a water conservancy project based on BIM according to claim 1, characterized in that, Determine the construction intensity correction value according to the area distribution and height design change of the operation area where the mutation points are located, including: Determine the area influence index according to the operation area where the mutation points are located, the initial grid, and the area comparison of the operation area within the initial grid; Analyze the height to be constructed in the operation area where the mutation points are located within the initial grid to determine the height influence index; Combine the height influence index and the area influence index to determine the construction intensity correction value.

5. The method for calculating the earthwork volume of a water conservancy project based on BIM according to claim 4, wherein, Determine the area influence index according to the operation area where the mutation points are located, the initial grid, and the area comparison of the operation area within the initial grid, including: Take the area of the operation area where the mutation points are located within the initial grid as the first area, the area of the operation area where the mutation points are located as the second area, and the area of the initial grid as the third area; Calculate the ratio of the first area to the second area as the first ratio, and the ratio of the first area to the third area as the second ratio respectively; Take the sum value of the first ratio and the second ratio as the area influence index.

6. The method for calculating the earthwork volume of a water conservancy project based on BIM according to claim 4, characterized in that, Analyze the height to be constructed in the operation area where the mutation points are located within the initial grid to determine the height influence index, including: Average and determine a preset number of height points in the operation area where the mutation points are located within the initial grid; Take the absolute value of the difference between the model height and the height of the designed elevation plane at the same point as the construction height; Take the average value of the construction heights of all height points as the height influence index.

7. A method for calculating the earthwork volume of a water conservancy project based on BIM according to claim 4, characterized in that, Combine the height influence index and the area influence index to determine the construction intensity correction value, including: Calculate the product value of the area influence index and the height influence index, and normalize it as the construction intensity correction value of the mutation points.

8. A method for calculating the earthwork volume of a water conservancy project based on BIM according to claim 1, characterized in that, Determine the reference value coefficient of the initial grid based on the number of mutation points on each side of the initial grid, the construction intensity correction value of the mutation points, and the number of zero points, including: Take the mean value of the construction intensity correction values of all mutation points on the same side as the edge intensity coefficient; Take the ratio of the number of zero points on the same side to the total number of zero points on the four sides of the initial grid as the zero point proportion; Calculate the product value of the edge intensity coefficient and the zero point proportion as the edge reference index on the initial grid side; Normalize the opposite number of the sum of the edge reference indexes of all sides of the initial grid as the reference value coefficient.

9. The method for calculating the earthwork volume of a water conservancy project based on BIM according to claim 1, characterized in that Combine all reference value coefficients with the same preset side length value to determine the optimal side length, including: Calculate the mean value of the reference value coefficients of all initial grids with the same preset side length value, and take the product value of the mean value and the preset side length value as the reference coefficient; Take the preset side length value with the largest reference coefficient as the optimal side length.

10. A method for calculating the earthwork volume of a water conservancy project based on BIM according to claim 1, characterized in that, The preset side length values include 1 meter, 3 meters, and 5 meters.

Images