CN116933535A - 一种地质灾害位移监测方法、装置、设备以及存储介质 - Google Patents

一种地质灾害位移监测方法、装置、设备以及存储介质 Download PDF

Info

Publication number
CN116933535A
CN116933535A CN202310912621.6A CN202310912621A CN116933535A CN 116933535 A CN116933535 A CN 116933535A CN 202310912621 A CN202310912621 A CN 202310912621A CN 116933535 A CN116933535 A CN 116933535A
Authority
CN
China
Prior art keywords
deformation
surface layer
layer
sliding
displacement
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
CN202310912621.6A
Other languages
English (en)
Other versions
CN116933535B (zh
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.)
Guangdong Nonferrous Mine Geological Disaster Prevention And Control Center
Original Assignee
Guangdong Nonferrous Mine Geological Disaster Prevention And Control Center
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 Guangdong Nonferrous Mine Geological Disaster Prevention And Control Center filed Critical Guangdong Nonferrous Mine Geological Disaster Prevention And Control Center
Priority to CN202310912621.6A priority Critical patent/CN116933535B/zh
Publication of CN116933535A publication Critical patent/CN116933535A/zh
Application granted granted Critical
Publication of CN116933535B publication Critical patent/CN116933535B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/32—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid
    • 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
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Tourism & Hospitality (AREA)
  • Pure & Applied Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Computing Systems (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

本申请涉及一种地质灾害位移监测方法、装置、设备以及存储介质，所述方法包括获取选定监测区域的地质条件参数和所述选定监测区域的地表层对应的地表层数据，根据所述地质条件参数和所述地表层数据，计算对应的地表层在当前地质条件下的滑动形变系数，获取所述地表层的倾斜角数据和对应的植被覆盖面积，根据所述倾斜角数据和所述植被覆盖面积计算所述地表层的最大滑动系数，获取所述地表层受到的拉伸应力数据，根据所述最大滑动系数和所述拉伸应力数据监测所述选定监测区域的形变位移。本申请具有提高当前灾害区域的位移检测准确性的效果。

Description

一种地质灾害位移监测方法、装置、设备以及存储介质
技术领域
本发明涉及地质灾害监测的技术领域，尤其是涉及一种地质灾害位移监测方法、装置、设备以及存储介质。
背景技术
目前，地质灾害在综合了人为因素和自然因素的作用下，在对环境造成巨大破坏的同时也对人类的生产生活造成了严重的损失，如坍塌、滑坡、泥石流、沉降以及地裂缝等，因此，在发生或者预防地质灾害时，对地质灾害的位移进行监测的精确度也提出了更高的要求。
现有的地质灾害位移监测方法通常为在可能发生地质灾害的区域埋设位移检测装置，在土方出现位移时，第一时间感知滑方位移从而快速得到对应的位移检测数据，根据位移大小预测当前的地质灾害的风险性，进而进行对应的预警处理，以便于降低地质灾害的土方位移滑动对人类生产生活的影响，但是，同一灾害区域的不同位置的土方滑移对位移检测装置所造成的滑移影响不同，影响对当前灾害区域的位移检测准确性。
发明内容
为了提高当前灾害区域的位移检测准确性，本申请提供一种地质灾害位移监测方法、装置、设备以及存储介质。
本申请的上述发明目的一是通过以下技术方案得以实现的：
一种地质灾害位移监测方法，包括：
获取选定监测区域的地质条件参数和所述选定监测区域的地表层对应的地表层数据；
根据所述地质条件参数和所述地表层数据，计算对应的地表层在当前地质条件下的滑动形变系数，以便于根据所述滑动形变系数分析所述地表层的滑动风险；
获取所述地表层的倾斜角数据和对应的植被覆盖面积，根据所述倾斜角数据和所述植被覆盖面积计算所述地表层达到位移滑动阈值的的最大滑动系数；
获取所述地表层受到的拉伸应力数据，根据所述最大滑动系数和所述拉伸应力数据监测所述选定监测区域的形变位移。
通过采用上述技术方案，在当前地质条件下选定发生灾害风险最大的的监测区域，获取选定监测区域的地址条件参数如构造层发育情况、地下水位、岩石类型和沟道形态等，对选定监测区域进行综合分析，并根据选定监测区域的地表层对应的地表层数据，对地表层是否发生滑动位移进行针对性监测，通过对当前地质条件下的地表层的滑动形变系数的计算，便于分析地表层达到地质灾害所需要的大幅度滑动风险的形变程度，并通过地表层的倾斜角数据和对应的植被覆盖面积，对地表层的当前土壤的吸水性以及与相邻土层之间的附着性进行多维度分析，并计算对应地表层的最大滑动系数，有助于对地表层趋近于塌方滑动的滑动阈值进行准确计算，并通过土层之间的相互挤压变动，获取地表层受到的拉伸应力数据，根据拉伸应力数据对地表层造成的滑动影响是否达到最大滑动系数的阈值范畴，来准确监测对应的选定监测区域是否发生形变位移，达到对地表层的形变位移进行预先监测的效果，从而提高对当前灾害区域的位移检测准确性和预防及时性。
本申请在一较佳示例中可以进一步配置为：所述获取所述地表层受到的拉伸应力数据，根据所述最大滑动系数和所述拉伸应力数据监测所述选定监测区域的形变位移，具体包括：
获取所述选定监测区域的稳定层发生地面塌陷时的重力下沉数据；
根据所述重力下沉数据，计算与所述稳定层相邻的形变层所受到的拉伸应力作为所述地表层对应的拉伸应力数据；
根据所述最大滑动系数与所述拉伸应力数据，计算所述稳定层坍塌的土壤挤压对所述地表层造成的滑动位移；
获取所述滑动位移对残坡积层造成的积层挤压重力，当所述积层挤压重力达到预设的挤压受力饱和状态时，生成对所述选定监测区域进行防控的灾害监测指令。
通过采用上述技术方案，在选定监测区域的稳定层发生地面塌陷时，根据地面下层空间分析稳定层的重力下沉所引起的重力势能，从而得到对应的重力下沉数据，有助于对重力下沉数据对相邻斜坡的形变层造成的拉伸应力变化进行受力变化监测，并根据重力下沉数据计算相邻的形变层所受到的拉伸应力，对重力下沉数据给不同位置的形变层造成的位移形变风险进行实时监测，从而得到地表层在地面塌陷时所受到的拉伸应力数据，有助于对地表层的受力变化情况进行及时监测，根据拉伸应力变化及时地计算对应地表层的形变程度，根据最大滑动系数和拉伸应力数据，分析稳定层地面坍塌对地表层造成的滑动位移，包括横向水平位移和改变倾斜角的竖向位移，对地表层的位移情况进行准确分析，并获取地表层发生滑动位移对残坡积层造成的积层挤压重力，有助于对残坡积层的塌方危险进行及时的监测，并在积层挤压重力达到预设的挤压受力饱和状态时，即残坡积层趋近于塌方时，能够通过灾害监测指令对残坡积层进行及时的防控，减少残坡积层发生大幅度滑动位移对生产生活造成的灾害影响。
本申请在一较佳示例中可以进一步配置为：在所述根据所述重力下沉数据，计算与所述稳定层相邻的形变层所受到的拉伸应力，得到所述地表层对应的拉伸应力数据之后，所述方法还包括：
获取所述形变层中每个土壤层与所述地表层之间的距离关系，以及每个土壤层对应的土壤含水率；
根据所述土壤含水率计算相邻土壤层之间的含水率差值，结合所述距离关系生成所述形变层的梯度土壤含水率；
根据所述梯度土壤含水率和对应土壤含水率的每个土壤层受到的拉伸应力，计算相邻土壤层之间的土壤层滑动位移。
通过采用上述技术方案，通过土壤层厚度测量工具测量形变层的每个土壤层的厚度，从而得到每个土壤层与地表层之间的距离关系，并通过地下水位测量工具获取每个土壤层的土壤含水率，有助于对不同土壤含水率的土壤层之间的相互作用的拉伸应力进行准确测量，并通过对相邻土壤层之间的含水率差值进行计算，通过含水率差值计算对应的拉伸应力差，并按照土壤层与地表层之间的距离关系，来构建形变层的梯度土壤含水率，有助于对形变层进行层级形变分析，将形变位移变化监测细化至每一个土壤层，并根据梯度土壤含水率和对应的拉伸应力，计算相邻土壤层之间的滑动位移，有助于对每个相邻土壤层之间的相互位移变化进行细分监测，提高形变层滑动位移的监测细致化程度。
本申请在一较佳示例中可以进一步配置为：所述根据所述最大滑动系数与所述拉伸应力数据，计算所述稳定层地面坍塌对所述地表层造成的滑动位移，具体包括：
获取所述地表层在当前拉伸应力作用下产生的形变倾角参数；
根据所述形变倾角参数对所述最大滑动系数进行形变纠正处理，得到与当前拉伸应力作用相对应的纠偏滑动系数；
计算所述纠偏滑动系数与所述拉伸应力数据的乘积，得到所述地表层在当前滑动纠偏作用下的综合受力参数；
当所述综合受力参数与所述地表层自身重力所引起的滑动阻力相适配时，计算所述综合受力参数对应的综合受力状态下对所述地表层造成的滑动位移。
通过采用上述技术方案，在当前拉伸应力作用下获取地表层的倾斜角变化情况，从而得到对应的形变倾角参数，有助于对地表层的倾斜程度进行实时监测，并根据形变倾角参数对最大滑动系数进行形变纠正，有助于通过在地表层施加形变干预措施对当前拉伸应力进行抵抗，从而得到对应的纠偏滑动系数，有助于根据纠偏滑动系数设置干预措施的反向拉伸应力值，从而降低当前拉伸应力对地表层所造成的形变影响，并通过纠偏滑动系数与拉伸应力数据的乘积，计算反向干预拉力的阻力值，从而得到形变层的综合受力参数，有助于通过对形变层的综合受力分析来判断当前形变层是否发生形变位移，提高形变层的形变位移监测准确性，从而在形变层的综合受力参数与当前自身重力相适配时，即说明形变层趋近于塌方的饱和状态，风险系数较大，通过综合受力参数，计算当前综合受力情况下的形变层可产生的滑动位移，有助于对形变层可产生的滑动位移，提前在对应的滑动位移范围内部署相应的防护措施，有助于降低滑动位移所带来的地质灾害风险。
本申请在一较佳示例中可以进一步配置为：所述获取所述滑动位移对残坡积层造成的积层挤压重力，当所述积层挤压重力达到预设的饱和状态时，生成对所述选定监测区域进行防控的灾害监测指令，还包括：
获取所述残坡积层的积层自身重力和对应的积层倾斜角度，根据所述积层自身重力和所述积层倾斜角度计算所述残坡积层的残坡拉伸应力；
根据所述残坡拉伸应力，分析当前土壤含水率的所述残坡积层达到含水率饱和状态所需要的吸水时间；
获取所述残坡积层处的当前降雨量，根据所述当前降雨量和所述吸水时间综合分析所述残坡积层的积层应力变化，得到积层应力分析结果；
根据所述积层应力分析结果对所述残坡积层进行偏距校正处理，根据校正结果输出所述残坡积层的分层级防护处理信息。
通过采用上述技术方案，通过残坡积层的积层自身重力和对应的积层倾斜角度等参数的获取，对残坡积层的当前状态进行实时监测，并进一步计算对应的残坡积层的残坡拉伸应力，来对当前残坡积层的当前受力情况进行分析，并通过残坡拉伸应力来分析当前土壤含水率下，残坡积层达到饱和状态所需要的吸水时间，有助于在当前吸水时间范围内及时地对残坡积层的坍塌风险进行提前防御，并根据当前降雨量和对应的吸水时间，分析残坡积层在持续性吸水情况下的积层应力变化，从而得到残坡积层的积层应力分析结果，有助于对残坡积层的应力变化情况进行实时监测，在残坡积层存在坍塌风险时能够及时进行防控，并通过对残坡积层进行偏距校正处理，输出残坡积层的分层级防护处理信息便于对残坡积层进行分层级防护，达到针对性防护的目的，提高不同位置的残坡积层的防护性能。
本申请在一较佳示例中可以进一步配置为：所述获取所述地表层的倾斜角数据和对应的植被覆盖面积，根据所述倾斜角数据和所述植被覆盖面积计算所述地表层的最大滑动系数，具体包括：
获取所述地表层的倾斜角数据和对应的植被覆盖面积，根据所述倾斜角数据和对应的所述植被覆盖面积，计算所述地表层的土壤吸水系数；
根据所述当前吸水系数，分析所述地表层的当前储水量对相邻形变层造成的孔隙水压力；当所述孔隙水压力大于形变层自身重力时，结合所述倾斜角数据，计算所述孔隙水压力与所述形变层自身重力之间的滑动形变合力；
将所述滑动形变合力与地表层重力之间的比值作为当前储水量下的地表层的最大滑动系数。
通过采用上述技术方案，根据地表层的倾斜角数据和对应的植被覆盖面积，来分析当前植被覆盖情况对当前倾斜角度下的地表层的土壤保护力度，并进一步地计算地表层的土壤吸水系数，有助于对当前倾斜角度和植被覆盖情况下的地表层水土流失平衡情况进行分析，并根据当前吸水系数，分析地表层的当前储水量对相邻形变层造成的孔隙水压力，提高相邻形变层的受力分析准确性，当孔隙水压力大于形变层自身重力时，说明形变层的滑动风险增大，并结合倾斜角数据，对形变层所受到的孔隙水压力与形变层自身重力之间的滑动形变合理进行计算，有助于对当前形变条件下的形变层进行准确的受力分析，提高形变层的形变位移风险检测准确性，并计算滑动形变合理与地表层重力之间的比值，将计算得到的比值作为当前储水量下的地表层的最大滑动系数，有助于根据最大滑动系数分析当前情况下的地表层与形变层之间的滑动风险，提高当前储水量下的地表层的滑动位移监测及时性。
本申请在一较佳示例中可以进一步配置为：所述根据所述地质条件参数和所述地表层数据，计算对应的地表层在当前地质条件下的滑动形变系数，还包括：
根据所述地质条件参数和所述地表层数据，分析对应的选定监测区域的地表层的地表形变，得到地表形变分析结果；
获取所述地表层的地表形变时间，根据所述地表形变时间和对应的地表形变分析结果，计算所述地表层的形变位移；
计算所述形变位移对所述选定监测区域造成的形变系数，根据所述形变系数预测对应形变位置发生断层的风险；
根据预测得到的断层风险，对所述形变位置进行构造层发育分析，并输出与构造层发育程度相对应的形变防控预警指令。
通过采用上述技术方案，根据地质条件参数和地表层数据，对选定监测区域中的地表层的地表形变进行分析，从而得到当前岩石和植被覆盖情况下的地表形变分析结果，有助于提高地表形变分析结果与地表层的实际地质情况的适配性，并通过地表形变时间的获取，结合对应的地表形变分析结果，对地表层的形变位移进行准确分析，从而有助于对当前地表层的地表位移趋势进行准确预测，并通过形变位移对选定监测区域造成的形变系数进行计算，根据计算得到的形变系数来预测对应的形变位置发生断层的风险，有助于对当前形变位置的灾害风险进行及时的预判，提高形变防护的及时性，并根据对应的断层风险预测结果，对形变位置进行构造层发育分析，对当前构造层的发育程度进行抗灾害预测，从而得到与构造层发育程度相对应的形变防控预警指令，有助于对当前构造层进行及时的风险防控，降低形变灾害的危害性。
本申请的上述发明目的二是通过以下技术方案得以实现的：
一种地质灾害位移监测装置，包括：
获取选定监测区域的地质条件参数和所述选定监测区域的地表层对应的地表层数据；
根据所述地质条件参数和所述地表层数据，计算对应的地表层在当前地质条件下的滑动形变系数；
获取所述地表层的倾斜角数据和对应的植被覆盖面积，根据所述倾斜角数据和所述植被覆盖面积计算所述地表层的最大滑动系数；
获取所述地表层受到的拉伸应力数据，根据所述最大滑动系数和所述拉伸应力数据监测所述选定监测区域的形变位移。
数据获取模块，用于获取选定监测区域的地质条件参数和所述选定监测区域的地表层对应的地表层数据；
形变计算模块，用于根据所述地质条件参数和所述地表层数据，计算对应的地表层在当前地质条件下的滑动形变系数，以便于根据所述滑动形变系数分析所述地表层的滑动风险；
滑动计算模块，用于获取所述地表层的倾斜角数据和对应的植被覆盖面积，根据所述倾斜角数据和所述植被覆盖面积计算所述地表层达到位移滑动阈值的的最大滑动系数；
位移监测模块，用于获取所述地表层受到的拉伸应力数据，根据所述最大滑动系数和所述拉伸应力数据监测所述选定监测区域的形变位移。
通过采用上述技术方案，在当前地质条件下选定发生灾害风险最大的的监测区域，获取选定监测区域的地址条件参数如构造层发育情况、地下水位、岩石类型和沟道形态等，对选定监测区域进行综合分析，并根据选定监测区域的地表层对应的地表层数据，对地表层是否发生滑动位移进行针对性监测，通过对当前地质条件下的地表层的滑动形变系数的计算，便于分析地表层达到地质灾害所需要的大幅度滑动风险的形变程度，并通过地表层的倾斜角数据和对应的植被覆盖面积，对地表层的当前土壤的吸水性以及与相邻土层之间的附着性进行多维度分析，并计算对应地表层的最大滑动系数，有助于对地表层趋近于塌方滑动的滑动阈值进行准确计算，并通过土层之间的相互挤压变动，获取地表层受到的拉伸应力数据，根据拉伸应力数据对地表层造成的滑动影响是否达到最大滑动系数的阈值范畴，来准确监测对应的选定监测区域是否发生形变位移，达到对地表层的形变位移进行预先监测的效果，从而提高对当前灾害区域的位移检测准确性和预防及时性。
本申请的上述目的三是通过以下技术方案得以实现的：
一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述地质灾害位移监测方法的步骤。
本申请的上述目的四是通过以下技术方案得以实现的：
一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述地质灾害位移监测方法的步骤。
综上所述，本申请包括以下至少一种有益技术效果：
1、在当前地质条件下选定发生灾害风险最大的的监测区域，获取选定监测区域的地址条件参数如构造层发育情况、地下水位、岩石类型和沟道形态等，对选定监测区域进行综合分析，并根据选定监测区域的地表层对应的地表层数据，对地表层是否发生滑动位移进行针对性监测，通过对当前地质条件下的地表层的滑动形变系数的计算，便于分析地表层达到地质灾害所需要的大幅度滑动风险的形变程度，并通过地表层的倾斜角数据和对应的植被覆盖面积，对地表层的当前土壤的吸水性以及与相邻土层之间的附着性进行多维度分析，并计算对应地表层的最大滑动系数，有助于对地表层趋近于塌方滑动的滑动阈值进行准确计算，并通过土层之间的相互挤压变动，获取地表层受到的拉伸应力数据，根据拉伸应力数据对地表层造成的滑动影响是否达到最大滑动系数的阈值范畴，来准确监测对应的选定监测区域是否发生形变位移，达到对地表层的形变位移进行预先监测的效果，从而提高对当前灾害区域的位移检测准确性和预防及时性；
2、在选定监测区域的稳定层发生地面塌陷时，根据地面下层空间分析稳定层的重力下沉所引起的重力势能，从而得到对应的重力下沉数据，有助于对重力下沉数据对相邻斜坡的形变层造成的拉伸应力变化进行受力变化监测，并根据重力下沉数据计算相邻的形变层所受到的拉伸应力，对重力下沉数据给不同位置的形变层造成的位移形变风险进行实时监测，从而得到地表层在地面塌陷时所受到的拉伸应力数据，有助于对地表层的受力变化情况进行及时监测，根据拉伸应力变化及时地计算对应地表层的形变程度，根据最大滑动系数和拉伸应力数据，分析稳定层地面坍塌对地表层造成的滑动位移，包括横向水平位移和改变倾斜角的竖向位移，对地表层的位移情况进行准确分析，并获取地表层发生滑动位移对残坡积层造成的积层挤压重力，有助于对残坡积层的塌方危险进行及时的监测，并在积层挤压重力达到预设的饱和状态时，即残坡积层趋近于塌方时，能够通过灾害监测指令对残坡积层进行及时的防控，减少残坡积层发生大幅度滑动位移对生产生活造成的灾害影响；
3、通过土壤层厚度测量工具测量形变层的每个土壤层的厚度，从而得到每个土壤层与地表层之间的距离关系，并通过地下水位测量工具获取每个土壤层的土壤含水率，有助于对不同土壤含水率的土壤层之间的相互作用的拉伸应力进行准确测量，并通过对相邻土壤层之间的含水率差值进行计算，通过含水率差值计算对应的拉伸应力差，并按照土壤层与地表层之间的距离关系，来构建形变层的梯度土壤含水率，有助于对形变层进行层级形变分析，将形变位移变化监测细化至每一个土壤层，并根据梯度土壤含水率和对应的拉伸应力，计算相邻土壤层之间的滑动位移，有助于对每个相邻土壤层之间的相互位移变化进行细分监测，提高形变层滑动位移的监测细致化程度。
附图说明
图1是本申请一实施例一种地质灾害位移监测方法的实现流程图。
图2是本申请一实施例一种地质灾害位移监测方法进行形变位移防控的实现流程图。
图3是本申请一实施例一种地质灾害位移监测方法步骤S30的实现流程图。
图4是本申请一实施例一种地质灾害位移监测方法步骤S40的实现流程图。
图5是本申请一实施例一种地质灾害位移监测方法对土壤层进行监测的实现流程图。
图6是本申请一实施例一种地质灾害位移监测方法步骤S403的实现流程图。
图7是本申请一实施例一种地质灾害位移监测方法进行积层分级防护的实现流程图。
图8是本申请一实施例一种地质灾害位移监测装置的结构框图。
图9是实现一种地质灾害位移监测方法的计算机设备的内部结构示意图。
图10是本实施例监测区域的地质层级结构示意图。
图中，1、稳定层，2、残坡积层，3、地表层，4、形变层，41、第一土壤层，42、土壤层2，43、土壤层3。
具体实施方式
以下结合附图对本申请作进一步详细说明。
本实施例中，如图10所示，稳定层1为监测斜坡中地形相对平稳的坡顶位置，残坡积层2为日常风化或者滚落的岩石、土壤等积累在地表层底部位置形成的岩石层，地表层3为监测斜坡的最外层土壤，通常覆盖有草地或者植被，形变层4为地表层与地下水之间的内部地质层，通常由岩石、土壤材质的差异划分为多个土壤层，如图中的第一土壤层41、第二土壤层42、第三土壤层43，多个相邻土壤层共同形成形变层4，且竖直方向上，土壤含水率按照坡顶至地下水地层的顺序由高到低依次升高，水平方向上，按照地质层与中心地下水地层的距离由外到内依次升高。在稳定层1的坡顶位置出现地层塌陷时，塌陷的土壤对相邻的形变层4产生挤压，推动形变层4、地表层3和残坡积层2等发生水平方向上的位移。
在一实施例中，如图1所示，本申请公开了一种地质灾害位移监测方法，具体包括如下步骤：
S10：获取选定监测区域的地质条件参数和选定监测区域的地表层对应的地表层数据。
具体的，通过选定监测区域的遥感数据获取对应的地质条件参数包括区域外形特征、植被覆盖情况和区域沟道形态等，通过预设的图像采集装置获取坡面裸露松散物质、岩石类型和人为活动范畴等地表层数据。
S20：根据地质条件参数和地表层数据，计算对应的地表层在当前地质条件下的滑动形变系数，以便于根据滑动形变系数分析地表层的滑动风险。
具体的，根据地址条件参数和地表层数据，假设外部条件不发生改变，在当前岩石类型以及地表层的土壤条件下，通过对地表层不同位置进行土壤采样，并计算多个土壤样品在当前含水率下的土壤粘力平均值，根据土壤粘力平均值与当前地表层的自身重力的比值作为当前地质条件下的滑动形变系数，其中当前地表层的自身重力通过预埋于地表层中的压力计采集得到。
在一实施例中，为了对地表层的形变位移进行提前防控，如图2所示，步骤S20还包括：
S201：根据地质条件参数和地表层数据，分析对应的选定监测区域的地表层的地表形变，得到地表形变分析结果。
具体的，根据地质条件参数和当前地势条件下的地表层数据，对地表层的地表位移或者地表倾斜角等进行综合分析，从而根据地表位移或地表倾斜角变动情况，分析地表层的地表形变情况，如当地表层存在地表位移时说明地表层存在水平方向的形变，当地表层的地表倾斜角变动，则说明地表层存在竖直方向的形变，当地表层同时存在地表位移和地表倾斜角变动时，则说明地表层同事存在水平方向和竖直方向上的形变，从而得到地表分析结果。
S202：获取地表层的地表形变时间，根据地表形变时间和对应的地表形变分析结果，计算地表层的形变位移。
具体的，当地表层的内部土壤发生流动或挤压时开始计算地表层的地表形变时间，当检测到内部土壤处于相对静止时，将统计得到的监测时间作为地表形变时间，地表形变时间与地表位移之间的乘积作为水平方向上的形变位移，将地表形变时间、地表倾斜角变动角度与地表沉总长度之间的乘积作为竖直方向上的形变位移，从而将水平方向与竖直方向上的形变位移之和作为地表层的总形变位移。
S203：计算形变位移对选定监测区域造成的形变系数，根据形变系数预测对应形变位置发生断层的风险。
具体的，当地表层的内部土壤受到由内而外的挤压应力时，地表层开始发生形变，当形变达到地表层的最大形变阈值时，地表层开始发生形变位移，通过预埋于地表层的压力计采集发生形变位移期间的地表层内部压力变化值，将内部压力变化值与未发生位移前的原始压力值之间的比值作为形变位移对应的形变系数，当形变系数超过预设的断层阈值如原始压力值的1/2时，说明当前形变位置存在有断层风险；或者计算当前形变位移所造成的倾斜角度变化值，将倾斜角度变化值与未发生位移前的原始倾斜角度之间的比值作为形变系数，当形变系数超过或者趋近于预设的倾斜阈值时，则说明当前形变位置存在的断层风险越高。
S204：根据预测得到的断层风险，对形变位置进行构造层发育分析，并生成与构造层发育程度相对应的形变防控预警指令。
具体的，根据形变位置预测得到的断层风险，对形变位置的构造层发育程度进行分析，包括地表层的沉积分布情况、岩石分布对地表土壤的拉伸应力的影响和地表植被覆盖程度等，根据预设的断层风险分析地表土壤中在当前岩石分布、植被覆盖情况下的不同位置的形变程度，从而得到构造层与当前形变相对应的土壤发育程度，并根据对应的形变位置生成对应的形变防控预警指令，对地表层中形变程度不同的多个位置进行针对性防控。
S30：获取地表层的倾斜角数据和对应的植被覆盖面积，根据倾斜角数据和植被覆盖面积计算地表层达到位移滑动阈值的的最大滑动系数。
具体的，如图3所示，步骤S30具体包括以下步骤：
S301：获取地表层的倾斜角数据和对应的植被覆盖面积，根据倾斜角数据和对应的植被覆盖面积，计算地表层的土壤吸水系数。
具体的，通过在地表层的表面位置布设多个水准仪进行倾斜角测量，通过与形变位置相邻水准仪之间的倾斜角度得到对应形变位置的倾斜角数据，并通过选定监测区域的遥感数据分析当前地表层的植被覆盖情况，从而得到对应的植被覆盖面积，获取当前倾斜角下的植被覆盖区域的水流速度，当前水流速度乘以对应的降水时间得到总流失水量，通过公共数据获取选定检测区域的当前总降雨量，将当前总降雨量与总流失水量的差值作为当前植被的土壤吸水量，土壤吸水量与当前总降雨量之间的比值作为地表层的土壤吸水系数。
S302：根据当前吸水系数，分析地表层的当前储水量对相邻形变层造成的孔隙水压力。
具体的，根据当前吸水系数与总地表层土壤的孔隙水压力变化，获取当前吸水状态下的地表层对相邻形变层造成的压力与土层之间的滑动摩擦力变化情况，在当前倾斜角下计算地表层压力与滑动摩擦力之间的地表合力，并结合相邻形变层的渗压计采集到的渗压变化值，将地表合力与渗压变化值之间的压力差作为形变层对应的孔隙水压力，根据孔隙水压力判断地表层的吸水状态是否导致地表层发生滑动位移。
S303：当孔隙水压力大于形变层自身重力时，结合倾斜角数据，计算孔隙水压力与形变层自身重力之间的滑动形变合力。
具体的，通过预埋的土压力计获取每个土层的自身压力，从而得到形变层自身重力，将孔隙水压力与形变层自身重力进行求差运算，当孔隙水压力大于形变层自身重力时，说明此时的地表层存在的滑移风险较大，则结合当前倾斜角度，计算孔隙水压力与形变层自身重力之间的滑动形变合力，计算过程如公式(1)所示：
F合＝G形-F孔cosθ2 (1)
其中，F合表示滑动形变合力，G形表示形变层自身重力，F孔表示孔隙水压力，θ2表示形变层当前倾斜角度。
S304：计算滑动形变合力与地表层重力之间的比值，并结合倾斜角数据，计算当前储水量下的地表层的最大滑动系数。
具体的，通过预埋于地表层与形变层相邻位置的土压力计采集当前位置的地表层重力，计算滑动形变合力与地表层重力之间的比值，并结合地表层的当前倾斜角度，计算当前储水量下的地表层的最大滑动系数，其中，地表层的当前储水量对地表层重力造成的变化可通过植被覆盖面积与选定监测区域的总面积之间的面积占比进行表示，其中地表层的最大滑动系数通过公式(2)计算得到，公式(2)如下所示：
其中，μ表示地表层的当前最大滑动系数，F合表示滑动形变合力，G表表示地表层重力，S植表示植被覆盖面积，S区表示选定监测区域的总面积，θ1表示地表层当前倾斜角，本实施例中的θ1的取值范围为40°-75°，需要说明的是可以根据实际地址条件需要选择对应的倾斜角取值进行计算，不局限于本实施例中的一种。
S40：获取地表层受到的拉伸应力数据，根据最大滑动系数和拉伸应力数据监测选定监测区域的形变位移。
具体的，如图4所示，步骤S40具体包括以下步骤：
S401：获取选定监测区域的稳定层发生地面塌陷时的重力下沉数据。
具体的，通过预设于选定监测区域的稳定层的水平监测设备获取稳定层的地面塌陷情况，在稳定层发生地面塌陷时通过水平监测设备获取地面塌陷空间范围，通过对塌陷区域的土质进行采样判断塌方的土质密度，将塌方空间范围计算对应的塌方体积，将塌方体积与土质密度之间乘积作为塌方对应的重力下沉数据。
S402：根据重力下沉数据，计算与稳定层相邻的形变层所受到的拉伸应力作为地表层对应的拉伸应力数据。
具体的，根据倾斜角数据获取稳定层与相邻形变层之间的夹角，将重力下沉数据与当前夹角的余弦值之间的乘积作为形变层所受到的拉伸应力，根据力的传递关系，假设形变层的土壤倾斜情况不发生改变，则地表层对应的拉伸应力数据与形变层受到的拉伸应力相同，拉伸应力是相邻地质层之间的相互抵抗阻力，拉伸应力形变阈值是在相邻地质层的影响下，监测地质层开始发生位移滑动时的阻力值，拉伸应力是地质层抵抗其他相邻地质层的滑动作用于地质层，且在地质层达到滑动位移阈值之前的受力情况。
在一实施例中，为了根据土壤含水量的不同对形变层的每个土壤层进行位移精准监测，如图5所示，步骤S402还包括：
S4021：获取形变层中每个土壤层与地表层之间的距离关系，以及每个土壤层对应的土壤含水率。
具体的，通过预埋于选定监测区域的多点式位移计获取形变层的每个土壤层与地表层之间的垂直距离，从而得到对应的距离关系，并通过对应的渗压计采集每个土壤层的含水量，从而得到对应的土壤含水率。
S4022：根据土壤含水率计算相邻土壤层之间的含水率差值，结合距离关系生成形变层的梯度土壤含水率。
具体的，含水率不同的土壤所对应的滑动阻力不同，如含水率越低的地层对应的滑动阻力越大，含水率越高的地层对应的滑动阻力越小，当含水率达到100％时即泥沙之间的阻力趋近于0随着水流被带走，从而形成泥石流；通过相邻土壤层之间的土壤含水率的差值计算得到相邻土壤层对应的含水率差值，并按照每个土壤层与地表层之间的距离关系，依次对含水率差值进行排序，从而得到形变层的梯度土壤含水率。
S4023：根据梯度土壤含水率和对应土壤含水率的每个土壤层受到的拉伸应力，计算相邻土壤层之间的土壤层滑动位移。
具体的，根据梯度土壤含水率，计算相邻含水率差值与当前土壤层的含水率之间的比值，将计算得到的比值与对应的土壤层受到的拉伸应力的乘积作为相邻土壤层之间的滑动拉力，假设滑动过程中每个土壤层的自身重力不发生改变，则计算滑动过程中造成的土壤层倾角变化值对应的余弦值，将滑动拉力与对应倾角余弦值之间的乘积与当前土壤层的检测位置的重力变化值之间的比值作为滑移加速度，当土壤层发生滑移时统计对应的滑移时间，从而根据滑移加速度和对应的滑移时间计算得到对应的滑移位移。
S403：根据最大滑动系数与拉伸应力数据，计算稳定层坍塌的土壤挤压对地表层造成的滑动位移。
具体的，如图6所示，步骤S403具体包括以下步骤：
S4031：获取地表层在当前拉伸应力作用下产生的形变倾角参数。
具体的，通过预设于地表层的倾角检测装置获取地表层静止状态下的原始倾斜角度，在地表层受到其他土壤层的外力干扰时，通过倾角检测装置采集当前拉伸应力作用下的形变倾角，将形变倾角与原始倾斜角度之间的角度差作为地表层的形变倾角参数。
S4032：根据形变倾角参数对最大滑动系数进行形变纠正处理，得到与当前拉伸应力作用相对应的纠偏滑动系数。
具体的，根据形变倾角参数的余弦值，计算相邻观测点之间的水平距离在地表层观测线上的投影长度，将计算得到的投影长度值作为当前拉伸应力作用下的纠偏滑动系数。
S4033：计算纠偏滑动系数与拉伸应力数据的乘积，得到地表层在当前滑动纠偏作用下的综合受力参数。
具体的，将地表层倾斜角的余弦值与地表层所受到的拉伸应力数据之间的乘积作为地表层对形变层施加的外力，结合稳定层对形变层施加的积压外力与纠偏滑动系数之间的乘积作为反向积压拉力，将地表层施加的外力与反向积压拉力之间的合力作为形变层的综合受力参数。
S4034：当综合受力参数与地表层自身重力所引起的滑动阻力相适配时，计算综合受力参数对应的综合受力状态下对地表层造成的滑动位移。
具体的，通过预埋于形变层的土压力计获取形变层在采集位置处的自身重力，通过综合受力参数与形变层自身重力之间的合力差值，判断综合受力参数是否与形变层自身重力相适配，当综合受力参数趋近于形变层自身重力时，说明形变层的滑动位移结束，形变层的综合受力保持相对稳定状态，则通过合力差值与形变层的形变时间之间的乘积得到对应的滑动位移。
S404：获取滑动位移对残坡积层造成的积层挤压重力，当积层挤压重力达到预设的挤压受力饱和状态时，生成对选定监测区域进行防控的灾害监测指令。
具体的，获取当前滑动位移对残坡积层造成的倾斜角度变化值，假定地表层在发生滑动位移过程中不发生水土流失，则将地表层所受到的拉伸应力与倾斜角度变化值的余弦值之间的乘积作为残坡积层的积层挤压重力，当积层挤压重力达到预设的饱和状态时，说明当前倾斜状态下的残坡积层容易发生土层滑落的风险，则根据当前倾斜角度设置对应的防控灾害监测指令，通过灾害监测指令对残坡积层进行定点防护，其中，本实施例中的积层饱和状态设置为积层的倾斜角度阈值，本实施例中的倾斜角度阈值设置为75°，可以根据实际需要进行设置，不局限于本实施例中的一种。
在一实施例中，为了根据实际吸水情况对残坡积层进行分级防护，如图7所示，步骤S404还包括：
S4041：获取残坡积层的积层自身重力和对应的积层倾斜角度，根据积层自身重力和积层倾斜角度计算残坡积层的残坡拉伸应力。
具体的，通过预埋于残坡积层的土压力计对残坡积层的自身重力进行数据采集，并通过预设于残坡积层上的多个水平测量仪器获取相邻测量位置之间的积层倾斜角度，根据积层自身重力与积层倾斜角度的余弦值之间的乘积得到对应的残坡拉伸应力。
S4042：根据残坡拉伸应力，分析当前土壤含水率的残坡积层达到含水率饱和状态所需要的吸水时间。
具体的，获取单位时间内的土壤含水率的变化值，计算对应含水率变化值造成的残坡拉伸应力变化值，当残坡拉伸应力达到与残坡自身重力相持平时，计算当前残坡拉伸应力变化值达到应力持平所需要的总时间，从而得到残坡积层达到饱和状态所需要的吸水时间。
S4043：获取残坡积层处的当前降雨量，根据当前降雨量和吸水时间综合分析残坡积层的积层应力变化，得到积层应力分析结果。
具体的，通过选定监测区域的降水量统计得到对应残坡积层的当前降水量，根据吸水时间和当前降水量的乘积得到残坡积层的总吸水量，并根据残坡积层的当前倾斜角度计算残坡积层的总流失水量，从而得到残坡积层的实际吸水量，根据实际吸水量对残坡积层自身重力的影响情况，分析残坡积层的积层应力变化，从而得到积层应力分析结果。
S4044：根据积层应力分析结果对残坡积层进行偏距校正处理，根据校正结果输出残坡积层的分层级防护处理信息。
具体的，根据积层应力分析结果，通过预埋的水压力计获取实际吸水量下的积层压力变化值，根据残坡积层上相邻监测点的水平距离在积层观测线上的投影长度，将积层压力变化值与当前积层投影长度之间的乘积作为偏距校正结果，根据偏距校正结果，对相邻观测点之间的残坡积层进行分层级防护，如设置对应拉伸应力的防护围栏或者沙坝等。
应理解，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。
在一实施例中，提供一种地质灾害位移监测装置，该地质灾害位移监测装置与上述实施例中地质灾害位移监测方法一一对应。如图8所示，该地质灾害位移监测装置包括数据获取模块、形变计算模块、滑动计算模块和位移监测模块。各功能模块详细说明如下：数据获取模块，用于获取选定监测区域的地质条件参数和选定监测区域的地表层对应的地表层数据。
形变计算模块，用于根据地质条件参数和地表层数据，计算对应的地表层在当前地质条件下的滑动形变系数，以便于根据滑动形变系数分析地表层的滑动风险。
滑动计算模块，用于获取地表层的倾斜角数据和对应的植被覆盖面积，根据倾斜角数据和植被覆盖面积计算地表层达到位移滑动阈值的的最大滑动系数。
位移监测模块，用于获取地表层受到的拉伸应力数据，根据最大滑动系数和拉伸应力数据监测选定监测区域的形变位移。
优选的，位移监测模块具体包括：
稳定数据获取子模块，用于获取选定监测区域的稳定层发生地面塌陷时的重力下沉数据。
应力计算子模块，用于根据重力下沉数据，计算与稳定层相邻的形变层所受到的拉伸应力作为地表层对应的拉伸应力数据。
位移计算子模块，用于根据最大滑动系数与拉伸应力数据，计算稳定层坍塌的土壤挤压对地表层造成的滑动位移。
积层挤压分析子模块，用于获取滑动位移对残坡积层造成的积层挤压重力，当积层挤压重力达到预设的挤压受力饱和状态时，生成对选定监测区域进行防控的灾害监测指令。
优选的，应力计算子模块还包括：
土壤数据获取单元，用于获取形变层中每个土壤层与地表层之间的距离关系，以及每个土壤层对应的土壤含水率。
含水率计算单元，用于根据土壤含水率计算相邻土壤层之间的含水率差值，结合距离关系生成形变层的梯度土壤含水率。
土壤位移计算单元，用于根据梯度土壤含水率和对应土壤含水率的每个土壤层受到的拉伸应力，计算相邻土壤层之间的土壤层滑动位移。
优选的，位移计算子模块具体包括：
倾角数据获取单元，用于获取地表层在当前拉伸应力作用下产生的形变倾角参数。
纠偏单元，用于根据形变倾角参数对最大滑动系数进行形变纠正处理，得到与当前拉伸应力作用相对应的纠偏滑动系数。
综合受力计算单元，用于计算纠偏滑动系数与拉伸应力数据的乘积，得到地表层在当前滑动纠偏作用下的综合受力参数。
综合位移单元，用于当综合受力参数与地表层自身重力所引起的滑动阻力相适配时，计算综合受力参数对应的综合受力状态下对地表层造成的滑动位移。
优选的，积层挤压分析子模块还包括：
残坡应力计算单元，用于获取残坡积层的积层自身重力和对应的积层倾斜角度，根据积层自身重力和积层倾斜角度计算残坡积层的残坡拉伸应力。
吸水时间计算单元，用于根据残坡拉伸应力，分析当前土壤含水率的残坡积层达到含水率饱和状态所需要的吸水时间。
应力分析单元，用于获取残坡积层处的当前降雨量，根据当前降雨量和吸水时间综合分析残坡积层的积层应力变化，得到积层应力分析结果。
分级防护单元，用于根据积层应力分析结果对残坡积层进行偏距校正处理，根据校正结果输出残坡积层的分层级防护处理信息。
优选的，滑动计算模块具体包括：
吸水系数计算子模块，用于获取地表层的倾斜角数据和对应的植被覆盖面积，根据倾斜角数据和对应的植被覆盖面积，计算地表层的土壤吸水系数。
水压力计算子模块，用于根据当前吸水系数，分析地表层的当前储水量对相邻形变层造成的孔隙水压力。
形变合力计算子模块，用于当孔隙水压力大于形变层自身重力时，结合倾斜角数据，计算孔隙水压力与形变层自身重力之间的滑动形变合力。
地表滑动系数计算子模块，用于将滑动形变合力与地表层重力之间的比值作为当前储水量下的地表层的最大滑动系数。
优选的，形变计算模块还包括：
地表形变分析子模块，用于根据地质条件参数和地表层数据，分析对应的选定监测区域的地表层的地表形变，得到地表形变分析结果。
地表形变位移计算子模块，用于获取地表层的地表形变时间，根据地表形变时间和对应的地表形变分析结果，计算地表层的形变位移。
风险分析子模块，用于计算形变位移对选定监测区域造成的形变系数，根据形变系数预测对应形变位置发生断层的风险。
形变防控子模块，用于根据预测得到的断层风险，对形变位置进行构造层发育分析，并输出与构造层发育程度相对应的形变防控预警指令。
关于地质灾害位移监测装置的具体限定可以参见上文中对于地质灾害位移监测方法的限定，在此不再赘述。上述地质灾害位移监测装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是服务器，其内部结构图可以如图9所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的数据库用于存储选定监测区域的不同土壤层发生位移滑动所产生的数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种地质灾害位移监测方法。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现一种地质灾害位移监测方法的步骤。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，仅以上述各功能单元、模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能单元、模块完成，即将所述装置的内部结构划分成不同的功能单元或模块，以完成以上描述的全部或者部分功能。
以上所述实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围，均应包含在本申请的保护范围之内。

Claims (10)

1.一种地质灾害位移监测方法，其特征在于，包括：
获取选定监测区域的地质条件参数和所述选定监测区域的地表层对应的地表层数据；
根据所述地质条件参数和所述地表层数据，计算对应的地表层在当前地质条件下的滑动形变系数，以便于根据所述滑动形变系数分析所述地表层的滑动风险；
获取所述地表层的倾斜角数据和对应的植被覆盖面积，根据所述倾斜角数据和所述植被覆盖面积计算所述地表层达到位移滑动阈值的的最大滑动系数；
获取所述地表层受到的拉伸应力数据，根据所述最大滑动系数和所述拉伸应力数据监测所述选定监测区域的形变位移。
2.根据权利要求1所述的地质灾害位移监测方法，其特征在于，所述获取所述地表层受到的拉伸应力数据，根据所述最大滑动系数和所述拉伸应力数据监测所述选定监测区域的形变位移，具体包括：
获取所述选定监测区域的稳定层发生地面塌陷时的重力下沉数据；
根据所述重力下沉数据，计算与所述稳定层相邻的形变层所受到的拉伸应力作为所述地表层对应的拉伸应力数据；
根据所述最大滑动系数与所述拉伸应力数据，计算所述稳定层坍塌的土壤挤压对所述地表层造成的滑动位移；
获取所述滑动位移对残坡积层造成的积层挤压重力，当所述积层挤压重力达到预设的挤压受力饱和状态时，生成对所述选定监测区域进行防控的灾害监测指令。
3.根据权利要求2所述的地质灾害位移监测方法，其特征在于，在所述根据所述重力下沉数据，计算与所述稳定层相邻的形变层所受到的拉伸应力，得到所述地表层对应的拉伸应力数据之后，所述方法还包括：
获取所述形变层中每个土壤层与所述地表层之间的距离关系，以及每个土壤层对应的土壤含水率；
根据所述土壤含水率计算相邻土壤层之间的含水率差值，结合所述距离关系生成所述形变层的梯度土壤含水率；
根据所述梯度土壤含水率和对应土壤含水率的每个土壤层受到的拉伸应力，计算相邻土壤层之间的土壤层滑动位移，。
4.根据权利要求2所述的地质灾害位移监测方法，其特征在于，所述根据所述最大滑动系数与所述拉伸应力数据，计算所述稳定层地面坍塌对所述地表层造成的滑动位移，具体包括：
获取所述地表层在当前拉伸应力作用下产生的形变倾角参数；
根据所述形变倾角参数对所述最大滑动系数进行形变纠正处理，得到与当前拉伸应力作用相对应的纠偏滑动系数；
计算所述纠偏滑动系数与所述拉伸应力数据的乘积，得到所述地表层在当前滑动纠偏作用下的综合受力参数；
当所述综合受力参数与所述地表层自身重力所引起的滑动阻力相适配时，计算所述综合受力参数对应的综合受力状态下对所述地表层造成的滑动位移。
5.根据权利要求2所述的地质灾害位移监测方法，其特征在于，所述获取所述滑动位移对残坡积层造成的积层挤压重力，当所述积层挤压重力达到预设的饱和状态时，生成对所述选定监测区域进行防控的灾害监测指令，还包括：
获取所述残坡积层的积层自身重力和对应的积层倾斜角度，根据所述积层自身重力和所述积层倾斜角度计算所述残坡积层的残坡拉伸应力；
根据所述残坡拉伸应力，分析当前土壤含水率的所述残坡积层达到含水率饱和状态所需要的吸水时间；
获取所述残坡积层处的当前降雨量，根据所述当前降雨量和所述吸水时间综合分析所述残坡积层的积层应力变化，得到积层应力分析结果；
根据所述积层应力分析结果对所述残坡积层进行偏距校正处理，根据校正结果输出所述残坡积层的分层级防护处理信息。
6.根据权利要求1所述的地质灾害位移监测方法，其特征在于，所述获取所述地表层的倾斜角数据和对应的植被覆盖面积，根据所述倾斜角数据和所述植被覆盖面积计算所述地表层的最大滑动系数，具体包括：
获取所述地表层的倾斜角数据和对应的植被覆盖面积，根据所述倾斜角数据和对应的所述植被覆盖面积，计算所述地表层的土壤吸水系数；
根据所述当前吸水系数，分析所述地表层的当前储水量对相邻形变层造成的孔隙水压力；
当所述孔隙水压力大于形变层自身重力时，结合所述倾斜角数据，计算所述孔隙水压力与所述形变层自身重力之间的滑动形变合力；
将所述滑动形变合力与地表层重力之间的比值作为当前储水量下的地表层的最大滑动系数。
7.根据权利要求1所述的地质灾害位移监测方法，其特征在于，在所述根据所述地质条件参数和所述地表层数据，计算对应的地表层在当前地质条件下的滑动形变系数之后，还包括：
根据所述地质条件参数和所述地表层数据，分析对应的选定监测区域的地表层的地表形变，得到地表形变分析结果；
获取所述地表层的地表形变时间，根据所述地表形变时间和对应的地表形变分析结果，计算所述地表层的形变位移；
计算所述形变位移对所述选定监测区域造成的形变系数，根据所述形变系数预测对应形变位置发生断层的风险；
根据预测得到的断层风险，对所述形变位置进行构造层发育分析，并输出与所述构造层发育程度相对应的形变防控预警指令。
8.一种地质灾害位移监测装置，其特征在于，包括：
数据获取模块，用于获取选定监测区域的地质条件参数和所述选定监测区域的地表层对应的地表层数据；
形变计算模块，用于根据所述地质条件参数和所述地表层数据，计算对应的地表层在当前地质条件下的滑动形变系数，以便于根据所述滑动形变系数分析所述地表层的滑动风险；
滑动计算模块，用于获取所述地表层的倾斜角数据和对应的植被覆盖面积，根据所述倾斜角数据和所述植被覆盖面积计算所述地表层达到位移滑动阈值的的最大滑动系数；
位移监测模块，用于获取所述地表层受到的拉伸应力数据，根据所述最大滑动系数和所述拉伸应力数据监测所述选定监测区域的形变位移。
9.一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至7任一项所述地质灾害位移监测方法的步骤。
10.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至7任一项所述地质灾害位移监测方法的步骤。
CN202310912621.6A 2023-07-24 2023-07-24 一种地质灾害位移监测方法、装置、设备以及存储介质 Active CN116933535B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310912621.6A CN116933535B (zh) 2023-07-24 2023-07-24 一种地质灾害位移监测方法、装置、设备以及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310912621.6A CN116933535B (zh) 2023-07-24 2023-07-24 一种地质灾害位移监测方法、装置、设备以及存储介质

Publications (2)

Publication Number Publication Date
CN116933535A true CN116933535A (zh) 2023-10-24
CN116933535B CN116933535B (zh) 2024-03-19

Family

ID=88380212

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310912621.6A Active CN116933535B (zh) 2023-07-24 2023-07-24 一种地质灾害位移监测方法、装置、设备以及存储介质

Country Status (1)

Country Link
CN (1) CN116933535B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117764402A (zh) * 2024-02-22 2024-03-26 山东光合云谷大数据有限公司 一种地质灾害信息处理系统

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102013150A (zh) * 2010-09-28 2011-04-13 浙江工业大学 基于雨强、斜坡土壤含水量和形变量的地质灾害预测系统
JP5731700B1 (ja) * 2014-07-25 2015-06-10 エー・シー・エス株式会社 水位予測に基づいた土砂災害予測システム
CN104846848A (zh) * 2015-03-23 2015-08-19 安徽理工大学 一种控制地下层状矿产开采沉陷区地表水土流失的方法
CN109493569A (zh) * 2018-12-14 2019-03-19 深圳高速工程检测有限公司 滑坡预警方法、装置、计算机设备和存储介质
CN111160644A (zh) * 2019-12-27 2020-05-15 成都理工大学 一种基于地质灾害风险评估的铁路选线方法及其装置
CN112598881A (zh) * 2020-12-03 2021-04-02 中煤航测遥感集团有限公司 一种地质灾害监测方法、装置及计算机设备
CN113392500A (zh) * 2021-05-11 2021-09-14 青海省地质环境监测总站 地质灾害监测预警系统
CN113566894A (zh) * 2021-08-13 2021-10-29 长春工程学院 基于物联网的地质灾害监测方法、系统、终端及介质
US20210341647A1 (en) * 2020-04-30 2021-11-04 Institute Of Geochemistry, Chinese Academy Of Sciences Method for determining surface runoff yield in vegetation-covered area
CN114021487A (zh) * 2022-01-10 2022-02-08 西南交通大学 一种山坡崩塌的预警方法、装置、设备及可读存储介质
US20220392326A1 (en) * 2021-05-27 2022-12-08 Northwest University Multi-state triggered system and method for remote automatic monitoring and early warning of landslide mass
CN116030597A (zh) * 2022-10-31 2023-04-28 成都智慧企业发展研究院有限公司 地质灾害风险预警处理方法、系统及电子设备
CN116415481A (zh) * 2022-10-21 2023-07-11 杭州鲁尔物联科技有限公司 区域滑坡灾害风险预测方法、装置、计算机设备及存储介质

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102013150A (zh) * 2010-09-28 2011-04-13 浙江工业大学 基于雨强、斜坡土壤含水量和形变量的地质灾害预测系统
JP5731700B1 (ja) * 2014-07-25 2015-06-10 エー・シー・エス株式会社 水位予測に基づいた土砂災害予測システム
CN104846848A (zh) * 2015-03-23 2015-08-19 安徽理工大学 一种控制地下层状矿产开采沉陷区地表水土流失的方法
CN109493569A (zh) * 2018-12-14 2019-03-19 深圳高速工程检测有限公司 滑坡预警方法、装置、计算机设备和存储介质
CN111160644A (zh) * 2019-12-27 2020-05-15 成都理工大学 一种基于地质灾害风险评估的铁路选线方法及其装置
US20210341647A1 (en) * 2020-04-30 2021-11-04 Institute Of Geochemistry, Chinese Academy Of Sciences Method for determining surface runoff yield in vegetation-covered area
CN112598881A (zh) * 2020-12-03 2021-04-02 中煤航测遥感集团有限公司 一种地质灾害监测方法、装置及计算机设备
CN113392500A (zh) * 2021-05-11 2021-09-14 青海省地质环境监测总站 地质灾害监测预警系统
US20220392326A1 (en) * 2021-05-27 2022-12-08 Northwest University Multi-state triggered system and method for remote automatic monitoring and early warning of landslide mass
CN113566894A (zh) * 2021-08-13 2021-10-29 长春工程学院 基于物联网的地质灾害监测方法、系统、终端及介质
CN114021487A (zh) * 2022-01-10 2022-02-08 西南交通大学 一种山坡崩塌的预警方法、装置、设备及可读存储介质
CN116415481A (zh) * 2022-10-21 2023-07-11 杭州鲁尔物联科技有限公司 区域滑坡灾害风险预测方法、装置、计算机设备及存储介质
CN116030597A (zh) * 2022-10-31 2023-04-28 成都智慧企业发展研究院有限公司 地质灾害风险预警处理方法、系统及电子设备

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
倪忠云;何政伟;赵银兵;赵勇;南希;: "基于RS和GIS的丹巴县植被盖度与地质灾害关系研究", 物探化探计算技术, no. 02, 15 March 2011 (2011-03-15) *
邱宇: "乐昌市两江镇某滑坡自动监测系统实施与变形数据分析", 采矿技术, vol. 17, no. 3, 31 May 2017 (2017-05-31), pages 61 - 63 *
邱宇: "探析贝岭镇米贝村泥石流灾害成因与防治措施", 西部资源, no. 3, 15 June 2021 (2021-06-15), pages 120 - 122 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117764402A (zh) * 2024-02-22 2024-03-26 山东光合云谷大数据有限公司 一种地质灾害信息处理系统
CN117764402B (zh) * 2024-02-22 2024-05-14 山东光合云谷大数据有限公司 一种地质灾害信息处理系统

Also Published As

Publication number Publication date
CN116933535B (zh) 2024-03-19

Similar Documents

Publication Publication Date Title
Fan et al. Spatio-temporal evolution of mass wasting after the 2008 M w 7.9 Wenchuan earthquake revealed by a detailed multi-temporal inventory
Askarinejad et al. Precursors of instability in a natural slope due to rainfall: a full-scale experiment
Park et al. Assessment of rainfall-induced shallow landslide susceptibility using a GIS-based probabilistic approach
CN109493569B (zh) 滑坡预警方法、装置、计算机设备和存储介质
Maurer et al. Evaluation of the liquefaction potential index for assessing liquefaction hazard in Christchurch, New Zealand
Corsini et al. Field monitoring of the Corvara landslide (Dolomites, Italy) and its relevance for hazard assessment
Chai et al. Finite element simulation of an embankment on soft clay–Case study
CN116933535B (zh) 一种地质灾害位移监测方法、装置、设备以及存储介质
Ching-Chuan et al. Internal soil moisture and piezometric responses to rainfall-induced shallow slope failures
CN104318103A (zh) 一种滑坡灾害监测预警降雨阈值判定方法
CN106884442B (zh) 一种沟谷状高填方地基多场耦合监测系统的实施方法
Zhang et al. Study of a roof water inrush prediction model in shallow seam mining based on an analytic hierarchy process using a grey relational analysis method
Rahardjo et al. Sensing and monitoring for assessment of rainfall-induced slope failures in residual soil
Wei et al. The effectiveness of a drainage tunnel in increasing the rainfall threshold of a deep-seated landslide
Rashidian et al. Modification of the liquefaction potential index to consider the topography in Christchurch, New Zealand
Wu et al. Model test of deformation evolution and multi factor prediction of anchorage slope stability under rainfall condition
Kramer et al. The use of numerical analysis in the interpretation of liquefaction case histories
Wubalem Landslide inventory, susceptibility, hazard and risk mapping
Bao et al. Evolution of high-filling loess slope under long-term seasonal fluctuation of groundwater
Ali et al. Field behavior and numerical simulation of coastal bund on soft marine clay loaded to failure
Simeoni et al. Field measurements in river embankments: Validation and management with spatial database and webGIS
Elkortbawi Insights into the Liquefaction Hazards in Napier and Hastings Based on the Assessment of Data from the 1931 Hawke's Bay, New Zealand, Earthquake
Camera et al. A coupled distributed hydrological-stability analysis on a terraced slope of Valtellina (northern Italy)
Maurer Moving towards an improved liquefaction hazard framework: lessons resulting from the 2010-2011 Canterbury, New Zealand, earthquake sequence
Lester et al. A new method for finite element modelling of prefabricated vertical drains

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