A method, apparatus, and storage medium for calculating the volume of earthwork for site leveling.

Abstract

This application discloses a method, apparatus, and storage medium for calculating earthwork volume in site leveling. The method includes: dividing the site leveling area into cells and obtaining earthwork data for each cell, including information on the soil within the area of ​​each cell from the top elevation of different soil layers to the site leveling design elevation; constructing a three-dimensional geological model and correcting the excavation and filling volumes of each cell in the earthwork data based on the three-dimensional geological model; calculating the excavation and filling volumes of each soil layer based on the corrected excavation and filling volumes, and for each soil layer, accumulating the excavation volumes of all cells to obtain the total excavation volume of that soil layer; modifying the site leveling design elevation in the site leveling design data, recalculating, and obtaining multiple sets of earthwork volume data; calculating the project cost based on the multiple sets of earthwork volume data, and selecting the optimal earthwork volume data based on the project cost. This application improves the accuracy of earthwork volume calculation by interpolating existing layer points to construct a refined three-dimensional geological model and correcting the designed excavation and filling volumes.

Description

Technical Field

[0001] This invention belongs to the field of earthwork engineering technology, and specifically relates to a method, device and storage medium for calculating the volume of earthwork for site leveling. Background Technology

[0002] Site leveling is fundamental to the implementation of a construction project, and the cost of site leveling is directly related to the calculation of earthwork volume. According to current budget calculation methods, the cost of earthwork must be priced according to the different properties of the soil.

[0003] However, in the current earthwork volume calculation, due to the influence of engineering cost control, the number of geological boreholes is limited, the three-dimensional geological model is relatively discrete, there are many abrupt changes in the distribution of strata in different regions, and the accuracy of the geological model is low. Since the site leveling design data is determined based on the three-dimensional geological model, the insufficient accuracy of the three-dimensional geological model leads to the inaccuracy of the site leveling design data, and therefore the designed excavation and filling volumes are not accurate enough.

[0004] Moreover, construction is usually carried out based on the initial design of the excavation and filling volumes, without considering factors such as project cost to optimize the excavation and filling volumes. As a result, the excavation and filling scheme adopted may not be the optimal one. Summary of the Invention

[0005] To address the aforementioned issues, the present invention aims to develop an optimized earthwork calculation program for site leveling projects. This program combines 3D terrain data processing technology with geological interface objects to rapidly calculate earthwork volume for site leveling, determine the design elevation for cut-fill balance, and minimize excess soil transportation. It quickly calculates earthwork volume, visually displays the soil properties and bearing layer conditions of each foundation layer, improves the accuracy of earthwork calculation, and simultaneously enhances the efficiency of the design project.

[0006] A method for calculating the volume of earthwork for site leveling includes the following steps:

[0007] Step S1: Divide the site leveling area into cells and obtain the earthwork data for each cell based on the site leveling design data. The earthwork data is the information of the soil within the area of ​​each cell from the top elevation of different soil layers to the site leveling design elevation, including cell number, area, soil layer elevation, site leveling design elevation, and earthwork volume. The earthwork volume is either fill volume or excavation volume.

[0008] Step S3: Construct a three-dimensional geological model and correct the earthwork volume of each cell in the earthwork data based on the three-dimensional geological model.

[0009] Step S4: Calculate the earthwork volume of each soil layer based on the corrected earthwork volume. For fill areas, the fill volume is the fill volume above ground level. The total fill volume is obtained by accumulating the fill volume of all cells. For excavation areas, starting from the ground level downwards, the excavation volume of each layer is the excavation volume of that layer minus the excavation volume of the next layer. For each soil layer, the total excavation volume of that layer is obtained by accumulating the excavation volume of all cells.

[0010] Step S5: Modify the site leveling design elevation in the site leveling design data, return to step S1, and obtain multiple sets of earthwork volume data;

[0011] Step S6: Calculate the project cost based on multiple sets of earthwork volume data, and select the optimal earthwork volume data based on the project cost.

[0012] Optional, a. If transportation is inconvenient and costs are high, adopt the site leveling design elevation used during cut-fill balance;

[0013] b. If transportation is convenient but soil extraction is inconvenient due to terrain, adopt a site leveling design elevation that results in zero fill volume;

[0014] c. If transportation is convenient and soil extraction is unrestricted, adopt the site leveling design elevation with the lowest cost.

[0015] Optionally, for areas within the leveled area with a relative height of less than 200 meters, the side length of the cell is 0.5m; for areas within the leveled area with a relative height of more than 200 meters, the side length of the cell is 0.2m.

[0016] Optionally, step S2 is also included, which checks the validity of the earthwork data. The validity of the earthwork data means that for each row of data, when the site leveling design elevation is higher than the soil layer elevation, the fill volume should be positive and the excavation volume should be zero; when the site leveling design elevation is lower than the soil layer elevation, the fill volume is zero and the excavation volume is positive.

[0017] Optionally, the construction of the three-dimensional geological model and the correction of the earthwork volume in each cell of the earthwork data based on the three-dimensional geological model include:

[0018] Step S31: For each cell, use the interpolation function to perform three-dimensional surface interpolation on the existing top surface of the stratum, and calculate the estimated earthwork volume V within the range from the top surface of each stratum to the site leveling design elevation based on the interpolated data;

[0019] Step S32: In each cell, the estimated earthwork volume V obtained by interpolation is divided by the earthwork volume of each cell in the soil layer data file to obtain the determination coefficient M.

[0020] In step S33, if 0.95≦M≦1.05, the earthwork volume of the corresponding cell is replaced with the calculated estimated earthwork volume V value; otherwise, return to step S31 to re-interpolate with denser insertion points and recalculate the estimated earthwork volume V value until the calculated determination coefficient M satisfies 0.95≦M≦1.05, thereby obtaining the corrected earthwork volume data for each cell.

[0021] Optionally, if a soil layer is missing during the excavation volume calculation, the earthwork volume of that soil layer will be counted as zero.

[0022] Optionally, the griddata function in MATLAB software can be used to perform three-dimensional surface interpolation based on discrete geological layer points.

[0023] Optionally, the step of calculating the estimated earthwork volume V within the range from the top surface of each stratum to the site leveling design elevation based on the interpolated data includes:

[0024] L = sum(sum(abs(Z-Z0)));

[0025] S=abs((y0(2)-y0(1))*(x0(2)-x0(1)));

[0026] V = L * S;

[0027] Among them, Z0 is the site leveling design elevation;

[0028] Z is the height of a point on the top surface of the stratum;

[0029] L is the difference between the top surface height of the cell and the site leveling design elevation;

[0030] x0(2), x0(1), y0(2), and y0(1) are the position coordinates of the top edge of the stratum in the horizontal X-axis and Y-axis directions;

[0031] S is the area of ​​the cell;

[0032] V is the estimated earthwork volume within the range from the top surface of the stratum to the site leveling design elevation.

[0033] This application also provides an electronic device comprising: a memory and a processor, wherein the memory stores a site leveling earthwork volume calculation program, and when the site leveling earthwork volume calculation program is executed by the processor, it implements the site leveling earthwork volume calculation method as described above.

[0034] This application also provides a computer-readable storage medium storing a computer program, the computer program including program instructions, which, when executed by a processor, implement the site leveling earthwork calculation method as described above.

[0035] Compared with the prior art, this application has the following advantages:

[0036] (1) Interpolate the existing layer points using the interpolation function to construct a refined three-dimensional geological model. Based on the three-dimensional geological model, correct the designed excavation and filling volumes to improve the accuracy of earthwork volume calculation.

[0037] (2) Set the site leveling design elevation as a parameter, and obtain multiple sets of earthwork volume calculation results by assigning values ​​to this parameter multiple times. Then, comprehensively evaluate and select the optimal project scheme to improve the project's economy and the design company's industry competitiveness. Attached Figure Description

[0038] Figure 1 This is a partial flowchart of the site leveling earthwork volume calculation method described in an embodiment of the present invention.

[0039] Figure 2 This is a schematic diagram of earthwork data as described in an embodiment of the present invention.

[0040] Figure 3 This is a schematic diagram of the soil layer described in an embodiment of the present invention.

[0041] Figure 4 This is a schematic diagram of the electronic device described in an embodiment of the present invention.

[0042] Figure 5 This is a schematic diagram of the program module described in an embodiment of the present invention.

[0043] Figure 6 This is a schematic diagram of the earthwork data in Excel format as described in an embodiment of the present invention. Detailed Implementation

[0044] The technical solution of the present invention will now be clearly and completely described with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0045] Terminology Explanation:

[0046] Site leveling design elevation refers to the site design elevation at the project location to meet the needs of engineering construction. In engineering design, a certain point is set as ±0:00, and then the relative elevation value of the building is based on this point.

[0047] Soil elevation refers to the elevation of the top surface of different soil layers below the ground level at the project location, based on the properties of the soil.

[0048] Figure 1The main steps in the method for calculating the earthwork volume of site leveling in this embodiment are shown below. In combination with Figure 1 the following, the method for calculating the earthwork volume of site leveling in this embodiment is described, including the following steps:

[0049] Step S1: Divide the site leveling area into square cells, and obtain the earthwork data of each cell according to the site leveling design data. The earthwork data is the information of the soil within the area of each cell from the elevation of the top surface of different soil layers to the site leveling design elevation. It may include cell number, area, soil layer elevation, site leveling design elevation, earthwork volume, and the earthwork volume is the filling volume or excavation volume.

[0050] Among them, the smaller the side length of the cell, the relatively more accurate the calculation result, and the corresponding calculation amount is larger. For projects in plain areas or areas with relatively small undulations in terrain, the side length of the cell can be 0.5 m; for mountainous areas and areas with complex terrain, the side length of the cell can be 0.2 m.

[0051] Specifically, the earthwork data is stored row by row according to the cell number, and each data in each row is separated by ",". Different soil layers are stored in different notepad files, and each notepad is named according to the nature of the soil layer and saved as soil layer data in the corresponding notepad.

[0052] Now an example is given to illustrate the storage format of each data. As Figure 2 shown on the left side in, there are 7 soil layers, which are respectively named plain fill, clay, silt, silty sand, fine sand, medium sand and coarse sand according to different soil properties. The file naming can add serial numbers 1, 2, 3,..., 7 in sequence according to the vertical distribution of the soil layer from the ground downwards. Figure 2 Shown on the right side in is the data in the file named 1-plain fill.txt.

[0053] It can be seen from it that vertically is the cell number. Corresponding to each cell number, each row is the information of a cell. The data in the first column is the cell number, and the subsequent columns are area, site elevation, design elevation, filling volume and excavation volume in sequence. The data in each column is separated by "," and no symbol is added after the data in the last column.

[0054] Step S2: Check the validity of the earthwork data, that is, for each row of data, when the site leveling design elevation is higher than the soil layer elevation, the filling volume should be a positive value and the excavation volume should be zero; when the site leveling design elevation is lower than the soil layer elevation, the filling volume is zero and the excavation volume is a positive value; in other cases, output the prompt message "It is recommended to reduce the unit area. There is a semi-fill and semi-excavation area in the calculation area, which affects the calculation accuracy".

[0055] After all the data is confirmed to be correct, it is stored in a public type list variable in sequence according to the cell number for subsequent calculation calls.

[0056] Step S3: Construct a three-dimensional geological model and correct the excavation and filling volumes of each cell in the earthwork data based on the three-dimensional geological model.

[0057] Specifically, due to the influence of engineering cost control, the number of geological boreholes is limited, the three-dimensional geological model is relatively discrete, there are many abrupt changes in the distribution of strata in different regions, and the accuracy of the geological model is low. Since the site leveling design data is determined based on the three-dimensional geological model, the insufficient accuracy of the three-dimensional geological model leads to the inaccuracy of the site leveling design data. Therefore, it is necessary to perform statistical analysis, interpolation and fitting on the existing geological data to construct a refined three-dimensional geological model, and then correct the excavation and filling volumes based on the refined three-dimensional geological model.

[0058] Specifically, it includes the following steps:

[0059] Step S31: For each cell, use an interpolation function to perform three-dimensional surface interpolation on the existing top surface of the strata, and calculate the earthwork volume within the range from the top surface of each stratum to the site leveling design elevation based on the interpolated data. For example, the griddata function in MATLAB software can be used to achieve three-dimensional surface interpolation based on several discrete geological layer points. The main calculation process is as follows:

[0060] clear,clc

[0061] %%Standard Case of Discrete Point Surface Interpolation

[0062] data=readmatrix('data.xlsx');

[0063] x = data(:,1);

[0064] y = data(:,2);

[0065] z = data(:,3);

[0066] [X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x))',linspace(min(y),max(y)),'cubic'); % Triangle-based cubic interpolation only supports interpolation.

[0067] %

[0068] [X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x))',linspace(min(y),max(y)),'v4'); % Biharmonic spline interpolation supports extrapolation

[0069] figure,contourf(X,Y,Z)%Contour plot

[0070] figure,surf(X,Y,Z); % 3D surface

[0071] After interpolation, the estimated earthwork volume is calculated within each cell range from the top surface of each stratum to the site leveling design elevation. The method is as follows:

[0072] L = sum(sum(abs(Z-Z0)));

[0073] S=abs((y0(2)-y0(1))*(x0(2)-x0(1)));

[0074] V = L * S; % Volume

[0075] Where sum(sum(abs(Z-Z0))) is the integral of the difference between the top surface height of the stratum and the site leveling design elevation;

[0076] Z0 is the site leveling design elevation;

[0077] Z is the height of a point on the top surface of the stratum;

[0078] x0(2), x0(1), y0(2), and y0(1) are the position coordinates of the top edge of the stratum in the horizontal X-axis and Y-axis directions. Specifically, x0(1): the x-coordinate of the bottom left point of the cell; y0(1): the y-coordinate of the bottom left point of the cell; x0(2): the x-coordinate of the top right point of the cell; y0(2): the y-coordinate of the top right point of the cell.

[0079] L is the difference between the top surface height of the cell and the site leveling design elevation;

[0080] S is the area of ​​the cell;

[0081] V is the estimated earthwork volume within the range from the top surface of the stratum to the site leveling design elevation.

[0082] In step S32, the estimated earthwork volume V obtained through interpolation in each cell is divided by the earthwork volume of each cell in the soil layer data file to obtain the determination coefficient M. In other words, a determination coefficient M will be calculated for each cell.

[0083] In step S33, if 0.95≦M≦1.05, it means that the cell size accuracy meets the engineering requirements. Then, the fill or excavation volume of that cell in the soil layer data file is replaced with the calculated estimated earthwork volume V value. Otherwise, the cell size is reduced, that is, the process returns to step S31 to re-interpolate with denser insertion points and recalculate the estimated earthwork volume V value until the calculated determination coefficient M satisfies 0.95≦M≦1.05. In this way, the earthwork volume of all cells is replaced with the estimated earthwork volume V value. The fill and excavation volumes used in the subsequent calculations are all revised data.

[0084] Step S34: Sort the soil layers according to the elevation of each soil layer in the cell, and store the soil layer area, soil layer elevation, site leveling design elevation, earthwork volume and other information in a public list variable called CellModel.

[0085] Step S4: Calculate the excavation and fill volumes, where for the fill area, the fill volume is the fill volume of the first soil layer ( Figure 3 The F1 key on the left accumulates the fill volume of all cells to obtain the total fill volume for that soil layer. For excavation areas, starting from the ground surface downwards, the excavation volume of each layer is the excavation volume of that layer minus the excavation volume of the next layer, such as... Figure 3 The excavation volume of the first soil layer on the right is F1-F2; the excavation volume of the last soil layer is Fn. For each soil layer, the excavation volume of all cells is accumulated to obtain the total excavation volume of that soil layer.

[0086] For example, attached Figure 3 Assuming the strata are divided into 5 layers based on different soil properties, calculate the excavation volume of the i-th cell. The excavation volume of soil layer 1 in this cell is: Ti1 = (F1 - F2). Summing Ti1 in all cells within the project scope yields... T1 represents the excavation volume of soil layer 1. Similarly, the excavation volumes of soil layers 2, 3, ..., 5 can be obtained.

[0087] Specifically, during the excavation calculation, if a soil layer is missing (i.e., the data in that row is blank), the excavation and fill volume of that soil layer is counted as zero and stored in the CellModel.

[0088] Step S5: Modify the site leveling design elevation in the site leveling design data, return to step S1, and obtain multiple sets of earthwork volume data. Each set of earthwork volume data includes the total excavation or filling volume of each stratum. The multiple sets of earthwork data can be selected as needed; for example, repeating steps S1 to S5 to obtain 5 sets of earthwork data.

[0089] Step S6: Calculate the project cost based on multiple sets of earthwork volume data, and select the optimal earthwork volume data based on the project cost.

[0090] Specifically, the site leveling design elevation can be set as a parameter. By assigning values ​​to this parameter multiple times, multiple sets of earthwork volume calculation results can be obtained, and thus the corresponding project cost can be obtained. This can then be used for project scheme selection.

[0091] a. If transportation is inconvenient and costs are high, the site elevation can be adopted when the fill volume and cut volume are as close as possible (i.e., fill-cut balance);

[0092] b. If transportation is convenient but soil extraction is inconvenient due to terrain, a scheme with a fill volume as close to 0 as possible can be adopted.

[0093] c. If transportation is convenient and soil extraction is unrestricted, the lowest-cost option can be chosen.

[0094] This application also provides an electronic device, Figure 4 This is a schematic diagram of the hardware architecture of an embodiment of the electronic device of the present invention. In this embodiment, the electronic device 2 is a device capable of automatically performing numerical calculations and / or information processing according to pre-set or stored instructions. For example, it can be a smartphone, tablet computer, laptop computer, desktop computer, server, etc. Figure 4 As shown, the electronic device 2 includes, but is not limited to, a memory 21 and a processor 22 that are communicatively connected to each other via a system bus. The memory 21 includes at least one type of computer-readable storage medium, including flash memory, hard disk, multimedia card, card-type memory (e.g., SD or DX memory), random access memory (RAM), static random access memory (SRAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), programmable read-only memory (PROM), magnetic memory, magnetic disk, optical disk, etc. In some embodiments, the memory 21 may be an internal storage unit of the electronic device 2, such as a hard disk or RAM of the electronic device 2. In other embodiments, the memory 21 may also be an external storage device of the electronic device 2, such as a plug-in hard disk, smart media card (SMC), secure digital (SD) card, flash card, etc. Of course, the memory 21 may also include both internal storage units and external storage devices of the electronic device 2. In this embodiment, the memory 21 is typically used to store the operating system and various application software installed on the electronic device 2, such as the code for calculating the earthwork volume of site leveling. Furthermore, the memory 21 can also be used to temporarily store various types of data that have been output or will be output.

[0095] In some embodiments, the processor 22 may be a central processing unit (CPU), controller, microcontroller, microprocessor, or other data processing chip. The processor 22 is typically used to control the overall operation of the electronic device 2, such as performing control and processing related to data interaction or communication with the electronic device 2. In this embodiment, the processor 22 is used to run program code stored in the memory 21 or process data, such as running the site leveling earthwork calculation program.

[0096] Optionally, the electronic device 2 may further include a display, which may also be called a display screen or display unit. In the display screen, the program interface is divided into three areas: a data input area, a calculation area, and a result display area. The result display area is used to display the calculation results. Clicking the calculation button completes the calculation process, and the calculation result is automatically displayed in the display area. Figure 5 As shown, the display area lists in detail the fill and cut volumes for each soil layer throughout the entire project. This data can be directly used in the project's preliminary budget calculations.

[0097] It should be pointed out that, Figure 4 Only the electronic device 2 with components 21-22 is shown; however, it should be understood that it is not required to implement all of the components shown, and more or fewer components may be implemented instead.

[0098] The memory 21, containing a readable storage medium, may include an operating system, a site leveling earthwork volume calculation program 50, etc. When the processor 22 executes the site leveling earthwork volume calculation 50 stored in the memory 21, it implements the steps described in the site leveling earthwork volume calculation method above. In this embodiment, the site leveling earthwork volume calculation program stored in the memory 21 can be divided into one or more program modules. The one or more program modules are stored in the memory 21 and can be executed by one or more processors (processor 22 in this embodiment) to complete the present invention.

[0099] For example, Figure 5 A schematic diagram of the program modules for the site leveling earthwork volume calculation program is shown. In this embodiment, the site leveling earthwork volume calculation program 50 can be divided into an earthwork data input module 501, an earthwork volume correction module 501, an earthwork volume calculation module 503, and an optimal earthwork volume data determination module 504. Here, the program module, as referred to in this invention, is a series of computer program instruction segments capable of performing specific functions, and is more suitable than a program for describing the execution process of the site leveling earthwork volume calculation program in the electronic device 2. The following description will specifically introduce the specific functions of the program modules.

[0100] Among them, the earthwork data input module 501 is used to execute step S1, divide the site leveling area into square cells, and obtain the earthwork data of each cell according to the site leveling design data. The earthwork data is the information of the soil within the area of each cell from the elevation of the top surface of different soil layers to the site leveling design elevation. It may include cell number, area, soil layer elevation, site leveling design elevation, and earthwork volume, where the earthwork volume is the filling volume or the excavation volume.

[0101] Among them, the smaller the side length of the cell, the relatively more accurate the calculation result, and the corresponding larger the calculation amount. For projects in plain areas or areas with relatively small undulations in terrain, the side length of the cell can be 0.5 m; for mountainous areas and areas with complex terrain, the side length of the cell can be 0.2 m.

[0102] Specifically, the earthwork data is stored row by row according to the cell number, and each data in each row is separated by ",". Different soil layers are stored in different notepad files, and each notepad is named according to the nature of the soil layer and saved in the corresponding notepad as the soil layer data.

[0103] Now an example is given to illustrate the storage format of each data. As Figure 2 shown on the left side, there are 7 soil layers in total, which are named plain fill, clay, silt, silty sand, fine sand, medium sand, and coarse sand respectively according to different soil properties. The file naming can add serial numbers 1, 2, 3,..., 7 in sequence according to the vertical distribution of the soil layer from the ground downwards. Figure 2 Shown on the right side is the data in the file named 1-plain fill.txt.

[0104] It can be seen from it that vertically is the cell number. Corresponding to each cell number, each row is the information of a cell. The first column data is the cell number, and the subsequent columns are the area, site elevation, design elevation, filling volume, and excavation volume in sequence. Each column of data is separated by ",", and no symbol is added after the last column of data.

[0105] Step S2, check the validity of the earthwork data, that is, for each row of data, when the site leveling design elevation is higher than the soil layer elevation, the filling volume should be positive and the excavation volume should be zero; when the site leveling design elevation is lower than the soil layer elevation, the filling volume is zero and the excavation volume is positive; in other cases, output the prompt message "It is recommended to reduce the unit area. There is a half-fill and half-excavation area in the calculation area, which affects the calculation accuracy".

[0106] When all the data is confirmed to be correct, it is stored in a public type list variable in sequence according to the cell number for subsequent calculation calls.

[0107] Among them, the earthwork volume correction module 502 is used to execute step S3, construct a three-dimensional geological model, and correct the excavation volume and filling volume of each cell in the earthwork data according to the three-dimensional geological model.

[0108] Specifically, due to the influence of engineering cost control, the number of geological boreholes is limited, the three-dimensional geological model is relatively discrete, there are many abrupt changes in the distribution of strata in different regions, and the accuracy of the geological model is low. Since the site leveling design data is determined based on the three-dimensional geological model, the insufficient accuracy of the three-dimensional geological model leads to the inaccuracy of the site leveling design data. Therefore, it is necessary to perform statistical analysis, interpolation and fitting on the existing geological data to construct a refined three-dimensional geological model, and then correct the excavation and filling volumes based on the refined three-dimensional geological model.

[0109] Specifically, it includes the following steps:

[0110] Step S31: For each cell, use an interpolation function to perform three-dimensional surface interpolation on the existing top surface of the strata, and calculate the earthwork volume within the range from the top surface of each stratum to the site leveling design elevation based on the interpolated data. For example, the griddata function in MATLAB software can be used to achieve three-dimensional surface interpolation based on several discrete geological layer points. The main calculation process is as follows:

[0111] clear,clc

[0112] %%Standard Case of Discrete Point Surface Interpolation

[0113] data=readmatrix('data.xlsx');

[0114] x = data(:,1);

[0115] y = data(:,2);

[0116] z = data(:,3);

[0117] [X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x))',linspace(min(y),max(y)),'cubic'); % Triangle-based cubic interpolation only supports interpolation.

[0118] %

[0119] [X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x))',linspace(min(y),max(y)),'v4'); % Biharmonic spline interpolation supports extrapolation

[0120] figure,contourf(X,Y,Z)%Contour plot

[0121] figure,surf(X,Y,Z); % 3D surface

[0122] After interpolation, the estimated earthwork volume is calculated within each cell range from the top surface of each stratum to the site leveling design elevation. The method is as follows:

[0123] L = sum(sum(abs(Z-Z0)));

[0124] S=abs((y0(2)-y0(1))*(x0(2)-x0(1)));

[0125] V = L * S; % Volume

[0126] Where sum(sum(abs(Z-Z0))) is the integral of the difference between the top surface height of the stratum and the site leveling design elevation;

[0127] Z0 is the site leveling design elevation;

[0128] Z is the height of a point on the top surface of the stratum;

[0129] x0(2), x0(1), y0(2), and y0(1) are the position coordinates of the top edge of the stratum in the horizontal X-axis and Y-axis directions. Specifically, x0(1): the x-coordinate of the bottom left point of the cell; y0(1): the y-coordinate of the bottom left point of the cell; x0(2): the x-coordinate of the top right point of the cell; y0(2): the y-coordinate of the top right point of the cell.

[0130] L is the difference between the top surface height of the cell and the site leveling design elevation;

[0131] S is the area of ​​the cell;

[0132] V is the estimated earthwork volume within the range from the top surface of the stratum to the site leveling design elevation.

[0133] In step S32, the estimated earthwork volume V obtained through interpolation in each cell is divided by the earthwork volume of each cell in the soil layer data file to obtain the determination coefficient M. In other words, a determination coefficient M will be calculated for each cell.

[0134] In step S33, if 0.95≦M≦1.05, it means that the cell size accuracy meets the engineering requirements. Then, the fill or excavation volume of that cell in the soil layer data file is replaced with the calculated estimated earthwork volume V value. Otherwise, the cell size is reduced, that is, the process returns to step S31 to re-interpolate with denser insertion points and recalculate the estimated earthwork volume V value until the calculated determination coefficient M satisfies 0.95≦M≦1.05. In this way, the earthwork volume of all cells is replaced with the estimated earthwork volume V value. The fill and excavation volumes used in the subsequent calculations are all revised data.

[0135] Step S34: Sort the soil layers according to the elevation of each soil layer in the cell, and store the soil layer area, soil layer elevation, site leveling design elevation, earthwork volume and other information in a public list variable called CellModel.

[0136] Earthwork volume calculation module 503 is used to execute step S4, calculating the excavation and fill volumes, where for the fill area, the fill volume is the fill volume of the first soil layer ( Figure 3 The F1 key on the left accumulates the fill volume of all cells to obtain the total fill volume for that soil layer. For excavation areas, starting from the ground surface downwards, the excavation volume of each layer is the excavation volume of that layer minus the excavation volume of the next layer, such as... Figure 3 The excavation volume of the first soil layer on the right is F1-F2; the excavation volume of the last soil layer is Fn. For each soil layer, the excavation volume of all cells is accumulated to obtain the total excavation volume of that soil layer.

[0137] For example, attached Figure 3 Assuming the strata are divided into 5 layers based on different soil properties, for the first soil layer, calculate the excavation volume of the i-th cell. The excavation volume of soil layer 1 within this cell is: Ti1 = (F1 - F2). Summing Ti1 in all cells within the project scope yields... m indicates that the soil layer is divided into m cells. T1 is the excavation volume of the first soil layer, and similarly, the excavation volumes of the second, ..., fifth soil layers can be obtained. This is exactly the data required for the project.

[0138] Specifically, during the excavation calculation, if a soil layer is missing (i.e., the data in that row is blank), the excavation and fill volume of that soil layer is counted as zero and stored in the CellModel.

[0139] Step S5: Modify the site leveling design elevation in the site leveling design data, return to step S1, and obtain multiple sets of earthwork volume data. Each set of earthwork volume data includes the total excavation or filling volume of each stratum. The multiple sets of earthwork data can be selected as needed; for example, repeating steps S1 to S5 to obtain 5 sets of earthwork data.

[0140] The optimal earthwork volume data determination module 504 is used to execute step S6, calculate the project cost based on multiple sets of earthwork volume data, and select the optimal earthwork volume data based on the project cost.

[0141] Specifically, the site leveling design elevation can be set as a parameter. By assigning values ​​to this parameter multiple times, multiple sets of earthwork volume calculation results can be obtained, and thus the corresponding project cost can be obtained. This can then be used for project scheme selection.

[0142] a. If transportation is inconvenient and costs are high, the site elevation can be adopted when the fill volume and cut volume are as close as possible (i.e., fill-cut balance);

[0143] b. If transportation is convenient but soil extraction is inconvenient due to terrain, a scheme with a fill volume as close to 0 as possible can be adopted.

[0144] c. If transportation is convenient and soil extraction is unrestricted, the lowest-cost option can be chosen.

[0145] The following describes the calculation process using the electronic device described in this application.

[0146] Step 1: Prepare data

[0147] Step 1.1 If there is Figure 1 For files formatted on the right, rename them according to the format "serial number-soil layer type.txt" and save them, for example, "1-plain fill soil.txt";

[0148] If step 1.2 is not available Figure 1 The file format on the right can also be in Excel format (e.g., ...). Figure 6 Compile a data file, renaming it to "Serial Number-Soil Layer Type.txt" and saving it, for example, "1-Plain Fill.txt"; the file format is .xlsx; in each sheet of the file; A1 contains "Cut and Fill Calculation"; A2 contains "Project Name:"; A3 to K3 contain "Number", "Area (㎡)", "Original Ground Elevation (m)", "Design Elevation (m)", "Excavation (m3)", and "Cut and Fill (m3)" respectively. Starting from row 4, the vertical column represents cell numbers, corresponding to each cell number. Each row contains information for one cell. The first column of data is the cell number, and the following columns are, in order, the area, site elevation, design elevation, fill volume, and excavation volume. The data in each column is separated by commas, and no symbols are added after the last column of data.

[0149] Step 2: Click "Start Program", select the prepared data file, and drag the mouse to the data file list box;

[0150] Step 3: The program iterates through each input file in turn, reads the data from the file and saves it in double DataFile(i), then sorts the data according to the soil elevation and stores it in a public list variable.

[0151] Step 4: The program iterates through each calculation cell to determine whether the cell is a fill area or a cut area.

[0152] Step 4.1 For embankment sections, the fill volume is the fill volume from the top surface of the first soil layer to the ground surface. The fill type is an open data output for designers to select as needed. During the cut calculation, if a soil layer is missing, the cut and fill volume of that soil layer is counted as zero and stored in CellModel.

[0153] Step 4.2 For cut sections, starting from the ground surface downwards, the cut volume for each layer is the cut volume of that layer minus the cut volume of the next layer, with the last layer using its cut volume. For each layer, the cut volume of all cells is accumulated to obtain the total cut volume for that layer. During the cut calculation process, if a layer is missing, the cut and fill volume for that layer is counted as zero and stored in the CellModel.

[0154] Step 4.3 Output the calculated data to "Earthwork Calculation Process Results.txt".

[0155] Step 5: Traverse each soil layer in each calculation unit and accumulate the fill and cut volumes for each soil layer.

[0156] Step 6 Output Results

[0157] Step 6.1 Output the final calculation results to "Earthwork Final Result.txt";

[0158] Step 6.2 Output the final calculation result to the "Calculation Result" output area of ​​the main program interface.

[0159] To facilitate use by designers, a drag-and-drop control, List1_OLEDragDro, was designed. It supports drag-and-drop processing of multiple data files and can process up to 300 soil layer data files at a time. The file names are stored in Public NameofDwgFile(1To 300)As String.

[0160] Of course, the present invention may have other various embodiments. Without departing from the spirit and essence of the present invention, those skilled in the art can make various corresponding changes and modifications according to the present invention, but these corresponding changes and modifications are all within the protection scope of the claims of the present invention.

Claims

1. A method for calculating the volume of earthwork for site leveling, characterized in that, Includes the following steps: Step S1: Divide the site leveling area into cells and obtain the earthwork data for each cell based on the site leveling design data. The earthwork data is the information of the soil within the area of ​​each cell from the top elevation of different soil layers to the site leveling design elevation, including cell number, area, soil layer elevation, site leveling design elevation, and earthwork volume. The earthwork volume is either fill volume or excavation volume. Earthwork data is stored in rows according to cell number, with each data in each row separated by commas; different soil layers are stored in different notepad files, and each notepad is named according to the properties of the soil layer, with the soil layer data saved in the corresponding notepad. Step S2: Check the validity of the earthwork data. The validity of the earthwork data means that for each row of data, when the site leveling design elevation is higher than the soil layer elevation, the fill volume should be positive and the excavation volume should be zero; when the site leveling design elevation is lower than the soil layer elevation, the fill volume should be zero and the excavation volume should be positive. Step S3: Construct a three-dimensional geological model, and correct the earthwork volume of each cell in the earthwork data based on the three-dimensional geological model, including: Step S31: For each cell, use the interpolation function to perform three-dimensional surface interpolation on the existing top surface of the stratum, and calculate the estimated earthwork volume V within the range from the top surface of each stratum to the site leveling design elevation based on the interpolated data; Step S32: In each cell, the estimated earthwork volume V obtained by interpolation is divided by the earthwork volume of each cell in the soil layer data file to obtain the determination coefficient M. In step S33, if 0.95≤M≤1.05, the earthwork volume of the corresponding cell is replaced with the calculated estimated earthwork volume V value; otherwise, return to step S31 to re-interpolate with denser insertion points and recalculate the estimated earthwork volume V value until the calculated determination coefficient M satisfies 0.95≤M≤1.05, thereby obtaining the corrected earthwork volume data for each cell. Step S4: Calculate the earthwork volume of each soil layer based on the corrected earthwork volume. For fill areas, the fill volume is the fill volume above ground level. The total fill volume is obtained by accumulating the fill volume of all cells. For excavation areas, starting from the ground level downwards, the excavation volume of each layer is the excavation volume of that layer minus the excavation volume of the next layer. For each soil layer, the total excavation volume of that layer is obtained by accumulating the excavation volume of all cells. Step S5: Modify the site leveling design elevation in the site leveling design data, return to step S1, and obtain multiple sets of earthwork volume data; Step S6: Calculate the project cost based on multiple sets of earthwork volume data, and select the optimal earthwork volume data based on the project cost.

2. The method for calculating the earthwork volume for site leveling according to claim 1, characterized in that, a. If transportation is inconvenient and costs are high, adopt the site leveling design elevation used in cut-fill balance; b. If transportation is convenient but soil extraction is inconvenient due to terrain, adopt a site leveling design elevation that results in zero fill volume; c. If transportation is convenient and soil extraction is unrestricted, adopt the site leveling design elevation with the lowest cost.

3. The method for calculating the earthwork volume for site leveling according to claim 1, characterized in that, For areas within the leveled site with a relative height of less than 200 meters, the side length of the cell is 0.5m; for areas within the leveled site with a relative height of more than 200 meters, the side length of the cell is 0.2m.

4. The method for calculating the earthwork volume for site leveling according to claim 1, characterized in that, If a soil layer is missing during the excavation volume calculation process, the earthwork volume of that soil layer is counted as zero.

5. The method for calculating the earthwork volume for site leveling according to claim 1, characterized in that, Using the griddata function in MATLAB software, three-dimensional surface interpolation is achieved based on discrete geological layer points.

6. The method for calculating the earthwork volume for site leveling according to claim 1, characterized in that, The calculation of the estimated earthwork volume V within the range from the top surface of each stratum to the site leveling design elevation based on the interpolated data includes: L = sum(sum(abs(Z-Z0))); S=abs((y0(2)-y0(1))*(x0(2)-x0(1))); V = L * S; Among them, Z0 is the site leveling design elevation; Z is the height of a point on the top surface of the stratum; L is the difference between the top surface height of the cell and the site leveling design elevation; x0(2), x0(1), y0(2), and y0(1) are the position coordinates of the top edge of the stratum in the horizontal X-axis and Y-axis directions; S is the area of ​​the cell; V is the estimated earthwork volume within the range from the top surface of the stratum to the site leveling design elevation.

7. An electronic device, characterized in that, The electronic device includes a memory and a processor, wherein the memory stores a site leveling earthwork volume calculation program, and when the site leveling earthwork volume calculation program is executed by the processor, it implements the site leveling earthwork volume calculation method as described in any one of claims 1-6.

8. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program, the computer program including program instructions, which, when executed by a processor, implement the site leveling earthwork volume calculation method as described in any one of claims 1-6.

Images