Limiting discharge calculation method and device based on tunnel periphery ecological water demand

Abstract

The present application provides a kind of based on tunnel perimeter ecological water demand limit discharge calculation method and device, it is related to tunnel engineering construction technical field, including: according to the simulation analysis of vegetation water data determines actual evapotranspiration data and potential evapotranspiration data, according to actual evapotranspiration data and potential evapotranspiration data determine the water stress factor of plant;The relationship between the water stress factor of different plants and groundwater level line is obtained by curve fitting;By setting threshold to the water stress factor of plant, determine the plant groundwater demand data corresponding to groundwater level line;The depth of emergency water well and the depth of water level drop caused by construction drainage are obtained, and the resident water reduction depth is determined according to the resident water data;Determine the well depth, determine the resident water demand data based on well depth and tunnel site water level data;According to plant groundwater demand data and resident water demand data, the tunnel limit discharge that meets the water demand of tunnel perimeter ecology is obtained.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of tunnel engineering construction, in particular to a limited discharge calculation method and device based on water demand of tunnel periphery ecology. BACKGROUND

[0002] This section is intended to provide background or context to the embodiments of the application. The description herein does not constitute admission that the information provided herein is prior art.

[0003] In the early process of mountain tunnel construction, the treatment measure for groundwater is mainly drainage, which drains all the groundwater seeping behind the lining. Although such design concept can avoid the lining bearing large water pressure and reduce the cost of water plugging, it will also cause many problems: a large amount of tunnel drainage will cause waste and damage of groundwater resources, decrease of surface water level, reduction or even dry-up of spring flow in the tunnel site area, and influence on the growth of surface vegetation and domestic water use. A large amount of groundwater drainage will easily take away a large amount of filling material in the surrounding rock fissures, expand the surrounding rock fissures and pipe cavity, and enhance the stratum permeability, causing vicious cycle. Since all the groundwater behind the lining is drained, the water pressure is generally not considered in the design of the lining, which may lead to insufficient design of the lining thickness and cause lining rupture during tunnel operation.

[0004] Although the groundwater treatment measure mainly by drainage can avoid the lining bearing large water pressure and reduce the cost of water plugging, it will also cause many problems, among which the most important problem is that a large amount of tunnel drainage will easily cause waste and damage of groundwater resources, decrease of surface water level, reduction or even dry-up of spring flow in the tunnel site area, and serious influence on the periphery ecology.

[0005] In summary, a technical scheme is needed which can effectively analyze the demand of periphery vegetation and domestic water, and avoid the ecological problems caused by excessive tunnel construction drainage through limited discharge. SUMMARY

[0006] To solve the problems in the prior art, the present application provides a limited discharge calculation method and device based on water demand of tunnel periphery ecology, which can make the tunnel drainage meet the demand of periphery ecological environment during tunnel construction, provide guidance for limited discharge of tunnel engineering construction with complex periphery ecological environment, and avoid the ecological problems caused by excessive tunnel construction drainage.

[0007] In a first aspect of the embodiments of the present application, a limited discharge calculation method based on water demand of tunnel periphery ecology is provided, which comprises:

[0008] Collecting vegetation water data and domestic water data;

[0009] determine actual evapotranspiration data and potential evapotranspiration data according to simulation analysis based on the vegetation water data, and determine the water stress factor of the plants according to the actual evapotranspiration data and the potential evapotranspiration data;

[0010] curve fitting is performed according to the water stress factor of the plants, so as to obtain the relationship between the water stress factor of different plants and the underground water level line; the underground water demand data of the plants corresponding to the underground water level line is determined by setting a threshold value for the water stress factor of the plants;

[0011] an underground water level depth acquisition module is configured to acquire the well depth for emergency water and the water level depth caused by construction drainage, and determine the resident water reduction depth according to the resident water data;

[0012] a resident water demand determination module is configured to determine the well depth for drilling according to the resident water reduction depth, the well depth for emergency water and the water level depth caused by construction drainage, and determine the resident water demand data based on the well depth for drilling and the water level data of the tunnel site area;

[0013] the tunnel limited discharge amount satisfying the ecological water demand around the tunnel is obtained according to the plant underground water demand data and the resident water demand data.

[0014] In a second aspect of the embodiment of the present application, a limited discharge calculation device based on the ecological water demand around the tunnel is provided, which comprises:

[0015] a data acquisition module configured to acquire vegetation water data and resident water data;

[0016] a water stress factor determination module configured to determine actual evapotranspiration data and potential evapotranspiration data according to simulation analysis based on the vegetation water data, and determine the water stress factor of the plants according to the actual evapotranspiration data and the potential evapotranspiration data;

[0017] a plant underground water demand determination module configured to perform curve fitting according to the water stress factor of the plants, so as to obtain the relationship between the water stress factor of different plants and the underground water level line; the underground water demand data of the plants corresponding to the underground water level line is determined by setting a threshold value for the water stress factor of the plants;

[0018] an underground water level depth acquisition module is configured to acquire the well depth for emergency water and the water level depth caused by construction drainage, and determine the resident water reduction depth according to the resident water data;

[0019] a resident water demand determination module is configured to determine the well depth for drilling according to the resident water reduction depth, the well depth for emergency water and the water level depth caused by construction drainage, and determine the resident water demand data based on the well depth for drilling and the water level data of the tunnel site area;

[0020] The tunnel limited discharge calculation module is used for obtaining the tunnel limited discharge meeting the water demand of the tunnel surrounding ecological water quantity according to the plant groundwater demand data and the resident water demand data.

[0021] In a third aspect of the embodiments of the present application, a computer device is provided, which comprises a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the limited discharge calculation method based on the water demand of the tunnel surrounding ecological water quantity when executing the computer program.

[0022] In a fourth aspect of the embodiments of the present application, a computer readable storage medium is provided, which stores a computer program, and the computer program implements the limited discharge calculation method based on the water demand of the tunnel surrounding ecological water quantity when executed by a processor.

[0023] In a fifth aspect of the embodiments of the present application, a computer program product is provided, which comprises a computer program, and the computer program implements the limited discharge calculation method based on the water demand of the tunnel surrounding ecological water quantity when executed by a processor.

[0024] The limited discharge calculation method and device based on the water demand of the tunnel surrounding ecological water quantity provided by the present application can calculate the water demand of vegetation growth and the water demand of resident life by analyzing the vegetation and resident life water of the engineering surrounding, so as to determine the tunnel limited discharge meeting the water demand of the tunnel surrounding ecological water quantity, provide guidance for the tunnel engineering construction limited discharge in the surrounding ecological environment, avoid the ecological problems caused by the excessive tunnel construction drainage, and have a wide application prospect, and provide strong data support for protecting the ecology and construction and rationality. BRIEF DESCRIPTION OF DRAWINGS

[0025] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings needed in the embodiment description will be briefly introduced. Obviously, the drawings in the following description are some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor.

[0026] Figure 1 is a flowchart of the limited discharge calculation method based on the water demand of the tunnel surrounding ecological water quantity according to an embodiment of the present application.

[0027] Figure 2 is a schematic diagram of the water transfer relationship of atmospheric precipitation according to an embodiment of the present application.

[0028] Figure 3 is a model schematic diagram of determining the well depth by comprehensively considering the construction dewatering, resident life water induced dewatering, and resident emergency water factors according to an embodiment of the present application.

[0029] Figure 4 Figure 1 is a schematic diagram of a tunnel surrounding ecological water demand-based limited discharge calculation device according to an embodiment of the present application.

[0030] Figure 5 Figure 2 is a schematic diagram of a computer device structure according to an embodiment of the present application. DETAILED DESCRIPTION

[0031] The principles and spirits of the present application will be described below with reference to several exemplary embodiments. It should be understood that these embodiments are given only to enable those skilled in the art to better understand and implement the present application, and do not limit the scope of the present application in any way. On the contrary, these embodiments are provided to make the present disclosure more thorough and complete, and to fully convey the scope of the present disclosure to those skilled in the art.

[0032] Those skilled in the art know that the embodiments of the present application can be implemented as a system, device, apparatus, method or computer program product. Therefore, the present disclosure can be embodied in the form of a complete hardware, a complete software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.

[0033] According to the embodiments of the present application, a tunnel surrounding ecological water demand-based limited discharge calculation method and device are provided, which relate to the technical field of tunnel construction.

[0034] The principles and spirits of the present application will be described below with reference to several exemplary embodiments. It should be understood that these embodiments are given only to enable those skilled in the art to better understand and implement the present application, and do not limit the scope of the present application in any way. On the contrary, these embodiments are provided to make the present disclosure more thorough and complete, and to fully convey the scope of the present disclosure to those skilled in the art.

[0035] Figure 1 Figure 1 is a schematic diagram of a tunnel surrounding ecological water demand-based limited discharge calculation device according to an embodiment of the present application. Figure 1 As shown in the figure, the method comprises:

[0036] S101, collecting vegetation water data and resident water data;

[0037] S102, determining actual evapotranspiration data and potential evapotranspiration data according to the vegetation water data, and determining water stress factors of plants according to the actual evapotranspiration data and the potential evapotranspiration data;

[0038] S103, performing curve fitting according to the water stress factors of plants to obtain the relationship between the water stress factors of different plants and the groundwater level line, and determining plant groundwater demand data corresponding to the groundwater level line by setting a threshold value for the water stress factors of plants;

[0039] S104, obtaining the well depth of emergency water and the water level drop depth caused by construction drainage, determining the resident water drop depth according to the resident water data;

[0040] S105, determining the well depth according to the resident water drop depth, the well depth of emergency water and the water level drop depth caused by construction drainage, determining the resident water demand data based on the well depth and the water level data of the tunnel site area;

[0041] S106, obtaining the tunnel limited discharge amount satisfying the tunnel surrounding ecological water demand according to the plant underground water demand data and the resident water demand data.

[0042] In order to more clearly explain the above-mentioned calculation method of limited discharge based on the tunnel surrounding ecological water demand, the following will be described in detail in combination with embodiments.

[0043] In an embodiment, for S102, the actual evapotranspiration data and the potential evapotranspiration data are determined by simulation analysis according to the vegetation water data, and the water stress factor of the plant is determined according to the actual evapotranspiration data and the potential evapotranspiration data, and the specific calculation method includes:

[0044] The water stress factor of the plant is determined based on the following calculation formula:

[0045]

[0046] In the formula, PWSI is the water stress factor of the plant; E ta is the actual evapotranspiration data; E tp is the potential evapotranspiration data;

[0047]

[0048] In the formula, E tp is the potential evapotranspiration data; L v is the latent heat of evaporation; R n is the net radiation flux of the vegetation canopy; Δ is the slope of the saturated water vapor pressure curve at the average temperature; ρ a is the air density; C p is the specific heat capacity of air at constant pressure; e s is the theoretical saturated water vapor pressure; e a is the actual saturated water vapor pressure; r a is the aerodynamic resistance; r s is the canopy surface resistance; γ is the dry-wet surface constant;

[0049]

[0050] In the formula, E ta is the actual evapotranspiration data: f(z) = 1 - exp(-z) for transpiration without compensation for water uptake: f umov f(z) = 1 - exp(-z) for transpiration without compensation for water uptake: f z = f(z) for potential evapotranspiration data: E r f(z) = 1 - exp(-z) for transpiration without compensation for water uptake: f tp E = f(z) for potential evapotranspiration data: E ia e = f(z) for interception evaporation: E rat z = f(z) for the ratio of intercepted potential evaporation to potential transpiration.

[0051] For the response function of soil water potential, the calculation relationship is:

[0052]

[0053] In the formula, f(ψ(z)) is the response function of soil water potential; ψ c is an empirical coefficient, the critical soil water potential value; ψ(z) is the soil water potential at depth z; p1 is an empirical coefficient; p2 is an empirical coefficient; E tp E = f(z) for potential evapotranspiration data; f θ is a function related to soil moisture content θ;

[0054] For the response function of osmotic potential, the calculation relationship is:

[0055]

[0056] In the formula, f(π(z)) is the response function of osmotic potential; n r is the total number of root layers; π(z) is the osmotic potential at depth z; π c is the critical osmotic potential value, an empirical coefficient; r i (Δz) is the relative root density of the i-th depth layer; Δz is the soil depth increment; p x is the sensitivity of osmotic potential to plant water uptake;

[0057] For the response function of soil temperature, the calculation relationship is:

[0058]

[0059] In the formula, f(T(z)) is the response function of soil temperature; T(z) is the soil temperature at depth z; T trig is the trigger temperature value, indicating the threshold value of soil temperature starting to significantly affect plant water uptake; t wA is an empirical coefficient; t wB is an empirical coefficient.

[0060] In an embodiment, for S103, a curve fitting is performed according to the water stress factor of the plant to obtain the relationship between the water stress factor of different plants and the groundwater level line; a threshold value is set for the water stress factor of the plant to determine the plant groundwater demand data corresponding to the groundwater level line, including:

[0061] By adjusting the groundwater level boundary of the simulation analysis, a plurality of sets of actual evapotranspiration data and potential evapotranspiration data are determined;

[0062] According to the plurality of sets of actual evapotranspiration data and potential evapotranspiration data, the water stress factors of a plurality of plants are determined;

[0063] According to the water stress factors of the plurality of plants, a curve fitting is performed to obtain the relationship between the water stress factor of different plants and the groundwater level line;

[0064] A threshold value is set for the water stress factors of the plurality of plants to determine the groundwater depth corresponding to the groundwater level line, and the plant groundwater demand data meeting the growth of the vegetation is determined based on the groundwater depth.

[0065] In an embodiment, for S104, the well drilling depth is determined according to the resident water reduction depth, the well depth of emergency water, and the water level drop depth caused by construction drainage, and the resident water demand data is determined based on the well drilling depth and the tunnel site area water level data, including:

[0066] The resident water reduction depth, the well depth of emergency water, and the water level drop depth caused by construction drainage are added to obtain the well drilling depth, and the calculation relationship is:

[0067] H = h1 + h2 + h3;

[0068] In the formula, H is the well drilling depth; h1 is the water level drop depth caused by construction drainage; h2 is the resident water reduction depth; and h3 is the well depth of emergency water.

[0069] More specifically, the method further includes:

[0070] The groundwater reduction depth is determined by the underground dynamics unstable flow theory, and the calculation formula is:

[0071]

[0072] In the formula, s is the groundwater reduction depth; q is the drainage rate; T is the water conductivity coefficient of the aquifer; W(u) is the well function for describing the relationship between the water level drop and the drainage condition; -0.577216 is the Euler constant; u is the specific yield; r is the average distance from the well wall; μ' is the flow potential of the drainage well; and t is the drainage time.

[0073] According to the groundwater lowering depth, the surface water level drawdown condition of the area around the well is analyzed to determine the influence range.

[0074] In the preset well water funnel drawdown (B) and influence range, the rainfall recharge water quantity is calculated, the drainage quantity is compared with the rainfall recharge water quantity, and the underground water balance condition is determined, wherein the comparison relationship is:

[0075]

[0076] Q2=qt2;

[0077] In the formula, Q1 is the rainfall recharge water quantity; p is the precipitation; W is the influence range; F is the precipitation frequency; t1 is the precipitation time; Q2 is the drainage quantity; q is the drainage rate; and t2 is the drainage time.

[0078] When Q1 < Q2, the rainfall recharge water quantity is smaller than the drainage quantity, and the water level line continuously decreases; when Q1 >= Q2, the water level does not decrease.

[0079] It should be noted that although the operations of the method of the present application are described in a specific order in the above embodiments and drawings, this does not require or imply that the operations must be performed in this specific order, or that all of the shown operations must be performed to achieve the desired result. Additionally or alternatively, certain steps can be omitted, multiple steps can be combined into one step, and / or one step can be divided into multiple steps.

[0080] The tunnel peripheral ecological water demand-based limited discharge calculation method of the present application will be described below in combination with a specific embodiment.

[0081] The present application can realize that the tunnel drainage quantity meets the peripheral ecological environment demand in the tunnel construction process through the calculation of the limited discharge standard of the high osmotic pressure tunnel in the water-rich environment. The present application mainly uses an ecological hydrological simulation simulation model and a resident water analysis model to analyze water.

[0082] 1. Ecological hydrological simulation simulation model:

[0083] Specifically, the ecological hydrological simulation simulation model (Coupmodel) is a comprehensive model for exchanging heat and material transport of the soil-vegetation-atmosphere system. Two coupled water and heat flow partial differential equations are solved by using the finite difference method. Reference Figure 2 As shown in the figure, the water balance of each day can be expressed by the following formula:

[0084] P=I+E s +E' ta +q deep +Δs+q surf ; formula (1)

[0085] where P is rainfall, mm; I is canopy interception, mm; E s is soil evaporation, mm; E' ta is vegetation transpiration, mm; q surf is surface runoff, mm; As is change of soil water storage, mm; q deep is deep percolation below the root zone, mm.

[0086] The model includes a series of associated sub-models for soil water dynamics, soil evaporation dynamics and plant water dynamics.

[0087] The soil water dynamics includes three sub-models, i.e., soil water flow process, surface runoff and deep percolation sub-models; the soil evaporation model has only one sub-model, i.e., an evaporation calculation model based on the Penman-Monteith equation; the plant water movement model includes four sub-models, i.e., vegetation potential transpiration, root water uptake, actual transpiration and canopy interception sub-models.

[0088] For the process of soil water dynamics, the water flow in soil is assumed to be laminar flow, mainly based on Darcy's law and Richards' equation, and to satisfy the law of conservation of mass:

[0089]

[0090] where q w is water flux (m 3 / m 2 ·s) ; k w is unsaturated hydraulic conductivity (ms) ; ψ is soil water potential (kPa) ; z is soil depth (m) ; c v is water vapor concentration in soil air (m 3 / m 2 ·s) ; D v is water diffusion rate in soil (m 2 ·s -1 ) ; q bypass is macropore flow (m 3 / m 2 ·s) ; θ is soil volumetric water content (m 3 / m 3 ) ; S w is a source or sink term.

[0091] To simplify the calculation, it is assumed that when the model soil pores are small and there is no water vapor in the soil air, c v = 0, q bypass = 0:

[0092]

[0093] Soil evaporation process:

[0094] Soil evaporation is based on the Penman-Monteith equation, whose calculation equation is:

[0095] The constraint condition is: q in = i cap , i cap ≤ q th ; formula (5)

[0096] i cap represents the inflow rate limit value; q in represents the inflow rate; q th is the maximum allowable inflow rate;

[0097]

[0098] In the formula, E s is the soil evaporation (mm·s -1 ); R ns is the net surface radiation flux (J·m -2 s -1 ); q h is the soil heat flux (wm -2 s -1 ); ρ a is the air density (kg / m 3 ); C p is the air constant-pressure specific heat capacity (J·m -1 ·℃ -1 ); e a is the actual saturated water vapor pressure (kPa); e s is the theoretical saturated water vapor pressure (kPa); r as is the surface aerodynamic resistance (s·m -1 ); r ss is the surface resistance (s·m -1 ); L v is the latent heat of evaporation (J·kg -1 ); Δ is the saturated water vapor pressure curve slope at the average air temperature (kPa·℃ -1 ); and γ is the dry and wet surface constant (kPa·℃ -1 ).

[0099] Plant water movement process:

[0100] Vegetation potential transpiration refers to the maximum possible evapotranspiration under conditions of ample and unrestricted water supply. This is primarily influenced by meteorological conditions and can also be calculated using the Penman-Monteith equation. Compared to the soil evaporation equation, the potential transpiration equation differs in three parameters, as shown in the following formula:

[0101] The constraint is: q in =i cap i cap ≤q th Equation (7)

[0102] i cap Indicates the inflow rate limit; q in Indicates the inflow rate; q th The maximum permissible inflow rate;

[0103]

[0104] In the formula, E tp Potential evapotranspiration data (mm·s) -1 );L v Latent heat of vaporization (J·kg) -1 ); R n Net radiation flux of vegetation canopy (J·m -2 s -1 ); Δ is the slope of the saturated vapor pressure curve at the average temperature (kPa·℃⁻¹); ρ a air density (kg / m³) 3 );C p The specific heat capacity of air at constant pressure (J·m) -1 ·℃ -1 );e s The theoretical saturated water vapor pressure (kPa); e a r is the actual saturated vapor pressure (kPa); a Aerodynamic impedance (s·m) -1 );r s Canopy surface impedance (s·m) -1 ); γ is the wet / dry constant (kPa℃) -1 );

[0105] Canopy surface impedance r s Related to the leaf area index (LAI) of plant canopy:

[0106]

[0107] In the formula, r s For the surface impedance of the canopy; g i g is an intermediate computational quantity. ris g is the radiation coefficient; max For maximum conductivity; gvpd is the water vapor pressure deficit; R is is the total solar radiation. For the conditions of normal oxygen supply to the roots, the root water uptake function follows the equation:

[0108]

[0109] where f θ is a function related to the soil moisture content θ; p ox is an empirical parameter; S ox is an intermediate calculation; θ is the soil moisture content; θ s is the soil saturated water content; θ ox is the soil moisture threshold.

[0110] The relationship between the actual transpiration of vegetation and the potential transpiration of vegetation, and the root morphology:

[0111]

[0112] where E ta is the actual evapotranspiration data (mm s -1 ): is the transpiration amount (mm s -1 ) without considering the compensation of water uptake: f umov is the compensation degree of root water uptake: is the potential transpiration amount (mm s -1 ) after the interception evaporation reduction: z r is the root depth (m): f(ψ(z)), f(π(z)), f(T(z)) are the response functions of the soil water potential, the osmotic potential, and the soil temperature, respectively: r(z) is the relative root density distribution: E tp is the potential evapotranspiration data (mm s -1 ): E ia is the interception evaporation amount (mm s -1 ): e rat is the ratio of the interception potential evaporation amount to the potential transpiration amount.

[0113] For the response function of the soil water potential, the calculation relationship is:

[0114]

[0115] where f(ψ(z)) is the response function of the soil water potential; ψ c is an empirical coefficient, the critical soil water potential value; ψ(z) is the soil water potential at the depth z; p1 is an empirical coefficient; p2 is an empirical coefficient; E tp is the potential evapotranspiration data; f θ is a function related to the soil moisture content θ;

[0116] The response function for osmotic potential is calculated as:

[0117]

[0118] where f(π(z)) is the response function for osmotic potential; n r is the total number of root layers; π(z) is the osmotic potential at depth z; π c is the critical osmotic potential value, an empirical coefficient; r i (Δz) is the relative root density of the ith depth layer; Δz is the soil depth increment; p x is the sensitivity of osmotic potential to plant water uptake;

[0119] For the osmotic potential π(z) at depth z, the relationship is as follows:

[0120]

[0121] where π(z) is the osmotic potential at depth z; R is the gas constant (J·mol -1 ·k -1 ); T is the soil temperature; C Cl (z) is the initial salt concentration of each protrusion; M Cl is the molar mass of Cl (g / mol).

[0122] The response function for soil temperature is calculated as:

[0123]

[0124] where f(T(z)) is the response function for soil temperature; T(z) is the soil temperature at depth z; T trig is the triggering temperature value, representing the threshold at which soil temperature begins to significantly affect plant water uptake; t wA is an empirical coefficient; t wB is an empirical coefficient.

[0125] The plant water stress indicator (PWSI) is used as an ecological indicator to quantitatively describe the impact of soil water deficit on vegetation caused by changes in hydrological processes. By considering the water demand and available water of plants, the ratio of actual evapotranspiration to potential evapotranspiration is used to represent the plant water stress indicator (PWSI) through simulation analysis:

[0126]

[0127] where PWSI is the plant water stress indicator; E ta is the actual evapotranspiration data; E tp is the potential evapotranspiration data;

[0128] The relationship between water stress factor of different plants and groundwater level line is obtained by curve fitting, and the threshold of water stress of the plants is set, so that the corresponding groundwater depth can be obtained, i.e., the water quantity meeting the growth requirement of the plants.

[0129] 2. Analysis model of well water for residents using water source:

[0130] For some areas, the original water source may be affected by tunnel construction drainage, and the water yield of the water source is reduced, or even cut off. In order to solve this problem, a new water source can be added by drilling a well.

[0131] The control depth of well drilling needs to consider two factors: 1. The water supply of the well, i.e. to ensure the daily water use of villagers and livestock and the emergency water use of residents; 2. The water permeability of the stratum, i.e. under the condition of the stratum, whether the water supply of atmospheric precipitation to the well water will cause the water level line to drop.

[0132] The model for determining the well drilling depth considering the factors such as construction dewatering, resident domestic water induced dewatering, and resident emergency water is shown in Figure 3 The width B of the dewatering funnel is the total width of the funnel shape formed by the dewatering.

[0133] The well drilling depth is obtained by adding the resident domestic water lowering depth, the well depth for emergency water, and the water level lowering depth caused by construction dewatering:

[0134] H=h1+h2+h3; equation (21)

[0135] In the equation, H is the well drilling depth (m); h1 is the water level lowering depth caused by construction dewatering (m); h2 is the resident domestic water lowering depth (m); and h3 is the well depth for emergency water (m).

[0136] The depth of the underground water level lowering (h1) caused by construction dewatering is mainly controlled by the measures of tunnel construction excavation. The well depth for emergency water (h3) is considered to be able to extract the emergency water quantity in a short time, so the groundwater recharge is ignored, and the water quantity of the well depth is determined.

[0137] The resident domestic water lowering depth (h2) is quantified by the theory of unsteady flow of underground dynamics through the permeability coefficient of the stratum and the resident water supply of the well, and finally the minimum control value of the well drilling depth H of different strata along the tunnel is determined.

[0138] The calculation of the drawdown of the underground water level in the well flow is a classic problem in the dynamics of groundwater. The non-steady flow of groundwater can be directly applied to single well dewatering according to the Theis formula. The expression is:

[0139]

[0140] wherein s is the groundwater lowering depth (m) ; q is the drainage rate (m 3 / d) ; T is the aquifer transmissivity (m 2 / d) ; W (u) is the well function, which is used to describe the relationship between the water level drawdown and the drainage condition; -0.577216 is the Euler constant; u is the specific yield; r is the average distance from the well wall (m) ; mu' is the flow potential of the drainage well; t is the drainage time (d) ;

[0141] The drawdown equation can be used to determine the drawdown of the surrounding water level when the well is drained, and the influence range can be determined by analyzing the drawdown of the surface water level in the area near the well.

[0142] When considering atmospheric rainfall recharge, it is considered that the rainfall recharge in the drainage influence range, and the water level outside the drainage influence area is not affected, and the excess water is discharged outside through horizontal runoff.

[0143]

[0144] Q2=qt2; equation (26)

[0145] wherein Q1 is the rainfall recharge water quantity; p is the precipitation; W is the influence range; F is the precipitation frequency; t1 is the precipitation time; Q2 is the drainage quantity; q is the drainage rate; t2 is the drainage time;

[0146] When Q1 < Q2, the rainfall recharge water quantity is less than the drainage quantity, and the water level line continues to drop; when Q1 >= Q2, the water level does not drop.

[0147] After the well water control depth is calculated, the groundwater level control standard is determined by combining the common well depth of the residents in the tunnel site area, that is, the water demand for meeting the residents' life.

[0148] Combined with the water quantity required for meeting the vegetation growth demand and the water demand for meeting the residents' life, the tunnel limited drainage quantity required for meeting the water demand of the surrounding ecology of the tunnel is calculated.

[0149] The ecological hydrological simulation model, the resident water analysis model adopted in the present application realize the calculation of the water demand of the surrounding environment of the tunnel, and provide a calculation method for meeting the vegetation demand and the water demand for residents' life for the establishment of the limited drainage standard of the complex tunnel surrounding ecological environment, avoid the ecological problems caused by excessive drainage quantity of the tunnel construction, and have a wide application prospect.

[0150] The present application considers the water demand of plants and residents, and provides guidance for the limited drainage of the water demand of the tunnel construction, and the overall scheme can protect the ecology, effectively protect the vegetation and the water demand for residents' life in the surrounding environment of the tunnel by calculating the water demand of the surrounding environment of the tunnel, and avoid the ecological problems caused by excessive drainage quantity of the tunnel construction.

[0151] Meanwhile, through the ecological hydrological simulation model and the resident water analysis model, the technology can accurately calculate the limited drainage standard during tunnel construction, realize the reasonable allocation and utilization of water resources, ensure the ecological environment water demand, and help maintain the safety and stability of the tunnel structure, avoiding problems such as tunnel lining rupture caused by improper groundwater discharge. By limiting the drainage amount of tunnel construction, the waste and damage to groundwater resources are reduced, and the risk of surface water level drop and spring flow reduction is lowered. The technology is not only suitable for tunnel projects with complex ecological environment, but also has wide popularization value and can be applied to various engineering projects that need to consider ecology and water resource management. The technology adopts the ecological hydrological simulation model and the resident water analysis model, which is an innovation in the field of tunnel construction, providing a new calculation method. The overall scheme has practical application potential, and by reasonably controlling the drainage amount of tunnel construction, the technology helps reduce the cost of environmental governance and ecological restoration, while ensuring the quality of life of residents, and has significant social and economic benefits.

[0152] After introducing the method of the exemplary embodiment of the present application, next, with reference to Figure 4 A calculation device for limited discharge based on water demand of ecology around tunnel is introduced.

[0153] The implementation of the calculation device for limited discharge based on water demand of ecology around tunnel can refer to the implementation of the above method, and the repeated parts will not be described again. The term "module" or "unit" used below can be a combination of software and / or hardware that realizes a predetermined function. Although the device described in the following embodiments is preferably realized in software, hardware or a combination of software and hardware is also possible and conceived.

[0154] Based on the same inventive concept, the present application also proposes a calculation device for limited discharge based on water demand of ecology around tunnel, as shown in Figure 4 The device comprises:

[0155] A data acquisition module 410 is configured to acquire vegetation water data and resident water data.

[0156] A water stress factor determination module 420 is configured to determine actual evapotranspiration data and potential evapotranspiration data through simulation analysis based on the vegetation water data, and determine the water stress factor of plants based on the actual evapotranspiration data and the potential evapotranspiration data.

[0157] A plant groundwater demand determination module 430 is configured to perform curve fitting based on the water stress factor of plants to obtain the relationship between the water stress factor of different plants and the groundwater level line, and determine the plant groundwater demand data corresponding to the groundwater level line by setting a threshold value for the water stress factor of plants.

[0158] The water level depth acquisition module 440 is configured to acquire the well depth of the emergency water and the water level drop depth caused by the construction drainage, and determine the resident water drop depth according to the resident water data;

[0159] The resident water demand determination module 450 is configured to determine the well depth according to the resident water drop depth, the well depth of the emergency water and the water level drop depth caused by the construction drainage, and determine the resident water demand data based on the well depth and the water level data of the tunnel site area.

[0160] The tunnel limited drainage amount calculation module 460 is configured to obtain the tunnel limited drainage amount for meeting the ecological water demand of the tunnel periphery according to the plant groundwater demand data and the resident water demand data.

[0161] It should be noted that although several modules of the limited drainage amount calculation device based on the ecological water demand of the tunnel periphery are mentioned in the foregoing detailed description, such division is merely exemplary and not mandatory. In fact, according to the embodiments of the present application, the features and functions of two or more modules described above can be embodied in one module. Conversely, the features and functions of one module described above can be further divided into several modules.

[0162] The present application is applicable to tunnel projects with complex surrounding ecological environment, and the overall scheme considers the water demand for the growth of vegetation around the tunnel, calculates the tunnel limited drainage amount for meeting the growth demand of vegetation by using an ecological hydrological simulation model, and considers the water demand for residents around the tunnel, calculates the tunnel limited drainage amount for meeting the growth demand of vegetation by using a resident water analysis model with well water as the water source.

[0163] The present application can realize the calculation of the water demand of the surrounding environment of the tunnel, provide guidance for the limited drainage amount of tunnel projects with complex surrounding ecological environment, and avoid ecological problems caused by excessive drainage amount of tunnel construction. The calculation method of the water demand of the surrounding environment of the tunnel proposed by the present application realizes the calculation of the water amount for meeting the ecological demand of the surrounding environment of the tunnel during the tunnel construction process, guides the limited drainage amount of the tunnel construction, and has strong popularization and economic value.

[0164] Based on the foregoing inventive concept, as shown in Figure 5 The present application further proposes a computer device 500, which comprises a memory 510, a processor 520, and a computer program 530 stored in the memory 510 and capable of running on the processor 520, and the processor 520 implements the foregoing limited drainage amount calculation method based on the ecological water demand of the tunnel periphery when executing the computer program 530.

[0165] Based on the foregoing inventive concept, the present application provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the foregoing tunnel surrounding ecological water demand based limited discharge calculation method.

[0166] Based on the foregoing inventive concept, the present application provides a computer program product, which comprises a computer program, and the computer program is executed by a processor to implement the tunnel surrounding ecological water demand based limited discharge calculation method.

[0167] The tunnel surrounding ecological water demand based limited discharge calculation method and device provided by the present application can calculate the water demand for vegetation growth and the water demand for resident life by analyzing the vegetation and resident life water demand of the surrounding engineering, so as to determine the tunnel limited discharge amount for meeting the tunnel surrounding ecological water demand, provide guidance for the tunnel engineering construction limited discharge amount in the surrounding complex ecological environment, avoid the ecological problems caused by excessive tunnel construction drainage amount, and have a wide application prospect, and provide strong data support for ecological protection and construction and rationality.

[0168] The data acquisition, storage, use, processing and the like in the technical solution of the present application comply with the relevant provisions of laws and regulations.

[0169] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, device or computer program product. Therefore, the present application can adopt a complete hardware embodiment, a complete software embodiment or an embodiment combining software and hardware aspects. Moreover, the present application can adopt a computer program product implemented on one or more computer usable storage media (including but not limited to disk memory, CD-ROM, optical memory, etc.) containing computer usable program code.

[0170] The present application is described with reference to flowcharts and / or block diagrams of the method and computer program product according to the embodiments of the present application. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of the flows and / or blocks in the flowcharts and / or block diagrams can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device produce a device for implementing the functions specified in the flowcharts and / or block diagrams. Figure 1 The device for implementing the functions specified in one flow or multiple flows and / or blocks. Figure 1 The device for implementing the functions specified in one block or multiple blocks.

[0171] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instructions which implement the flow Figure 1 of the flow or flows and / or blocks Figure 1 of the block or blocks specified in the flow.

[0172] The computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the flow Figure 1 of the flow or flows and / or blocks Figure 1 of the block or blocks specified in the flow.

[0173] Finally, it should be noted that the above-described embodiments are merely intended to illustrate the technical solutions of the present application, not to limit the same. The protection scope of the present application is not limited to the above-described embodiments, and although the above-described embodiments have been described in detail, those skilled in the art should understand that any modification or easy-to-think change or equivalent replacement of the technical features recorded in the above-described embodiments within the technical scope disclosed by the present application can be made by those skilled in the art, and these modifications, changes or replacements do not make the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application, and should be covered within the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A method for calculating limited emissions based on ecological water demand around tunnels, characterized in that, The method includes: Collect data on vegetation water use and residential water use; Based on the vegetation water use data, simulation analysis was performed to determine the actual evapotranspiration data and potential evapotranspiration data. Based on the actual evapotranspiration data and potential evapotranspiration data, the water stress factors of the plants were determined. Curve fitting was performed on the water stress factors of plants to obtain the relationship between the water stress factors of different plants and the groundwater level; by setting thresholds for the water stress factors of plants, the groundwater demand data of plants corresponding to the groundwater level were determined. Obtain the well depth for emergency water supply and the water level drop depth caused by construction drainage, and determine the drop depth for residential water supply based on the residential water use data; The well depth is determined based on the depth of the reduction in residential water use, the well depth for emergency water use, and the depth of the water level drop caused by construction drainage. Based on the well depth and the water level data of the tunnel site area, the residential water demand data is determined. Based on the plant groundwater demand data and the residents' water demand data, the tunnel discharge limit is determined to meet the ecological water demand around the tunnel.

2. The method for calculating limited emissions based on ecological water demand around tunnels according to claim 1, characterized in that, Based on the vegetation water use data, simulation analysis was performed to determine the actual evapotranspiration data and potential evapotranspiration data. Based on the actual and potential evapotranspiration data, the water stress factors for plants were determined, including: The water stress factors of plants are determined based on the following formula: In the formula, PWSI represents the plant's water stress factor; E ta This is actual evaporation data; E tp For potential evapotranspiration data; In the formula, E tp For potential evapotranspiration data; L v The latent heat of vaporization; R n ρ is the net radiation flux of the vegetation canopy; Δ is the slope of the saturated water vapor pressure curve at the average temperature; a C is the density of air. p The specific heat capacity of air at constant pressure; e s e is the theoretical saturated water vapor pressure; a r is the actual saturated vapor pressure; a For aerodynamic impedance; r s γ is the surface impedance of the canopy; γ is the dry and wet surface constant. In the formula, E ta Actual evaporation data: To disregard evaporation to compensate for water absorption: f umov Root water absorption compensation: To retain the potential evaporation after evaporation reduction: z r For root depth: f(ψ(z)), f(π(z)), and f(T(z)) are the response functions of soil water potential, osmotic potential, and soil temperature, respectively; r(z) is the relative root density distribution; E tp For potential evapotranspiration data: E ia To retain evaporation: e rat The ratio of potential evaporation to potential transpiration.

3. The method for calculating limited emissions based on ecological water demand around tunnels according to claim 2, characterized in that, The response function for soil water potential is calculated as follows: In the formula, f(ψ(z)) is the response function of soil water potential; ψ c ψ(z) is the empirical coefficient, the critical soil water potential value; p1 is the empirical coefficient; p2 is the empirical coefficient; E tp For potential evapotranspiration data; f θ It is a function related to soil moisture content θ; The response function for the osmotic potential is calculated as follows: In the formula, f(π(z)) is the response function of the osmotic potential; n r The root system has a total number of layers; π(z) is the permeability potential at depth z; π c r is the critical osmotic potential value, an empirical coefficient; i (Δz) represents the relative root density at the i-th depth layer; Δz represents the soil depth increment; p x The sensitivity of osmotic potential to the effect of plant water absorption; The response function for soil temperature is calculated using the following formula: In the formula, f(T(z)) is the soil temperature response function; T(z) is the soil temperature at depth z; T trig The trigger temperature value represents the threshold at which soil temperature begins to significantly affect plant water absorption; t wA For empirical coefficients; t wB This is an empirical coefficient.

4. The method for calculating limited emissions based on ecological water demand around tunnels according to claim 1, characterized in that, Curve fitting was performed on plant water stress factors to obtain the relationship between water stress factors of different plants and groundwater level; by setting thresholds for plant water stress factors, the groundwater demand data of plants corresponding to groundwater level were determined, including: By adjusting the groundwater level boundary in the simulation analysis, multiple sets of actual evapotranspiration data and potential evapotranspiration data were determined. Based on the aforementioned sets of actual evapotranspiration data and potential evapotranspiration data, water stress factors for various plants were determined; Based on the water stress factors of various plants, curve fitting was used to obtain the relationship between the water stress factors of different plants and the groundwater level. Thresholds are set for water stress factors of various plants to determine the groundwater depth corresponding to the groundwater level, and the groundwater demand data of plants that meet the growth requirements of the vegetation are determined based on the groundwater depth.

5. The method for calculating limited emissions based on ecological water demand around tunnels according to claim 1, characterized in that, The well depth is determined based on the depth of water level reduction for residential use, the depth of wells for emergency water use, and the depth of water level drop caused by construction drainage. Based on the well depth and water level data in the tunnel site area, residential water demand data is determined, including: The well depth is obtained by adding the depth of the residential water supply, the well depth for emergency water supply, and the water level drop caused by construction drainage. The calculation formula is as follows: H = h1 + h2 + h3; In the formula, H is the drilling depth; h1 is the water level drop caused by construction drainage; h2 is the drop in water level for residential use; and h3 is the well depth for emergency water use.

6. The method for calculating limited emissions based on ecological water demand around tunnels according to claim 1, characterized in that, The method also includes: The groundwater drawdown depth is determined using the theory of unsteady flow in subsurface dynamics, and the calculation formula is as follows: Where, s is the depth of groundwater drawdown; q is the drainage rate; T is the hydraulic conductivity of the aquifer; W(u) is the well function, used to describe the relationship between the drawdown and the drainage conditions; -0.577216 is the Euler constant; u is the specific yield; r is the average distance from the wellbore; μ' is the flow potential of the drainage well; t is the drainage time; Analyze the surface water drawdown situation in the area around the well according to the depth of groundwater drawdown, and determine the influence range.

7. The method for calculating limited emissions based on ecological water demand around tunnels according to claim 6, characterized in that, This method further includes: Within the preset drawdown and influence range of the well water funnel, calculate the rainfall recharge amount, compare the drainage amount with the rainfall recharge amount, and determine the groundwater balance situation, where the comparison relationship is: Q2 = qt2; Where, Q1 is the rainfall recharge amount; p is the precipitation; W is the influence range; F is the precipitation frequency; t1 is the precipitation time; Q2 is the drainage amount; q is the drainage rate; t2 is the drainage time; When Q1 < Q2, the rainfall recharge amount is less than the drainage amount, and the water level line continues to drop; when Q1 ≥ Q2, the water level does not drop.

8. A limited emission calculation device based on the ecological water demand around a tunnel, characterized in that, This device includes: A data acquisition module, used to acquire vegetation water use data and domestic water use data; A water stress factor determination module, used to determine the actual evapotranspiration data and potential evapotranspiration data through simulation analysis based on the vegetation water use data, and determine the plant water stress factor according to the actual evapotranspiration data and potential evapotranspiration data; A plant groundwater demand determination module, used to perform curve fitting based on the plant water stress factor to obtain the relationship between the water stress factor of different plants and the groundwater level line; by setting a threshold for the plant water stress factor, determine the plant groundwater demand data corresponding to the groundwater level line; A water level depth acquisition module, used to acquire the well depth of emergency water use and the water level drop depth caused by construction drainage, and determine the domestic water use reduction depth according to the domestic water use data; A domestic water demand determination module, used to determine the well drilling depth according to the domestic water use reduction depth, the well depth of emergency water use and the water level drop depth caused by construction drainage, and determine the domestic water demand data based on the well drilling depth and the water level data in the tunnel site area; A tunnel limit discharge calculation module, used to obtain the tunnel limit discharge that meets the ecological water use demand around the tunnel according to the plant groundwater demand data and the domestic water demand data.

9. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program, and when the computer program is executed by the processor, it implements the method according to any one of claims 1 to 7.

Images