[go: up one dir, main page]

CN114254484B - A method for arranging monitoring wells and determining tracer injection boundaries for landfill leakage detection - Google Patents

A method for arranging monitoring wells and determining tracer injection boundaries for landfill leakage detection Download PDF

Info

Publication number
CN114254484B
CN114254484B CN202111407272.XA CN202111407272A CN114254484B CN 114254484 B CN114254484 B CN 114254484B CN 202111407272 A CN202111407272 A CN 202111407272A CN 114254484 B CN114254484 B CN 114254484B
Authority
CN
China
Prior art keywords
tracer
concentration
landfill
injection
migration
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111407272.XA
Other languages
Chinese (zh)
Other versions
CN114254484A (en
Inventor
王巧
谢海建
费爽珂
丁昊
詹良通
陈赟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202111407272.XA priority Critical patent/CN114254484B/en
Publication of CN114254484A publication Critical patent/CN114254484A/en
Application granted granted Critical
Publication of CN114254484B publication Critical patent/CN114254484B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Processing Of Solid Wastes (AREA)
  • Examining Or Testing Airtightness (AREA)

Abstract

本发明公开了一种用于填埋场渗漏探测的监测井布置和示踪剂注入边界确定方法。本发明旨在利用示踪剂在填埋体和阻隔屏障的运移规律确定现场示踪剂试验注入井中示踪剂的注入量及监测井的布置。考虑示踪剂在填埋体和阻隔屏障中的扩散和对流作用,结合菲克定律、达西定律和质量守恒方程建立示踪剂运移的二维控制方程,利用拉普拉斯‑傅里叶变换和Talbot数值反演得到方程的解。根据填埋场覆盖层、垃圾体和底部衬垫的尺寸、孔隙率和渗透特性等不同特征,确定示踪剂的注入量和阻隔屏障渗漏监测井的布置方法。本发明提出的方法计算高效、可靠,可用于在役填埋场阻隔屏障的渗漏探测和高效定位,对于填埋场的治理具有重要的工程价值。

The present invention discloses a method for arranging monitoring wells and determining tracer injection boundaries for landfill leakage detection. The present invention aims to determine the injection amount of tracers in field tracer test injection wells and the arrangement of monitoring wells by using the migration law of tracers in landfill bodies and barriers. Considering the diffusion and convection of tracers in landfill bodies and barriers, a two-dimensional control equation for tracer migration is established in combination with Fick's law, Darcy's law and mass conservation equation, and the solution of the equation is obtained by Laplace-Fourier transform and Talbot numerical inversion. According to different characteristics of the landfill cover layer, garbage body and bottom liner, such as size, porosity and permeability characteristics, the injection amount of tracers and the arrangement method of barrier barrier leakage monitoring wells are determined. The method proposed by the present invention is efficient and reliable in calculation, can be used for leakage detection and efficient positioning of barrier barriers in service landfills, and has important engineering value for the management of landfills.

Description

Monitoring well arrangement and tracer injection boundary determination method for landfill leakage detection
Technical Field
The invention belongs to the field of landfill leakage monitoring, and particularly relates to a monitoring well arrangement and tracer injection boundary determination method for landfill leakage detection.
Background
The landfill is composed of a covering layer, a garbage body and a bottom liner, and meanwhile, the pollution of the landfill to underground water, soil and atmosphere is reduced through the bottom liner and the covering layer. However, due to factors such as construction, dry and wet circulation, climate change and the like, the defects such as tearing, loophole and clay cracking of the geomembrane forming the covering layer and the bottom liner of the landfill are easily caused. This results in landfill gas escaping from the overburden and leachate escaping from the bottom pad, contaminating the atmosphere and groundwater. Assuming that the performance of other facilities in the landfill is kept good, the leakage of a 25cm diameter leak in the bottom pad can reach 15 kg/day, about 5.4 tons/year. Leakage of landfill cover layer allows landfill gas to escape, resulting in greenhouse effect and odor disturbing the citizens. There are hundreds of malodorous gases in landfill gas, which once leaked, can severely pollute the environment and affect the health of surrounding residents.
The gas tracer test method can be used for determining various properties of the gas-coated belt system, such as water content, gas flow rate, curvature and the like, and can also be used for testing a gas flow mode in a landfill site and quantifying the residence time of volatile pollutants. The main idea of using tracers to detect leaks is to inject a harmless, environmentally low content, inert and easily monitored tracer gas (such as SF 6) through the injection port into the barrier system (cover system and bottom liner system) creating a large area tracer source. When there is a gap in the barrier, the tracer gas will diffuse through the gap into the adjacent medium, and the concentration of the tracer gas will be measured at different times in the soil gas of the sampling hole outside the barrier. The resulting trace gas concentration versus time data (with the position of the sampling holes, the nature of the medium and the trace gas source concentration) will be provided to the global optimization code. Based on the input relevant parameters, the code iteration will find an optimal solution to accurately locate the leak location. However, the migration law of the tracer in the landfill system is not clear at present, and in order to determine the injection amount of the tracer in the injection well and determine the arrangement method of the monitoring well, it is very necessary to know the migration law of the tracer in the landfill body.
Disclosure of Invention
The invention aims to provide a tracer injection amount and leakage position determining method for determining leakage detection of a landfill, aims at detecting leakage in the landfill in a later operation stage of the landfill in the case of serious leakage problem of the landfill in China, provides a calculation method for detecting the tracer injection amount by using a tracer and a placement method for monitoring a well, and has important significance for preliminary judgment of leakage detection of the landfill.
The invention aims at the leakage of barrier barriers such as a bottom liner, a covering layer and the like in a landfill site by utilizing a tracer, and determines the injection quantity and the leakage position of the tracer in an on-site tracer test injection well by utilizing the migration rule of the tracer in the landfill body and the barrier barriers. And (3) taking the diffusion and convection effects of the tracer in the landfill and the barrier into consideration, establishing a two-dimensional control equation of tracer migration by combining the Phak law, the Darcy law and the mass conservation equation, and obtaining a solution of the equation by utilizing Laplace-Fourier transform and Talbot numerical inversion. The injection amount of the tracer and the arrangement method of the barrier leakage monitoring well are determined according to different characteristics such as the size, the porosity, the permeability and the like of the covering layer, the garbage body and the bottom liner of the landfill. The method provided by the invention can be used for leakage detection and positioning of the blocking barrier of the landfill, and has important engineering value for the treatment of the landfill.
The invention is realized by the following technical scheme:
The invention firstly provides a monitoring well arrangement for landfill leak detection and a tracer injection boundary determination method, which are characterized by comprising the following steps:
step 1), constructing a physical model of the gas phase tracer injected into the landfill;
Step 2), establishing a control equation and solving the control equation aiming at the physical model;
Step 3.1), calculating the concentration profile and release flux of the tracer agent moving in the landfill system, and drawing the time and space change curves of the tracer agent;
and 3.2) determining the position coordinates of the monitoring points according to the concentration profile and the detection limit of the tracer drawn in the step 3.1), and finding out the points with the concentration equal to the detection limit of the tracer on the calculated concentration profile, wherein the positions of the points and the points with the concentration larger than the points are positions where the monitoring wells can be arranged.
And 4) determining the release flux and concentration of the tracer at the monitoring point, comparing the release flux and the detection limit, determining the injection concentration boundary of the injection well in the step 1) as the actual injection concentration if the concentration or release flux of the tracer at the monitoring point is larger than the detection limit, and repeating the steps 3.1) -4) if the concentration or release flux of the tracer at the monitoring point is smaller than the detection limit, the tracer cannot be detected at the monitoring point, increasing the injection boundary concentration or flux of the injection well.
As a preferred embodiment of the present invention, the step 1) specifically includes the following steps:
Determining the thickness of the garbage body as l 1, the thickness of the covering layer as l 2 and the thickness of the bottom liner as l 3 based on the construction scale of the sealed field landfill, measuring the water content of the covering layer, the garbage body and the bottom liner of the landfill on site, obtaining the diffusion coefficient and the convection rate of the gas phase tracer in the covering layer, the garbage body and the bottom liner of the landfill according to the diffusion coefficient calculation model (formula 1) and the Darcy law (formula 2),
Wherein D and D 0 are the diffusion coefficients of the tracer in the soil body and the air respectively, theta a is the gas-containing porosity of the soil body, n is the porosity of the soil body, v is the convection rate of the tracer in the soil body, k a is the permeability coefficient of the tracer in the soil body, and dP/dz is the air pressure gradient at two ends of the soil body.
As a preferred embodiment of the present invention, the step 2) specifically includes the following steps:
step 2.1) determining control equations and boundary conditions
Establishing a two-dimensional radial axisymmetric coordinate system by taking a garbage body middle point as an origin, establishing a two-dimensional control equation of tracer migration, injecting the tracer under the condition of constant concentration C 0, and determining boundary conditions used by a model;
Step 2.2) obtaining the migration rule of the tracer in the landfill body, which specifically comprises the following steps:
Step 2.2.1) calculating the migration concentration profile of the tracer on the Laplace domain;
Step 2.2.2) calculating the release flux of the surface tracer on the covering layer on the Laplace domain, and calculating the concentration of the tracer at the bottom of the liner on the Laplace domain;
Step 2.2.3) obtaining the concentration profile and the release flux of the tracer by utilizing Talbot numerical inversion.
As a preferred embodiment of the present invention, the step 2.2.1) specifically includes:
Injecting the tracer agent into the landfill body at a constant concentration boundary
Wherein C 0 is the concentration of the tracer in the injection well, p is the time factor on the Laplace domain, mu n and omega n are characteristic roots, r is the radial coordinate, r w is the radius of the injection well, z is the vertical coordinate, and m 2 is the migration parameter of the tracer in the garbage body, and can be calculated by the following formula:
Wherein v is the convection rate of the tracer in the garbage body, theta is the pore rate of the garbage body, D z is the diffusion coefficient of the tracer in the garbage body along the axial direction, and D r is the diffusion coefficient of the tracer in the garbage body along the radial direction;
The concentration of the tracer in the coating is
Where m 3 and m 4 are migration parameters of the tracer in the overburden, and can be calculated from the following formulas, respectively:
Wherein v u is the convection rate of the tracer in the coating, D u is the diffusion coefficient of the tracer in the coating along the axial direction;
Concentration profile of tracer in bottom pad
Where m 5 and m 6 are migration parameters of the tracer in the bottom liner, and are solved by the following equations, respectively:
Where v l is the convection velocity of the tracer in the bottom liner and D l is the diffusion coefficient of the tracer in the bottom liner in the axial direction.
As a preferred embodiment of the present invention, the step 2.2.2) specifically includes:
calculating the surface tracer release flux of the cover layer on the Laplace domain:
the liner bottom tracer concentration on the laplace domain was calculated:
Step 2.2.3), specifically:
the concentration profile and the release flux of the tracer are obtained by utilizing Talbot numerical inversion, and the calculation formula is as follows:
Wherein M is a superposition constant, taken as 64, re is a real part of a complex number, gamma ii is a Laplace inversion constant, and is defined by the following formula:
Wherein a is an imaginary factor
Compared with the prior art, the invention has the beneficial effects that:
According to the method, the injection quantity of the tracer in the injection well in the on-site tracer experiment can be calculated on the premise of considering the diffusion and convection of the tracer, and meanwhile the arrangement of the monitoring well is determined.
The invention can calculate the migration concentration profile and the release flux of the tracer according to the different water contents of the covering layer, the garbage body and the bottom liner of the landfill. So that the migration law of the tracer can be obtained under different weather conditions. By comparison with the detection limit of the tracer, the placement of the monitoring well can be determined, while the injection boundary of the injection well tracer is determined. The approximate position of the leakage point can be obtained preliminarily by comparing with the field actual measurement data.
By comparing with the numerical simulation theory calculation, the method has good reliability.
Drawings
FIG. 1 is a flow chart of the calculation of landfill leak detection tracer migration calculation of the present invention.
Fig. 2 is a computational model of the present invention.
Figure 3 is a cloud of calculated tracer concentration profiles of the invention.
FIG. 4 is a graph of tracer release flux over time for different gas-containing porosities of the computational model of the invention.
Detailed description of the preferred embodiments
The invention is further illustrated in the following figures and examples, which are not intended to limit the scope of the invention.
The invention mainly determines the size of a landfill according to a landfill design drawing, and mainly comprises the depth of the landfill, the thickness l 3 of a bottom liner, the thickness l 1 of a garbage body and the thickness l 2 of a covering layer. The water content of landfill garbage, the water content of the upper liner, the types of the covering layer and the bottom liner materials are obtained through on-site detection, and the diffusion coefficient and the convection rate of the tracer in a landfill system are determined. Injection boundary conditions for the tracer in the injection well are established and the placement of the well is monitored.
The invention is further illustrated by the following calculation of migration law of SF 6 tracer in simplified model of landfill as an example:
As shown in figure 1, the method comprises 4 steps, namely, constructing a physical model of the migration of the tracer in the landfill (injecting the gas-phase tracer into the physical model of the landfill), mainly determining the types and parameters of each part of the landfill, the migration parameters of the tracer and the boundary conditions of an injection well, namely, establishing a control equation and solving the equation to obtain a tracer migration concentration profile and a surface release flux expression, and step 3, drawing a tracer migration concentration profile and a release flux curve, simultaneously determining the positions of monitoring points and the distances between the monitoring points and the injection well (namely, calculating the concentration profile and the release flux of the tracer in a landfill system, drawing a time-space change curve, and determining the position coordinates of the monitoring points according to the drawn concentration profile and the detection limit of the tracer). And step 4, determining the concentration and the release flux of the tracer at the monitoring point, comparing the concentration or the release flux of the tracer at the monitoring point with the detection limit, determining the injection concentration boundary of the injection well in the step 1 as the actual injection concentration if the concentration or the release flux of the tracer at the monitoring point is larger than the detection limit, and increasing the injection boundary concentration or the flux of the injection well if the concentration or the release flux of the tracer at the monitoring point is smaller than the detection limit, and repeating the steps 3 and 4.
The invention can obtain the injection quantity of the tracer and the arrangement method of the monitoring well. In the subsequent monitoring process, the actual concentration of the tracer obtained by the monitoring well is compared with the calculated result, and the approximate position of the leakage point can be obtained preliminarily.
The invention is further illustrated by the following examples. Step 1) physical model construction by tracer migration in landfill
As shown in fig. 2, the physical model of the tracer migration in a landfill is mainly composed of a cover ①, a refuse receptacle ②, a bottom pad ③, and an injection well ④. Wherein the thickness of the cover layer is assumed to be 1m, the thickness of the garbage body is assumed to be 2m, and the thickness of the bottom liner is assumed to be 1m. According to the materials used for the covering layer and the bottom liner of the typical landfill in China, clay is selected as a main covering material, and the soil porosity n is selected to be 0.46, so that the specific working conditions and parameters of the covering layer and the bottom liner take the values as shown in the following table 1:
TABLE 1 migration parameters of tracers in soil
The migration parameters of the tracers in the garbage body are as follows in table 2:
TABLE 2 migration parameters of tracers in refuse bodies
Assuming that the injection boundary of the injection well is a constant concentration injection, i.e., the concentration of the tracer in the injection well remains constant, the tracer enters the landfill system in the form of diffusion and convective migration due to concentration differences and pressure differences. The implantation boundaries are as follows:
For ease of field operations and calculations, the present example selects an initial concentration of injection well tracer set to C 0=10g/m3.
Step 2) model establishment and solution
Concentration profiles of the tracer in landfill cover, refuse body and bottom pad on the laplace domain, and release flux expressions of the tracer on the cover surface are obtained by formulas (1) - (10). Numerical inversion is then performed by equations (13) - (17).
Step 3.1) tracer migration concentration profile and Release flux Curve plotting
A cloud of concentration profiles of the tracer in the landfill system is shown in fig. 3. The working conditions are that the gas-containing porosities of the covering layer, the garbage body and the bottom liner are respectively 0.1,0.2 and 0.1. Table 1 and table 2 give the diffusion coefficient and convection rate of the tracer in the cover layer, the waste body and the bottom liner, respectively. Assuming that the detection limit of the SF 6 detector adopted by the invention is 1ppm, a curve with the concentration of 1ppm can be obtained according to the concentration profile cloud chart of fig. 3, and the radial coordinates of monitoring points at different depths can be determined at t=0.5 year. Monitoring points may be placed 0.5m below the overburden surface layer, i.e. z=1.5m, at a radial distance r=5m from the injection well.
Step 4) injection well injection amount determination and parameter impact analysis
Assuming a detection limit of 1ppm for the SF 6 detector employed in the present invention, the concentration profile plotted according to fig. 3 yields a tracer concentration at the monitoring point location of less than 1ppm at t=0.5 years. The method is characterized in that when the tracer can be detected at the monitoring point for t=0.5 years, if the concentration of the tracer in the injection well needs to be increased, the concentration of the tracer in the injection well is increased by 10 times, the tracer is substituted into the steps 2 and 3 to be calculated, and the calculation is repeated until the concentration of the tracer obtained at the monitoring point is larger than the detection limit, and finally the value of the concentration of the tracer injected into the injection well is the concentration of the tracer injected in the field test. The resulting tracer injection concentration ensures that the tracer is detected at the monitoring point location (concentration greater than the limit of detection) at t=0.5 years. In the subsequent monitoring process, for example, 0.5 year, the rough position of the leakage point can be obtained initially by detecting the tracer concentration data of each monitoring point. The graph of the tracer release flux from the surface of the coating under different conditions of air void content (water content) is shown in fig. 4. From the figure it is seen that the gas porosity of the cover layer has the greatest effect on the tracer release flux, while the gas porosity of the bottom liner has little effect on the tracer release flux. From fig. 4 it can be seen that when θ=0.2, θ l =0.1, t=100 day, the gas-containing porosity of the cover layer increases from 0.1 to 0.3, the tracer release flux increases by 5 orders of magnitude at r=6m. When the gas-containing porosity of the garbage body increases from 0.2 to 0.4 for θ u=0.3,θl = 0.1,100 days, the release flux of the tracer at r=6m increases by a factor of 2000. Thus, it is desirable to utilize the method of the present invention to monitor well placement and tracer injection boundary determination depending on the specifics of each landfill.
The above embodiments are illustrative of the present invention and not limiting, and any simple modifications of the present invention are within the scope of the present invention.

Claims (4)

1.一种用于填埋场渗漏探测的监测井布置和示踪剂注入边界确定方法,其特征在于,包括如下步骤:1. A method for arranging monitoring wells and determining tracer injection boundaries for landfill leakage detection, comprising the following steps: 步骤1):基于已封场填埋场的建设规模,确定垃圾体厚度为l1、覆盖层厚度为l2、底部衬垫厚为l3;现场测定填埋场覆盖层、垃圾体及底部衬垫的含水率,根据扩散系数计算模型(式1)和达西定律(式2)得到气相示踪剂在填埋场覆盖层、垃圾体及底部衬垫中的扩散系数和对流速率,Step 1): Based on the construction scale of the closed landfill, determine the thickness of the garbage body as l 1 , the thickness of the cover layer as l 2 , and the thickness of the bottom liner as l 3 ; measure the moisture content of the landfill cover layer, garbage body and bottom liner on site, and obtain the diffusion coefficient and convection rate of the gas phase tracer in the landfill cover layer, garbage body and bottom liner according to the diffusion coefficient calculation model (Formula 1) and Darcy's law (Formula 2). D=D0a 10/3/n2) (1)D=D 0a 10/3 /n 2 ) (1) 式中D和D0分别为示踪剂在土体和空气中的扩散系数,θa为土体的含气孔隙率,n为土体孔隙度,v为示踪剂在土体中的对流速率,ka为示踪剂在土体的渗透系数,dP/dz为土体两端的气压梯度;Where D and D0 are the diffusion coefficients of the tracer in the soil and air respectively, θa is the air porosity of the soil, n is the porosity of the soil, v is the convection velocity of the tracer in the soil, k a is the permeability coefficient of the tracer in the soil, and dP/dz is the pressure gradient at both ends of the soil; 步骤2.1):确定控制方程和边界条件Step 2.1): Determine the governing equations and boundary conditions 以垃圾体中间点为原点建立二维径向轴对称坐标系,建立示踪剂运移的二维控制方程,示踪剂注入条件为定浓度C0注入,确定模型所用边界条件;A two-dimensional radial axisymmetric coordinate system is established with the middle point of the garbage body as the origin, and the two-dimensional control equation of tracer migration is established. The tracer injection condition is a constant concentration C 0 injection, and the boundary conditions used in the model are determined; 步骤2.2):得到示踪剂在填埋体中的运移规律,具体包括:Step 2.2): Obtain the migration law of the tracer in the landfill, specifically including: 步骤2.2.1):计算拉普拉斯域上示踪剂运移浓度剖面;Step 2.2.1): Calculate the tracer migration concentration profile in the Laplace domain; 步骤2.2.2):计算拉普拉斯域上覆盖层表面示踪剂释放通量,计算拉普拉斯域上衬垫底部示踪剂浓度;Step 2.2.2): Calculate the tracer release flux on the overburden surface in the Laplace domain and calculate the tracer concentration on the bottom of the liner in the Laplace domain; 步骤2.2.3):利用Talbot数值反演得到示踪剂的浓度分布剖面和释放通量;Step 2.2.3): Use Talbot numerical inversion to obtain the tracer concentration distribution profile and release flux; 步骤3.1):计算示踪剂在填埋场系统中运移的浓度剖面和释放通量,并绘制其随时间和空间变化曲线;Step 3.1): Calculate the concentration profile and release flux of the tracer in the landfill system and plot its variation curve over time and space; 步骤3.2):根据步骤3.1)所绘制的浓度剖面及示踪剂的检出限确定监测点的位置坐标;Step 3.2): Determine the position coordinates of the monitoring point based on the concentration profile drawn in step 3.1) and the detection limit of the tracer; 步骤4):确定监测点示踪剂的释放通量及浓度,并将其与检出限进行对比,若监测点示踪剂浓度或释放通量大于检出限,则确定步骤1)的注入井注入浓度边界为实际注入浓度;若监测点示踪剂浓度或释放通量小于检出限,则监测点不能检测到示踪剂,增大注入井注入边界浓度或通量,重复步骤3.1)-步骤4)。Step 4): Determine the release flux and concentration of the tracer at the monitoring point, and compare them with the detection limit. If the tracer concentration or release flux at the monitoring point is greater than the detection limit, the injection concentration boundary of the injection well in step 1) is determined to be the actual injection concentration; if the tracer concentration or release flux at the monitoring point is less than the detection limit, the tracer cannot be detected at the monitoring point, increase the injection boundary concentration or flux of the injection well, and repeat steps 3.1)-step 4). 2.如权利要求1所述的一种用于填埋场渗漏探测的监测井布置和示踪剂注入边界确定方法,其特征在于,所述的步骤2.2.1),具体为:2. The method for arranging monitoring wells and determining tracer injection boundaries for landfill leakage detection according to claim 1, wherein the step 2.2.1) is specifically: 注入井定浓度边界注入,示踪剂在填埋体中运移浓度值为The injection well is injected at the fixed concentration boundary, and the tracer migration concentration in the landfill is 式中C0为注入井中示踪剂的浓度,p为拉普拉斯域上的时间因子,μn和ωn为特征根,r为径向坐标,rw为注入井的半径,z为竖向坐标,m2为示踪剂在垃圾体中的运移参数,可用下式进行计算:Where C0 is the concentration of the tracer in the injection well, p is the time factor on the Laplace domain, μn and ωn are characteristic roots, r is the radial coordinate, rw is the radius of the injection well, z is the vertical coordinate, and m2 is the migration parameter of the tracer in the garbage body, which can be calculated using the following formula: 式中v为示踪剂在垃圾体中的对流速率;θ为垃圾体的含气孔隙率;Dz为示踪剂在垃圾体中沿轴向的扩散系数;Dr为示踪剂在垃圾体中沿径向的扩散系数;Where v is the convection velocity of the tracer in the garbage body; θ is the air-containing porosity of the garbage body; Dz is the diffusion coefficient of the tracer in the garbage body along the axial direction; Dr is the diffusion coefficient of the tracer in the garbage body along the radial direction; 示踪剂在覆盖层中的浓度为The concentration of the tracer in the overburden is 式中,l1、l2为垃圾体厚度和覆盖层厚度,m3和m4为示踪剂在覆盖层中的运移参数,可分别由下式进行计算:Where l 1 and l 2 are the thickness of the garbage body and the thickness of the cover layer, m 3 and m 4 are the migration parameters of the tracer in the cover layer, which can be calculated by the following formulas: 式中vu为示踪剂在覆盖层中的对流速率;Du为示踪剂在覆盖层中沿轴向的扩散系数;Where v u is the convection velocity of the tracer in the covering layer; Du is the diffusion coefficient of the tracer in the covering layer along the axial direction; 示踪剂在底部衬垫中浓度剖面为The tracer concentration profile in the bottom liner is 式中m5和m6为示踪剂在底部衬垫中的运移参数,分别由下式进行求解:Where m5 and m6 are the migration parameters of the tracer in the bottom liner, which are solved by the following equations: 式中vl为示踪剂在底部衬垫中的对流速率;Dl为示踪剂在底部衬垫中沿轴向的扩散系数。Where v l is the convection velocity of the tracer in the bottom liner; D l is the diffusion coefficient of the tracer in the bottom liner along the axial direction. 3.如权利要求2所述的一种用于填埋场渗漏探测的监测井布置和示踪剂注入边界确定方法,其特征在于,所述的步骤2.2.2),具体为:3. A method for arranging monitoring wells and determining tracer injection boundaries for landfill leakage detection according to claim 2, characterized in that said step 2.2.2) is specifically: 计算拉普拉斯域上覆盖层表面示踪剂释放通量:Calculate the tracer release flux on the overburden surface in the Laplace domain: 计算拉普拉斯域上衬垫底部示踪剂浓度:Calculate the tracer concentration at the bottom of the pad on the Laplace domain: 式中l3为底部衬垫厚度。Where l 3 is the bottom pad thickness. 4.如权利要求3所述的一种用于填埋场渗漏探测的监测井布置和示踪剂注入边界确定方法,其特征在于,所述的步骤2.2.3),具体为:4. A method for arranging monitoring wells and determining tracer injection boundaries for landfill leakage detection according to claim 3, characterized in that said step 2.2.3) is specifically: 利用Talbot数值反演得到示踪剂的浓度分布剖面和释放通量,计算公式如下:The concentration distribution profile and release flux of the tracer are obtained by Talbot numerical inversion. The calculation formula is as follows: 式中M为叠加常数,Re为取复数的实部,γi,δi为拉普拉斯反演常数,通过下述公式定义:Where M is the superposition constant, Re is the real part of the complex number, γ i , δ i are Laplace inversion constants, which are defined by the following formula: 式中a为虚数因子。Where a is an imaginary factor.
CN202111407272.XA 2021-11-24 2021-11-24 A method for arranging monitoring wells and determining tracer injection boundaries for landfill leakage detection Active CN114254484B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111407272.XA CN114254484B (en) 2021-11-24 2021-11-24 A method for arranging monitoring wells and determining tracer injection boundaries for landfill leakage detection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111407272.XA CN114254484B (en) 2021-11-24 2021-11-24 A method for arranging monitoring wells and determining tracer injection boundaries for landfill leakage detection

Publications (2)

Publication Number Publication Date
CN114254484A CN114254484A (en) 2022-03-29
CN114254484B true CN114254484B (en) 2025-02-11

Family

ID=80791129

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111407272.XA Active CN114254484B (en) 2021-11-24 2021-11-24 A method for arranging monitoring wells and determining tracer injection boundaries for landfill leakage detection

Country Status (1)

Country Link
CN (1) CN114254484B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114912158B (en) * 2022-05-13 2023-02-14 华中科技大学 Design method of soil covering layer of solid waste landfill with consideration of plant root system effect
CN115828528B (en) * 2022-11-09 2025-01-21 浙江大学 A landfill leakage channel detection method based on tracer migration
CN117911639B (en) * 2024-03-19 2024-06-04 江苏华东八一四地球物理勘查有限公司(江苏省有色金属华东地质勘查局八一四队) Three-dimensional geological modeling method based on well-electricity combined exploration landfill environment investigation

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105971026A (en) * 2016-05-31 2016-09-28 河海大学 Method for forecasting service life of landfill anti-seepage system through indication pollutants

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8949096B2 (en) * 2012-05-09 2015-02-03 The United States Of America, As Represented By The Secretary Of The Navy Three-dimensional tracer dispersion model

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105971026A (en) * 2016-05-31 2016-09-28 河海大学 Method for forecasting service life of landfill anti-seepage system through indication pollutants

Also Published As

Publication number Publication date
CN114254484A (en) 2022-03-29

Similar Documents

Publication Publication Date Title
CN114254484B (en) A method for arranging monitoring wells and determining tracer injection boundaries for landfill leakage detection
US6035701A (en) Method and system to locate leaks in subsurface containment structures using tracer gases
WO2016074510A1 (en) Method for determining earth surface interpenetrated crack distribution and air leakage characteristics in shallow burial coal mining
CN103207212B (en) Method for simulation of pollutant vertical isolation barrier
Cho et al. Estimating natural gas emissions from underground pipelines using surface concentration measurements☆
CN104865165B (en) The on-site soil solute transfer dispersion coefficient method for measuring of overall process flux conservation
US6611760B2 (en) Method and system for estimating gas production by a landfill or other subsurface source
CN111872095B (en) Flexible bushing and method for monitoring contaminated sites combined with flexible bushing
Oldenburg et al. Modeling gas transport in the shallow subsurface during the ZERT CO 2 release test
Xie et al. Analytical model for methane migration through fractured unsaturated landfill cover soil
Morales et al. Solute transport characterization in karst aquifers by tracer injection tests for a sustainable water resource management
CN108732068B (en) Soil gas invasion simulation device and indoor air TCA concentration evaluation method thereof
CN105628892B (en) For studying the model equipment that dusty gas is migrated when contaminated site pile foundation is set
Hou et al. Leakage locating and sampling optimization of small-rate leakage on medium-and-low-pressure unground natural gas pipelines
CN203275342U (en) Device for simulating vertical isolation barrier of pollutant
Jakubick et al. Vacuum testing of the permeability of the excavation damaged zone
Yazdani et al. Quantifying capture efficiency of gas collection wells with gas tracers
Joun et al. Reproducing natural variations in CO2 concentration in vadose zone wells with observed atmospheric pressure and groundwater data
CN206573100U (en) Settling column test instrument for pore pressure distribution measuring
Collignan et al. Methodology for the in situ characterisation of soil vapor contaminants and their impact on the indoor air quality of buildings
Xie et al. A two-dimensional semi-analytical model for tracer gas injection and transport in a landfill system
Wang et al. Leak detection method for bottom plate of oil tank based on oil/gas leak detection
Jayarathne et al. Understanding mid-to large underground leaks from buried pipelines as affected by soil and atmospheric conditions–Field scale experimental study
Sun et al. Analysis of pressure response at an observation well against pressure build-up by early stage of CO2 geological storage project
Esposito* et al. Migration of natural gas through heterogeneous sandy soils affected by atmospheric boundary conditions

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant