CN106370802A - Method for determining water loss of reservoir by using hydrogen and oxygen isotopes - Google Patents

Method for determining water loss of reservoir by using hydrogen and oxygen isotopes Download PDF

Info

Publication number
CN106370802A
CN106370802A CN201610673761.2A CN201610673761A CN106370802A CN 106370802 A CN106370802 A CN 106370802A CN 201610673761 A CN201610673761 A CN 201610673761A CN 106370802 A CN106370802 A CN 106370802A
Authority
CN
China
Prior art keywords
reservoir
water
isotope
hydrogen
oxygen
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201610673761.2A
Other languages
Chinese (zh)
Inventor
胡海英
黄国如
余水灵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN201610673761.2A priority Critical patent/CN106370802A/en
Publication of CN106370802A publication Critical patent/CN106370802A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/18—Water

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

The invention discloses a method for determining the water loss of a reservoir by using hydrogen and oxygen isotopes, belonging to the field of regional water balance and water resource management. The method comprises the following steps: (1) collecting hydro-meteorological, hydraulic engineering and other basic data of the reservoir; (2) investigating and analyzing the items of water bodies of a variety of water sources flowing in and out of the reservoir; (3) carrying out field monitoring experiments and acquiring water samples of the variety of water sources in the reservoir; (4) determining the composition of oxygen and hydrogen isotopes in every incoming and outgoing items of the reservoir in a laboratory; (5) establishing a calculation model; and (6) determining the water loss of the reservoir. A stable-isotope method is a useful tool in modern hydrological research; isotope tracer technology is introduced into the calculation of water balance in the reservoir, so the water loss of the reservoir can be effectively analyzed. In addition, the method provided by the invention has few measuring factors, and is simple in calculation method and reasonable in results; and the method is worth reference when applied as a novel method.

Description

A kind of method that application Hydrogen-oxygen Isotope determines reservoir water loss
Technical field
The present invention relates to regional water balance, the research field of water resources management, particularly to a kind of application same position of hydrogen-oxygen The method that element determines reservoir water loss.
Background technology
Hydrogen-oxygen Isotope is the important component part in natural water, its shared ratio very little, but very sensitively rings Answer the change of surrounding.Hydrogen-oxygen Isotope can spike water circulation effectively, inquire into many water resources, water environmental problems.? In isotope hydrology circulating research, the report of existing many researchs in terms of isotope.Surface water body (such as river, lake, water Storehouse) to be dried out by evaporation, during Water Evaporation, the Hydrogen-oxygen Isotope content of composition hydrone can occur small change, this Plant to change and caused due to isotope fractionation.In hydrologic cycle, there are many processes can cause stable isotope fractionation in water body As evaporation, condense, infiltration etc., typically based on the mixed influence at different water sources.Stable isotope in water body, is equivalent to human body " dna ", the water body (atmospheric vapour, sea water, river, subsoil water, lake etc.) of separate sources, its isotopics is different.According to The difference of separate sources water body stable isotope can be explained and analyze many phenomenons, and in water body water balanced calculation.
After reservoir dam water-retention, dam body, the dam foundation and dam abutment and all easily occur oozing to adjacent paddy in the presence of Ku Shui Leakage phenomenon, the In Limestone Area particularly developed in karst cave, reservoir leakage meeting ratio is more serious, to dam safety and reservoir Benefit is all unfavorable.For water reservoir management, set up the water yield prediction of reservoir, within certain period, quantify reservoir loss amount over the ground Area's Water Resource Balance and management are all significant.In carrying out reservoir water yield balance studies, only apply conventional hydrology side Method estimation reservoir water loss is generally relatively difficult or even is difficult to solve.
Content of the invention
Present invention is primarily targeted at overcoming shortcoming and the deficiency of prior art, provide a kind of application Hydrogen-oxygen Isotope true The method determining reservoir water loss, using isotope tracer technique, different according to the water body Hydrogen-oxygen Isotope content of separate sources, By measuring reservoir each revenue and expenditure item water body Hydrogen-oxygen Isotope composition, set up the reservoir water yield and isotopic mass conservation model, inquire into Reservoir evaporation amount, leakage etc. lose item.
In order to achieve the above object, the present invention employs the following technical solutions:
The method that a kind of present invention application Hydrogen-oxygen Isotope determines reservoir water loss, comprises the steps:
(1) collect reservoir hydrologic meteorology and physical geography data, specifically include reservoir watershed hydrometeorology data, water conservancy work Journey basic data；
(2) according to reservoir watershed hydrometeorology and hydraulic engineering basic data, analyze and determine and flow in and out reservoir Each source water item number；
(3) carry out field monitoring experiment, gather reservoir each revenue and expenditure item water sample, specifically include rainfall, river, subsoil water and water Storehouse water body etc.；
(4) measure reservoir each revenue and expenditure item stable isotope data, measure each water source using lgr aqueous water isotope analyses instrument Water body Hydrogen-oxygen Isotope concentration, measuring accuracy is respectively within 0.3 ‰ and 0.1 ‰；
(5) set up computation model；
(6) inquire into reservoir water loss item.
As preferred technical scheme, in step (3), during sample collecting and storage, to strictly avoid generation same The plain fractional distillation in position.
It is as preferred technical scheme, in step (5), different according to reservoir separate sources water body Hydrogen-oxygen Isotope content, Set up the reservoir water yield and stable isotope Mass balance model:
d v d t = q i s + q i g + p - q o l - q o x - e - - - ( 1 )
d ( c r v ) d t = c i s q i s + c i g q i g + c p p - c o l q o l - c o x q o x - c e e - - - ( 2 )
In formula, v is reservoir water body volume, qis、qig、qol、qoxRepresent reservoir surface water inbound traffics, the becoming a mandarin of subsoil water respectively Amount, leakage and vent flow, p is the rainfall in the range of reservoir, and e is reservoir evaporation amount, cr、cis、cig、cp、col、cox、ce After respectively storehouse water body, reservoir surface water inbound traffics, subsoil water inbound traffics, rainfall, leakage, vent flow, reservoir water body evaporate The concentration of Hydrogen-oxygen Isotope.
As preferred technical scheme, in step (6), substitute into reservoir each revenue and expenditure item stable isotope data, solution procedure (5) the reservoir water yield set up and stable isotope Mass balance model, determine reservoir water loss amount.
As preferred technical scheme, in step (6), within certain period, if water reservoir system hydrographic features or same position Plain concentration composition is in steady statue, then have dcr/ dt=0, dv/dt=0, then equation (1) and (2) be converted into:
qis+qig+p-qol-qox- e=0 (3)
cisqis+cigqig+cpp-colqol-coxqox-ceE=0 (4)
Joint equation (3) and equation (4) can solve any two variable, in the case of reservoir is sufficiently mixed, equation (4) Middle reservoir leakage amount and the isotopes concentration of vent flow are equal i.e. c with the isotopes concentration of reservoir water bodyol=cr, cox= cr, when other variables it is known that reservoir evaporation amount e and reservoir leakage amount q are obtained by below equationol,
e = q i s ( c i s - c r c e - c r ) + p ( c p - c r c e - c r ) + q i g ( c i g - c r c e - c r ) - - - ( 5 )
qol=qis+qig+p-e-qox(6)
Isotope composition c as preferred technical scheme, in step (6), in the net evaporation of described reservoir water bodyeCan by under Formula is given:
c e = αc r - c a h - ϵ 1 - h + δ ϵ - - - ( 7 )
In formula: h is the relative humidity of air；caIt is the isotopic content of free air above reservoir；δ ε, ε are respectively dynamic Power fractionation factor and total fractionation factor, wherein ε=1- α+δ ε, α is the Equilibrium fractionation factor between air and the water surface, and δ ε is power Fractionation factor, δ ε=k (1-h), k value is 14.3 ‰；
When temperature changes from 273.15k to 373.15k, fractionation factor α can be described with following formula with the relation of temperature t:
Right18O:
To d:
The present invention compared with prior art, has the advantage that and beneficial effect:
1st, the present invention inquires into reservoir water loss using Steady Isotope, can not directly measure the water of continuous item Amount, and only need to measure the content of its Hydrogen-oxygen Isotope, calculated by reservoir water balance and stable isotope mass balance equation.
2nd, this method is applied to the natural lake that natural depression is formed, and is also applied for lake and the reservoir of hand excavation.Surely Determine the useful tool that isotope method is modern hydrology research, isotope technology is introduced the meter of reservoir water balance Calculate, as a kind of application of new method, be worth using for reference.
Brief description
Fig. 1 is method of the present invention flow chart.
Specific embodiment
With reference to embodiment and accompanying drawing, the present invention is described in further detail, but embodiments of the present invention do not limit In this.
Embodiment
As shown in figure 1, the method that a kind of the present embodiment application Hydrogen-oxygen Isotope determines reservoir water loss, walk including following Rapid:
(1) collect reservoir hydrologic meteorology and physical geography data, specifically include reservoir watershed hydrometeorology data, water conservancy work Journey basic data；
(2) according to reservoir watershed hydrometeorology and hydraulic engineering basic data, analyze and determine and flow in and out reservoir Each source water item number；
(3) carry out field monitoring experiment, gather reservoir each revenue and expenditure item water sample, specifically include rainfall, river, subsoil water and water Storehouse water body etc.；
(4) measure reservoir each revenue and expenditure item stable isotope data；
(5) set up computation model；
(6) inquire into reservoir water loss amount.
The present embodiment is selected certain lake and is located at Northeastern Tibetan Plateau occasionally, is in east China monsoon region, northwestward arid biogeographic zone With the area that crosses of west and south High-cold regions, lake region mean annual temperature be 1.2 DEG C, relative humidity be 55%, Lake level change is about 4400km2.The stable isotope of lake water varies less with depth and horizontal range, myopia can regard as and be sufficiently mixed.Open first Exhibition field monitoring experiment, collection flows in and out the Different Waters water sample in lake.Lake Different Waters oxygen isotope concentration is average Value is respectively as follows: rainfall (- 15.19 ‰), subsoil water (- 8.19 ‰), river (- 7.39 ‰), lake water (1.97 ‰), lake surface air (- 21.88‰).
This lake is the inland lake of a closing, and evaporation from water surface is unique expenditure item of water loss, does not contemplate stream Impact, set up the water yield and the oxygen isotope concentration mass balance equation in lake:
d v d t = q i s + q i g + p - e - - - ( 10 )
d ( c r v ) d t = c i s q i s + c i g q i g + c p p - c e e - - - ( 11 )
In formula, v is lake water volume, qisAnd qigRepresent the inbound traffics of lake surface water inbound traffics and subsoil water respectively, p is Rainfall in the range of lake, e is lake evaporation capacity, cr、cis、cig、cp、ceIt is respectively water body in lake, lake surface water becomes a mandarin Amount, subsoil water inbound traffics, rainfall, water body in lake evaporation after (net evaporation capacity) oxygen isotope concentration.
Connection solve equation (10) and (11) formula obtain lake evaporate item:
e = ( c r - c i g ) d v d t + ( c i g - c p ) p + ( c i g - c i s ) i s c i g - c e - - - ( 12 )
Known lake mean annual temperature t=1.2 DEG C, the relative humidity h=55% of air, lake surface air oxygen isotope is dense Degree ca=-21.88 ‰, lake water oxygen isotope concentration cr=1.97 ‰.Substitute into formula (7) and (8), can try to achieve in the net evaporation of lake water Isotopes concentration ce.By lake Different Waters (rainfall, subsoil water, river, lake water) oxygen isotope concentration value, lake annual Rainfall, surface water inbound traffics and year lake water water yield variable quantity substitute into formula (12), can evaporate in the hope of the many annuals in this lake Measure as 40.9 hundred million m3, and many annual mean ground temperature are surveyed for 40.5 hundred million m in lake3, with Steady Isotope, finally counted year puts down Compared with measured value, absolute error and relative error are respectively 0.4 hundred million m to all evaporation capacity3With 1%, therefore apply Hydrogen-oxygen Isotope Method estimation Lake water quantity loss item is feasible.
Above-described embodiment is the present invention preferably embodiment, but embodiments of the present invention are not subject to above-described embodiment Limit, other any spirit without departing from the present invention and the change made under principle, modification, replacement, combine, simplify, All should be equivalent substitute mode, be included within protection scope of the present invention.

Claims (5)

1. a kind of application Hydrogen-oxygen Isotope determines the method for reservoir water loss it is characterised in that comprising the steps:
(1) collect reservoir hydrologic meteorology and physical geography data, specifically include reservoir watershed hydrometeorology data, hydraulic engineering base Plinth data；
(2) according to reservoir watershed hydrometeorology and hydraulic engineering basic data, analyze and determine each water flowing in and out reservoir Source water body item number；
(3) carry out field monitoring experiment, gather reservoir each revenue and expenditure item water sample, specifically include rainfall, river, subsoil water and reservoir water Body；
(4) measure reservoir each revenue and expenditure item stable isotope data, measure each source water using lgr aqueous water isotope analyses instrument Hydrogen-oxygen Isotope concentration, measuring accuracy is respectively within 0.3 ‰ and 0.1 ‰；
(5) set up computation model；
(6) inquire into reservoir water loss.
2. according to claim 1 a kind of application Hydrogen-oxygen Isotope determine reservoir water loss method it is characterised in that In step (3), during sample collecting and storage, strictly avoid isotope fractionation.
3. according to claim 1 a kind of application Hydrogen-oxygen Isotope determine reservoir water loss method it is characterised in that In step (5), different according to reservoir separate sources water body Hydrogen-oxygen Isotope content, set up the reservoir water yield and stable isotope quality Balance model:
In formula, v is reservoir water body volume, qis、qig、qol、qoxRepresent respectively reservoir surface water inbound traffics, the inbound traffics of subsoil water, Leakage and vent flow, p is the rainfall in the range of reservoir, and e is reservoir evaporation amount, cr、cis、cig、cp、col、cox、ceRespectively For hydrogen-oxygen after storehouse water body, reservoir surface water inbound traffics, subsoil water inbound traffics, rainfall, leakage, vent flow, reservoir water body evaporation Isotopic concentration.
4. according to claim 3 a kind of application Hydrogen-oxygen Isotope determine reservoir water loss method it is characterised in that In step (6), substitute into reservoir each revenue and expenditure item stable isotope data, the reservoir water yield that solution procedure (5) is set up is same with stable The plain Mass balance model in position, inquires into reservoir water loss item.
5. according to claim 4 a kind of application Hydrogen-oxygen Isotope determine reservoir water loss method it is characterised in that In step (6), within certain period, if water reservoir system hydrographic features or isotopes concentration composition are in steady statue, have dcr/ dt=0, dv/dt=0, then equation (1) and (2) be converted into:
qis+qig+p-qol-qox- e=0 (3)
cisqis+cigqig+cpp-colqol-coxqox-ceE=0 (4)
Joint equation (3) and equation (4) can solve any two variable, in the case of reservoir is sufficiently mixed, water in equation (4) The isotopes concentration of storehouse leakage and vent flow is equal with the isotopes concentration of reservoir water body, i.e. col=cr, cox=cr, when Other variables are it is known that obtain reservoir evaporation amount e and reservoir leakage amount q by below equationol,
qol=qis+qig+p-e-qox(6) .
CN201610673761.2A 2016-08-16 2016-08-16 Method for determining water loss of reservoir by using hydrogen and oxygen isotopes Pending CN106370802A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610673761.2A CN106370802A (en) 2016-08-16 2016-08-16 Method for determining water loss of reservoir by using hydrogen and oxygen isotopes

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610673761.2A CN106370802A (en) 2016-08-16 2016-08-16 Method for determining water loss of reservoir by using hydrogen and oxygen isotopes

Publications (1)

Publication Number Publication Date
CN106370802A true CN106370802A (en) 2017-02-01

Family

ID=57878691

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610673761.2A Pending CN106370802A (en) 2016-08-16 2016-08-16 Method for determining water loss of reservoir by using hydrogen and oxygen isotopes

Country Status (1)

Country Link
CN (1) CN106370802A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107991377A (en) * 2017-11-15 2018-05-04 东华理工大学 For analyzing the isotope extracting method of River Basin Hydrology forming process
CN109598082A (en) * 2018-12-17 2019-04-09 河海大学 The calculation method of lake evaporation capacity and crucial hydrographic information based on Hydrogen-oxygen Isotope
CN109767027A (en) * 2018-11-30 2019-05-17 深圳春沐源控股有限公司 Reservoir water percolating capacity prediction technique and device
CN111504277A (en) * 2020-03-20 2020-08-07 河海大学 Lake water supply tracing method utilizing hydrogen and oxygen isotopes
CN111539090A (en) * 2020-03-20 2020-08-14 河海大学 Watershed evapotranspiration estimation method combining water body hydrogen and oxygen isotopes and chloride ions

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105021732A (en) * 2014-11-03 2015-11-04 中国食品发酵工业研究院 Method for fast determination of oxygen isotope composition of water in beverage wine
CN105651913A (en) * 2016-02-29 2016-06-08 浙江大学 Device and method for determining amino acid adsorption capacity of soil

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105021732A (en) * 2014-11-03 2015-11-04 中国食品发酵工业研究院 Method for fast determination of oxygen isotope composition of water in beverage wine
CN105651913A (en) * 2016-02-29 2016-06-08 浙江大学 Device and method for determining amino acid adsorption capacity of soil

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
包为民等: "同位素技术在湖泊水文学中的应用研究", 《水电能源科学》 *
包为民等: "稳定同位素方法在湖泊水量平衡研究中的应用", 《人民黄河》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107991377A (en) * 2017-11-15 2018-05-04 东华理工大学 For analyzing the isotope extracting method of River Basin Hydrology forming process
CN109767027A (en) * 2018-11-30 2019-05-17 深圳春沐源控股有限公司 Reservoir water percolating capacity prediction technique and device
CN109598082A (en) * 2018-12-17 2019-04-09 河海大学 The calculation method of lake evaporation capacity and crucial hydrographic information based on Hydrogen-oxygen Isotope
CN111504277A (en) * 2020-03-20 2020-08-07 河海大学 Lake water supply tracing method utilizing hydrogen and oxygen isotopes
CN111539090A (en) * 2020-03-20 2020-08-14 河海大学 Watershed evapotranspiration estimation method combining water body hydrogen and oxygen isotopes and chloride ions
CN111504277B (en) * 2020-03-20 2021-05-04 河海大学 Lake water supply tracing method utilizing hydrogen and oxygen isotopes
CN111539090B (en) * 2020-03-20 2023-04-07 河海大学 Watershed evapotranspiration estimation method combining water body hydrogen and oxygen isotopes and chloride ions

Similar Documents

Publication Publication Date Title
Xiao et al. Eutrophic Lake Taihu as a significant CO2 source during 2000–2015
Yi et al. A coupled isotope tracer method to characterize input water to lakes
CN106370802A (en) Method for determining water loss of reservoir by using hydrogen and oxygen isotopes
Qian et al. Isotopic characteristics of precipitation, surface and ground waters in the Yinchuan plain, Northwest China
Yakir et al. The use of stable isotopes to study ecosystem gas exchange
Fabre et al. Assessment of sediment and organic carbon exports into the Arctic ocean: The case of the Yenisei River basin
Kumar et al. Tracing isotopic signatures (δD and δ18O) in precipitation and glacier melt over Chorabari Glacier–Hydroclimatic inferences for the Upper Ganga Basin (UGB), Garhwal Himalaya
Sun et al. Stable isotopes of atmospheric precipitation and its environmental drivers in the Eastern Chinese Loess Plateau, China
Wu et al. Spatial distribution and controlling factors of surface water stable isotope values (δ18O and δ2H) across Kazakhstan, Central Asia
Yang et al. Isotope analysis of the nutrient supply in Xiangxi Bay of the Three Gorges Reservoir
Qiu et al. Determining the origin, circulation path and residence time of geothermal groundwater using multiple isotopic techniques in the Heyuan Fault Zone of Southern China
Gibson et al. Regional trends in evaporation loss and water yield based on stable isotope mass balance of lakes: The Ontario Precambrian Shield surveys
Liu et al. Combined use of tracer approach and numerical simulation to estimate groundwater recharge in an alluvial aquifer system: A case study of Nasunogahara area, central Japan
Timsic et al. Spatial variability in stable isotope values of surface waters of Eastern Canada and New England
Ozyurt et al. Characterization of the Gacka River basin karst aquifer (Croatia): Hydrochemistry, stable isotopes and tritium-based mean residence times
Miller et al. Methane ebullition in temperate hydropower reservoirs and implications for US policy on greenhouse gas emissions
Looman et al. The carbon dioxide evasion cycle of an intermittent first-order stream: contrasting water–air and soil–air exchange
Liu et al. Nutrient transports in the Baltic Sea–results from a 30-year physical–biogeochemical reanalysis
Bai et al. Evaluation of evapotranspiration for exorheic basins in China using an improved estimate of terrestrial water storage change
CN111504277B (en) Lake water supply tracing method utilizing hydrogen and oxygen isotopes
Wang et al. Impact of groundwater extraction on hydrological process over the Beijing-Tianjin-Hebei region, China
Su et al. Using 222Rn temporal and spatial distributions to estimate the groundwater discharge rate and associated nutrient fluxes into high salinity lakes in Badain Jaran Desert, Northwest China
Jurikova et al. Variations in triple isotope composition of dissolved oxygen and primary production in a subtropical reservoir
Wang et al. Quantifying moisture recycling of a leeward oasis in arid central Asia using a Bayesian isotopic mixing model
Cai et al. Stable water isotope and surface heat flux simulation using ISOLSM: Evaluation against in-situ measurements

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20170201