Limited emission calculation method and device based on ecological water consumption demand around tunnel

Abstract

The invention provides a limited emission calculation method and device based on the ecological water consumption demand around a tunnel, and relates to the technical field of tunnel engineering construction.The method comprises the steps that simulation analysis is conducted according to vegetation water consumption data to determine actual evapotranspiration data and potential evapotranspiration data; determining a water stress factor of the plant according to the actual evapotranspiration data and the potential evapotranspiration data; the relation between the water stress factors of different plants and the underground water level line is obtained through curve fitting; determining plant groundwater demand data corresponding to the groundwater level line by setting a threshold value for the water stress factor of the plant; the well depth of emergency water use and the water level lowering depth caused by construction drainage are obtained, and the resident water use lowering depth is determined according to the resident water use data; determining well digging depth, and determining residential water demand data based on the well digging depth and the water level data of the tunnel site area; and according to the plant groundwater demand data and the resident water demand data, the tunnel discharge limit meeting the ecological water demand around the tunnel is obtained.

Description

Technical Field

[0001] The invention relates to the technical field of tunnel engineering construction, and in particular to a limited emission calculation method and device based on ecological water demand around a tunnel. Background Art

[0002] This section is intended to provide a background or context for embodiments of the present invention. No description herein is admitted to be prior art by virtue of its inclusion in this section.

[0003] In the early construction of mountain tunnels, the treatment measures for groundwater were mainly drainage, which was to drain all the groundwater that had infiltrated behind the lining. Although this design concept can prevent the lining from being subjected to large water pressure and reduce the cost of water blocking, it will also cause the following problems: large-scale drainage of the tunnel will cause waste and damage to groundwater resources, a drop in surface water levels, a reduction in the flow of springs in the tunnel site, or even dryness, which will affect the growth of surface vegetation and the living water use of residents. Large-scale discharge of groundwater can easily carry away a large amount of filling materials in the cracks of the surrounding rock, expand the cracks and cavities of the surrounding rock, enhance the permeability of the formation, and cause a vicious cycle. Since all the groundwater behind the lining is discharged, water pressure is generally not considered in the lining design, which may lead to insufficient design of the lining thickness and lining rupture during tunnel operation.

[0004] Although drainage-based groundwater treatment measures can prevent the lining from being subjected to high water pressure and reduce water blocking costs, they can also cause many problems. The most important problem is that large amounts of drainage from tunnels can easily lead to waste and damage of groundwater resources, causing surface water levels to drop, and spring water flow in the tunnel area to decrease or even dry up, which has a serious impact on the surrounding ecology.

[0005] To sum up, there is a technical solution that can effectively analyze the needs of surrounding vegetation and residents' domestic water use, and avoid ecological problems caused by excessive drainage during tunnel construction by limiting discharge. Summary of the invention

[0006] In order to solve the problems existing in the prior art, the present invention proposes a limited emission calculation method and device based on the ecological water demand in the surrounding area of ​​the tunnel, which can ensure that the tunnel drainage volume meets the surrounding ecological environment needs during the tunnel construction process, provide guidance for the limited drainage construction of tunnel projects with complex surrounding ecological environment, and avoid ecological problems caused by excessive drainage during tunnel construction.

[0007] In a first aspect of an embodiment of the present invention, a limited discharge calculation method based on the ecological water demand around a tunnel is proposed, the method comprising:

[0008] Collect vegetation water use data and residential water use data;

[0009] Performing simulation analysis based on the vegetation water use data to determine actual evapotranspiration data and potential evapotranspiration data, and determining a water stress factor of the plant based on the actual evapotranspiration data and the potential evapotranspiration data;

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

[0011] Obtain the depth of the well for emergency water use and the depth of the water level drop caused by construction drainage, and determine the depth of the residential water use reduction based on the residential water use data;

[0012] Determine the depth of the well according to the depth of the reduction of residential water use, the depth of the well for emergency water use, and the depth of the water level drop caused by construction drainage, and determine the residential water demand data based on the well depth and the water level data of the tunnel site;

[0013] According to the plant groundwater demand data and the residents' water demand data, the tunnel discharge limit that meets the ecological water demand around the tunnel is obtained.

[0014] In a second aspect of an embodiment of the present invention, a limited emission calculation device based on the ecological water demand around a tunnel is proposed, the device comprising:

[0015] Data collection module, used to collect vegetation water use data and residents' water use data;

[0016] A water stress factor determination module, used to determine actual evapotranspiration data and potential evapotranspiration data by simulation analysis based on the vegetation water use data, and determine the water stress factor of the plant based on the actual evapotranspiration data and the potential evapotranspiration data;

[0017] The plant groundwater demand determination module is used to perform curve fitting based on the water stress factor of the plant to obtain the relationship between the water stress factor of different plants and the groundwater level; by setting a threshold value for the water stress factor of the plant, the plant groundwater demand data corresponding to the groundwater level is determined;

[0018] A water level depth acquisition module is used to obtain the well depth of emergency water use and the depth of water level drop caused by construction drainage, and determine the depth of residential water use reduction based on the residential water use data;

[0019] A residential water demand determination module is used to determine the well depth according to the residential water reduction depth, the emergency water well depth and the water level drop depth caused by construction drainage, and determine the residential water demand data based on the well depth and the water level data of the tunnel site;

[0020] The tunnel discharge limit calculation module is used to obtain the tunnel discharge limit that meets the ecological water demand demand around the tunnel based on the plant groundwater demand data and the residents' water demand data.

[0021] In a third aspect of an embodiment of the present invention, a computer device is proposed, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, a limited emission calculation method based on the ecological water demand around a tunnel is implemented.

[0022] In a fourth aspect of an embodiment of the present invention, a computer-readable storage medium is proposed, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, a limited emission calculation method based on the ecological water demand around a tunnel is implemented.

[0023] In a fifth aspect of an embodiment of the present invention, a computer program product is proposed, which includes a computer program, and when the computer program is executed by a processor, a limited emission calculation method based on the ecological water demand around a tunnel is implemented.

[0024] The limited discharge calculation method and device based on the ecological water demand in the tunnel periphery proposed in the present invention analyzes the vegetation and the water use conditions of the residents in the vicinity of the project, calculates the water demand to meet the vegetation growth demand and the water demand to meet the residents' lives, thereby determining the tunnel limited discharge volume that meets the ecological water demand in the tunnel periphery, providing guidance for the limited discharge volume of tunnel engineering construction with complex surrounding ecological environment, avoiding ecological problems caused by excessive drainage during tunnel construction, having broad application prospects, and providing strong data support for protecting ecology and construction and rationality. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0026] Figure 1 It is a flow chart of a limited emission calculation method based on the ecological water demand around a tunnel according to an embodiment of the present invention.

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

[0028] Figure 3 It is a schematic diagram of a model for determining the depth of a well by comprehensively considering construction precipitation, precipitation caused by residents' domestic water use, and residents' emergency water use factors according to an embodiment of the present invention.

[0029] Figure 4 It is a schematic diagram of the architecture of a limited emission calculation device based on the ecological water demand around a tunnel according to an embodiment of the present invention.

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

[0031] The principles and spirit of the present invention will be described below with reference to several exemplary embodiments. It should be understood that these embodiments are provided only to enable those skilled in the art to better understand and implement the present invention, and are not intended to limit the scope of the present invention 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 will appreciate that the embodiments of the present invention may be implemented as a system, device, apparatus, method or computer program product. Therefore, the present disclosure may be specifically implemented in the following forms, namely: complete hardware, complete software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.

[0033] According to an embodiment of the present invention, a limited emission calculation method and device based on the ecological water demand around a tunnel are proposed, which relate to the technical field of tunnel engineering construction.

[0034] The principle and spirit of the present invention are explained in detail below with reference to several representative embodiments of the present invention.

[0035] Figure 1 FIG. 1 is a flow chart of a method for calculating limited discharge based on ecological water demand around a tunnel according to an embodiment of the present invention. Figure 1 As shown, the method includes:

[0036] S101, collect vegetation water use data and residents’ water use data;

[0037] S102, performing simulation analysis according to the vegetation water use data to determine actual evapotranspiration data and potential evapotranspiration data, and determining a water stress factor of the plant according to the actual evapotranspiration data and the potential evapotranspiration data;

[0038] S103, performing curve fitting 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; by setting a threshold value for the water stress factor of the plant, determining the groundwater demand data of the plant corresponding to the groundwater level;

[0039] S104, obtaining the well depth of emergency water use and the depth of water level drop caused by construction drainage, and determining the depth of residential water use reduction based on the residential water use data;

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

[0041] S106, obtaining a tunnel discharge limit that meets the ecological water demand demand around the tunnel based on the plant groundwater demand data and the residents' water demand data.

[0042] In order to explain more clearly the above-mentioned limited emission calculation method based on the ecological water demand around the tunnel, it is described in detail below in conjunction with an embodiment.

[0043] In one embodiment, for S102, actual evapotranspiration data and potential evapotranspiration data are determined by simulation analysis based on the vegetation water use data, and the water stress factor of the plant is determined based on the actual evapotranspiration data and the potential evapotranspiration data. The specific calculation method includes:

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

[0045]

[0046] Where PWSI is the water stress factor of plants; 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 impedance; r s is the canopy surface impedance; γ is the psychrometric constant;

[0049]

[0050] In the formula, E ta The actual evapotranspiration data is: The transpiration without considering compensating water absorption: f umov The root water absorption compensation degree is: is the potential evaporation after interception evaporation is reduced: r is the root depth: f(ψ(z)), f(π(z)), 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 E is the potential evapotranspiration data: ia is the intercepted evaporation capacity: e rat It is the ratio of intercepted potential evaporation to potential transpiration.

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

[0052]

[0053] Where f(ψ(z)) is the response function of soil water potential; ψ c is the empirical coefficient, critical soil water potential value; ψ(z) is the soil water potential at depth z; p1 is the empirical coefficient; p2 is the empirical coefficient; E tp is the 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] Where 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, 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 the osmotic potential to the water absorption of plants;

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

[0058]

[0059] Where 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 at which soil temperature begins to significantly affect plant water absorption; t wA is the empirical coefficient; t wB is the empirical coefficient.

[0060] In one embodiment, for S103, 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; by setting a threshold value for the water stress factor of the plant, the groundwater demand data of the plant corresponding to the groundwater level is determined, including:

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

[0062] determining water stress factors of a plurality of plants according to the plurality of sets of actual evapotranspiration data and potential evapotranspiration data;

[0063] According to the water stress factors of the various plants, curve fitting is used to obtain the relationship between the water stress factors of different plants and the groundwater level;

[0064] Thresholds are set for water stress factors of various plants, the groundwater depth corresponding to the groundwater level is determined, and plant groundwater demand data for satisfying vegetation growth is determined based on the groundwater depth.

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

[0066] The depth of the well for residential water use, the depth of the well for emergency water use, and the depth of the water level drop caused by construction drainage are added together to obtain the well depth. The calculation formula is:

[0067] H=h1+h2+h3;

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

[0069] More specifically, the method further includes:

[0070] The groundwater drawdown depth is determined by the underground dynamic instability flow theory, and the calculation formula is:

[0071]

[0072] 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, which is used to describe the relationship between water level drawdown and drainage conditions; -0.577216 is the Euler constant; u is the water supply degree; r is the average distance from the well wall; μ' is the flow potential of the drainage well; t is the drainage time;

[0073] According to the depth of groundwater drawdown, analyze the drawdown of the water table in the area around the well to determine the affected range.

[0074] Within the preset drawdown (B) of the well water funnel and the affected range, calculate the rainfall recharge volume, compare the drainage volume with the rainfall recharge volume to determine the groundwater balance situation, where the comparison relationship is:

[0075]

[0076] Q2 = qt2;

[0077] In the formula, Q1 is the rainfall recharge volume; p is the precipitation; W is the affected range; F is the precipitation frequency; t1 is the precipitation time; Q2 is the drainage volume; q is the drainage rate; t2 is the drainage time;

[0078] When Q1 < Q2, the rainfall recharge volume is less than the drainage volume, and the water level line continues to drop; when Q1 ≥ Q2, the water level does not drop.

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

[0080] The following describes the limited discharge calculation method of the present invention based on the ecological water consumption demand around the tunnel in combination with a specific embodiment.

[0081] Through the calculation of the limited discharge standard of the high-osmotic pressure tunnel in the rich water environment, the present invention can ensure that the tunnel drainage volume meets the ecological environment requirements during the tunnel construction process. The present invention mainly uses the ecological hydrological simulation model and the domestic water analysis model for water use analysis.

[0082] 1. Ecological hydrological simulation model:

[0083] Specifically, the ecological hydrological simulation model (Coupmodel) refers to a comprehensive model for the exchange of heat and mass transfer between the soil-vegetation and the atmosphere system. Through two coupled partial differential equations of water and heat flow, the equations are solved using the finite difference method. Refer to Figure 2 As shown, it is a schematic diagram of the water transfer relationship of atmospheric precipitation in an embodiment of the present invention. The daily water balance can be expressed by the following formula:

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

[0085] Where P is the rainfall, mm; I is the canopy interception, mm; E s is soil evaporation, mm; E' ta is the transpiration of vegetation, mm; q surf is the surface runoff, mm; Δs is the change in soil water storage, mm; q deep Deep penetration below the root layer, mm.

[0086] The model includes a series of related sub-models including soil moisture dynamics simulation, soil evaporation dynamics simulation and plant moisture dynamics simulation.

[0087] Among them, soil moisture dynamic simulation includes three sub-models, namely, soil moisture flow process, surface runoff and deep infiltration sub-models; soil evaporation simulation has only one sub-model, namely, the evaporation calculation model based on the Penman equation (Penman-Monteith equation); plant water movement simulation includes four sub-models, namely, vegetation potential transpiration, root water absorption, actual transpiration and canopy interception sub-models.

[0088] For the soil moisture dynamics simulation process, it is assumed that the water flow in the soil is laminar, mainly based on Darcy's law and Richards equation, and the law of conservation of mass must be satisfied:

[0089]

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

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

[0092]

[0093] Soil evaporation process:

[0094] Soil evaporation is based on the Penman-Monteith equation, which is calculated as:

[0095] The constraints are: in =i cap ,i cap ≤q th ; Formula (5)

[0096] i cap Indicates the inflow rate limit value; q in represents the inflow rate; q th is the maximum permissible 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 specific heat capacity of air at constant pressure (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 impedance (s·m -1 );r ss is the surface impedance (s·m -1 );L v is the latent heat of evaporation (J·kg -1 );Δ is the slope of the saturated water vapor pressure curve at the average temperature (kPa·℃ -1 ); γ is the psychrometric constant (kPa·℃ -1 ).

[0099] Plant water movement process:

[0100] Vegetation potential evapotranspiration refers to the maximum possible evapotranspiration under the condition of sufficient water supply and no restrictions. This is mainly affected by meteorological conditions and can also be calculated using the Penman-Monteith equation. Compared with the soil evaporation equation, there are three different parameters in the potential evapotranspiration equation, as shown in the following formula:

[0101] The constraints are: in =i cap ,i cap ≤q th ; Formula (7)

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

[0103]

[0104] In the formula, E tp is the potential evapotranspiration data (mm·s -1 );L v is the latent heat of evaporation (J·kg -1 );R n is the net radiation flux of the vegetation canopy (J·m -2 s -1 );Δ is the slope of the saturated water vapor pressure curve at the average temperature (kPa·℃-1); ρ a is the air density (kg / m 3 );C p is the specific heat capacity of air at constant pressure (J·m -1 ℃ -1 );e s is the theoretical saturated water vapor pressure (kPa); e a is the actual saturated water vapor pressure (kPa); r a is the aerodynamic impedance (s·m -1 );r s is the canopy surface impedance (s·m -1 ); γ is the psychrometric constant (kPa℃ -1 );

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

[0106]

[0107] In the formula, r s is the canopy surface impedance; g i is the intermediate calculation amount; g ris is the radiation coefficient; g max is the maximum conductivity; gvpd is the vapor pressure deficit; R is is the total solar radiation. Under the condition of normal oxygen supply to the roots, the root water absorption function follows the following equation:

[0108]

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

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

[0111]

[0112] In the formula, E ta is the actual evapotranspiration data (mm·s -1 ): is the transpiration without considering the compensation of water absorption (mm·s -1 ):f umov The root water absorption compensation degree is: is the potential evaporation after interception evaporation reduction (mm·s -1 ):z r is the root depth (m); 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 is the potential evapotranspiration data (mm·s -1 ): E ia is the intercepted evaporation capacity (mm·s -1 ):e rat It is the ratio of intercepted potential evaporation to potential transpiration.

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

[0114]

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

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

[0117]

[0118] Where 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, 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 the osmotic potential to the water absorption of plants;

[0119] For the penetration 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 soil temperature; C Cl (z) is the initial salt concentration of each protrusion; M Cl is the molar mass of Cl (g / mol).

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

[0123]

[0124] Where 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 at which soil temperature begins to significantly affect plant water absorption; t wA is the empirical coefficient; t wB is the empirical coefficient.

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

[0126]

[0127] Where PWSI is the water stress factor of plants; E ta is the actual evapotranspiration data; E tp is the potential evapotranspiration data;

[0128] Curve fitting is used to obtain the relationship between the water stress factors of different plants and the groundwater level. Then a threshold is set for the water stress of plants, which can correspond to the groundwater depth, that is, the amount of water required to meet the growth of vegetation.

[0129] 2. Analysis model of water use for residents with well water as water source:

[0130] In some areas, the original water source may be affected by the drainage of tunnel construction, and the water output may be reduced or even cut off. To solve this problem, additional wells can be used as new water sources.

[0131] There are two major factors to consider when controlling the depth of a well: 1. The water supply of the well, that is, to ensure daily water use for villagers and livestock and emergency water use for residents; 2. The permeability characteristics of the stratum, that is, the recharge effect of atmospheric precipitation on well water under the conditions of the stratum - whether the water supply will cause the water level to drop.

[0132] The model for determining the depth of wells by comprehensively considering factors such as construction precipitation, precipitation caused by residents' domestic water use, and residents' emergency water use is as follows: Figure 3 As shown, the precipitation funnel width B is the total width of the funnel shape formed by the entire precipitation.

[0133] The depth of the well is obtained by adding the depth of the residential water reduction, the depth of the emergency water well and the depth of the water level drop caused by construction drainage:

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

[0135] In the formula, H is the depth of the well (m); h1 is the depth of the water level drop caused by construction drainage (m); h2 is the depth of the reduction in residents' water use (m); and h3 is the depth of the well for emergency water use (m).

[0136] The depth of groundwater level drop caused by drainage during construction (h1) is mainly controlled by tunnel excavation measures. The depth of the well for emergency water use by residents (h3) is determined by the water volume of the well depth, considering that the emergency water volume can be extracted in a relatively short time. Therefore, the groundwater recharge is ignored.

[0137] The depth (h2) of water reduction for residents is quantified based on the permeability coefficient of the stratum and the water supply of the well through the theory of underground dynamic instability flow, and finally the minimum control value of the well drilling depth H in different strata along the tunnel is determined.

[0138] The calculation of groundwater level drawdown in well flow is a classic problem in groundwater dynamics. The Theis formula is used to calculate groundwater unsteady flow, which can be directly applied to single well dewatering. Its expression is:

[0139]

[0140] In the formula, s is the groundwater lowering depth (m); q is the drainage rate (m 3 / d); T is the hydraulic conductivity of the aquifer (m 2 / d); W(u) is the well function, which is used to describe the relationship between the water level drawdown and the drainage conditions; -0.577216 is the Euler constant; u is the specific yield; r is the average distance from the wellbore (m); μ' 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 during well drainage, 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 the recharge of atmospheric rainfall, it is considered that rainfall recharge occurs within the drainage influence range, and the water level is considered unaffected outside the drainage influence area, and the excess water is discharged through horizontal runoff.

[0143]

[0144] Q2 = qt2; Equation (26)

[0145] In the formula, 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;

[0146] When Q1 < Q2, that is, the rainfall recharge amount is less than the drainage amount, the water level line continues to drop; when Q1 ≥ Q2, the water level does not drop.

[0147] After calculating the control depth of the well water, combined with the common well depth of the residents in the tunnel site area, the groundwater level control standard is determined, that is, to meet the water use needs of residents' lives.

[0148] Combined with the water volume required to meet the vegetation growth needs and the water use needs of residents' lives, the limited drainage volume of the tunnel that meets the ecological water use requirements around the tunnel is calculated.

[0149] The ecological hydrological simulation and simulation model and the residential water use analysis model adopted by the present invention realize the calculation of the water demand of the surrounding environment of the tunnel, and provide a calculation method that meets the needs of the surrounding vegetation and the residential water use for formulating the limited drainage standard of the tunnel with a complex surrounding ecological environment, avoiding ecological problems caused by excessive drainage volume during tunnel construction, and having a wide application prospect.

[0150] The present invention considers the water use of plants and residents, provides guidance for the limited discharge of water used in tunnel construction, and the overall scheme can protect the ecology. By calculating the water demand of the surrounding ecological environment of the tunnel, it effectively protects the vegetation and the residential water use around the tunnel, and avoids ecological problems caused by excessive drainage volume during tunnel construction.

[0151] At the same time, through the eco-hydrological simulation model and the residential water use analysis model, the technology can accurately calculate the limited drainage standards during the tunnel construction process and realize the rational allocation and utilization of water resources. While ensuring the water demand of the ecological environment, the technology also helps to maintain the safety and stability of the tunnel structure and avoid problems such as tunnel lining rupture caused by improper groundwater discharge. By limiting the drainage volume of tunnel construction, the waste and damage to groundwater resources are reduced, and the risk of surface water level decline and spring flow reduction is reduced. This technology is not only suitable for tunnel projects with complex ecological environment, but also has a wide range of promotion value. It can be applied to various engineering projects that need to consider ecological and water resource management. This technology adopts the eco-hydrological simulation model and the residential water use analysis model. The combination of these models is an innovation in the field of tunnel engineering construction and provides a new calculation method. The overall solution has the potential for practical application. By reasonably controlling the drainage volume of tunnel construction, the technology helps to 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 invention, next, refer to Figure 4 A limited emission calculation device based on ecological water demand around a tunnel according to an exemplary embodiment of the present invention is introduced.

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

[0154] Based on the same inventive concept, the present invention also proposes a limited discharge calculation device based on the ecological water demand around the tunnel, such as Figure 4 As shown, the device comprises:

[0155] The data collection module 410 is used to collect vegetation water use data and residential water use data;

[0156] A water stress factor determination module 420 is used to perform simulation analysis based on the vegetation water use data to determine actual evapotranspiration data and potential evapotranspiration data, and determine the water stress factor of the plant based on the actual evapotranspiration data and the potential evapotranspiration data;

[0157] The plant groundwater demand determination module 430 is used to perform curve fitting 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; by setting a threshold value for the water stress factor of the plant, the plant groundwater demand data corresponding to the groundwater level is determined;

[0158] The water level depth acquisition module 440 is used to obtain the well depth of emergency water use and the water level drop depth caused by construction drainage, and determine the residential water use reduction depth according to the residential water use data;

[0159] The residential water demand determination module 450 is used to determine the well depth according to the residential water reduction depth, the emergency water well depth and the water level drop depth caused by construction drainage, and determine the residential water demand data based on the well depth and the water level data of the tunnel site;

[0160] The tunnel discharge limit calculation module 460 is used to obtain the tunnel discharge limit that meets the ecological water demand demand around the tunnel based on the plant groundwater demand data and the residents' water demand data.

[0161] It should be noted that although several modules of the limited emission calculation device based on the ecological water demand around the tunnel are mentioned in the above detailed description, this division is only exemplary and not mandatory. In fact, according to an embodiment of the present invention, the features and functions of two or more modules described above can be concretized in one module. Conversely, the features and functions of one module described above can be further divided into multiple modules for concretization.

[0162] The present invention is suitable for tunnel projects with complex surrounding ecological environments. The overall solution takes into account the water demand for vegetation growth around the tunnel, and uses an eco-hydrological simulation model to calculate the tunnel discharge limit that meets the vegetation growth needs; and takes into account the water consumption of residents around the tunnel, and uses a residents water use analysis model with well water as the water source to calculate the tunnel discharge limit that meets the vegetation growth needs.

[0163] The present invention can realize the calculation of the water demand of the surrounding environment of the tunnel, provide guidance for the limited discharge of tunnel engineering construction with complex surrounding ecological environment, and avoid ecological problems caused by excessive drainage during tunnel construction. The present invention proposes a calculation method for the water demand of the surrounding environment of the tunnel, realizes the calculation of the water consumption that meets the ecological needs of the surrounding area of ​​the tunnel during the tunnel construction process, guides the limited discharge of tunnel construction, and has strong promotion and economic value.

[0164] Based on the above invention concept, Figure 5 As shown, the present invention also proposes a computer device 500, including a memory 510, a processor 520, and a computer program 530 stored in the memory 510 and executable on the processor 520, wherein the processor 520 implements the aforementioned limited emission calculation method based on the ecological water demand around the tunnel when executing the computer program 530.

[0165] Based on the aforementioned inventive concept, the present invention proposes a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, it implements the aforementioned limited emission calculation method based on the ecological water demand around the tunnel.

[0166] Based on the aforementioned inventive concept, the present invention proposes a computer program product, which includes a computer program. When the computer program is executed by a processor, a limited emission calculation method based on the ecological water demand around the tunnel is implemented.

[0167] The limited discharge calculation method and device based on the ecological water demand in the tunnel periphery proposed in the present invention analyzes the vegetation and the water use conditions of the residents in the vicinity of the project, calculates the water demand to meet the vegetation growth demand and the water demand to meet the residents' lives, thereby determining the tunnel limited discharge volume that meets the ecological water demand in the tunnel periphery, providing guidance for the limited discharge volume of tunnel engineering construction with complex surrounding ecological environment, avoiding ecological problems caused by excessive drainage during tunnel construction, having broad application prospects, and providing strong data support for protecting ecology and construction and rationality.

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

[0169] It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as methods, devices, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0170] The present invention is described with reference to flowcharts and / or block diagrams of methods and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, 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 generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0171] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0172] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0173] Finally, it should be noted that the above-described embodiments are only specific implementations of the present invention, which are used to illustrate the technical solutions of the present invention, rather than to limit them. The protection scope of the present invention is not limited thereto. Although the present invention is described in detail with reference to the above-described embodiments, ordinary technicians in the field should understand that any technician familiar with the technical field can still modify the technical solutions recorded in the above-described embodiments within the technical scope disclosed by the present invention, or can easily think of changes, or make equivalent replacements for some of the technical features therein; and these modifications, changes or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims

1. A limited discharge calculation method based on the ecological water demand around the tunnel, characterized in that: The method includes: Collect vegetation water use data and residential water use data; Performing simulation analysis based on the vegetation water use data to determine actual evapotranspiration data and potential evapotranspiration data, and determining a water stress factor of the plant based on the actual evapotranspiration data and the potential evapotranspiration data; According to the water stress factor of plants, curve fitting is performed to obtain the relationship between the water stress factor of different plants and the groundwater level; by setting a threshold value for the water stress factor of plants, the groundwater demand data of plants corresponding to the groundwater level is determined; Obtain the depth of the well for emergency water use and the depth of the water level drop caused by construction drainage, and determine the depth of the residential water use reduction based on the residential water use data; Determine the depth of the well according to the depth of the reduction of residential water use, the depth of the well for emergency water use, and the depth of the water level drop caused by construction drainage, and determine the residential water demand data based on the well depth and the water level data of the tunnel site; According to the plant groundwater demand data and the residents' water demand data, the tunnel discharge limit that meets the ecological water demand around the tunnel is obtained.

2. The limited discharge calculation method based on the ecological water demand around the tunnel according to claim 1 is characterized in that: Performing simulation analysis based on the vegetation water use data to determine actual evapotranspiration data and potential evapotranspiration data, and determining the water stress factor of the plant based on the actual evapotranspiration data and the potential evapotranspiration data, including: The water stress factor of plants was determined based on the following calculation formula: Where PWSI is the water stress factor of plants; E ta is the actual evapotranspiration data; E tp is the potential evapotranspiration data; 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 impedance; r s is the canopy surface impedance; γ is the psychrometric constant; In the formula, E ta The actual evapotranspiration data is: The transpiration without considering compensating water absorption: f umov The root water absorption compensation degree is: is the potential evaporation after interception evaporation is reduced: r is the root depth: f(ψ(z)), f(π(z)), 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 E is the potential evapotranspiration data: ia is the intercepted evaporation capacity: e rat It is the ratio of intercepted potential evaporation to potential transpiration.

3. The limited discharge calculation method based on the ecological water demand around the tunnel according to claim 2 is characterized in that: For the response function of soil water potential, the calculation relationship is: Where f(ψ(z)) is the response function of soil water potential; ψ c is the empirical coefficient, critical soil water potential value; ψ(z) is the soil water potential at depth z; p1 is the empirical coefficient; p2 is the empirical coefficient; E tp is the potential evapotranspiration data; f θ is a function related to soil moisture content θ; For the response function of osmotic potential, the calculation relationship is: Where 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, 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 the osmotic potential to the water absorption of plants; For the response function of soil temperature, the calculation relationship is: Where 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 at which soil temperature begins to significantly affect plant water absorption; t wA is the empirical coefficient; t wB is the empirical coefficient.

4. The limited discharge calculation method based on the ecological water demand around the tunnel according to claim 1 is characterized in that: According to the water stress factor of plants, curve fitting is performed to obtain the relationship between the water stress factor of different plants and the groundwater level; by setting a threshold for the water stress factor of plants, the groundwater demand data of plants corresponding to the groundwater level is determined, including: By adjusting the groundwater level boundary of the simulation analysis, multiple sets of actual evapotranspiration data and potential evapotranspiration data are determined; determining water stress factors of a plurality of plants according to the plurality of sets of actual evapotranspiration data and potential evapotranspiration data; According to the water stress factors of the various plants, curve fitting is 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, the groundwater depth corresponding to the groundwater level is determined, and plant groundwater demand data for satisfying vegetation growth is determined based on the groundwater depth.

5. The limited discharge calculation method based on the ecological water demand around the tunnel according to claim 1 is characterized in that: The depth of the well is determined according to the depth of the reduction of residential water use, the depth of the well 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, the water demand data of residents is determined, including: The depth of the well for residential water use, the depth of the well for emergency water use, and the depth of the water level drop caused by construction drainage are added together to obtain the well depth. The calculation formula is: H=h1+h2+h3; In the formula, H is the depth of the well; h1 is the depth of the water level drop caused by construction drainage; h2 is the depth of the reduction in residential water use; and h3 is the depth of the well for emergency water use.

6. The limited discharge calculation method based on the ecological water demand around the tunnel according to claim 1 is characterized in that: The method further includes: The groundwater drawdown depth is determined by the underground dynamic instability flow theory, and the calculation formula is: In the formula, s is the groundwater lowering depth; 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; Based on the groundwater lowering depth, analyze the drawdown of the surface water level in the area around the well and determine the influence range.

7. The limited discharge calculation method based on the ecological water demand around the tunnel according to claim 6 is 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; In the formula, 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 discharge calculation device based on the ecological water demand around the tunnel, characterized in that: This device includes: A data acquisition module for acquiring vegetation water use data and domestic water use data; A water stress factor determination module for determining actual evapotranspiration data and potential evapotranspiration data through simulation analysis based on the vegetation water use data, and determining the plant water stress factor based on the actual evapotranspiration data and the potential evapotranspiration data; A plant groundwater demand determination module for performing 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, determining the plant groundwater demand data corresponding to the groundwater level line; A water level depth acquisition module for acquiring the well depth of emergency water use and the water level drop depth caused by construction drainage, and determining the domestic water use reduction depth based on the domestic water use data; A domestic water demand determination module for determining the well drilling depth based on the domestic water use reduction depth, the well depth of emergency water use, and the water level drop depth caused by construction drainage, and determining the domestic water demand data based on the well drilling depth and the water level data in the tunnel site area; A tunnel limited discharge calculation module for obtaining the tunnel limited discharge amount that meets the ecological water use demand around the tunnel based on 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