本发明公开了一种基于InSAR技术高填方压实度检测与沉降补偿方法，包括以下步骤：步骤一、获取所施工高填方的SAR图像；步骤二、获取所施工高填方的DEM数据；三、所施工高填方压实度的检测；四、所成型高填方的初始体积和初始压实度计算；五、获取沉降过程中所成型高填方的SAR图像；六、获取沉降过程中所成型高填方的DEM数据；七、沉降过程中所成型高填方的损失体积和压实度计算；八、沉降后所成型高填方的补偿。本发明方法步骤简单，对所施工高填方压实度以及对沉降过程中所成型高填方的损失体积和压实度准确检测，并获取所成型高填方的损失体积和所施工高填方的压实度之间关系，为沉降后所成型高填方的补偿提供依据。
1.一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于，该方法包括以下步骤：
步骤一、获取所施工高填方的SAR图像：在所施工高填方填筑至成型过程中，采用机载干涉合成孔径雷达InSAR对所施工高填方进行测量，获取每一个测量时间内的两幅所施工高填方单视SAR图像，并将所获取的每一个测量时间内的两幅所施工高填方单视SAR图像存储至与数据处理设备相接的存储器内；其中，两幅所述所施工高填方单视SAR图像均为正射影像，且两幅所述所施工高填方单视SAR图像分别为所施工高填方SAR主图像和所施工高填方SAR辅图像；
步骤二、获取所施工高填方的DEM数据：采用所述数据处理设备按照时间先后顺序对步骤一中各个测量时间内存储至所述存储器中的两幅所述所施工高填方单视SAR图像分别进行雷达数据处理，获取各个测量时间内的所施工高填方DEM数据，并存储至所述存储器内；
其中，所施工高填方DEM数据包括原始坡地或沟谷高程Z0(x,y)、当前填筑层的顶面高程Zi(x,y)、上一填筑层的顶面高程Zi-1(x,y)和所成型高填方的顶面高程Zd(x,y)；
步骤三、所施工高填方压实度的检测：所施工高填方由下至上分多层填筑层，每个所述填筑层的压实度检测方法均相同，其中对任一所述填筑层的压实度进行检测，具体过程如下：步骤301、当前填筑层的体积计算：采用所述数据处理设备且根据公式Vi＝∫∫∑(Zi(x,y)-Zi-1(x,y))ds，并结合步骤二中所存储的当前填筑层的顶面高程Zi(x,y)和步骤二中所存储的上一填筑层的顶面高程Zi-1(x,y)，得到当前填筑层的体积Vi，并存储至所述存储器内；
其中，ds为积分单元面积，i为正整数且i＝1、2、…、M，M为所述填筑层的数量；
步骤302、当前填筑层的干密度计算：采用所述数据处理设备且根据公式 并结合步骤301中所存储的当前填筑层的体积Vi，得到当前填筑层的干密度ρi，并存储至所述存储器内；
步骤303、当前填筑层的压实度计算：采用所述数据处理设备且根据公式
并结合步骤302中所存储的当前填筑层的干密度ρi，得到当前填筑层的压
实度λi，并存储至所述存储器内，所述数据处理设备还传输至显示器进行输出显示；其中，mi为当前填筑层的质量，ρimax为当前填筑层所用土的最大干密度，ωi为当前填筑层的含水量；
步骤304、多次重复步骤301至步骤303，完成多个所述填筑层的压实度检测，直至所施工高填方的压实度满足设定的压实度，所施工高填方施工成型；
步骤四、所成型高填方的初始体积和初始压实度计算：所施工高填方施工成型后，对所成型高填方的初始体积和初始压实度进行计算的具体过程为：首先，采用所述数据处理设备且根据公式V1＝∫∫∑(Zd(x,y)-Z0(x,y))ds，并结合步骤二中所存储的所成型高填方的顶面高程Zd(x,y)和步骤二中所存储的原始坡地或沟谷高程Z0(x,y)，得到所成型高填方的初始体积V1，并存储至所述存储器内；然后，采用所述数据处理设备且根据公式 得到所成型高填方的干密度ρ1，并存储至所述存储器内；再采用所述数据处理设备且根据公式 得到所成型高填方的初始压实度λ1，并存储至所述存储器内；其中，m1为所成型高填方的质量，ρmax为所成型高填方所用土的最大干密度，ω1为所成型高填方的含水量；
步骤五、获取沉降过程中所成型高填方的SAR图像：所施工高填方施工成型后，所成型高填方发生沉降，在所成型高填方沉降的过程中，采用所述机载干涉合成孔径雷达InSAR对所成型高填方进行测量，获取每一个测量时间内的两幅所成型高填方单视SAR图像，并将所获取的每一个测量时间内的两幅所成型高填方单视SAR图像存储至所述存储器内；其中，两幅所述所成型高填方单视SAR图像均为正射影像，且两幅所述所成型高填方单视SAR图像分别为所成型高填方SAR主图像和所成型高填方SAR辅图像；
步骤六、获取沉降过程中所成型高填方的DEM数据：采用所述数据处理设备按照时间先后顺序对步骤五中各个测量时间内存储至所述存储器中的两幅所述所成型高填方单视SAR图像分别进行雷达数据处理，获取各个测量时间内的沉降过程中所成型高填方DEM数据，并存储至所述存储器内；其中，所述沉降过程中所成型高填方DEM数据包括沉降过程中所成型高填方顶面高程Zj(x,y)；
步骤七、沉降过程中所成型高填方的损失体积和压实度计算：所施工高填方施工成型后，在所成型高填方沉降的过程中，对所成型高填方的损失体积和压实度进行计算，并获取所成型高填方的损失体积和所施工高填方的压实度之间关系，具体过程为：步骤701、所成型高填方的损失体积和压实度计算：首先，采用所述数据处理设备且根据公式Vj＝∫∫∑(Zd(x,y)-Zj(x,y))ds，并结合步骤二中所存储的所成型高填方的顶面高程Zd(x,y)和步骤六中所存储的沉降过程中所成型高填方的顶面高程Zj(x,y)，得到沉降过程中所成型高填方的损失体积Vj；然后，采用所述数据处理设备且根据公式 并结合步骤四中所成型高填方的质量m1，得到沉降过程中所成型高填方的干密度ρj，并存储至所述存储器内；再采用所述数据处理设备且根据公式 得到沉降过程中所成型高填方的压实度λj，并存储至所述存储器内；
步骤702、沉降过程中所成型高填方的损失体积和所成型高填方的压实度之间关系的获取：采用所述数据处理设备且根据公式ρ1V1＝ρj(V1-Vj)，并结合步骤四中所存储的所成型高填方的初始压实度λ1和步骤四中所存储的所成型高填方的初始体积V1，得到沉降过程中所成型高填方的损失体积和所成型高填方的压实度之间的关系式 并存储至所述存储器内；
步骤八、沉降后所成型高填方的补偿：在所成型高填方发生沉降的过程中，获取沉降后所成型高填方的损失体积和压实度，并对沉降后的所成型高填方进行补偿，具体过程为：步骤801：获取沉降后所成型高填方的损失体积和压实度：当所成型高填达到沉降稳定状态时，采用所述数据处理设备将沉降过程中所成型高填方的最终损失体积记录为沉降后所成型高填方的损失体积，沉降过程中所成型高填方的最终压实度记录为沉降后所成型高填方的压实度，并存储至所述存储器内；
步骤802、沉降后所成型高填方的补偿：根据步骤801中得到的沉降后所成型高填方的损失体积和沉降后所成型高填方的压实度，对沉降后的所成型高填方进行补偿，使补偿后的所施工高填方的高程达到沉降前所施工高填方顶面高程Zd(x,y)，同时沉降后所成型高填方中补偿部分的压实度达到沉降后所成型高填方的压实度。
2.按照权利要求1所述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤二中对各个测量时间内的两幅所述所施工高填方单视SAR图像的雷达数据处理方法均相同，其中，对任一测量时间内的两幅所述所施工高填方单视SAR图像进行雷达数据处理时，具体包括以下步骤：步骤201、图像匹配：采用所述数据处理设备调用图像匹配模块对所施工高填方SAR主图像和所施工高填方SAR辅图像进行匹配，获取匹配后的所施工高填方SAR主图像和匹配后的所施工高填方SAR辅图像；
步骤202、相位干涉图生成：采用所述数据处理设备调用相位干涉图模块，将步骤201中得到的匹配后的所施工高填方SAR主图像和匹配后的所施工高填方SAR辅图像进行相位干涉图生成处理，得到所施工高填方的相位干涉图；
步骤203、相位解缠：采用所述数据处理设备调用相位解缠模块对步骤202中得到的所施工高填方的相位干涉图进行相位解缠处理，得到所施工高填方的解缠图像；
步骤204、获取所施工高填方DEM数据：采用所述数据处理设备调用相位高程转换模块对步骤203中得到的所施工高填方的解缠图像进行相位高程转换处理，得到所施工高填方DEM数据；其中，所述所施工高填方DEM数据包括原始坡地或沟谷高程Z0(x,y)、当前填筑层的顶面高程Zi(x,y)、上一填筑层的顶面高程Zi-1(x,y)和所成型高填方顶面高程Zd(x,y)。
3.按照权利要求1或2所述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤六中对各个测量时间内的两幅所述所成型高填方单视SAR图像的雷达数据处理方法均相同，其中，对任一测量时间内的两幅所述所成型高填方单视SAR图像进行雷达数据处理时，具体包括以下步骤：步骤601、图像匹配：采用所述数据处理设备调用图像匹配模块对所成型高填方SAR主图像和所成型高填方SAR辅图像进行匹配，获取匹配后的所成型高填方SAR主图像和匹配后的所成型高填方SAR辅图像；
步骤602、相位干涉图生成：采用所述数据处理设备调用相位干涉图模块，将步骤601中得到的匹配后的所成型高填方SAR主图像和匹配后的所成型高填方SAR辅图像进行相位干涉图生成处理，得到所成型高填方的相位干涉图；
步骤603、相位解缠：采用所述数据处理设备调用相位解缠模块对步骤602中得到的所成型高填方的相位干涉图进行相位解缠处理，得到所成型高填方的解缠图像；
步骤604、获取沉降过程中所成型高填方DEM数据：采用所述数据处理设备调用相位高程转换模块对步骤603中得到的所成型高填方的解缠图像进行相位高程转换处理，得到沉降过程中所成型高填方DEM数据；其中，所述沉降过程中所成型高填方DEM数据包括沉降过程中所成型高填方顶面高程Zj(x,y)。
4.按照权利要求1或2所述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤304中设定的压实度的取值范围为0.83～0.93。
5.按照权利要求1或2所述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤801中判断所成型高填达到沉降稳定状态，具体过程为：步骤Ⅰ、采用所述数据处理设备根据判定条件 进行判定，当判定条件
不成立时，执行步骤Ⅱ；当判定条件 成立时，执行步骤Ⅲ；其中，α的
步骤Ⅱ：当判定条件 不成立时，说明所成型高填方未达到沉降稳定状态，采用所述机载干涉合成孔径雷达InSAR对所成型高填方继续进行测量，获取沉降过程中所成型高填方的损失体积和所成型高填方的压实度；
步骤Ⅲ：当判定条件 成立时，说明所成型高填达到沉降稳定状态，采用所述数据处理设备将沉降过程中所成型高填方的最终损失体积Vj+1记录为沉降后所成型高填方的损失体积Vj+1，沉降过程中所成型高填方的最终压实度λj+1记录为沉降后所成型高填方的压实度λj+1，并存储至所述存储器内。
6.按照权利要求1或2所述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤三和步骤四中ρimax＝ρmax；
所述当前填筑层所用土的最大干密度ρimax和步骤四中所成型高填方所用土的最大干密度ρmax均为采用击实试验获取。
一种基于InSAR技术高填方压实度检测与沉降补偿方法
技术领域
[0001] 本发明属于高填方工程检测技术领域，尤其是涉及一种基于InSAR技术高填方压实度检测与沉降补偿方法。
背景技术
[0002] 高填方工程即在原始坡地或沟谷中通过分层碾压回填土将所要求的构筑物达到比周围地势明显高出的施工工程，高填方工程的变形监测与压实度检测是工程建设中确保工程质量的两个必备环节。
[0003] 目前高填方的沉降变形监测一般采用实地水准测量和GPS测量，实地水准测量需要野外人工作业，并且周期长、花费大，GPS测量虽然解决了水准测量周期长的问题，然而仍然无法满足大面积、大体量高填方工程的监测要求。因此，传统的沉降变形监测方法，工作量大、效率低、系统误差大、监测点位离散、易受天气变化影响，无法满足大体量、高精度高填方工程的监测要求。另外，高填方填筑过程中，高填方中填筑层的压实度检测一般采用环刀现场取样，一般以填筑层中一个点的压实度来评估一定面积或范围内填筑层面的压实度，再以各个填筑层面的压实度去等效整个高填方的压实度，是“以点代面，以面代体”的检测方法，其检测结果能否真实反应高填方工程的整体压实度水平，往往取决于检测的密度和频率。因此，难以客观评估填料多样性和岩土参数离散性对高填方工程整体压实质量带来的影响；其次，高填方中填筑层的压实度检测需要现场进行测试，用单点或单体填筑检测标准进行大面积、大体量高填方工程压实度检测，势必会造成工程成本大幅飙升；同时，由于人工检测的速度限制，会延长和耽误施工工期；最后，大面积、大体量高填方施工完成后，高填方在长期自重作用下，高填方通常会发生沉降变形，沉降过程较长且沉降不均匀，另外沉降过程中高填方的压实度也发生改变需要进一步提高，尤其，对于某些关键、重大高填方工程(如机场高填方)需要对高填方的沉降变形进行补偿以达到工程设计要求是非常重要的。
[0004] 由上述内容可知，利用传统的压实度检测和沉降变形监测方法，由于无法规避自身缺陷，使得压实度检测结果和沉降变形监测与实际工程有所偏离，从而造成高填方的补偿存在较大误差，更难以建立沉降变形与压实度之间的关系，不能满足施工准确性需求。
发明内容
[0005] 本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种基于InSAR技术高填方压实度检测与沉降补偿方法，其方法步骤简单、设计合理且成本低，使用操作简便，能准确地对所施工高填方的压实度进行检测，同时对沉降过程中所成型高填方的损失体积和压实度进行检测，检测结果准确，并获取所成型高填方的损失体积和所施工高填方的压实度之间关系，为沉降后所成型高填方的补偿提供依据，使沉降后所成型高填方得到有效补偿，且沉降后所成型高填方的补偿效果好，保证所成型高填方补偿后的稳定性，实用性强。
[0006] 为解决上述技术问题，本发明采用的技术方案是：一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于该方法包括以下步骤：
[0007] 步骤一、获取所施工高填方的SAR图像：在所施工高填方填筑至成型过程中，采用机载干涉合成孔径雷达InSAR对所施工高填方进行测量，获取每一个测量时间内的两幅所施工高填方单视SAR图像，并将所获取的每一个测量时间内的两幅所施工高填方单视SAR图像存储至与数据处理设备相接的存储器内；其中，两幅所述所施工高填方单视SAR图像均为正射影像，且两幅所述所施工高填方单视SAR图像分别为所施工高填方SAR主图像和所施工高填方SAR辅图像；
[0008] 步骤二、获取所施工高填方的DEM数据：采用所述数据处理设备按照时间先后顺序对步骤一中各个测量时间内存储至所述存储器中的两幅所述所施工高填方单视SAR图像分别进行雷达数据处理，获取各个测量时间内的所施工高填方DEM数据，并存储至所述存储器内；其中，所施工高填方DEM数据包括原始坡地或沟谷高程Z0(x,y)、当前填筑层的顶面高程Zi(x,y)、上一填筑层的顶面高程Zi-1(x,y)和所成型高填方的顶面高程Zd(x,y)；
[0009] 步骤三、所施工高填方压实度的检测：所施工高填方由下至上分多层填筑层，每个所述填筑层的压实度检测方法均相同，其中对任一所述填筑层的压实度进行检测，具体过程如下：
[0010] 步骤301、当前填筑层的体积计算：采用所述数据处理设备且根据公式Vi＝∫∫∑(Zi(x,y)-Zi-1(x,y))ds，并结合步骤二中所存储的当前填筑层的顶面高程Zi(x,y)和步骤二中所存储的上一填筑层的顶面高程Zi-1(x,y)，得到当前填筑层的体积Vi，并存储至所述存储器内；其中，ds为积分单元面积，i为正整数且i＝1、2、…、M，M为所述填筑层的数量；
[0011] 步骤302、当前填筑层的干密度计算：采用所述数据处理设备且根据公式并结合步骤301中所存储的当前填筑层的体积Vi，得到当前填筑层的干密度ρi，并存储至所述存储器内；
[0012] 步骤303、当前填筑层的压实度计算：采用所述数据处理设备且根据公式并结合步骤302中所存储的当前填筑层的干密度ρi，得到当前填筑层的压实度λi，并存储至所述存储器内，所述数据处理设备还传输至显示器进行输出显示；其中，mi为当前填筑层的质量，ρi max为当前填筑层所用土的最大干密度，ωi为当前填筑层的含水量；
[0013] 步骤304、多次重复步骤301至步骤303，完成多个所述填筑层的压实度检测，直至所施工高填方的压实度满足设定的压实度，所施工高填方施工成型；
[0014] 步骤四、所成型高填方的初始体积和初始压实度计算：所施工高填方施工成型后，对所成型高填方的初始体积和初始压实度进行计算的具体过程为：
[0015] 首先，采用所述数据处理设备且根据公式V1＝∫∫∑(Zd(x,y)-Z0(x,y))ds，并结合步骤二中所存储的所成型高填方的顶面高程Zd(x,y)和步骤二中所存储的原始坡地或沟谷高程Z0(x,y)，得到所成型高填方的初始体积V1，并存储至所述存储器内；然后，采用所述数据处理设备且根据公式 得到所成型高填方的干密度ρ1，并存储至所述存储器内；再采用所述数据处理设备且根据公式 得到所成型高填方的初始压实度λ1，并存储至所述存储器内；其中，m1为所成型高填方的质量，ρmax为所成型高填方所用土的最大干密度，ω1为所成型高填方的含水量；
[0016] 步骤五、获取沉降过程中所成型高填方的SAR图像：所施工高填方施工成型后，所成型高填方发生沉降，在所成型高填方沉降的过程中，采用所述机载干涉合成孔径雷达InSAR对所成型高填方进行测量，获取每一个测量时间内的两幅所成型高填方单视SAR图像，并将所获取的每一个测量时间内的两幅所成型高填方单视SAR图像存储至所述存储器内；其中，两幅所述所成型高填方单视SAR图像均为正射影像，且两幅所述所成型高填方单视SAR图像分别为所成型高填方SAR主图像和所成型高填方SAR辅图像；
[0017] 步骤六、获取沉降过程中所成型高填方的DEM数据：采用所述数据处理设备按照时间先后顺序对步骤五中各个测量时间内存储至所述存储器中的两幅所述所成型高填方单视SAR图像分别进行雷达数据处理，获取各个测量时间内的沉降过程中所成型高填方DEM数据，并存储至所述存储器内；其中，所述沉降过程中所成型高填方DEM数据包括沉降过程中所成型高填方顶面高程Zj(x,y)；
[0018] 步骤七、沉降过程中所成型高填方的损失体积和压实度计算：所施工高填方施工成型后，在所成型高填方沉降的过程中，对所成型高填方的损失体积和压实度进行计算，并获取所成型高填方的损失体积和所施工高填方的压实度之间关系，具体过程为：
[0019] 步骤701、所成型高填方的损失体积和压实度计算：首先，采用所述数据处理设备且根据公式Vj＝∫∫∑(Zd(x,y)-Zj(x,y))ds，并结合步骤二中所存储的所成型高填方的顶面高程Zd(x,y)和步骤六中所存储的沉降过程中所成型高填方的顶面高程Zj(x,y)，得到沉降过程中所成型高填方的损失体积Vj；然后，采用所述数据处理设备且根据公式并结合步骤四中所成型高填方的质量m1，得到沉降过程中所成型高填方的干密度ρj，并存储至所述存储器内；再采用所述数据处理设备且根据公式 得到沉降过程中所成型高填方的压实度λj，并存储至所述存储器内；
[0020] 步骤702、沉降过程中所成型高填方的损失体积和所成型高填方的压实度之间关系的获取：采用所述数据处理设备且根据公式ρ1V1＝ρj(V1-Vj)，并结合步骤四中所存储的所成型高填方的初始压实度λ1和步骤四中所存储的所成型高填方的初始体积V1，得到沉降过程中所成型高填方的损失体积和所成型高填方的压实度之间的关系式 并存储至所述存储器内；
[0021] 步骤八、沉降后所成型高填方的补偿：在所成型高填方发生沉降的过程中，获取沉降后所成型高填方的损失体积和压实度，并对沉降后的所成型高填方进行补偿，具体过程为：
[0022] 步骤801：获取沉降后所成型高填方的损失体积和压实度：当所成型高填达到沉降稳定状态时，采用所述数据处理设备将沉降过程中所成型高填方的最终损失体积记录为沉降后所成型高填方的损失体积，沉降过程中所成型高填方的最终压实度记录为沉降后所成型高填方的压实度，并存储至所述存储器内；
[0023] 步骤802、沉降后所成型高填方的补偿：根据步骤801中得到的沉降后所成型高填方的损失体积和沉降后所成型高填方的压实度，对沉降后的所成型高填方进行补偿，使补偿后的所施工高填方的高程达到沉降前所施工高填方顶面高程Zd(x,y)，同时沉降后所成型高填方中补偿部分的压实度达到沉降后所成型高填方的压实度。
[0024] 上述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤二中对各个测量时间内的两幅所述所施工高填方单视SAR图像的雷达数据处理方法均相同，其中，对任一测量时间内的两幅所述所施工高填方单视SAR图像进行雷达数据处理时，具体包括以下步骤：
[0025] 步骤201、图像匹配：采用所述数据处理设备调用图像匹配模块对所施工高填方SAR主图像和所施工高填方SAR辅图像进行匹配，获取匹配后的所施工高填方SAR主图像和匹配后的所施工高填方SAR辅图像；
[0026] 步骤202、相位干涉图生成：采用所述数据处理设备调用相位干涉图模块，将步骤201中得到的匹配后的所施工高填方SAR主图像和匹配后的所施工高填方SAR辅图像进行相位干涉图生成处理，得到所施工高填方的相位干涉图；
[0027] 步骤203、相位解缠：采用所述数据处理设备调用相位解缠模块对步骤202中得到的所施工高填方的相位干涉图进行相位解缠处理，得到所施工高填方的解缠图像；
[0028] 步骤204、获取所施工高填方DEM数据：采用所述数据处理设备调用相位高程转换模块对步骤203中得到的所施工高填方的解缠图像进行相位高程转换处理，得到所施工高填方DEM数据；其中，所述所施工高填方DEM数据包括原始坡地或沟谷高程Z0(x,y)、当前填筑层的顶面高程Zi(x,y)、上一填筑层的顶面高程Zi-1(x,y)和所成型高填方顶面高程Zd(x,y)。
[0029] 上述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤六中对各个测量时间内的两幅所述所成型高填方单视SAR图像的雷达数据处理方法均相同，其中，对任一测量时间内的两幅所述所成型高填方单视SAR图像进行雷达数据处理时，具体包括以下步骤：
[0030] 步骤601、图像匹配：采用所述数据处理设备调用图像匹配模块对所成型高填方SAR主图像和所成型高填方SAR辅图像进行匹配，获取匹配后的所成型高填方SAR主图像和匹配后的所成型高填方SAR辅图像；
[0031] 步骤602、相位干涉图生成：采用所述数据处理设备调用相位干涉图模块，将步骤601中得到的匹配后的所成型高填方SAR主图像和匹配后的所成型高填方SAR辅图像进行相位干涉图生成处理，得到所成型高填方的相位干涉图；
[0032] 步骤603、相位解缠：采用所述数据处理设备调用相位解缠模块对步骤602中得到的所成型高填方的相位干涉图进行相位解缠处理，得到所成型高填方的解缠图像；
[0033] 步骤604、获取沉降过程中所成型高填方DEM数据：采用所述数据处理设备调用相位高程转换模块对步骤603中得到的所成型高填方的解缠图像进行相位高程转换处理，得到沉降过程中所成型高填方DEM数据；其中，所述沉降过程中所成型高填方DEM数据包括沉降过程中所成型高填方顶面高程Zj(x,y)。
[0034] 上述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤304中设定的压实度的取值范围为0.83～0.93。
[0035] 上述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤801中判断所成型高填达到沉降稳定状态，具体过程为：
[0036] 步骤Ⅰ、采用所述数据处理设备根据判定条件 进行判定，当判定条件不成立时，执行步骤Ⅱ；当判定条件 成立时，执行步骤Ⅲ；
[0037] 步骤Ⅱ：当判定条件 不成立时，说明所成型高填方未达到沉降稳定状态，采用所述机载干涉合成孔径雷达InSAR对所成型高填方继续进行测量，获取沉降过程中所成型高填方的损失体积和所成型高填方的压实度；
[0038] 步骤Ⅲ：当判定条件 成立时，说明所成型高填达到沉降稳定状态，采用所述数据处理设备将沉降过程中所成型高填方的最终损失体积Vj+1记录为沉降后所成型高填方的损失体积Vj+1，沉降过程中所成型高填方的最终压实度λj+1记录为沉降后所成型高填方的压实度λj+1，并存储至所述存储器内。
[0039] 上述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤Ⅰ中α的取值范围为0.05％～0.15％。
[0040] 上述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤三和步骤四中ρi max＝ρmax；
[0041] 所述当前填筑层所用土的最大干密度ρi max和步骤四中所成型高填方所用土的最大干密度ρmax均为采用击实试验获取。
[0042] 本发明与现有技术相比具有以下优点：
[0043] 1、本发明方法步骤简单且实现方便，检测效率和准确度高，使用效果好。
[0044] 2、本发明压实度检测方法设计合理且检测简便、快速、准确，所施工高填方由下至上分多层填筑层，通过对每一填筑层的压实度进行检测，保证所施工高填方的压实度满足设定的压实度，满足所施工高填方施工要求。
[0045] 3、本发明对任一所述填筑层的压实度进行检测时，首先采用数据处理设备计算得到当前填筑层的体积，再采用数据处理设备计算当前填筑层的干密度，然后，采用数据处理设备计算当前填筑层的压实度；当当前填筑层的压实度满足设定的压实度时，开始进行下一层填筑层的填筑并计算下一层填筑层的压实度；这样，实现每一填筑层的压实度检测，直至完成所施工高填方的填筑，这样一方面既避免工作人员现场取样检测的复杂度，又避免人力物力和时间消耗多，填筑层的压实度检测效率高；另一方面填筑层的压实度检测由点改为面，不再以每一填筑层中一个点的压实度来评估每一填筑层的压实度，使填筑层的压实度检测准确度高，满足大面积、大体量高填方工程压实度检测需求。
[0046] 4、所施工高填方施工成型后，首先采用数据处理设备计算得到所成型高填方的体积，再采用数据处理设备计算所成型高填方的干密度，然后，采用数据处理设备计算所成型高填方的压实度，使所成型高填方的检测方式由面改为体，增加了测量的准确性，消除了因各个填筑层面的压实度去等效整个所成型高填方的压实度而引起的检测误差，能满足高填方压实度检测需求。
[0047] 5、所成型高填方沉降的过程中，首先采用数据处理设备计算得到沉降过程中所成型高填方的损失体积，再采用数据处理设备计算沉降过程中所成型高填方的干密度，然后，采用数据处理设备计算沉降过程中所成型高填方的压实度，且采用所述数据处理设备根据沉降过程中所成型高填方的质量守恒，获取降过程中所成型高填方的损失体积和所成型高填方的压实度之间的关系式，为沉降后所成型高填方的补偿提供依据。
[0048] 6、在所成型高填方发生沉降的过程中，当所成型高填达到沉降稳定状态，采用数据处理设备获取沉降后所成型高填方的损失体积和压实度，并且根据沉降后所成型高填方的损失体积和压实度，并对沉降后的所成型高填方进行补偿，使补偿后的所施工高填方的高程达到沉降前所施工高填方顶面高程，同时沉降后所成型高填方中补偿部分的压实度达到沉降后所成型高填方的压实度，使沉降后所成型高填方得到有效补偿，且沉降后所成型高填方的补偿效果好，保证所成型高填方补偿后的稳定性，实用性强。
[0049] 7、本发明操作简便且使用效果好，为提高检测数据准确性，通过检测每一填筑层的压实度，进而保证每一填筑层的质量，从而最终保证了所施工高填方体的质量，消除了对所施工高填方体最后检测而存在的误差，并且检测效率高、检测结果准确；另外，在所成型高填方过程沉降的过程中，能对所成型高填方的损失体积和压实度进行检测，获取降过程中所成型高填方的损失体积和所成型高填方的压实度之间的关系式，为沉降后所成型高填方的补偿提供依据，使沉降后所成型高填方得到有效补偿，提高所成型高填方补偿后的稳定性。因而，采用本发明能有效提高检测数据的准确性，并提高沉降补偿效率和质量。
[0050] 综上所述，本发明方法步骤简单、设计合理且成本低，使用操作简便，能准确地对所施工高填方的压实度进行检测，同时对沉降过程中所成型高填方的损失体积和压实度进行检测，检测结果准确，并获取所成型高填方的损失体积和所施工高填方的压实度之间关系，为沉降后所成型高填方的补偿提供依据，使沉降后所成型高填方得到有效补偿，且沉降后所成型高填方的补偿效果好，保证所成型高填方补偿后的稳定性，实用性强。
[0051] 下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
[0052] 图1为本发明一种基于InSAR技术高填方压实度检测与沉降补偿方法的流程框图。
具体实施方式
[0053] 如图1所示，一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：
[0054] 步骤一、获取所施工高填方的SAR图像：在所施工高填方填筑至成型过程中，采用机载干涉合成孔径雷达InSAR对所施工高填方进行测量，获取每一个测量时间内的两幅所施工高填方单视SAR图像，并将所获取的每一个测量时间内的两幅所施工高填方单视SAR图像存储至与数据处理设备相接的存储器内；其中，两幅所述所施工高填方单视SAR图像均为正射影像，且两幅所述所施工高填方单视SAR图像分别为所施工高填方SAR主图像和所施工高填方SAR辅图像；
[0055] 步骤二、获取所施工高填方的DEM数据：采用所述数据处理设备按照时间先后顺序对步骤一中各个测量时间内存储至所述存储器中的两幅所述所施工高填方单视SAR图像分别进行雷达数据处理，获取各个测量时间内的所施工高填方DEM数据，并存储至所述存储器内；其中，所施工高填方DEM数据包括原始坡地或沟谷高程Z0(x,y)、当前填筑层的顶面高程Zi(x,y)、上一填筑层的顶面高程Zi-1(x,y)和所成型高填方的顶面高程Zd(x,y)；
[0056] 步骤三、所施工高填方压实度的检测：所施工高填方由下至上分多层填筑层，每个所述填筑层的压实度检测方法均相同，其中对任一所述填筑层的压实度进行检测，具体过程如下：
[0057] 步骤301、当前填筑层的体积计算：采用所述数据处理设备且根据公式Vi＝∫∫∑(Zi(x,y)-Zi-1(x,y))ds，并结合步骤二中所存储的当前填筑层的顶面高程Zi(x,y)和步骤二中所存储的上一填筑层的顶面高程Zi-1(x,y)，得到当前填筑层的体积Vi，并存储至所述存储器内；其中，ds为积分单元面积，i为正整数且i＝1、2、…、M，M为所述填筑层的数量；
[0058] 步骤302、当前填筑层的干密度计算：采用所述数据处理设备且根据公式并结合步骤301中所存储的当前填筑层的体积Vi，得到当前填筑层的干密度ρi，并存储至所述存储器内；
[0059] 步骤303、当前填筑层的压实度计算：采用所述数据处理设备且根据公式并结合步骤302中所存储的当前填筑层的干密度ρi，得到当前填筑层的压实度λi，并存储至所述存储器内，所述数据处理设备还传输至显示器进行输出显示；其中，mi为当前填筑层的质量，ρi max为当前填筑层所用土的最大干密度，ωi为当前填筑层的含水量；
[0060] 步骤304、多次重复步骤301至步骤303，完成多个所述填筑层的压实度检测，直至所施工高填方的压实度满足设定的压实度，所施工高填方施工成型；
[0061] 本实施例中，对任一所述填筑层的压实度进行检测时，当当前填筑层的压实度满足设定的压实度时，开始进行下一层填筑层的填筑并计算下一层填筑层的压实度；这样，实现每一填筑层的压实度检测，直至完成所施工高填方的填筑，这样一方面既避免工作人员现场取样检测的复杂度，又避免人力物力和时间消耗多，填筑层的压实度检测效率高；另一方面填筑层的压实度检测由点改为面，不再以每一填筑层中一个点的压实度来评估每一填筑层的压实度，使填筑层的压实度检测准确度高，满足大面积、大体量高填方工程压实度检测需求。
[0062] 步骤四、所成型高填方的初始体积和初始压实度计算：所施工高填方施工成型后，对所成型高填方的初始体积和初始压实度进行计算的具体过程为：
[0063] 首先，采用所述数据处理设备且根据公式V1＝∫∫∑(Zd(x,y)-Z0(x,y))ds，并结合步骤二中所存储的所成型高填方的顶面高程Zd(x,y)和步骤二中所存储的原始坡地或沟谷高程Z0(x,y)，得到所成型高填方的初始体积V1，并存储至所述存储器内；然后，采用所述数据处理设备且根据公式 得到所成型高填方的干密度ρ1，并存储至所述存储器内；再采用所述数据处理设备且根据公式 得到所成型高填方的初始压实度λ1，并存储至所述存储器内；其中，m1为所成型高填方的质量，ρmax为所成型高填方所用土的最大干密度，ω1为所成型高填方的含水量；
[0064] 本实施例中，所施工高填方施工成型后计算所成型高填方的压实度，使所成型高填方的检测方式由面改为体，增加了测量的准确性，消除了因各个填筑层面的压实度去等效整个所成型高填方的压实度而引起的检测误差，能满足高填方压实度检测需求。
[0065] 步骤五、获取沉降过程中所成型高填方的SAR图像：所施工高填方施工成型后，所成型高填方发生沉降，在所成型高填方沉降的过程中，采用所述机载干涉合成孔径雷达InSAR对所成型高填方进行测量，获取每一个测量时间内的两幅所成型高填方单视SAR图像，并将所获取的每一个测量时间内的两幅所成型高填方单视SAR图像存储至所述存储器内；其中，两幅所述所成型高填方单视SAR图像均为正射影像，且两幅所述所成型高填方单视SAR图像分别为所成型高填方SAR主图像和所成型高填方SAR辅图像；
[0066] 步骤六、获取沉降过程中所成型高填方的DEM数据：采用所述数据处理设备按照时间先后顺序对步骤五中各个测量时间内存储至所述存储器中的两幅所述所成型高填方单视SAR图像分别进行雷达数据处理，获取各个测量时间内的沉降过程中所成型高填方DEM数据，并存储至所述存储器内；其中，所述沉降过程中所成型高填方DEM数据包括沉降过程中所成型高填方顶面高程Zj(x,y)；
[0067] 步骤七、沉降过程中所成型高填方的损失体积和压实度计算：所施工高填方施工成型后，在所成型高填方沉降的过程中，对所成型高填方的损失体积和压实度进行计算，并获取所成型高填方的损失体积和所施工高填方的压实度之间关系，具体过程为：
[0068] 步骤701、所成型高填方的损失体积和压实度计算：首先，采用所述数据处理设备且根据公式Vj＝∫∫∑(Zd(x,y)-Zj(x,y))ds，并结合步骤二中所存储的所成型高填方的顶面高程Zd(x,y)和步骤六中所存储的沉降过程中所成型高填方的顶面高程Zj(x,y)，得到沉降过程中所成型高填方的损失体积Vj；然后，采用所述数据处理设备且根据公式并结合步骤四中所成型高填方的质量m1，得到沉降过程中所成型高填方的干密度ρj，并存储至所述存储器内；再采用所述数据处理设备且根据公式 得到沉降过程中所成型高填方的压实度λj，并存储至所述存储器内；
[0069] 步骤702、沉降过程中所成型高填方的损失体积和所成型高填方的压实度之间关系的获取：采用所述数据处理设备且根据公式ρ1V1＝ρj(V1-Vj)，并结合步骤四中所存储的所成型高填方的初始压实度λ1和步骤四中所存储的所成型高填方的初始体积V1，得到沉降过程中所成型高填方的损失体积和所成型高填方的压实度之间的关系式 并存储至所述存储器内；
[0070] 本实施例中，采用所述数据处理设备计算沉降过程中所成型高填方的压实度λj，且采用所述数据处理设备根据沉降过程中所成型高填方的质量守恒ρ1V1＝ρj(V1-Vj)，获取降过程中所成型高填方的损失体积和所成型高填方的压实度之间的关系式，为沉降后所成型高填方的补偿提供依据。
[0071] 本实施例中，采用所述数据处理设备且根据公式ρ1V1＝ρj(V1-Vj)，得到沉降过程中所成型高填方的损失体积和所成型高填方的压实度之间的关系式 时，首先根据公式 和公式 得到
[0072] 步骤八、沉降后所成型高填方的补偿：在所成型高填方发生沉降的过程中，获取沉降后所成型高填方的损失体积和压实度，并对沉降后的所成型高填方进行补偿，具体过程为：
[0073] 步骤801：获取沉降后所成型高填方的损失体积和压实度：当所成型高填达到沉降稳定状态时，采用所述数据处理设备将沉降过程中所成型高填方的最终损失体积记录为沉降后所成型高填方的损失体积，沉降过程中所成型高填方的最终压实度记录为沉降后所成型高填方的压实度，并存储至所述存储器内；
[0074] 步骤802、沉降后所成型高填方的补偿：根据步骤801中得到的沉降后所成型高填方的损失体积和沉降后所成型高填方的压实度，对沉降后的所成型高填方进行补偿，使补偿后的所施工高填方的高程达到沉降前所施工高填方顶面高程Zd(x,y)，同时沉降后所成型高填方中补偿部分的压实度达到沉降后所成型高填方的压实度。
[0075] 本实施例中，步骤一中采用机载干涉合成孔径雷达InSAR对所施工高填方进行测量，不需要工作人员现场检测采集数据，便能定时地对所施工高填方中每一填筑层进行检测，从而获取原始坡地或沟谷高程Z0(x,y)、当前填筑层的顶面高程Zi(x,y)、上一填筑层的顶面高程Zi-1(x,y)和所成型高填方的顶面高程Zd(x,y)等所施工高填方DEM数据。
[0076] 本实施例中，步骤五中采用机载干涉合成孔径雷达InSAR对所成型高填方进行测量，不需要工作人员现场检测采集数据，便能定时地对沉降过程中所成型高填方进行检测，从而获取沉降过程中所成型高填方顶面高程Zj(x,y)等沉降过程中所成型高填方DEM数据。
[0077] 实际使用过程中，可采用星载干涉合成孔径雷达InSAR对施工高填方和所成型高填方进行测量，获取检测数据，方便快捷。
[0078] 本实施例中，步骤四中所述所成型高填方的质量m1＝M×mi。
[0079] 本实施例中，步骤二中对各个测量时间内的两幅所述所施工高填方单视SAR图像的雷达数据处理方法均相同，其中，对任一测量时间内的两幅所述所施工高填方单视SAR图像进行雷达数据处理时，具体包括以下步骤：
[0080] 步骤201、图像匹配：采用所述数据处理设备调用图像匹配模块对所施工高填方SAR主图像和所施工高填方SAR辅图像进行匹配，获取匹配后的所施工高填方SAR主图像和匹配后的所施工高填方SAR辅图像；
[0081] 步骤202、相位干涉图生成：采用所述数据处理设备调用相位干涉图模块，将步骤201中得到的匹配后的所施工高填方SAR主图像和匹配后的所施工高填方SAR辅图像进行相位干涉图生成处理，得到所施工高填方的相位干涉图；
[0082] 步骤203、相位解缠：采用所述数据处理设备调用相位解缠模块对步骤202中得到的所施工高填方的相位干涉图进行相位解缠处理，得到所施工高填方的解缠图像；
[0083] 步骤204、获取所施工高填方DEM数据：采用所述数据处理设备调用相位高程转换模块对步骤203中得到的所施工高填方的解缠图像进行相位高程转换处理，得到所施工高填方DEM数据；其中，所述所施工高填方DEM数据包括原始坡地或沟谷高程Z0(x,y)、当前填筑层的顶面高程Zi(x,y)、上一填筑层的顶面高程Zi-1(x,y)和所成型高填方顶面高程Zd(x,y)。
[0084] 本实施例中，步骤201图像匹配是为了将所施工高填方SAR主图像和所施工高填方SAR辅图像进行滤波，除掉所施工高填方SAR主图像和所施工高填方SAR辅图像中干斑噪声，提高所施工高填方SAR主图像和所施工高填方SAR辅图像的精度，同时也为了进行步骤202做基础处理。
[0085] 本实施例中，步骤202相位干涉图生成是为了将步骤201中得到的匹配后的所施工高填方SAR主图像和匹配后的所施工高填方SAR辅图像进行复共轭相乘，生成所施工高填方的相位干涉图，为后续根据所施工高填方的相位干涉图中干涉相位图得到相位差反演所施工高填方中填筑层的相对高程的必要处理步骤。
[0086] 本实施例中，步骤203中相位解缠因为所施工高填方的相位干涉图中干涉相位差值是真实相位差的卷叠，所以通过步骤203相位解缠对所施工高填方的相位干涉图进行解缠，即给每一个相位测量值加上整数倍的相位周期，相位解缠最后得到一副真正反映所施工高填方至两天线之间的斜距信息的所施工高填方二维绝对相位图像。
[0087] 本实施例中，步骤204中相位高程转换模块对所述所施工高填方二维绝对相位图像处理得到相位干涉图中每个像素点对应的所施工高填方高程值，最后将雷达影像的坐标转换到大地坐标系中，则获取所述所施工高填方DEM数据包括原始坡地或沟谷高程Z0(x,y)、当前填筑层的顶面高程Zi(x,y)、上一填筑层的顶面高程Zi-1(x,y)和所成型高填方顶面高程Zd(x,y)等所施工高填方DEM数据。
[0088] 本实施例中，步骤六中对各个测量时间内的两幅所述所成型高填方单视SAR图像的雷达数据处理方法均相同，其中，对任一测量时间内的两幅所述所成型高填方单视SAR图像进行雷达数据处理时，具体包括以下步骤：
[0089] 步骤601、图像匹配：采用所述数据处理设备调用图像匹配模块对所成型高填方SAR主图像和所成型高填方SAR辅图像进行匹配，获取匹配后的所成型高填方SAR主图像和匹配后的所成型高填方SAR辅图像；
[0090] 步骤602、相位干涉图生成：采用所述数据处理设备调用相位干涉图模块，将步骤601中得到的匹配后的所成型高填方SAR主图像和匹配后的所成型高填方SAR辅图像进行相位干涉图生成处理，得到所成型高填方的相位干涉图；
[0091] 步骤603、相位解缠：采用所述数据处理设备调用相位解缠模块对步骤602中得到的所成型高填方的相位干涉图进行相位解缠处理，得到所成型高填方的解缠图像；
[0092] 步骤604、获取沉降过程中所成型高填方DEM数据：采用所述数据处理设备调用相位高程转换模块对步骤603中得到的所成型高填方的解缠图像进行相位高程转换处理，得到沉降过程中所成型高填方DEM数据；其中，所述沉降过程中所成型高填方DEM数据包括沉降过程中所成型高填方顶面高程Zj(x,y)。
[0093] 本实施例中，步骤601图像匹配是为了将所成型高填方SAR主图像和所成型高填方SAR辅图像进行滤波，除掉所成型高填方SAR主图像和所成型高填方SAR辅图像中干斑噪声，提高所成型高填方SAR主图像和所成型高填方SAR辅图像的精度，同时也为了进行步骤602做基础处理。
[0094] 本实施例中，步骤602相位干涉图生成是为了将步骤601中得到的匹配后的所成型高填方SAR主图像和匹配后的所成型高填方SAR辅图像进行复共轭相乘，生成所成型高填方的相位干涉图，为后续根据所成型高填方的相位干涉图中干涉相位图得到相位差反演所成型高填方中填筑层的相对高程的必要处理步骤。
[0095] 本实施例中，步骤603中相位解缠因为所成型高填方的相位干涉图中干涉相位差值是真实相位差的卷叠，所以通过步骤603相位解缠对所成型高填方的相位干涉图进行解缠，即给每一个相位测量值加上整数倍的相位周期，相位解缠最后得到一副真正反映所成型高填方至两天线之间的斜距信息的所成型高填方二维绝对相位图像。
[0096] 本实施例中，步骤604中相位高程转换模块对所述所成型高填方二维绝对相位图像处理得到相位干涉图中每个像素点对应的所成型高填方高程值，最后将雷达影像的坐标转换到大地坐标系中，则获取沉降过程中所成型高填方顶面高程Zj(x,y)等沉降过程中所成型高填方DEM数据。
[0097] 本实施例中，步骤304中设定的压实度λs的取值范围为0.83～0.93。
[0098] 上述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤801中判断所成型高填达到沉降稳定状态，具体过程为：
[0099] 步骤Ⅰ、采用所述数据处理设备根据判定条件 进行判定，当判定条件不成立时，执行步骤Ⅱ；当判定条件 成立时，执行步骤Ⅲ；
[0100] 步骤Ⅱ：当判定条件 不成立时，说明所成型高填方未达到沉降稳定状态，采用所述机载干涉合成孔径雷达InSAR对所成型高填方继续进行测量，获取沉降过程中所成型高填方的损失体积和所成型高填方的压实度；
[0101] 步骤Ⅲ：当判定条件 成立时，说明所成型高填达到沉降稳定状态，采用所述数据处理设备将沉降过程中所成型高填方的最终损失体积Vj+1记录为沉降后所成型高填方的损失体积，沉降过程中所成型高填方的最终压实度λj+1记录为沉降后所成型高填方的压实度λj+1，并存储至所述存储器内。
[0102] 本实施例中，在所成型高填方发生沉降的过程中，用数据处理设备判断所成型高填是否达到沉降稳定状态，当所成型高填达到沉降稳定状态，采用数据处理设备获取沉降后所成型高填方的损失体积和压实度，并且根据沉降后所成型高填方的损失体积和压实度，并对沉降后的所成型高填方进行补偿，使补偿后的所施工高填方的高程达到沉降前所施工高填方顶面高程，同时沉降后所成型高填方中补偿部分的压实度达到沉降后所成型高填方的压实度，使沉降后所成型高填方得到有效补偿，且沉降后所成型高填方的补偿效果好，保证所成型高填方补偿后的稳定性，实用性强。
[0103] 本实施例中，步骤Ⅰ中α的取值范围为0.05％～0.15％。所述α判稳标准0.05％～0.15％，需要根据高填方工程的建筑目的而确定。
[0104] 上述的一种基于InSAR技术高填方压实度检测与沉降补偿方法，其特征在于：步骤三和步骤四中ρi max＝ρmax；
[0105] 所述当前填筑层所用土的最大干密度ρi max和步骤四中所成型高填方所用土的最大干密度ρmax均为采用击实试验获取。
[0106] 本实施例中，所成型高填方的含水量ω1和当前填筑层的含水量ωi可采取烘干法获取。
[0107] 由上述内容可知，本发明方法简单、设计合理且使用效果好且检测结果准确，能对所施工高填方中的每一填筑层进行检测，保证每一填筑层质量满足施工要求，以此类推，直至完成所施工高填方的检测，能准确地对所施工高填方的压实度进行检测，从而保证所施工高填方施工完成后的所成型高填方体满足施工质量要求；同时对沉降过程中所成型高填方的损失体积和压实度进行检测，检测结果准确，并获取所成型高填方的损失体积和所施工高填方的压实度之间关系，为沉降后所成型高填方的补偿提供依据，使沉降后所成型高填方得到有效补偿，且沉降后所成型高填方的补偿效果好，提高沉降补偿效率和质量，保证所成型高填方补偿后的稳定性，实用性强。
[0108] 以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。
