1.本发明涉及地层变形监测领域，具体说是涉及一种区域地层岩土层压缩变形量分布的计算方法。
背景技术：
2.随着人类对地下水资源的大量攫取及采矿过程中地下水资源的散失，含水层失水后孔隙水压力减小造成地层骨架所承担的有效应力显著增大，从而造成含水层的压缩。特别是对于松散含水层来说，其孔隙率大，抗压性能弱，可压缩性好，层内水位降低将会造成含水层快速压缩，压缩量沿地层传递到地表后造成地表沉陷，威胁到地表及地下建构筑物的安全使用。因此，研究承压含水松散层区域水位的变化特征，从而计算承压含水松散层区域失水压缩量分布，对后期地表沉陷推演及地表、地下建构筑物变形控制方案的制定，有着重要的理论和现实意义。
技术实现要素：
3.本发明旨在提供一种简单高效的承压含水松散层区域失水压缩量分布的计算方法，实现对区域地层变形的准确监测。
4.为了实现上述技术目的，首先在计算区域内布设水文孔，并进行全孔取样，获取各岩土层的层厚数据，通过力学实验获取目标土层的密度、厚度、应力应变特征，上覆岩土层的密度、厚度等力学参数，持续进行区域内水位监测，获取目标含水层水位变化。通过目标含水层水位的变化量计算孔隙水压力的减少量，获取土体骨架所承担的应力增量。根据目标含水层不同部分的土地弹性模量计算该含水层的总压缩量，获取区域内不同孔位目标含水层的压缩量，然后绘制等值线形成承压含水松散层区域失水压缩量分布图。
5.具体的，为了实现以上目的，本发明提供了一种承压含水松散层区域失水压缩量分布的计算方法，包括以下步骤:
6.步骤1、获取目标含水层水位；
7.在监测区域地表均匀施工多个向下穿透目标含水层的水文观测孔，对所有水文观测孔进行编号，并进行观测，获取第i个水文观测孔第j次观测时的水头高度
8.步骤2、获取第i个水文观测孔所在位置对应目标含水层的厚度hi，并将该处目标含水层自上而下分为n个分量层，对应第k个分量层的厚度为h
ki
，且并获取每个分量层的应力应变关系，ε
ik
＝f(σ
ik
)，ε
ik
为第i个水文观测孔处第k个分量层的应变量，σ
ik
为第i个水文观测孔处第k个分量层中土骨架承担的有效应力；
9.步骤3、计量参数的获取
10.令任意一次观测时第i个水文观测孔对应的n个分量层中任意一个分量层所承受的上覆岩土体总应力均为σi，且其中，l为第i个水文观测孔上方的上覆岩
土体总层数，并对每层上覆岩土体编号，其中，第o层上覆岩土体的厚度为hso且饱和密度为ρso，g为重力加速度；
11.令第j次观测时第i个水文观测孔对应的n个分量层中任意一个分量层所承受的孔隙水压力均为且其中，ρw为目标含水层水体平均密度；
12.令第j次观测时第i个水文观测孔对应的n个分量层中第k个分量层中土骨架承担的有效应力为
13.步骤4、压缩量的计算
14.令第j次观测到第j+m次观测期间第i个水文观测孔对应的n个分量层中第k个分量层的应变变化量为对应目标含水层在第j次观测到第j+m次观测期间于第i个水文观测孔处的压缩量为
15.进一步，还包括步骤5、计算目标含水层于第j次观测到第j+m次观测期间各个水文观测孔处的压缩量，并绘制压缩量等值线，形成承压含水松散层区域失水压缩量分布图。
16.进一步，目标含水层中各水文观测孔对应的分量层总数为均相同、均不同或部分相同中任意一情况。
17.进一步，同一水文观测孔对应的所有分量层中各层的厚度为均相同、均不同或部分相同中任意一情况。
18.进一步，步骤1中还包括在钻孔内下放用于隔绝水力沟通的套管。
19.进一步，步骤1中水文数据通过水位传感器。
20.本发明相比现有技术具有以下优点：
21.本发明能够较为精确的确定监测区域含水层的实际压缩量，为区域地表沉降的反演预计提供重要的数据支撑，为区域地表建构筑物的破坏预警及保护提供重要参考。
具体实施方式
22.下面将结合本发明实施例，对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
23.一种承压含水松散层区域失水压缩量分布的计算方法，包括以下步骤:
24.步骤1、获取目标含水层水位；
25.在监测区域地表均匀施工多个向下穿透目标含水层的水文观测孔，对所有水文观测孔进行编号，同时进行全孔取芯，在水文观测孔内下放套管至目标含水层，使其与其他含水层隔绝不会产生水力沟通，放入水位传感器，对水文观测孔内水位进行实时监测，获取第i个水文观测孔第j次观测时的水头高度
26.步骤2、根据取芯结果，测定其物理力学特性，获取第i个水文观测孔所在位置对应目标含水层的厚度hi，并将该处目标含水层自上而下分为n个分量层，第k个分量层的厚度
为h
ki
，且分析每个分量层排水条件下的应力应变曲线，并获取每个分量层的应力应变关系，ε
ik
＝f(σ
ik
)，ε
ik
为第i个水文观测孔处第k个分量层的应变量，σ
ik
为第i个水文观测孔处第k个分量层中土骨架承担的有效应力；
27.步骤3、计量参数的获取
28.根据饱和土有效应力原理土骨架所承担的有效应力σ'＝σ-μw，进行以下设置(其中σ为土体所承受的总应力，μw为孔隙水压力，σ'为土骨架所承受的有效应力)：
29.考虑到目标含水层上覆岩土体在观测过程中基本不会发生变化则其产生的土压力在观测期间也不会发生变化，因此在本实施例中，令任意一次观测时第i个水文观测孔对应的n个分量层中任意一个分量层所承受的上覆岩土体总应力均为σi，且其中，l为第i个水文观测孔上方的上覆岩土体总层数，并对每层上覆岩土体编号，第o层上覆岩土体的厚度为hso而饱和密度为ρso，g为重力加速度；
30.考虑到承压水在不同观测时间点其含量可能会发生变化(也是产生失水压缩的主要因素)，发生变化时会导致同一位置产生的水文观测数据不同，孔隙水压力也会发生变化，而对应目标含水层水体密度变化则可忽略不计，因此在本实施例中，令第j次观测时第i个水文观测孔对应的n个分量层中任意一个分量层所承受的孔隙水压力均为且其中，ρw为目标含水层水体平均密度；
31.令第j次观测时第i个水文观测孔对应的n个分量层中第k个分量层内土骨架承担的有效应力为
32.步骤4、压缩量的计算
33.令第j次观测到第j+m次观测期间第i个水文观测孔对应的n个分量层中第k个分量层的应变变化量为则目标含水层在第j次观测到第j+m次观测期间于第i个水文观测孔处的压缩量为
34.步骤5、计算目标含水层于第j次观测到第j+m次观测期间各个水文观测孔处的压缩量，并绘制压缩量等值线，形成承压含水松散层区域失水压缩量分布图。
35.特别的，考虑到目标含水层不同位置的深度可能不同，则根据实际需要目标含水层中各水文观测孔对应的分量层总数为均相同、均不同或部分相同中任意一情况。
36.特别的，考虑到目标含水层同一位置不同高度土骨架强度不同，则根据实际需要同一水文观测孔对应的所有分量层中各层的厚度为均相同、均不同或部分相同中任意一情况。
37.以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。