高精度土方工程量计算及其软件开发

高精度土方工程量计算及其软件开发刘仁钊1，吕永洪2，张玉堂3（1，3.湖北国土资源职业学院，湖北武汉，430090；2.杭州阵列科技有限公司，浙江杭州，310000）摘要：土方工程是各项工程建设前期工作的主要项目之一，其主要目的是通过填挖土方使平整地块的平面位置和高程满足设计要求。

由于受原始地形复杂多样、外业数据采集方法、以及土方计算模型等因素影响，要得到精确的土方计算结果很困难。

本文首先评述了各类土方计算原理和方法，分析了土方计算中存在的问题和当前土方计算软件存在的不足，然后提出了高精度土方计算的各类措施和方法，并介绍了作者以此为基础开发的高精度专业土方计算软件HTCAD。

关键词：土方计算，方格网法，TIN法，土方精度，软件HTCAD0.引言各项工程建设中，土方工程是工程建设前期工作的一项重要工作内容，其主要目的是通过填挖场区土方使平整地块的平面位置和高程满足设计要求。

土方工程量的大小，不仅直接影响投资方的资金投入，同时也影响到工程建设进度和工期。

由于受原始地形复杂多样、外业数据采集方法、以及土方计算模型等因素影响，导致土方工程量的不精确性，容量引起业主和施工方之间的经济纠纷，出现通过虚报土方工程量套取国家建设项目资金等经济问题。

在国土行业土地整治中，土方工程量在整个工程量中占有比例高达40%，土方工程量占用资金往往是重点审计内容。

因此，土方工程量计算是一项涉及到项目各方利益的工作，高精度土方量的确定是项目各方的共同要求，因而研究提高土方量的计算方法是非常必要的。

本文对工程土方量计算中常用的断面法、方格网法、等高线法及基于数字高程模型的不规则三角网（TIN）法的基本原理、方法和优缺点进行比较分析，重点探讨了当前土方计算软件存在的不足，提出了高精度土方计算的各类措施和方法，并以此介绍了作者开发的专业土方计算软件HTCAD。

1.计算方法分析土方计算的目的主要是计算同一场地土方工程施工前后的填、挖方量，实际上就是计算移动土体的体积。

根据不同情况，土方计算可选用不同的方法。

常用的土方计算方法有方格网法、断面法、等高线法及基于DEM的计算法等方法。

在实际应用中，由于地形表面的复杂性和工程项目的特殊性，土方计算的准确性和精确性是人们最关心的问题。

1.1 断面法断面计算法，又称为截面法。

断面法的工作原理是在线状工程建设中，以线状工程中桩线为基准，按一定的间距沿垂直中线方向划分为若干个相互平行的横截面，量出各横断面之间的距离，按照设计高程与地面线所组成的断面图，计算每条断面线所围成的填、挖面积，再由两端横断面的平均填、挖面积乘以两端横断面之间的距离，求出每相邻两断面间的填、挖体积；最后将各相邻断面的填、挖体积加起来，求出填、挖总体积，这种方法称为断面法。

如下图1所示。

断面法是与其它方法完全不同的一种土方计算方法，它是通过取得地形的剖面来计算剖面之间的体积。

实践证明，对于线状建设工程，例如道路、沟渠、堤坝等这类设计地形有突变的情形，采用断面法是十分适合的，比其它方法更具优势。

计算相对简单方便，而且具有一定的精度保证。

目前有不少测量专业软件和土方计算软件都有断面法计算土方功能。

实践表明，要提高断面法的土方计算精度，需要着重考虑以下方面因素：（1）断面间距。

断面间距对于不同的工程，规范有相应的要求。

例如水利工程断面间距多为50—100米，公道工程断面间距一般为20-50米等。

容易理解，间距越短，剖面越能近似反映真实地表面，精度愈高，但外业工作量显著增大。

对于地形变化较大和中线方向发生变化的地段要缩短断面间距。

（2）实测断面。

对断面线上的每个变坡点用测量仪器实测，代替用地形散点内插，其精度要高很多。

（3）采用精度更高的计算公式计算，土方精度更高一些。

1.2等高线法等高线法就是将相邻两等高线截面间的形体看作是几何台体，如图2。

根据几何台体公式计算其体积，最后把所有相邻台体体积相加，得到总体积。

利用现成的绘有等高线的地形图，计算相邻两等高线所围的面积，然后乘以等高距求得台体体积。

当两相邻等高线形状相近，且面积大小相差不大时，可采取平均值的计算方法求取台体体积。

等高线法计算土方是教科书中常提及的一种方法，容易理解和实施，主要用于山地或高山地形。

在实际中多用在水库库容量计算中。

等高线是经过地形散点拟合后得到的，由于散点高程精度不高，加之测量位置可能不当，造成等高线局部可能不合理，因此计算精度不高，在目前土方计算中较少采用。

有时由于工程投资预算上的需要，同时为了充分利用已有的地形图，等高线往往和其他方法，如方格网法、TIN 法等结合进行。

因为等高线具有坐标和高程信息，具体做法是将等高线按一定的间距离散化后得到地形高程点，没有等高线的地方用测量仪器测量一些地形散点，然后用格网法、TIN 法等方法计算土方。

1.3 方格网法方格网法是根据场地情况和土方计算精度要求，将其划分为边长5—50m 的正方形方格网，然后通过测量手段测量或通过周围高程点内插得出每个方格四个角点的标高，再根据四个角点的设计标高计算出设计标高与自然地面标高的差值，即为各角点的施工高度（挖或填），将施工高度标注于角点上，然后分别计算每一方格地填挖土方量，并算出场地边坡的土方量。

将挖方区（或填方区）所有方格计算的土方量和边坡土方量汇总，即得场地挖方量和填方量的总土方量。

图1断面线与断面图2 等高线及山体水平投影图方格网法是专业教材中介绍最多的方法，目前都是通过地形散点高程采用三点线性内插得出方格角点的地面标高和设计标高。

由于方格网法土方计算特点，在实际应用中多限于地势变化平缓的低丘陵、平原地和起伏不大的场地，对于地形起伏变化较大的丘陵和山地地形误差较大。

此外，方格网边长是影响土方计算的一个重技术参数，下表1是在一个标准的地形情况下（没有坡、坎等突变地形），用南方CASS 软件计算的结果。

表1 方格网法边长和土方计算表实际结果：填：28294，挖：68343总和：96637由上表可见，即使是同一场地，都是方格网法，其结果差别相差很大。

按图下方的实际土方进行比较，边长取20米时，计算土方误差达到6.1%。

，取边长4米最为合适，误差仅0.05%。

根据经验，下面给出了提高方格网法土方计算精度的一些建议。

（1）地形高程点的密度不能太小。

目前规范还没有详细规定和要求，在大量的土方计算验证基础上，建议按不低于1：500比例尺地形测图要求，即地形高程点间距在5-15米较为合适，这样既能减少外业工作量，又能保证土方精度要求。

（2）地形高程点的高程精度尽可能高一些。

建议地形高程点测量精度要求：其高程中误差小于10cm ，目前使用的全站仪和GPS RTK 都能达到，经纬仪和平板仪测量的地形点高程只能达高分米级精度。

（3）方格网边长选取要适当。

边长不同，其计算结果可能相差很大，这可以通过试算进行确定，当计算结果变化趋于稳定就可以了。

也可以采用经验统计值，若边长选择为场区内散点高程间距的1/2左右时，其结果比较理想，可以保证计算误差小于0.1%。

（4）尽量避免陡坎横穿方格中间，如图3，造成内插高程实真。

在这种情况下可调整或加密方格网，并在外业时对坎上坎下同时采集高程点，以保证内插高程的正确性。

显然，对于梯田坎较多的土地整理项目来说，涉及多个田块合并时，采用方格网法计算土方量其精度不高。

（5）采用更好的拟合模型，提高方格角点内插高程精度，例如多点加权线性模型、二次曲面模型和加权最小二乘拟合曲面模型等。

目前市面上的土方软件都是采用三点线性内插模型。

线性内插是使用最靠近插值点的三个已知数据点，确定一个平面，继而求出内插点的高程，这种模型严重依赖三个已知点，因此其内插精度不是很高。

1.4 TIN 法数字地面模型（Digital Terrain Model ，DTM ）是用一群地面点的平面坐标和高程描述地表形状的一种方式。

当模型中以高程为特征值的DTM 也称为数字高程模型（Didital Elevation Model ，DEM ）。

DEM 是用数字形式x ，y ，z 坐标来表达区域内的地貌形态，以缩微图3 陡坎穿过方格中间的形式再现了地表形态起伏变化特征，如图4,具有形象、直观、精确等特点，是建立GIS 最重要的内容，广泛应用于各种工程规划和地形分析，并被用于土方工程量的计算。

按照数据结构划分数字地面模型可分为栅格形式的规则格网模型（GRID ）和矢量形式的不规则三角网（Triangulated Irregular Network ，TIN ）模型。

TIN 模型主要是通过有限的地面测量点构建三角网的方式模拟地表现状。

TIN 模型有多种建立方法，目前GIS 中广泛应用带约束的狄洛尼（Delaunay ）方法，基于散点的狄洛尼构网算法理论严密、唯一性好。

然后以每个三角形为基础，根据设计面计算每个不规则三角形所构成的立柱体的体积，再把区域内所有三角形立柱体的体积相加得到场区总体积。

无论是理论还是实践都表明，除了线状工程采用断面法外，基于数字地面模型的TIN 法目前是土方计算精度相对最高的一种方法，适合于各种不同地形类型。

TIN 法直接采用外业地形高程点构建三角网，只要地形高程点布置合理，其土方计算精度是较高的，特别适用于设计面不规则的场地，如大型土方工程的每月土方统计，露天采矿区每月采矿量统计等，其开采面和开挖面是不规则的。

不难看出，要提高TIN 法的土方计算精度，要注意以下几方面：（1）按照狄洛尼自动生成的三角网，一定要与实际地形相吻合。

在实际地形中，存在诸多地形变化特征线，如：山脊线、山凹线、坡顶坡脚线、陡坎线等，当出现与实际地形不相吻合，需要调整三角网，使任何一个三角形都不穿过地形变化特征线，才能建立起与现状地面相吻合的数字地面模型，如图5。

（2）当设计地面是不可表达的数学函数时，也涉及到内插设计高程。

为了提高地形点内插设计高程精度，采用多点加权线性模型、二次曲面模型和加权最小二乘拟合曲面模型等。

（3）按照如前所述的建议，按不低于1：500比例尺地形测图要求，保持地形高程点间距在5-15米较为合适，并特别注意增加变坡点、坎上坎下高程点采集（如图5），确保数字地面模型与实际地形的一致性。

1.5 各类土方计算方法精度统计作者综合已有精度分析资料，并根据市场上几种主要的商用土方计算软件，结合多种地形进行模拟计算，统计出各类方法计算土方量的适用范围和可能达到的精度，供工程技术人员计算参考。

表2 各类土方计算方法精度统计（单位：百分比）从上表可以看出，方格网法和TIN 法适用地形最广。

当场地为地形变化平缓低丘陵、平原时，方格网法和TIN 法精度基本相当，在其它地形场地TIN 法明显要比方格网法精度高。

图5 调整前、后的三角形图图4 由TIN 生成的DEM 模型2.高精度土方计算软件HTCAD如何进一步提高软件对各类不同地形的适应性，在计算时充分考虑各种坡、坎、沟、坝、塘等特殊线状和面状地物，并提高模型计算精度，是目前土方计算软件的两大发展方向。