CN114139467A - Method for calculating distribution of water loss compression amount of pressure-bearing water-containing loose layer region - Google Patents

Method for calculating distribution of water loss compression amount of pressure-bearing water-containing loose layer region Download PDF

Info

Publication number
CN114139467A
CN114139467A CN202111297430.0A CN202111297430A CN114139467A CN 114139467 A CN114139467 A CN 114139467A CN 202111297430 A CN202111297430 A CN 202111297430A CN 114139467 A CN114139467 A CN 114139467A
Authority
CN
China
Prior art keywords
hydrological
ith
observation hole
water
hydrological observation
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.)
Pending
Application number
CN202111297430.0A
Other languages
Chinese (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.)
Huaibei Mining Group Co Ltd
Original Assignee
Huaibei Mining Group 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 Huaibei Mining Group Co Ltd filed Critical Huaibei Mining Group Co Ltd
Priority to CN202111297430.0A priority Critical patent/CN114139467A/en
Publication of CN114139467A publication Critical patent/CN114139467A/en
Pending legal-status Critical Current

Links

Classifications

    • 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
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention discloses a method for calculating the distribution of water loss compression amount of a confined water-containing unconsolidated layer region, which relates to the field of stratum deformation monitoring and comprises the steps of constructing a plurality of hydrological observation holes from the ground surface downwards, acquiring relevant parameters of a target aquifer and an overlying rock-soil body thereof through the hydrological observation holes, carrying out multilayer division on the position of each hydrological observation hole in the target aquifer, acquiring the stress-strain relationship of each layer, then acquiring the total stress of the overlying rock-soil body corresponding to each hydrological observation hole during current observation and the pore water pressure corresponding to each layer of component layer, and acquiring the strain variation corresponding to the observation time, thereby calculating the compression amount of the soil body in a target time period; the method can effectively change the soil compression amount of the target aquifer and provide important data support for inversion prediction of regional ground surface settlement.

Description

Method for calculating distribution of water loss compression amount of pressure-bearing water-containing loose layer region
Technical Field
The invention relates to the field of stratum deformation monitoring, in particular to a method for calculating the distribution of compression deformation of regional stratum rock-soil layers.
Background
Along with the capture of a large amount of underground water resources by human beings and the loss of underground water resources in the mining process, the effective stress borne by the stratum framework is obviously increased due to the reduction of the pore water pressure after the water loss of the aquifer, thereby causing the compression of the aquifer. Particularly for loose aquifers, the porosity is large, the compression resistance is weak, the compressibility is good, the water level in the aquifer is reduced to cause the rapid compression of the aquifer, the compression quantity is transmitted to the ground surface along the stratum to cause the ground surface to sink, and the safe use of the ground surface and underground building structures is threatened. Therefore, the change characteristics of the water level of the pressure-bearing water-containing loose layer region are researched, so that the distribution of the water loss compression amount of the pressure-bearing water-containing loose layer region is calculated, and the method has important theoretical and practical significance for later-stage earth surface subsidence deduction and formulation of earth surface and underground structure deformation control schemes.
Disclosure of Invention
The invention aims to provide a simple and efficient calculation method for the distribution of the water loss compression amount of a region of a pressure-bearing water-containing loose layer, so that the deformation of the region stratum can be accurately monitored.
In order to achieve the technical purpose, hydrological holes are distributed in a calculation area, full-hole sampling is carried out, layer thickness data of each geotechnical layer is obtained, the density, thickness and stress strain characteristics of a target soil layer are obtained through a mechanical experiment, mechanical parameters such as the density and the thickness of the geotechnical layer are covered, water level monitoring in the area is carried out continuously, and water level change of the target aquifer is obtained. And calculating the reduction amount of the pore water pressure according to the variable quantity of the water level of the target aquifer to obtain the stress increment born by the soil framework. And calculating the total compression amount of the aquifer according to the elastic modulus of the land of different parts of the target aquifer, acquiring the compression amount of the target aquifer of different hole sites in the region, and then drawing an isoline to form a distribution map of the water loss compression amount of the region of the confined aquifer.
Specifically, in order to achieve the above object, the invention provides a method for calculating the distribution of the water loss compression amount of a region of a pressure-bearing water-containing loose layer, which comprises the following steps:
step 1, acquiring a target aquifer water level;
uniformly constructing a plurality of hydrological observation holes penetrating through a target aquifer downwards on the earth surface of the monitored area, numbering all the hydrological observation holes, observing, and obtaining the water head height of the ith hydrological observation hole in the jth observation
Figure BDA0003337039210000026
Step 2, obtaining the thickness H of the target aquifer corresponding to the position of the ith hydrological observation holeiDividing the target aquifer into n component layers from top to bottom, wherein the thickness of the corresponding k component layer is Hk iAnd is and
Figure BDA0003337039210000021
and obtaining the stress-strain relationship, epsilon, of each component layerik＝f(σik)，εikFor the strain quantity, σ, of the kth component layer at the ith hydrological observation holeikEffective stress born by a soil framework in the kth component layer at the ith hydrological observation hole;
step 3, obtaining the measurement parameter
The total stress of the overlying rock-soil mass borne by any one of the n component layers corresponding to the ith hydrological observation hole in any one-time observation is sigmaiAnd is and
Figure BDA0003337039210000022
wherein l is the total number of the overlying rock-soil mass layers above the ith hydrological observation hole, and each overlying rock-soil mass layer is numbered, wherein the ithThe thickness of the overlying rock-soil body of the o layer is HsoAnd a saturation density of ρ soG is the acceleration of gravity;
the pore water pressure borne by any one of the n component layers corresponding to the ith hydrological observation hole in the jth observation is all
Figure BDA0003337039210000023
And is
Figure BDA0003337039210000024
Where ρ iswThe average density of the target aquifer water body is obtained;
the effective stress born by the soil framework in the kth component layer in the n component layers corresponding to the ith hydrological observation hole in the jth observation is set as
Figure BDA0003337039210000025
Step 4, calculating the compression amount
The strain variation of the kth component layer in the n component layers corresponding to the ith hydrological observation hole in the j + m observation period is set as
Figure BDA0003337039210000031
The compression amount of the corresponding target aquifer at the ith hydrological observation hole during the j + m observations
Figure BDA0003337039210000032
And further, step 5, calculating the compression amount of each hydrological observation hole of the target water-bearing stratum during the j + m observation period in the j th time, and drawing a compression amount contour line to form a distribution map of the water loss compression amount of the area of the confined water-bearing loose layer.
Further, the total number of the component layers corresponding to each hydrological observation hole in the target aquifer is the same, different or partially the same.
Further, the thicknesses of all the component layers corresponding to the same hydrological observation hole are the same, different or partially the same.
Further, step 1 includes lowering a casing in the borehole to isolate hydraulic communication.
Further, the hydrological data in step 1 pass through a water level sensor.
Compared with the prior art, the invention has the following advantages:
the method can accurately determine the actual compression amount of the aquifer in the monitored area, provide important data support for inversion prediction of the surface settlement of the area, and provide important reference for damage early warning and protection of the surface building structure of the area.
Detailed Description
The technical solutions of the present invention will be described clearly and completely with reference to the following embodiments of the present invention, and it should be understood that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
A calculation method for distribution of water loss compression amount of a region of a pressure-bearing water-containing loose layer comprises the following steps:
step 1, acquiring a target aquifer water level;
uniformly constructing a plurality of hydrological observation holes penetrating through a target aquifer downwards on the ground surface of a monitored area, numbering all the hydrological observation holes, simultaneously carrying out whole-hole coring, putting a casing into the hydrological observation holes to isolate the hydrological observation holes from other aquifers without hydraulic communication, putting a water level sensor into the hydrological observation holes, monitoring the water level in the hydrological observation holes in real time, and obtaining the water head height of the ith hydrological observation hole in the jth observation
Figure BDA0003337039210000041
Step 2, measuring the physical and mechanical properties of the core according to the core taking result, and obtaining the thickness H of the target aquifer corresponding to the position of the ith hydrological observation holeiDividing the target aquifer into n component layers from top to bottom and the thickness of the kth component layerDegree of Hk iAnd is and
Figure BDA0003337039210000042
analyzing the stress-strain curve of each component layer under the condition of water drainage, and obtaining the stress-strain relation, epsilon, of each component layerik＝f(σik)，εikFor the strain quantity, σ, of the kth component layer at the ith hydrological observation holeikEffective stress born by a soil framework in the kth component layer at the ith hydrological observation hole;
step 3, obtaining the measurement parameter
According to the effective stress principle of saturated soil, the effective stress sigma' borne by the soil framework is sigma-muwThe following settings are made (where σ is the total stress experienced by the soil, μwPore water pressure, σ' is the effective stress borne by the soil skeleton):
in consideration of the fact that the overburden rock-soil mass of the target aquifer does not change basically in the observation process and the generated soil pressure does not change in the observation period, in the embodiment, the total stress of the overburden rock-soil mass borne by any one of the n component layers corresponding to the ith hydrological observation hole in any one observation is sigmaiAnd is and
Figure BDA0003337039210000043
wherein l is the total number of the overlying rock-soil mass layers above the ith hydrological observation hole, each overlying rock-soil mass layer is numbered, and the thickness of the overlying rock-soil mass layer o is HsoAnd the saturation density is rho soG is the acceleration of gravity;
considering that the content of the confined water may change (which is also a main factor for generating dehydration compression) at different observation time points, the hydrographic observation data generated at the same position are different when the content of the confined water changes, the pore water pressure also changes, and the change of the water density of the corresponding target aquifer is negligible, so in this embodiment, the pore water pressure borne by any one component layer of the n component layers corresponding to the ith hydrographic observation hole is made to be the same as the pore water pressure borne by any one component layer of the n component layers corresponding to the ith hydrographic observation hole during the jth observation
Figure BDA0003337039210000051
And is
Figure BDA0003337039210000052
Where ρ iswThe average density of the target aquifer water body is obtained;
the effective stress born by the soil framework in the kth component layer in the n component layers corresponding to the ith hydrological observation hole in the jth observation is made to be
Figure BDA0003337039210000053
Step 4, calculating the compression amount
The strain variation of the kth component layer in the n component layers corresponding to the ith hydrological observation hole in the j + m observation period is set as
Figure BDA0003337039210000054
The compression of the target aquifer at the ith hydrological observation hole during the jth observation of the jth + mth observation is
Figure BDA0003337039210000055
And 5, calculating the compression amount of each hydrological observation hole of the target water-bearing stratum during the j + m observation period in the j th time, and drawing a compression amount contour line to form a distribution map of the water loss compression amount of the pressure-bearing water-bearing loose layer region.
Particularly, considering that the depths of different positions of the target aquifer may be different, the total number of the component layers corresponding to each hydrological observation hole in the target aquifer is the same, different or partially the same according to actual needs.
Particularly, considering that the soil frameworks of the target aquifer at the same position and different heights have different strengths, the thicknesses of all the component layers corresponding to the same hydrological observation hole are the same, different or partially the same according to actual needs.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.

Claims (6)

1. A calculation method for distribution of water loss compression amount of a region of a pressure-bearing water-containing loose layer is characterized by comprising the following steps:
step 1, acquiring a target aquifer water level;
uniformly constructing a plurality of hydrological observation holes penetrating through a target aquifer downwards on the earth surface of the monitored area, numbering all the hydrological observation holes, observing, and obtaining the water head height of the ith hydrological observation hole in the jth observation
Figure RE-FDA0003488490300000016
Step 2, obtaining the thickness H of the target aquifer corresponding to the position of the ith hydrological observation holeiDividing the target aquifer into n component layers from top to bottom, wherein the thickness of the corresponding k component layer is Hk iAnd is and
Figure RE-FDA0003488490300000011
and obtaining the stress-strain relationship, epsilon, of each component layerik＝f(σik)，εikFor the strain quantity, σ, of the kth component layer at the ith hydrological observation holeikEffective stress born by a soil framework in the kth component layer at the ith hydrological observation hole;
step 3, obtaining the measurement parameter
The total stress of the overlying rock-soil mass borne by any one of the n component layers corresponding to the ith hydrological observation hole in any one-time observation is sigmaiAnd is and
Figure RE-FDA0003488490300000012
wherein l is the total number of the overlying rock-soil mass layers above the ith hydrological observation hole, and each overlying rock-soil mass layer is numbered, wherein the thickness of the overlying rock-soil mass layer o is HsoAnd a saturation density of ρ soG is gravity plusSpeed;
the pore water pressure borne by any one of the n component layers corresponding to the ith hydrological observation hole in the jth observation is all
Figure RE-FDA0003488490300000013
And is
Figure RE-FDA0003488490300000014
Where ρ iswThe average density of the target aquifer water body is obtained;
the effective stress born by the soil framework in the kth component layer in the n component layers corresponding to the ith hydrological observation hole in the jth observation is set as
Figure RE-FDA0003488490300000015
Step 4, calculating the compression amount
The strain variation of the kth component layer in the n component layers corresponding to the ith hydrological observation hole in the j + m observation period is set as
Figure RE-FDA0003488490300000021
The compression amount of the corresponding target aquifer at the ith hydrological observation hole during the j + m observations
Figure RE-FDA0003488490300000022
2. The method of claim 1, further comprising the step 5 of calculating the compression of the target aquifer at each hydrographic observation hole during the j + m observation period, and drawing a compression contour line to form a distribution map of the loss of water compression of the confined aquifer region.
3. The method of claim 2, wherein the total number of component layers corresponding to each hydrological observation hole in the target aquifer is any one of the same, different and partially the same.
4. The method of claim 3, wherein the thicknesses of the layers in all component layers corresponding to the same hydrological observation hole are the same, different or partially the same.
5. The method of claim 4, wherein step 1 further comprises lowering a casing within the borehole for isolating hydraulic communication.
6. The method of claim 5, wherein the hydrological data in step 1 is passed through a water level sensor.
CN202111297430.0A 2021-11-04 2021-11-04 Method for calculating distribution of water loss compression amount of pressure-bearing water-containing loose layer region Pending CN114139467A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111297430.0A CN114139467A (en) 2021-11-04 2021-11-04 Method for calculating distribution of water loss compression amount of pressure-bearing water-containing loose layer region

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111297430.0A CN114139467A (en) 2021-11-04 2021-11-04 Method for calculating distribution of water loss compression amount of pressure-bearing water-containing loose layer region

Publications (1)

Publication Number Publication Date
CN114139467A true CN114139467A (en) 2022-03-04

Family

ID=80392603

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111297430.0A Pending CN114139467A (en) 2021-11-04 2021-11-04 Method for calculating distribution of water loss compression amount of pressure-bearing water-containing loose layer region

Country Status (1)

Country Link
CN (1) CN114139467A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117456239A (en) * 2023-10-25 2024-01-26 江苏省无锡探矿机械总厂有限公司 Drilling machine equipment monitoring system and method for soil remediation

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117456239A (en) * 2023-10-25 2024-01-26 江苏省无锡探矿机械总厂有限公司 Drilling machine equipment monitoring system and method for soil remediation
CN117456239B (en) * 2023-10-25 2024-05-14 江苏省无锡探矿机械总厂有限公司 Drilling machine equipment monitoring system and method for soil remediation

Similar Documents

Publication Publication Date Title
CN104831762B (en) Deep basal pit artesian water precipitation dynamic monitoring system and monitoring method
Scherrer et al. Formation of runoff at the hillslope scale during intense precipitation
Rahardjo et al. Effect of antecedent rainfall on pore‐water pressure distribution characteristics in residual soil slopes under tropical rainfall
Wood et al. Quantifying macropore recharge: Examples from a semi‐arid area
Chae et al. A method for predicting the factor of safety of an infinite slope based on the depth ratio of the wetting front induced by rainfall infiltration
Askrinejad et al. Rainfall induced instabilities: a field experiment on a silty sand slope in northern Switzerland
CN102912780B (en) Sandy soil water-immersion testing method for loess collapsible deformation
Cheng et al. Is annual recharge coefficient a valid concept in arid and semi-arid regions?
Springman et al. Lesson learnt from field tests in some potentially unstable slopes in Switzerland
Cuthbert et al. Linking soil moisture balance and source-responsive models to estimate diffuse and preferential components of groundwater recharge
Rahardjo et al. Numerical analyses and monitoring performance of residual soil slopes
Cotecchia et al. Analysis of climate-driven processes in clayey slopes for early warning system design
Rezaur et al. Hydrologic behavior of residual soil slopes in Singapore
CN114139467A (en) Method for calculating distribution of water loss compression amount of pressure-bearing water-containing loose layer region
CN110377935B (en) Method for predicting change of excess pore water pressure along depth of soil body around pile in pile sinking process
CN114297950A (en) Method for calculating distribution of water loss compression amount of pressure-bearing water-containing loose layer region
Pisani et al. Hydrogeological model and hydraulic behaviour of a large landslide in the Italian Western Alps
Pirone et al. Cyclical suction characteristics in unsaturated slopes
Baiocchi et al. Sustainable yield of fractured rock aquifers: The case of crystalline rocks of Serre Massif (Calabria, Southern Italy)
Rahardjo et al. Monitoring and modeling of slope response to climate changes
Olivares et al. The role of suction and its changes on stability of steep slopes in unsaturated granular soils
Roman Quintero et al. Understanding hydrologic controls of slope response to precipitations through Machine Learning analysis applied to synthetic data
CN113075385A (en) Method for measuring transient saturation region range of deep slip mass
Askarinejad et al. Water Exfiltration from Bedrock: A Drastic Landslide Triggering Mechanism
Mojtahedi et al. Measurement of moisture and temperature profiles in different layers of soil

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