CN112861222A - Earthwork calculation method and device capable of distinguishing soil quality - Google Patents

Earthwork calculation method and device capable of distinguishing soil quality Download PDF

Info

Publication number
CN112861222A
CN112861222A CN202110100004.7A CN202110100004A CN112861222A CN 112861222 A CN112861222 A CN 112861222A CN 202110100004 A CN202110100004 A CN 202110100004A CN 112861222 A CN112861222 A CN 112861222A
Authority
CN
China
Prior art keywords
tensor
current
excavation
point cloud
soil
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110100004.7A
Other languages
Chinese (zh)
Other versions
CN112861222B (en
Inventor
殷达
廖曾平
黎江
黎忠豪
詹锐彪
杨秀武
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CCCC Guangzhou Dredging Co Ltd.
Original Assignee
CCCC Guangzhou Dredging Co Ltd.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CCCC Guangzhou Dredging Co Ltd. filed Critical CCCC Guangzhou Dredging Co Ltd.
Priority to CN202110100004.7A priority Critical patent/CN112861222B/en
Publication of CN112861222A publication Critical patent/CN112861222A/en
Application granted granted Critical
Publication of CN112861222B publication Critical patent/CN112861222B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05—Geographic models
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Computer Graphics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Structural Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Civil Engineering (AREA)
  • Remote Sensing (AREA)
  • Architecture (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention relates to an earthwork calculation method and device capable of distinguishing soil texture. Meanwhile, a mapping rule that the rasterized model is converted into a tensor model is set to calculate the earthwork excavation amount at the current period, so that quantitative analysis of engineering geology is facilitated, related personnel can establish application programs for executing discrete tasks in various fields according to the mapping rule, and the earthwork calculation amount can be accurately and efficiently completed in various projects and terrains.

Description

Earthwork calculation method and device capable of distinguishing soil quality
Technical Field
The invention relates to the technical field of earthwork measurement and calculation, in particular to an earthwork calculation method and device capable of distinguishing soil texture.
Background
Earthwork engineering is one of the main engineering in construction engineering, and includes the steps of cutting earth (stone) side, digging, filling, transportation, drainage, precipitation and the like. The construction cost of the square project is controlled by soil quality categories, and the influence is obvious. Therefore, in the development of earthwork, the excavation work amount needs to be calculated according to various soil qualities, and a reference foundation is provided for calculation of the earthwork amount.
The traditional method for developing the earthwork calculation amount mainly depends on-site observation and sampling, and the method cannot directly realize the classified output of the earthwork engineering amount of each soil texture and is difficult to be applied to hidden engineering (such as tunnel engineering and dredging engineering) which has urgent calculation requirements on the engineering amount data of each soil texture. The calculated earthwork calculation amount can be based on geological exploration drilling information, and a corresponding three-dimensional geological digital model can be established according to relevant specifications and standards. However, the three-dimensional geological digital model still has difficulty in carrying out earth calculation for distinguishing the soil texture, and the method has low calculation efficiency or poor robustness and is not suitable for iterative calculation of large projects or multi-stage landforms.
Therefore, the conventional earthwork calculation method is difficult to perform earthwork calculation according to soil property classification, and the technical difficulty of calculating the earthwork calculation according to the soil property classification is that when the conventional three-dimensional graph method performs Boolean operation on multi-period terrain, the grid quality of a three-dimensional model of each soil property is gradually reduced, so that the technical problems of tolerance and the like are caused. One solution is to reconstruct a three-dimensional model, such as extracting a model point cloud, but it is difficult to ensure the conformity of the common plane of the three-dimensional models of the respective soil properties after reconstruction and the consistency of the positions of the key points. The rasterization method can solve the consistency problem, but the rasterization method has the problems of space complexity and time complexity, and when the rasterization method is applied to actual engineering, the calculation efficiency is difficult to guarantee while the model precision is guaranteed.
In summary, the conventional earth mass calculation method has various drawbacks as described above.
Disclosure of Invention
Therefore, it is necessary to provide a method and an apparatus for earthwork calculation capable of distinguishing the soil quality, aiming at the defects of the conventional earthwork calculation method.
An earthwork calculation method capable of distinguishing soil texture, comprising the steps of:
acquiring a three-dimensional geological digital model, and rasterizing the three-dimensional geological digital model to obtain a discrete grid array;
determining function mapping of each field according to the grid array; wherein the form of the function mapping is a three-dimensional tensor;
extracting sub-tensors corresponding to various soil textures from the three-dimensional tensors in various fields;
respectively obtaining a current topographic survey point cloud and a current topographic survey point cloud according to a time sequence, calculating a current accumulated excavation tensor according to the current topographic survey point cloud fitting, and obtaining a current excavation tensor according to the current accumulated excavation tensor and the current accumulated excavation tensor;
performing point multiplication operation on the sub-tensors of all kinds of soil textures in each field and the current excavation tensor of the corresponding field to obtain the current excavation tensor of each soil texture, and updating the sub-tensors according to the current excavation tensor of each soil texture;
and acquiring the sum of elements of the current-period excavation tensor of various types of soil in each field, and converting the sum of elements into the space occupancy rate of an Euclidean space according to the grid scale to be used as the current-period excavation amount of various types of soil.
The earthwork calculation method capable of distinguishing the soil texture constructs an intermediate data standard through tensor based on the field concept, and solves the problems of tolerance and grid quality reduction existing in the traditional three-dimensional graphic method by rasterizing the three-dimensional geological digital model. Meanwhile, a mapping rule that the rasterized model is converted into a tensor model is set to calculate the earthwork excavation amount at the current period, so that quantitative analysis of engineering geology is facilitated, related personnel can establish application programs for executing discrete tasks in various fields according to the mapping rule, and earthwork calculation amount of various types of soil can be accurately and efficiently completed in various engineering and terrains.
In one embodiment, the process of calculating the cumulative excavation tensor according to the point cloud fitting of the current topographic survey comprises the following steps:
acquiring point cloud original data of a last-stage topographic measurement point cloud;
based on the point cloud original data, generating a curved surface function by applying a three-dimensional curved surface interpolation method, and reconstructing by taking a specific grid scale to generate a new point cloud array; wherein the grid dimension is a grid dimension of the grid array;
and outputting the up-term accumulated excavation tensor based on the new point cloud array.
In one embodiment, the process of tensor operation comprises the steps of:
tensor operations are performed by a sparse matrix compression algorithm.
In one embodiment, the process of tensor operation comprises the steps of:
tensor operations are performed by parallel operations.
In one embodiment, the process of converting the element sum into the space occupancy of euclidean space according to the grid scale as the current earthwork volume of various types of soil includes the steps of:
the current amount of earth excavation is obtained from the product of the element and the cube of the grid dimension of the grid array.
In one embodiment, the grid array has a grid dimension of 1 m.
In one embodiment, the three-dimensional geosynthetic model comprises a BIM three-dimensional geological model.
An earth mass calculating apparatus capable of distinguishing the quality of earth, comprising:
the rasterization processing module is used for acquiring the three-dimensional geological digital model and rasterizing the three-dimensional geological digital model to obtain a discrete grid array;
the three-dimensional tensor determining module is used for determining function mapping of each field according to the grid array; wherein the form of the function mapping is a three-dimensional tensor;
the sub-tensor extraction module is used for extracting sub-tensors corresponding to various soil textures from the three-dimensional tensors in various fields;
the current excavation tensor determining module is used for respectively obtaining a current topographic survey point cloud and a current topographic survey point cloud according to the time sequence, calculating a current accumulated excavation tensor according to the current topographic survey point cloud fitting, and obtaining a current excavation tensor according to the current accumulated excavation tensor and the current accumulated excavation tensor;
the soil texture current excavation tensor determining module is used for performing point multiplication operation on the sub tensors of various soil textures in each field and the current excavation tensor of the corresponding field to obtain the current excavation tensor of each soil texture, and updating the sub tensors according to the current excavation tensor of each soil texture;
and the current-period earthwork excavation amount determining module is used for acquiring the sum of elements of current-period excavation tensor of various types of soil in each field, and converting the sum of elements into the space occupancy rate of an Euclidean space according to the grid scale to be used as the current-period earthwork excavation amount of various types of soil.
The earthwork calculation amount device capable of distinguishing the soil texture constructs the intermediate data standard through tensor based on the field concept, and solves the problems of tolerance and grid quality reduction existing in the traditional three-dimensional graphic method through rasterization processing of the three-dimensional geological digital model. Meanwhile, a mapping rule that the rasterized model is converted into a tensor model is set to calculate the earthwork excavation amount at the current period, so that quantitative analysis of engineering geology is facilitated, related personnel can establish application programs for executing discrete tasks in various fields according to the mapping rule, and the earthwork calculation amount can be accurately and efficiently completed in various projects and terrains.
A computer storage medium having stored thereon computer instructions which, when executed by a processor, implement the method of any of the above embodiments for earth mass distinguishable earthmoving computation.
The computer storage medium constructs an intermediate data standard through tensor based on the field concept, and solves the problems of tolerance and grid quality reduction existing in the traditional three-dimensional graphic method through rasterization processing of the three-dimensional geological digital model. Meanwhile, a mapping rule that the rasterized model is converted into a tensor model is set to calculate the earthwork excavation amount at the current period, so that quantitative analysis of engineering geology is facilitated, related personnel can establish application programs for executing discrete tasks in various fields according to the mapping rule, and the earthwork calculation amount can be accurately and efficiently completed in various projects and terrains.
A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the method of distinguishable earthiness computation of any of the above embodiments when executing the program.
The computer device constructs the intermediate data standard through tensor based on the field concept, and solves the problems of tolerance and grid quality reduction existing in the traditional three-dimensional graphic method through rasterization processing of the three-dimensional geological digital model. Meanwhile, a mapping rule that the rasterized model is converted into a tensor model is set to calculate the earthwork excavation amount at the current period, so that quantitative analysis of engineering geology is facilitated, related personnel can establish application programs for executing discrete tasks in various fields according to the mapping rule, and the earthwork calculation amount can be accurately and efficiently completed in various projects and terrains.
Drawings
FIG. 1 is a flow chart of an earthwork calculation method for differentiating soil properties according to an embodiment;
FIG. 2 is a schematic diagram of tensor mapping for a rasterized geological model;
FIG. 3 is a schematic diagram of the principle of excavation work volume for each soil mass obtained based on the tensor method;
FIG. 4 is a flow chart of an earthwork calculation method for differentiating soil properties according to another embodiment;
FIG. 5 is a schematic diagram of the original distribution of the measured point cloud;
FIG. 6 is a schematic projection diagram of the original distribution of the measured point cloud on the xoy plane;
FIG. 7 is a schematic projection diagram of the reconstructed measurement point cloud on the xoy plane;
FIG. 8 is a schematic diagram of three-dimensional distribution of a reconstructed point cloud;
FIG. 9 is a schematic logic diagram of a sparse matrix compression algorithm;
FIG. 10 is a flowchart of an earthwork calculation method for differentiating soil properties according to yet another embodiment;
fig. 11 is a block diagram of an earth mass calculating device capable of distinguishing the quality of earth according to an embodiment.
Detailed Description
For better understanding of the objects, technical solutions and effects of the present invention, the present invention will be further explained with reference to the accompanying drawings and examples. Meanwhile, the following described examples are only for explaining the present invention, and are not intended to limit the present invention.
The embodiment of the invention provides an earthwork calculation method capable of distinguishing soil texture.
Fig. 1 is a flowchart of an earth calculation method for distinguishing an earth mass according to an embodiment, and as shown in fig. 1, the earth calculation method for distinguishing an earth mass according to an embodiment includes steps S100 to S105:
s100, acquiring a three-dimensional geological digital model, and rasterizing the three-dimensional geological digital model to obtain a discrete grid array;
and (3) creating a geological attribute scalar field by introducing a three-dimensional geological digital model, and rasterizing the scalar field with certain precision. The mathematical book-continuous geological model is replaced by a discrete grid array.
In one embodiment, the grid dimensions of the grid array are determined based on the accuracy of the three-dimensional geosgital model itself. As a preferred embodiment, the three-dimensional geological digital model comprises a BIM (building Information model) three-dimensional geological model or a GIS-based three-dimensional geological model. Based on this, the grid dimension of the grid array is 0.5m to 2 m. As a preferred embodiment, the grid array has a grid dimension of 1 m.
In order to better explain the technical solution of the embodiment of the present invention, the following is explained by using an application example. It should be noted that the application examples are only for convenience of explanation and do not represent the only limitations of the embodiments of the present invention. In the application example, the grid dimension l is 1 m. And setting an independent numerical value code j for each soil texture, and establishing function mapping of the address attribute at the central point of each grid in the Euclidean space.
S101, determining function mapping of each field according to the grid array; wherein the form of the function mapping is a three-dimensional tensor;
in one embodiment, the research domain is divided into domains, and the function mapping of the grid array in each domain is output respectively. The form of the function mapping is a three-dimensional tensor, so that the memory can be written conveniently.
In an application example, fig. 2 is a schematic diagram of tensor mapping of a rasterized geological model, as shown in fig. 2, a three-dimensional tensor TiAnd writing into a memory, wherein i represents a field number.
S102, extracting sub-tensors corresponding to various soil textures from the three-dimensional tensors in various fields;
in an application example, fig. 3 is a schematic diagram illustrating the principle of obtaining excavation work volume of each soil texture based on a tensor method, and as shown in fig. 3, a sub-tensor T corresponding to each type of soil texture is extracted from a three-dimensional tensor of each fieldijSetting the sub-tensor TijIs the initial residual tensor.
S103, respectively obtaining a current topographic survey point cloud and a current topographic survey point cloud according to a time sequence, fitting and calculating a current accumulated excavation tensor according to the current topographic survey point cloud, and obtaining a current excavation tensor according to the current accumulated excavation tensor and the current accumulated excavation tensor;
the method comprises the steps of forming topographic measurement point cloud data of continuous periods in a time sequence, and when determining a topographic measurement point cloud of a current period, taking the topographic measurement point cloud of the previous period of the current period as the topographic measurement point cloud of the previous period. In an application example, after the current topographic measurement point cloud and the current topographic measurement point cloud are obtained and guided in, fitting calculation is carried out, and the current accumulated excavation tensor E is calculated according to the current topographic measurement point cloud fittingkCalculating the cumulative excavation tensor E of the previous period according to the point cloud fitting of the topographic survey of the previous periodk-1. Where k denotes the number of periods, and when k is 0, it denotes the initial topography.
In one embodiment, the non-0 element in the neutron tensor in step S103 is returned to 1 to indicate that the state of the soil property at the corresponding spatial position is present.
In one example, fig. 4 is a flowchart of an earthwork calculation method for differentiating soil texture according to another embodiment, and as shown in fig. 4, the process of calculating the current accumulated excavation tensor according to the current topographic survey point cloud fitting in step S103 includes steps S200 to S202:
s200, acquiring point cloud original data of a last-stage topographic measurement point cloud;
s201, based on the point cloud original data, generating a curved function by applying a three-dimensional curved surface interpolation method, and reconstructing by taking a specific grid scale to generate a new point cloud array; wherein the grid dimension is a grid dimension of the grid array;
and S202, outputting the up-term accumulated excavation tensor based on the new point cloud array.
In an application example, fig. 5 is a schematic diagram of original distribution of measurement point cloud, fig. 6 is a schematic diagram of projection of the original distribution of the measurement point cloud on a xoy plane, as shown in fig. 5 and fig. 6, a three-dimensional curved surface interpolation method is applied to generate a curved surface function, and point cloud original data are reconstructed by taking the same grid scale l. Fig. 7 is a schematic projection diagram of the reconstructed measurement point cloud on the xoy plane, and the reconstructed point cloud raw data is as shown in fig. 7. And generating a new point cloud array based on the reconstruction shown in fig. 7, wherein fig. 8 is a schematic diagram of three-dimensional distribution of the reconstructed point cloud, and the new point cloud array is shown in fig. 8. The accumulated excavation tensor of the geological tensor is output by the new point cloud array, and the values of excavation tensor elements above each point (in the z + direction) in the point cloud are set to be 1, and the values of elements below are set to be 0.
S104, performing point multiplication operation on the sub tensors of all the soil textures in each field and the current excavation tensor of the corresponding field to obtain the current excavation tensor of each soil texture, and updating the sub tensors according to the current excavation tensor of each soil texture;
in the application example, the number i indicates the field, the number j indicates the soil texture, and the number ij indicates the sub-tensor T of each soil texture in each field of the soil texture j of the field iijCurrent excavation tensor F corresponding to corresponding domainijkPerforming dot product operation, i.e. O ═ T × F, and outputting current excavation tensor O of each soil texture based on the operationijk. According to the nature of the soilCurrent excavation tensor OijkUpdating the sub-tensor TijI.e. the updated sub-tensor T' ═ T-O.
In one embodiment, the above steps S100 to S104 relate to tensor operation, and include the steps of: tensor operations are performed by a sparse matrix compression algorithm.
The storage and calculation efficiency of the large tensor is optimized through the sparse matrix compression algorithm, so that the storage and calculation of the large tensor are facilitated. Fig. 9 is a schematic logic diagram of a sparse matrix compression algorithm, and as shown in fig. 9, in the operation of the large sparse tensor, the sparse row compression format (CSR) or the sparse column compression format (CSC) is used for storage and operation.
In one embodiment, the above steps S100 to S104 relate to tensor operation, and include the steps of: tensor operations are performed by parallel operations.
The operation of the large sparse tensor can be set into multi-path parallel operation, and the parallel number of the parallel operation can be set according to the size of the three-dimensional geological digital model, so that the tensor operation efficiency is improved.
And S105, acquiring the sum of the elements of the current-period excavation tensor of each type of soil texture in each field, and converting the sum of the elements into the space occupancy rate of an Euclidean space according to the grid scale to be used as the current-period excavation amount of each type of soil texture.
The element sum of the current excavation tensor of various soil qualities is fed back through the execution of the algorithm, and the space occupancy rate of the Euclidean space is converted according to the element sum.
In one embodiment, fig. 10 is a flowchart of an earthwork calculation method for differentiating soil texture according to a further embodiment, and as shown in fig. 10, the step S105 of converting the element sum into the space occupancy rate of euclidean space according to the grid scale as the current earthwork excavation amount of each type of soil texture includes the step S300:
and S300, obtaining the current earth excavation volume according to the product of the element and the cube of the grid scale of the grid array.
In the application example, the current excavation tensor O of each soil texture in each field is calculated respectivelyijkAnd a andijkthe voxel sums are converted into space occupation of the corresponding euclidean space according to the grid scale. The volume V of the earth mass j in the k-th period in the area i in the step S300ijk＝l3aijk。
The earthwork calculation method capable of distinguishing the soil texture according to any one of the embodiments described above constructs an intermediate data standard by tensor based on the concept of a field, and solves the problems of tolerance and grid quality degradation in the conventional three-dimensional figure method by rasterizing the three-dimensional geological digital model. Meanwhile, a mapping rule that the rasterized model is converted into a tensor model is set to calculate the earthwork excavation amount at the current period, so that quantitative analysis of engineering geology is facilitated, related personnel can establish application programs for executing discrete tasks in various fields according to the mapping rule, and earthwork calculation amount of various types of soil can be accurately and efficiently completed in various engineering and terrains.
The embodiment of the invention also provides an earthwork calculation device capable of distinguishing the soil texture.
Fig. 11 is a block diagram illustrating an earth mass calculating device capable of distinguishing the soil mass according to an embodiment, and as shown in fig. 11, the earth mass calculating device capable of distinguishing the soil mass according to an embodiment includes a block 100, a block 101, a block 102, a block 103, a block 104, and a block 105:
the rasterization processing module 100 is configured to obtain a three-dimensional geological digital model, and perform rasterization processing on the three-dimensional geological digital model to obtain a discrete grid array;
a three-dimensional tensor determination module 101, configured to determine function mappings of the respective fields according to the grid array; wherein the form of the function mapping is a three-dimensional tensor;
the sub-tensor extraction module 102 is configured to extract sub-tensors corresponding to various soil textures from three-dimensional tensors in various fields;
the current excavation tensor determining module 103 is used for respectively obtaining a current topographic survey point cloud and a current topographic survey point cloud according to a time sequence, calculating a current accumulated excavation tensor according to the current topographic survey point cloud fitting, and obtaining a current excavation tensor according to the current accumulated excavation tensor and the current accumulated excavation tensor;
the soil texture current excavation tensor determining module 104 is used for performing point multiplication operation on the sub tensors of various soil textures in each field and the current excavation tensor of the corresponding field to obtain the current excavation tensor of each soil texture, and updating the sub tensors according to the current excavation tensor of each soil texture;
and the current-period earthwork excavation amount determining module 105 is used for acquiring the sum of elements of current-period excavation tensor of various types of soil in each field, and converting the sum of elements into the space occupancy rate of an Euclidean space according to the grid scale to be used as the current-period earthwork excavation amount of various types of soil.
The earthwork calculation amount device capable of distinguishing the soil texture constructs the intermediate data standard through tensor based on the field concept, and solves the problems of tolerance and grid quality reduction existing in the traditional three-dimensional graphic method through rasterization processing of the three-dimensional geological digital model. Meanwhile, a mapping rule that the rasterized model is converted into a tensor model is set to calculate the earthwork excavation amount at the current period, so that quantitative analysis of engineering geology is facilitated, related personnel can establish application programs for executing discrete tasks in various fields according to the mapping rule, and the earthwork calculation amount can be accurately and efficiently completed in various projects and terrains.
Embodiments of the present invention further provide a computer storage medium having stored thereon computer instructions, which when executed by a processor, implement the method for earth computation capable of distinguishing soil quality according to any of the above embodiments.
Those skilled in the art will understand that: all or part of the steps for implementing the method embodiments may be implemented by hardware related to program instructions, and the program may be stored in a computer readable storage medium, and when executed, the program performs the steps including the method embodiments; and the aforementioned storage medium includes: various media capable of storing program codes, such as a removable Memory device, a Random Access Memory (RAM), a Read-Only Memory (ROM), a magnetic disk, and an optical disk.
Alternatively, the integrated unit of the present invention may be stored in a computer-readable storage medium if it is implemented in the form of a software functional module and sold or used as a separate product. Based on such understanding, the technical solutions of the embodiments of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a terminal, or a network device) to execute all or part of the methods of the embodiments of the present invention. And the aforementioned storage medium includes: a removable storage device, a RAM, a ROM, a magnetic or optical disk, or various other media that can store program code.
Corresponding to the computer storage medium, in one embodiment, there is also provided a computer device including a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor executes the computer program to implement any one of the above-described soil mass-distinguishable earthwork calculation methods.
The computer device constructs an intermediate data standard through tensor based on the field concept, and solves the problems of tolerance and grid quality reduction existing in the traditional three-dimensional graphic method through rasterization processing of the three-dimensional geological digital model. Meanwhile, a mapping rule that the rasterized model is converted into a tensor model is set to calculate the earthwork excavation amount at the current period, so that quantitative analysis of engineering geology is facilitated, related personnel can establish application programs for executing discrete tasks in various fields according to the mapping rule, and the earthwork calculation amount can be accurately and efficiently completed in various projects and terrains.
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above examples only show some embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent shall be subject to the appended claims.

Claims (10)

1. An earthwork calculation method capable of distinguishing soil texture, comprising the steps of:
acquiring a three-dimensional geological digital model, and rasterizing the three-dimensional geological digital model to obtain a discrete grid array;
determining function mapping of each field according to the grid array; wherein the functional mapping is in the form of a three-dimensional tensor;
extracting sub-tensors corresponding to various soil textures from the three-dimensional tensors in various fields;
respectively acquiring a current-period topographic measurement point cloud and a current-period topographic measurement point cloud according to a time sequence, fitting and calculating a current-period accumulated excavation tensor according to the current-period topographic measurement point cloud, and acquiring a current-period excavation tensor according to the current-period accumulated excavation tensor and the current-period accumulated excavation tensor;
performing point multiplication operation on the sub-tensors of various soil textures in each field and the current excavation tensor of the corresponding field to obtain the current excavation tensor of each soil texture, and updating the sub-tensors according to the current excavation tensor of each soil texture;
the method comprises the steps of obtaining element sums of current-period excavation tensor of various types of soil in various fields, and converting the element sums into space occupancy rates of Euclidean spaces according to grid scales to serve as current-period excavation amounts of the various types of soil.
2. The method for calculating earth mass according to claim 1, wherein the calculating of the cumulative excavation tensor of the upper period according to the point cloud fitting of the upper period topographic survey comprises the steps of:
acquiring point cloud original data of the current topographic measurement point cloud;
based on the point cloud original data, generating a curved surface function by applying a three-dimensional curved surface interpolation method, and reconstructing a specific grid scale to generate a new point cloud array; wherein the grid dimension is a grid dimension of the grid array;
and outputting the up accumulated excavation tensor based on the new point cloud array.
3. The method of claim 1, wherein each tensor operation procedure comprises the steps of:
performing the tensor operation by a sparse matrix compression algorithm.
4. The method of claim 1, wherein each tensor operation procedure comprises the steps of:
the tensor operation is performed by a parallel operation.
5. The method for earthwork computation of distinguishable earthiness according to claim 1, wherein the process of converting the element sum into space occupancy of euclidean space according to grid dimension as the current earthwork volume of each kind of earthiness comprises the steps of:
obtaining the current amount of earth excavation from a product of the element and a cube of a grid scale of the grid array.
6. The method of earthwork computation of distinguishable earthiness according to any one of claims 1 to 5, wherein the grid array has a grid dimension of 1 m.
7. The method of earthwork differential earthiness according to any one of claims 1 to 5, wherein the three-dimensional geological digital model comprises a BIM three-dimensional geological model.
8. An earth mass calculating device capable of distinguishing the quality of earth, comprising:
the rasterization processing module is used for acquiring a three-dimensional geological digital model and rasterizing the three-dimensional geological digital model to obtain a discrete grid array;
the three-dimensional tensor determining module is used for determining function mapping of each field according to the grid array; wherein the functional mapping is in the form of a three-dimensional tensor;
the sub-tensor extraction module is used for extracting sub-tensors corresponding to various soil textures from the three-dimensional tensors in various fields;
the current excavation tensor determining module is used for respectively obtaining a current topographic survey point cloud and a current topographic survey point cloud according to a time sequence, calculating a current accumulated excavation tensor according to the current topographic survey point cloud fitting, and obtaining the current excavation tensor according to the current accumulated excavation tensor and the current accumulated excavation tensor;
the soil texture current excavation tensor determining module is used for performing point multiplication operation on the sub-tensors of various soil textures in each field and the current excavation tensor of the corresponding field to obtain the current excavation tensor of each soil texture, and updating the sub-tensors according to the current excavation tensor of each soil texture;
and the current-period earthwork excavation amount determining module is used for acquiring the sum of elements of current-period excavation tensor of various types of soil in each field, and converting the sum of elements into the space occupancy rate of an Euclidean space according to the grid scale to be used as the current-period earthwork excavation amount of various types of soil.
9. A computer storage medium having computer instructions stored thereon, wherein the computer instructions, when executed by a processor, implement the method of earth mass differentiable earth mass computation of any one of claims 1 to 7.
10. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the method of earth mass discrimination calculation according to any one of claims 1 to 7 when executing the program.
CN202110100004.7A 2021-01-25 2021-01-25 Earthwork calculation method and device capable of distinguishing soil texture Active CN112861222B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110100004.7A CN112861222B (en) 2021-01-25 2021-01-25 Earthwork calculation method and device capable of distinguishing soil texture

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110100004.7A CN112861222B (en) 2021-01-25 2021-01-25 Earthwork calculation method and device capable of distinguishing soil texture

Publications (2)

Publication Number Publication Date
CN112861222A true CN112861222A (en) 2021-05-28
CN112861222B CN112861222B (en) 2023-07-14

Family

ID=76008868

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110100004.7A Active CN112861222B (en) 2021-01-25 2021-01-25 Earthwork calculation method and device capable of distinguishing soil texture

Country Status (1)

Country Link
CN (1) CN112861222B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113919025A (en) * 2021-09-06 2022-01-11 广东电网有限责任公司广州供电局 Method and device for determining pit engineering quantity, computer equipment and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002081008A (en) * 2000-09-11 2002-03-22 Kokusai Kogyo Co Ltd Earthwork design support system using digital map data
CN106777823A (en) * 2017-01-25 2017-05-31 天津大学 Construction schemes of underground structure group progress simulation optimization method based on ventilation numerical simulation
CN110532719A (en) * 2019-09-04 2019-12-03 上海宝冶集团有限公司 A kind of construction site Method of Earthwok Calculation and device
CN111595403A (en) * 2020-05-15 2020-08-28 中交广州航道局有限公司 Engineering earthwork measuring method based on point cloud measuring technology
CN112100715A (en) * 2020-08-20 2020-12-18 中国建筑第八工程局有限公司 Three-dimensional oblique photography technology-based earthwork optimization method and system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002081008A (en) * 2000-09-11 2002-03-22 Kokusai Kogyo Co Ltd Earthwork design support system using digital map data
CN106777823A (en) * 2017-01-25 2017-05-31 天津大学 Construction schemes of underground structure group progress simulation optimization method based on ventilation numerical simulation
CN110532719A (en) * 2019-09-04 2019-12-03 上海宝冶集团有限公司 A kind of construction site Method of Earthwok Calculation and device
CN111595403A (en) * 2020-05-15 2020-08-28 中交广州航道局有限公司 Engineering earthwork measuring method based on point cloud measuring technology
CN112100715A (en) * 2020-08-20 2020-12-18 中国建筑第八工程局有限公司 Three-dimensional oblique photography technology-based earthwork optimization method and system

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
廖伟华 等: "基于空间内插的场地平整土方工程量计算方法研究", 《黑龙江工程学院学报（自然科学版）》 *
廖伟华 等: "基于空间内插的场地平整土方工程量计算方法研究", 《黑龙江工程学院学报（自然科学版）》, vol. 21, no. 4, 31 December 2007 (2007-12-31), pages 14 - 17 *
殷达 等: "基于微结构张量的岩石各向异性弹塑性本构及其应用", 《岩土工程学报》 *
殷达 等: "基于微结构张量的岩石各向异性弹塑性本构及其应用", 《岩土工程学报》, vol. 42, no. 9, 30 September 2020 (2020-09-30), pages 1751 - 1754 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113919025A (en) * 2021-09-06 2022-01-11 广东电网有限责任公司广州供电局 Method and device for determining pit engineering quantity, computer equipment and storage medium

Also Published As

Publication number Publication date
CN112861222B (en) 2023-07-14

Similar Documents

Publication Publication Date Title
AU2005277378B2 (en) Multiple-point statistics simulation with enhanced computational efficiency
Arpat et al. Conditional simulation with patterns
Mostafavi et al. Delete and insert operations in Voronoi/Delaunay methods and applications
Shahzad et al. TecDEM: A MATLAB based toolbox for tectonic geomorphology, Part 1: Drainage network preprocessing and stream profile analysis
Zovi et al. Identification of high-permeability subsurface structures with multiple point geostatistics and normal score ensemble Kalman filter
CN111339691A (en) Intelligent geotechnical engineering parameter three-dimensional analysis and evaluation system and method based on voxler software
Conolly Geographical information systems and landscape archaeology
Wu et al. An effective parallelization algorithm for DEM generalization based on CUDA
CN111383336A (en) Three-dimensional geological model construction method
CN107221028A (en) A kind of geologic body closed surface three-dimensional rebuilding method based on seismic interpretation data
CN112861222B (en) Earthwork calculation method and device capable of distinguishing soil texture
Marschallinger et al. A visual LISP program for voxelizing AutoCAD solid models
Ferreira et al. Computed-tomography-based discrete fracture-matrix modeling: An integrated framework for deriving fracture networks
CN115758792A (en) Geological disaster assessment method and device based on digital numerical integration
Rengers et al. Information Technology applied to engineering geology
Wang et al. Towards automatic and rapid 3D geological modelling of urban sedimentary strata from a large amount of borehole data using a parallel solution of implicit equations
Allgeier et al. Surrogate‐Model Assisted Plausibility‐Check, Calibration, and Posterior‐Distribution Evaluation of Subsurface‐Flow Models
Tonini et al. Three-dimensional reconstruction of the Carrara Syncline (Apuane Alps, Italy): An approach to reconstruct and control a geological model using only field survey data
Fan et al. Automatic reconstruction of geological reservoir models based on conditioning data constraints and BicycleGAN
Hebeler et al. 14 Modeling DEM Data Uncertainties for Monte Carlo Simulations of Ice Sheet Models
Zhang et al. Conditional stochastic simulation of fluvial reservoirs using multi-scale concurrent generative adversarial networks
Liu et al. Study of a high-precision complex 3D geological modelling method based on a fine KNN and kriging coupling algorithm: a case study for Jiangsu, China
Hamidi et al. Developing an ArcGIS-based toolbox for automatic analysis of instrumentation data for engineering structures: case of embankment dam
Loudon Digital spatial models and geological maps
Zhuang et al. Estimation of the volume of a rock layer using kriging method for energy efficient excavation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant