CN115407049B - Method for measuring and calculating carbon sink potential of soil in water and soil conservation engineering measures - Google Patents

Method for measuring and calculating carbon sink potential of soil in water and soil conservation engineering measures Download PDF

Info

Publication number
CN115407049B
CN115407049B CN202211070781.2A CN202211070781A CN115407049B CN 115407049 B CN115407049 B CN 115407049B CN 202211070781 A CN202211070781 A CN 202211070781A CN 115407049 B CN115407049 B CN 115407049B
Authority
CN
China
Prior art keywords
sediment
soil
engineering
type
area
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
CN202211070781.2A
Other languages
Chinese (zh)
Other versions
CN115407049A (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.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
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 China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN202211070781.2A priority Critical patent/CN115407049B/en
Publication of CN115407049A publication Critical patent/CN115407049A/en
Application granted granted Critical
Publication of CN115407049B publication Critical patent/CN115407049B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00Sampling; Preparing specimens for investigation
    • G01N1/02Devices for withdrawing samples
    • G01N1/04Devices for withdrawing samples in the solid state, e.g. by cutting
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00Sampling; Preparing specimens for investigation
    • G01N1/02Devices for withdrawing samples
    • G01N1/10Devices for withdrawing samples in the liquid or fluent state
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N5/00Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
    • 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00Sampling; Preparing specimens for investigation
    • G01N1/02Devices for withdrawing samples
    • G01N1/10Devices for withdrawing samples in the liquid or fluent state
    • G01N2001/1006Dispersed solids
    • G01N2001/1012Suspensions
    • G01N2001/1025Liquid suspensions; Slurries; Mud; Sludge

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Geology (AREA)
  • Mathematical Analysis (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Mathematical Optimization (AREA)
  • Environmental & Geological Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Hydrology & Water Resources (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention discloses a method for measuring and calculating carbon sink potential of soil in water and soil conservation engineering measures, which comprises the steps of dividing the water and soil conservation engineering measures into a point type and a surface type; for spot-type tooling Cheng Cuoshi: the sediment accumulation and carbon fixation amount of the spot type tool Cheng Cuoshi is finally calculated by measuring the maximum storage capacity, the layered sediment accumulation of the blocked sediment amount, calculating sediment accumulation potential of the part of the spot type engineering which is not full, and the like; for face-type tooling Cheng Cuoshi: the carbon fixation amount of the surface engineering measures is finally calculated by calculating parameters such as the service life of the engineering measures, the area of a control area of the surface engineering measures, the soil erosion modulus of the surface engineering measures, and the like. And finally, calculating the total carbon sequestration amount of the engineering measures in the area by combining the results of the point-type engineering measures and the surface-type engineering measures. The invention scientifically and comprehensively considers the category of water and soil conservation engineering measures and the loss of carbon in the sediment transportation process, and the measuring and calculating result is more true and accurate.

Description

一种测算水土保持工程措施土壤碳汇潜力的方法A method for measuring soil carbon sink potential of soil and water conservation engineering measures

技术领域technical field

本发明属于水土保持碳汇技术领域,尤其涉及水土保持工程措施土壤碳汇测算的方法。The invention belongs to the technical field of soil and water conservation carbon sinks, in particular to a method for measuring and calculating soil carbon sinks for water and soil conservation engineering measures.

背景技术Background technique

水土保持工程措施可以有效减少水土流失,固定和增加土壤有机碳含量,是土壤碳汇的重要途径。李勇等(2003)选取典型黄土高原丘陵沟壑区进行研究,以位于陕西省延安市宝塔区东北部的碾庄沟流域12个淤地坝为研究对象,结果表明:1957~2000年碾庄沟流域12个淤地坝共储存有机碳1.73×105t,流域碳储存强度提高了0.13~5.03t·hm2·a-1,平均为1.28t·hm2·a-1。符素华等(2009)发现鱼鳞坑、水平条等水土保持工程措施的多年平均减沙效益超过90%,同时还有研究发现沉沙池、小型水库的修建可以强化泥沙淤积掩埋,形成不易矿化的深层土壤碳从而实现碳埋藏,增加土壤碳汇。淤地坝是水土流失沟道治理工程的骨干力量,设置淤地坝不仅可以提高土壤有机碳含量,同时拦泥淤地,使荒沟变成了人造小平原,增加耕地面积,大大提高地表植物生产力。袁立敏等(2014)研究发现沙障可以减少风蚀,结合种植植被来增加土壤有机质含量,起到固碳作用。岳遥等(2016)计算了小流域尺度土壤侵蚀/淤积引起的碳通量特征,同时计算了全国尺度由侵蚀引起的土壤有机质和CO2通量变化。Soil and water conservation engineering measures can effectively reduce soil erosion, fix and increase soil organic carbon content, and are an important way for soil carbon sinks. Li Yong et al. (2003) selected typical hilly and gully areas of the Loess Plateau for research, and took 12 check dams in the Nianzhuanggou watershed located in the northeast of Baota District, Yan'an City, Shaanxi Province as the research object. The results showed that: 1957~2000 The 12 check dams in the watershed have stored a total of 1.73×105 t of organic carbon, and the carbon storage intensity in the watershed has increased by 0.13-5.03 t·hm 2 ·a -1 , with an average of 1.28 t·hm 2 ·a -1 . Fu Suhua et al. (2009) found that the multi-year average sediment reduction benefits of water and soil conservation engineering measures such as fish scale pits and horizontal strips exceeded 90%. At the same time, other studies found that the construction of sedimentation tanks and small reservoirs can strengthen sedimentation and burial, forming a mineralized environment that is difficult to mineralize. The deep soil carbon can realize carbon burial and increase soil carbon sink. The check dam is the backbone of the water and soil erosion ditch control project. The installation of the check dam can not only increase the organic carbon content of the soil, but also block the mud and silt, turning the barren ditch into a small artificial plain, increasing the area of cultivated land, and greatly improving the surface vegetation. productive forces. Yuan Limin et al. (2014) found that sand barriers can reduce wind erosion, combined with planting vegetation to increase soil organic matter content and play a role in carbon sequestration. Yue et al. (2016) calculated the carbon flux characteristics caused by soil erosion/siltation at the small watershed scale, and at the same time calculated the changes in soil organic matter and CO 2 fluxes caused by erosion at the national scale.

目前计算土壤侵蚀对土壤碳影响时未考虑水土保持工程措施对土壤碳的作用,而工程措施的研究一般为实测土壤碳的变化且为单一工程措施的研究,均未对水土保持工程措施进行分类研究,目前尚无水土保持工程措施土壤碳汇的计算方法。At present, when calculating the impact of soil erosion on soil carbon, the effect of soil and water conservation engineering measures on soil carbon is not considered, and the research on engineering measures is generally based on the actual measurement of soil carbon changes and is a study of a single engineering measure, and none of the soil and water conservation engineering measures are classified. However, there is currently no calculation method for soil carbon sinks in soil and water conservation engineering measures.

发明内容Contents of the invention

针对现有对水土保持工程措施碳汇计算的需求,本发明的目的是提供一种普遍适用的针对不同水土保持工程措施的土壤碳汇潜力测算方法。In view of the existing demand for the calculation of carbon sinks of water and soil conservation engineering measures, the purpose of the present invention is to provide a universally applicable method for measuring and calculating soil carbon sink potentials for different water and soil conservation engineering measures.

本发明的目的是通过以下技术方案实现的:The purpose of the present invention is achieved by the following technical solutions:

一种测算水土保持工程措施土壤碳汇潜力的方法,包括以下步骤:A method for calculating soil carbon sink potential of soil and water conservation engineering measures, comprising the following steps:

步骤一,将所采用水土保持工程措施进行分类:分为面型措施和点型措施;所述面型措施为对面尺度土壤侵蚀控制作用的工程措施;所述点型措施为对于已侵蚀的泥沙进行拦截存储作用的措施;Step 1, classifying the adopted soil and water conservation engineering measures: divided into surface measures and point measures; the surface measures are engineering measures for controlling soil erosion on the surface scale; Measures for the interception and storage of sand;

步骤二,确定点型工程措施的最大库容Q最大和已拦蓄泥沙量的分层淤积量,分层淤积量的总和为点型工程措施已拦蓄泥沙量QStep 2, determine the maximum storage capacity Q of the point-type engineering measures and the layered deposition amount of the deposited sediment amount, the sum of the layered deposition amount is the point-type engineering measure and the accumulated sediment amount Q ;

步骤三,计算点型工程未淤满部分泥沙拦蓄潜力:采用点型工程措施的最大库容减去点型工程措施已拦蓄泥沙量,即:Q最大-QStep 3: Calculating the sediment storage potential of the unsilted part of the point-type project: using the maximum storage capacity of the point-type project measures minus the amount of sediment stored by the point-type project measures, that is: Q max - Q already ;

步骤四,计算点型工程达到最大库容后可发挥的其他减蚀作用的系数α;Step 4, calculate the coefficient α of other corrosion reduction effects that can be exerted after the point-type project reaches the maximum storage capacity;

步骤五,测定点型工程淤积层泥沙分层有机碳含量和侵蚀区土壤有机碳含量;Step 5, measuring the organic carbon content of the silt layer in the point-type engineering sedimentation layer and the soil organic carbon content in the erosion area;

步骤六,计算点型工程措施拦蓄泥沙固碳量:采用公式(1)~(5)进行计算:Step 6: Calculating the amount of sediment and carbon fixation in point-type engineering measures: use formulas (1) to (5) to calculate:

C=C点已+C点未           (1)C point = C point already + C point not (1)

Figure GDA0003887022500000021
Figure GDA0003887022500000021

Figure GDA0003887022500000022
Figure GDA0003887022500000022

C点未=Q潜力G侵蚀/1000                   (4) Point C = Q Potential G Erosion /1000 (4)

Q潜力=(Q最大-Q)+Q最大α                   (5)Q Potential = ( QMax - QHave ) + QMax α (5)

式中:C为点型工程措施拦蓄泥沙固碳量,单位为t;C点已为点型工程措施已淤积泥沙固碳量,单位为t;C点未为点型工程措施未来固碳潜力,单位为t;Qi为淤积泥沙第i层淤积量,单位为t;Gi为淤积泥沙第i层有机碳含量,单位为kg/t;Q为点型工程措施已拦蓄泥沙量,单位为t;Q潜力为点型工程泥沙拦蓄总潜力,单位为t;G侵蚀为点型工程措施侵蚀区土壤有机碳含量,单位为kg/t;Q最大为点型工程措施的最大库容,单位为t;α为点型工程达到最大库容后可发挥的其他减蚀作用的系数,取值范围0~1,无单位;In the formula: point C is the amount of sediment carbon sequestered by point-type engineering measures, and the unit is t; point C is the amount of deposited sediment carbon-fixed carbon by point-type engineering measures, and the unit is t; point C is the future of point-type engineering measures Carbon sequestration potential, unit is t; Q i is deposition amount of sediment layer i, unit is t; G i is organic carbon content of sediment layer i, unit is kg/t; Q is point engineering measure The amount of sediment that has been stored is in t; the Q potential is the total potential of sediment storage in point- type projects, and the unit is t; The maximum storage capacity of type engineering measures, the unit is t; α is the coefficient of other corrosion reduction effects that point type projects can play after reaching the maximum storage capacity, the value range is 0~1, no unit;

步骤七,计算工程措施使用年Y(例如工程措施建成时间为1990年,计算当年为2022年,则Y等于2022年-1990年=32年);Step 7, calculating the use year Y of engineering measures (for example, the construction time of engineering measures is 1990, and the calculation year is 2022, then Y is equal to 2022-1990=32 years);

步骤八,测定面型工程措施内部区域表层土壤有机碳含量G面内Step 8, measuring the surface soil organic carbon content G in the internal area of the surface engineering measures;

步骤九,测量面型工程措施的管控区域面积A(通过实地或在地形图上测量得到);Step 9, measure the area A of the control area of the surface-type engineering measures (obtained by measuring on the spot or on the topographic map);

步骤十,分别建立无面型工程措施和有面型工程措施径流小区并测定小区年总产沙量;无面型工程措施径流小区年总产沙量除以小区面积,得到无面型工程措施的土壤侵蚀模数X0;有面型工程措施径流小区年总产沙量除以小区面积,得到有面型工程措施的土壤侵蚀模数X;Step 10: Establish the runoff plots of the non-surface engineering measures and the surface engineering measures respectively and measure the total annual sediment production of the plots; divide the annual total sediment production of the runoff plots of the non-surface engineering measures by the area of the plot to obtain the non-surface engineering measures The soil erosion modulus X 0 of the surface engineering measures; the annual total sediment production of runoff plots with surface engineering measures divided by the area of the plot, to obtain the soil erosion modulus X of surface engineering measures;

步骤十一,测量输沙量,计算泥沙输移比SDR;Step 11, measure the amount of sediment transported, and calculate the sediment transport ratio SDR;

步骤十二,计算面型工程措施拦蓄泥沙固碳量,采用公式(6)-(8)进行计算:Step 12: Calculating the amount of sediment carbon fixation retained by surface-type engineering measures, using formulas (6)-(8) for calculation:

C=QG面内Y/1000                               (6)C plane = Q plane G plane Y/1000 (6)

Q=A(X0-X)                                  (7)Q surface = A surface (X 0 -X) (7)

X=X0E                                     (8)X=X 0 E (8)

式中:C为面型工程措施拦蓄泥沙固碳量,单位为t;Q为面型工程措施区年土壤侵蚀减沙/固沙量,单位为t/a;G面内为面型工程措施内部区域表层土壤有机碳含量,单位为kg/t;Y为工程措施使用年,单位为a;A为面型工程措施的管控区域面积,单位为km2;X0为无面型工程措施的土壤侵蚀模数,单位为t/(km2·a);X为有面型工程措施的土壤侵蚀模数,单位为t/(km2·a);E为工程措施因子(取值参考《中国水土流失方程(CSLE)》);In the formula: surface C is the amount of carbon sequestration and carbon sequestration of surface engineering measures, in t; surface Q is the annual soil erosion and sand reduction/sedimentation amount in the area of surface engineering measures, and the unit is t/a; surface G is surface The organic carbon content of the surface soil in the internal area of the engineering measure is in kg/t; Y is the service year of the engineering measure in a; the surface A is the control area of the surface engineering measure in km 2 ; X 0 is the non-face type Soil erosion modulus of engineering measures, unit is t/(km 2 ·a); X is soil erosion modulus of surface engineering measures, unit is t/(km 2 ·a); E is engineering measure factor (take Value refer to "China Soil Erosion Equation (CSLE)");

步骤十三,在面型工程措施控制区测定侵蚀区土壤有机碳含量G侵蚀和淤积区泥沙有机碳含量G淤积Step 13, measuring the soil organic carbon content G in the erosion area and the sediment organic carbon content G deposition in the deposition area in the surface engineering measure control area;

步骤十四,计算区域内工程措施总固碳量:Step 14, calculate the total carbon sequestration of engineering measures in the area:

不考虑泥沙输移比时采用公式(9)进行计算:When the sediment transport ratio is not considered, formula (9) is used for calculation:

C=C+C             (9)C=C surface +C point (9)

考虑泥沙输移比时采用公式(10)~(11)进行计算:When considering the sediment transport ratio, formulas (10) to (11) are used for calculation:

C=CSDR+C+C(1-SDR)b         (10)C=C plane SDR+C point +C plane (1-SDR)b (10)

b=1-G淤积/G控侵                (11)b=1-G deposition /G control invasion (11)

式中:C为区域内固碳量,单位为t;SDR为工程措施管控区的泥沙输移比;b为土壤被侵蚀后,搬运-淤积过程中损失的碳比例,在0~1之间;G控侵为面型工程措施控制区被侵蚀区土壤有机碳含量,单位为kg/t;G淤积为淤积区泥沙有机碳含量,单位为kg/t。In the formula: C is the amount of carbon sequestration in the region, and the unit is t; SDR is the sediment transport ratio in the control area of engineering measures; G control encroachment is the soil organic carbon content in the eroded area of the control area of surface engineering measures, the unit is kg/t; G sedimentation is the sediment organic carbon content in the deposition area, the unit is kg/t.

进一步的,步骤一中,所述面型措施包括:梯田、水平阶、水平沟、竹节沟、鱼鳞坑、大型果树坑、沙障固沙、工程护路,所述点型措施包括:坡面小型蓄排工程、沟头防护、谷坊、於地坝。Further, in step 1, the surface measures include: terraced fields, horizontal steps, horizontal ditches, bamboo ditch, fish scale pits, large fruit tree pits, sand barriers for sand fixation, engineering road protection, and the point measures include: slope surface Small-scale storage and drainage projects, ditch head protection, Gufang, Yudiba.

进一步的,步骤二中,所述已拦蓄泥沙量的分层淤积量,选取20cm为一层,淤积泥沙第i层淤积量表示为QiFurther, in step 2, for the layered deposition volume of the deposited sediment volume, 20 cm is selected as a layer, and the deposition volume of the i-th layer of sedimentation sediment is expressed as Q i ,

Qi=Ai·hi·BDi/1000Q i =A i ·h i ·BD i /1000

其中:Qi为淤积泥沙第i层淤积量,单位为t;Ai为淤积泥沙第i层的底面积,单位为平方米;hi为淤积泥沙第i层的厚度,单位为m;BDi为淤积泥沙第i层的淤积泥沙容重,单位为kg/m3Among them: Q i is the silting volume of the i-th layer of silt, in t; A i is the bottom area of the i-th layer of silt, in square meters; h i is the thickness of the i-th layer of silt, in t m; BD i is the sediment bulk density of the i-th layer of sediment, in kg/m 3 .

进一步的,步骤四中,α=1-(S/S),其中:S为点型工程措施淤满后点型工程措施出口产沙量,单位为t/yr;S为未建点型工程措施时出口产沙量,单位为t/yr;yr表示年。Further, in step 4, α=1-(S is full /S is not ), wherein: S is full is the point-type engineering measure outlet sediment yield after the point-type engineering measure is silted up, and the unit is t/yr; Export sediment yield when building point-type engineering measures, the unit is t/yr; yr means year.

进一步的,步骤五中,具体的操作为:点型工程采集整个淤积层泥沙样品,分层使用有机碳分析仪测定淤积泥沙第i层有机碳含量Gi,取20cm为一层;侵蚀区土壤有机碳含量根据采集的侵蚀区土壤表层5cm土样,使用有机碳分析仪测定,点型工程措施侵蚀区土壤有机碳含量表示为G控侵Further, in step 5, the specific operation is: collect the sediment samples of the entire alluvial layer in a point-type project, use an organic carbon analyzer to measure the organic carbon content G i of the i-th layer of the alluvial sediment, and take 20cm as a layer; The soil organic carbon content in the area is based on the 5cm soil samples collected from the soil surface of the erosion area, and is measured by an organic carbon analyzer. The soil organic carbon content in the erosion area of point-type engineering measures is expressed as G control erosion .

进一步的,步骤八中,具体操作为:采集面型工程措施内部区域表层5cm土壤样品,使用有机碳分析仪测定土壤有机碳含量G面内Further, in Step 8, the specific operation is: collect a 5 cm soil sample in the surface layer of the internal area of the surface-shaped engineering measures, and use an organic carbon analyzer to measure the soil organic carbon content G inside the surface .

进一步的,步骤十中,建立无面型工程措施和有面型工程措施径流小区并测定小区年总产沙量的具体操作为:无面型工程措施和有面型工程措施径流小区除工程措施外其他参数均一样,小区下方放置泥沙集蓄设备,收集径流中的泥沙,每次降雨后对泥沙进行采集和称重,获取每次产沙量,一年中每次产沙量之和为年总产沙量。Further, in step 10, the specific operation of establishing runoff plots with surfaceless engineering measures and surfaced engineering measures and measuring the total annual sediment production in the plots is: removing engineering measures for runoff plots with surfaceless engineering measures and surfaced engineering measures The other parameters are the same, and the sediment storage equipment is placed under the community to collect the sediment in the runoff. After each rainfall, the sediment is collected and weighed to obtain the amount of each sediment production, and the annual sediment production amount The sum is the total annual sediment production.

进一步的,步骤十一中,测量输沙量,计算泥沙输移比SDR的具体操作为:在区域出口设置卡口站,放置泥沙集蓄设备,收集径流中的泥沙,每次降雨后对泥沙进行采集和称重,获取每次输沙量,一年中每次输沙量之和为年总输沙量,总输沙量与总产沙量的比值为泥沙输移比SDR。Further, in step eleven, the specific operation of measuring the amount of sediment transport and calculating the sediment transport ratio SDR is as follows: set up checkpoint stations at the exit of the area, place sediment storage equipment, collect sediment in the runoff, and Afterwards, the sediment is collected and weighed to obtain the amount of sediment transported each time. The sum of each sediment transport amount in a year is the total annual sediment transport amount, and the ratio of the total sediment transport volume to the total sediment production is the sediment transport volume. than SDR.

进一步的,步骤十三中,具体的操作为:在面型工程措施控制区的侵蚀区采集的侵蚀区土壤表层5cm土样,在淤积区采集整个淤积层的泥沙样品,分别使用有机碳分析仪测定侵蚀区土壤有机碳含量G控侵和淤积区泥沙有机碳含量G积。Further, in step thirteen, the specific operation is: collect 5 cm soil samples in the soil surface layer of the erosion area in the erosion area of the surface engineering measure control area, and collect the sediment samples of the entire alluvial layer in the deposition area, and use organic carbon analysis respectively. Determination of the soil organic carbon content G in the erosion area to control erosion and the sediment organic carbon content G deposition in the sedimentation area.

本发明的优点和有益效果是:本发明首次提出了针对面型和点型水土保持工程措施,在考虑了土壤侵蚀-输移-淤积过程土壤有机碳变化的情况下,根据不同资料情况,普遍适用的水土保持工程措施碳汇计算方法。本发明提供不同水土保持工程措施的土壤碳汇潜力计算方法,相比现有技术中其他计算方法更加科学全面,考虑了水土保持工程措施的类别和泥沙在输移过程中碳的损失,结果更准确。The advantages and beneficial effects of the present invention are: the present invention proposes for the first time the project measures for surface and point water and soil conservation, and in consideration of the change of soil organic carbon in the process of soil erosion-transportation-siltation, according to different data situations, generally Applicable methods for calculating carbon sequestration of water and soil conservation engineering measures. The present invention provides a soil carbon sink potential calculation method for different soil and water conservation engineering measures, which is more scientific and comprehensive than other calculation methods in the prior art, considering the types of soil and water conservation engineering measures and the carbon loss of sediment during transport, and the result more acurrate.

附图说明Description of drawings

下面结合附图和实施例对本发明作进一步说明。The present invention will be further described below in conjunction with drawings and embodiments.

图1是本发明实施例方法流程示意图;Fig. 1 is a schematic flow chart of the method of the embodiment of the present invention;

图2是本发明实施例1中2梯田和卡口站位置图。Fig. 2 is a location diagram of 2 terraced fields and bayonet stations in Embodiment 1 of the present invention.

具体实施方式Detailed ways

实施例1Example 1

如图1所示,一种测算水土保持工程措施土壤碳汇潜力的方法,包括以下步骤:As shown in Figure 1, a method for measuring soil carbon sink potential of soil and water conservation engineering measures includes the following steps:

步骤一,对水土保持工程措施进行分类,对面尺度土壤侵蚀控制作用的工程措施划分为面型,对于已侵蚀的泥沙进行拦截存储作用的措施划分为点型措施。The first step is to classify the engineering measures of water and soil conservation. The engineering measures for surface-scale soil erosion control are divided into surface-type measures, and the measures for intercepting and storing eroded sediment are divided into point-type measures.

表1水土保持工程措施分类Table 1 Classification of water and soil conservation engineering measures

本专利适用的工程措施Engineering measures applicable to this patent 分类Classification 梯田terraced fields 面型face shape 水平阶horizontal step 面型face shape 水平沟horizontal groove 面型face shape 竹节沟Zhujiegou 面型face shape 鱼鳞坑fish scale pit 面型face shape 大型果树坑large fruit tree pit 面型face shape 沙障固沙Sand barrier 面型face shape 工程护路Engineering road protection 面型face shape 坡面小型蓄排工程Small slope storage and drainage project 点型point type 沟头防护ditch head protection 点型point type 谷坊Gufang 点型point type 淤地坝check dam 点型point type

步骤二,确定点型工程措施的最大库容Q最大和已拦蓄泥沙量的分层淤积量,分层淤积量的总和为点型工程措施已拦蓄泥沙量QStep 2, determine the maximum storage capacity Q of the point-type engineering measures and the layered deposition amount of the deposited sediment amount, the sum of the layered deposition amount is the point-type engineering measure and the accumulated sediment amount Q ;

所述已拦蓄泥沙量的分层淤积量,选取20cm为一层,淤积泥沙第i层淤积量表示为QiFor the layered deposition amount of the deposited sediment amount, 20cm is selected as a layer, and the deposition amount of the i layer of the deposited sediment is expressed as Q i ,

Qi=Ai·hi·BDi/1000Q i =A i ·h i ·BD i /1000

其中:Qi为淤积泥沙第i层淤积量,单位为t;Ai为淤积泥沙第i层的底面积,单位为平方米;hi为淤积泥沙第i层的厚度,单位为m;BDi为淤积泥沙第i层的淤积泥沙容重,单位为kg/m3Among them: Q i is the silting volume of the i-th layer of silt, in t; A i is the bottom area of the i-th layer of silt, in square meters; h i is the thickness of the i-th layer of silt, in t m; BD i is the sediment bulk density of the i-th layer of sediment, in kg/m 3 .

步骤三,计算点型工程未淤满部分泥沙拦蓄潜力:采用点型工程措施的最大库容减去点型工程措施已拦蓄泥沙量,即:Q最大-QStep 3: Calculating the sediment storage potential of the unsilted part of the point-type project: using the maximum storage capacity of the point-type project measures minus the amount of sediment stored by the point-type project measures, that is : Qmax - Qha.

步骤四,计算点型工程达到最大库容后可发挥的其他减蚀作用的系数α,为1减去点型工程措施淤满后点型工程措施出口产沙量与未建点型工程措施时出口产沙量的比值;Step 4: Calculate the coefficient α of other corrosion reduction effects that can be exerted by point-type engineering measures after reaching the maximum storage capacity, which is 1 minus the amount of sediment produced at the outlet of point-type engineering measures after the point-type engineering measures are silted and the outlet when no point-type engineering measures are built The ratio of sediment yield;

α=1-(S/S),其中:S为点型工程措施淤满后点型工程措施出口产沙量,单位为t/yr;S为未建点型工程措施时出口产沙量,单位为t/yr。α=1-(S full /S unfinished ), where: S full is the sediment yield at the outlet of the point-type engineering measures after the point-type engineering measures are silted, and the unit is t/yr; S is not the outlet when no point-type engineering measures are built Sediment production, unit is t/yr.

步骤五,测定点型工程淤积层泥沙分层有机碳含量和侵蚀区土壤有机碳含量:点型工程采集整个淤积层泥沙样品,分层使用有机碳分析仪测定淤积泥沙第i层有机碳含量Gi,取20cm为一层;侵蚀区土壤有机碳含量根据采集的侵蚀区土壤表层5cm土样,使用有机碳分析仪测定点型工程措施侵蚀区土壤有机碳含量G蚀。Step 5: Determining the organic carbon content of sediment layers in point-type projects and soil organic carbon content in erosion areas: collect sediment samples from the entire sediment layer in point-type projects, and use an organic carbon analyzer to measure the organic carbon content of the i-th layer of sediment layers. For the carbon content G i , 20 cm is taken as a layer; the soil organic carbon content in the erosion area is based on the 5 cm soil samples collected from the soil surface in the erosion area, and the organic carbon analyzer is used to measure the soil organic carbon content G erosion in the erosion area of point-type engineering measures.

步骤六,根据公式(1)-(5)计算点型工程措施拦蓄泥沙固碳量;Step 6, according to the formula (1)-(5), calculate the amount of sediment and carbon fixation in point-type engineering measures;

C=C点已+C点未          (1)C point = C point already + C point not (1)

Figure GDA0003887022500000071
Figure GDA0003887022500000071

Figure GDA0003887022500000072
Figure GDA0003887022500000072

C点未=Q潜力G侵蚀/1000        (4) Point C = Q Potential G Erosion /1000 (4)

Q潜力=(Q最大-Q)+Q最大α       (5)Q Potential = ( QMax - QHave ) + QMax α (5)

式中:C为点型工程措施拦蓄泥沙固碳量,单位为t;C点已为点型工程措施已淤积泥沙固碳量,单位为t;C点未为点型工程措施未来固碳潜力,单位为t;Qi为淤积泥沙第i层淤积量,单位为t;Gi为淤积泥沙第i层有机碳含量,单位为kg/t;Q为点型工程措施已拦蓄泥沙量,单位为t;Q潜力为点型工程总泥沙拦蓄总潜力,单位为t;G侵蚀为点型工程措施侵蚀区土壤有机碳含量,单位为kg/t;Q最大为点型工程措施的最大库容,单位为t;α为点型工程达到最大库容后可发挥的其他减蚀作用的系数,取值范围0~1,无单位。In the formula: point C is the amount of sediment carbon sequestered by point-type engineering measures, and the unit is t; point C is the amount of deposited sediment carbon-fixed carbon by point-type engineering measures, and the unit is t; point C is the future of point-type engineering measures Carbon sequestration potential, unit is t; Q i is deposition amount of sediment layer i, unit is t; G i is organic carbon content of sediment layer i, unit is kg/t; Q is point engineering measure The amount of sediment that has been stored is in t; the Q potential is the total sediment storage potential of the point-type project, and the unit is t; G erosion is the soil organic carbon content in the eroded area of the point-type project measure, and the unit is kg/t; the maximum Q is The maximum storage capacity of point-type engineering measures, the unit is t; α is the coefficient of other corrosion reduction effects that point-type projects can play after reaching the maximum storage capacity, and the value ranges from 0 to 1, with no unit.

步骤七,计算工程措施使用年Y,例如工程措施建成时间为1990年,计算当年为2022年,则Y等于2022年-1990年=32年;Step 7, calculate the use year Y of the engineering measures, for example, the construction time of the engineering measures is 1990, and the calculation year is 2022, then Y is equal to 2022-1990=32 years;

步骤八,测定面型工程措施内部区域表层土壤有机碳含量G面内;具体操作为:采集面型工程措施内部区域表层5cm土壤样品,使用有机碳分析仪测定土壤有机碳含量G面内Step 8, measure the surface soil organic carbon content G in the inner area of the surface-shaped engineering measures; the specific operation is: collect 5 cm soil samples in the surface layer of the inner area of the surface-shaped engineering measures, and use an organic carbon analyzer to measure the soil organic carbon content G in the surface .

步骤九,测量面型工程措施的管控区域面积A(通过实地或在地形图上测量得到)。Step 9: Measure the area A of the control area of the surface-type engineering measures (obtained by on-site or topographical map measurement).

步骤十,分别建立无面型工程措施和有面型工程措施径流小区并测定小区年总产沙量;无面型工程措施径流小区年总产沙量除以小区面积,得到无面型工程措施的土壤侵蚀模数X0;有面型工程措施径流小区年总产沙量除以小区面积,得到有面型工程措施的土壤侵蚀模数X;建立无面型工程措施和有面型工程措施径流小区并测定小区年总产沙量的具体操作为:无面型工程措施和有面型工程措施径流小区除工程措施外其他参数均一样,小区下方放置泥沙集蓄设备,收集径流中的泥沙,每次降雨后对泥沙进行采集和称重,获取每次产沙量,一年中每次产沙量之和为年总产沙量。Step 10: Establish the runoff plots of the non-surface engineering measures and the surface engineering measures respectively and measure the total annual sediment production of the plots; divide the annual total sediment production of the runoff plots of the non-surface engineering measures by the area of the plot to obtain the non-surface engineering measures The soil erosion modulus X 0 of the surface engineering measure; the annual total sediment production of the runoff area of the surface engineering measure is divided by the area of the plot, and the soil erosion modulus X of the surface engineering measure is obtained; the establishment of the non-surface engineering measure and the surface engineering measure The specific operation of the runoff plot and the determination of the total annual sediment production in the plot is as follows: the runoff plots with no surface engineering measures and the runoff plots with surface engineering measures have the same parameters except engineering measures, and sediment collection equipment is placed under the plots to collect the sediment in the runoff. For sediment, the sediment is collected and weighed after each rainfall to obtain the amount of each sediment production, and the sum of each sediment production in a year is the total annual sediment production.

步骤十一,测量输沙量,计算泥沙输移比SDR的具体操作为:在区域出口设置卡口站,放置泥沙集蓄设备,收集径流中的泥沙,每次降雨后对泥沙进行采集和称重,获取每次输沙量,一年中每次输沙量之和为年总输沙量,总输沙量与总产沙量的比值为泥沙输移比SDR。Step 11: Measure the amount of sediment transported and calculate the sediment transport ratio SDR. The specific operations are as follows: set up a bayonet station at the exit of the area, place sediment storage equipment, collect the sediment in the runoff, and remove the sediment after each rainfall. Collecting and weighing are carried out to obtain the amount of sediment transported each time. The sum of each sediment transport amount in a year is the total annual sediment transport amount, and the ratio of the total sediment transport amount to the total sediment production is the sediment transport ratio SDR.

步骤十二,计算面型工程措施拦蓄泥沙固碳量,采用公式(6)-(8)进行计算:Step 12: Calculating the amount of sediment carbon fixation retained by surface-type engineering measures, using formulas (6)-(8) for calculation:

C=QG面内Y/1000                               (6)C plane = Q plane G plane Y/1000 (6)

Q=A(X0-X)                                 (7)Q surface = A surface (X 0 -X) (7)

X=X0E                                    (8)X=X 0 E (8)

式中:C为面型工程措施拦蓄泥沙固碳量,单位为t;Q为面型工程措施区年土壤侵蚀减沙/固沙量,单位为t/a;G面内为面型工程措施内部区域表层土壤有机碳含量,单位为kg/t;Y为工程措施使用年,单位为a;A为面型工程措施的管控区域面积,单位为km2;X0为无面型工程措施的土壤侵蚀模数,单位为t/(km2·a);X为有面型工程措施的土壤侵蚀模数,单位为t/(km2·a);E为工程措施因子(取值参考《中国水土流失方程(CSLE)》)。In the formula: surface C is the amount of carbon sequestration and carbon sequestration of surface engineering measures, in t; surface Q is the annual soil erosion and sand reduction/sedimentation amount in the area of surface engineering measures, and the unit is t/a; surface G is surface The organic carbon content of the surface soil in the internal area of the engineering measure is in kg/t; Y is the service year of the engineering measure in a; the surface A is the control area of the surface engineering measure in km 2 ; X 0 is the non-face type Soil erosion modulus of engineering measures, unit is t/(km 2 ·a); X is soil erosion modulus of surface engineering measures, unit is t/(km 2 ·a); E is engineering measure factor (take The values refer to "China Soil Erosion Equation (CSLE)").

步骤十三,在面型工程措施控制区测定侵蚀区土壤有机碳含量G控侵和淤积区泥沙有机碳含量G淤积;具体的操作为:在面型工程措施控制区的侵蚀区采集的侵蚀区土壤表层5cm土样,在淤积区采集整个淤积层的泥沙样品,分别使用有机碳分析仪测定侵蚀区土壤有机碳含量G控侵和淤积区泥沙有机碳含量G淤积Step thirteen, measure the soil organic carbon content G in the erosion zone in the surface engineering measure control area, control erosion and sediment organic carbon content G deposition in the sedimentation area; the specific operation is: the erosion collected in the erosion area of the surface engineering measure control area The 5cm soil samples from the surface layer of the soil in the siltation area were collected, and the sediment samples of the entire silt layer were collected in the siltation area, and the organic carbon analyzer was used to measure the soil organic carbon content G in the erosion area and the sediment organic carbon content in the siltation area G deposition .

步骤十四,计算区域内工程措施总固碳量:Step 14, calculate the total carbon sequestration of engineering measures in the area:

不考虑泥沙输移比时采用公式(9)进行计算:When the sediment transport ratio is not considered, formula (9) is used for calculation:

C=C+C        (9)C=C surface +C point (9)

考虑泥沙输移比时采用公式(10)~(11)进行计算:When considering the sediment transport ratio, formulas (10) to (11) are used for calculation:

C=CSDR+C+C(1-SDR)b                (10)C=C plane SDR+C point +C plane (1-SDR)b (10)

b=1-G淤积/G控侵                      (11)b=1-G deposition /G control invasion (11)

式中:C为区域内固碳量,单位为t;SDR为工程措施管控区的泥沙输移比;b为土壤被侵蚀后,搬运-淤积过程中损失的碳比例,在0~1之间;G控侵为面型工程措施控制区被侵蚀区土壤有机碳含量,单位为kg/t;G淤积为淤积区泥沙有机碳含量,单位为kg/t。In the formula: C is the amount of carbon sequestration in the region, and the unit is t; SDR is the sediment transport ratio in the control area of engineering measures; G control encroachment is the soil organic carbon content in the eroded area of the control area of surface engineering measures, the unit is kg/t; G sedimentation is the sediment organic carbon content in the deposition area, the unit is kg/t.

本实施例中:In this example:

点型工程措施计算数据如下:The calculation data of point-type engineering measures are as follows:

点工程措施选取淤地坝,淤地坝建坝设计最大库容Q最大500t,侵蚀区实测土壤有机碳含量为G侵蚀2.3kg/t,淤地坝淤满后实测出口产沙量S为400t/a,未建淤地坝时实测出口产沙量S为500t/a,根据实测的分层淤积量和有机碳含量(表1)和公式(1)-(5)可以得到:已淤积泥沙Q420t,已淤积部分的固碳量C点已为0.825t,α为1-400t/a÷500t/a=0.2,未来淤积潜力Q最大-Q为80t,总潜力Q潜力=(Q最大-Q)+Q最大α,计算得Q潜力为180t,点型工程措施未来固碳潜力C点未计算得0.414,淤地坝拦蓄泥沙固碳量为C=C点已+C点未为1.239t。The point engineering measure is to select the check dam, the maximum storage capacity Q of the check dam is 500t, the measured soil organic carbon content in the erosion area is G erosion 2.3kg/t, and the measured outlet sediment production S is 400t after the check dam is full . /a, when the check dam is not built, the measured outlet sediment yield S is 500t/a, according to the measured layered sedimentation and organic carbon content (Table 1) and formulas (1)-(5), it can be obtained: Sediment Q has been 420t, the amount of carbon sequestration in the silted part C has been 0.825t, α is 1-400t/a÷500t/a=0.2, the future silting potential Q is the largest - Q is 80t, the total potential Q potential = ( Qmax - Qha )+ Qmaxα , the calculated Q potential is 180t, the future carbon sequestration potential of point-type engineering measures is not calculated at point C, which is 0.414, and the carbon sequestration amount of the check dam is C point = C point has been +C point is not 1.239t.

表2已淤积泥沙固碳量Table 2 Carbon sequestration amount of deposited sediment

深度(cm)Depth (cm) <![CDATA[分层淤积量Q<sub>i</sub>(t)]]><![CDATA[Layered deposition volume Q<sub>i</sub>(t)]]> <![CDATA[有机碳含量G<sub>i</sub>(kg/t)]]><![CDATA[Organic carbon content G<sub>i</sub>(kg/t)]]> <![CDATA[Q<sub>i</sub>*G<sub>i</sub>固碳量(t)]]><![CDATA[Q<sub>i</sub>*G<sub>i</sub>Carbon Sequestration(t)]]> 0-200-20 5050 2.12.1 0.1050.105 20-4020-40 7070 2.22.2 0.1540.154 0-600-60 8080 1.51.5 0.120.12 60-8060-80 100100 1.71.7 0.170.17 80-10080-100 120120 2.32.3 0.2760.276

面型工程计算数据如下:Surface engineering calculation data are as follows:

本实施例面型工程措施选取梯田,根据遥感影像(图2)确定梯田范围,测量梯田面积A为0.1km2,梯田建设期为1995-1997年,1998年开始使用,工程措施使用年Y为24年(2022-1998年),根据土壤样品,梯田多年平均有机碳含量G面内为3.4kg/t,在区域出口建设卡口站,测量区域产沙量,无梯田措施时产沙量为20t/a。根据CSLE,计算该无梯田时土壤侵蚀模数X0为500t/(km2·a),石坎水平梯田工程因子E为0.124,根据公式(6)-(8),计算得到:泥沙输移比SDR为0.4,工程措施建立后土壤侵蚀模数X为62t/(km2·a),梯田年土壤侵蚀减沙量Q为43.8t/a,梯田拦蓄泥沙固碳量C为3.57t。Terraces are selected for surface engineering measures in this embodiment, and the range of terraces is determined according to remote sensing images (Figure 2). The measured terrace area A is 0.1km 2 . For 24 years (2022-1998), according to the soil samples, the multi-year average organic carbon content of the terraced fields in the G plane is 3.4kg/t. A bayonet station will be built at the regional exit to measure the regional sand production. The sand production without terrace measures It is 20t/a. According to CSLE, the soil erosion modulus X 0 is calculated as 500t/(km 2 ·a) when there are no terraces, and the engineering factor E of stone ridge level terraces is 0.124. According to the formulas (6)-(8), the calculation results are as follows: sediment transport The specific SDR is 0.4, the soil erosion modulus X is 62t/(km 2 ·a) after the engineering measures are established, the annual soil erosion and sediment reduction amount Q of terraced fields is 43.8t/a, and the amount of carbon sequestration by terraced fields is 3.57 t.

工程措施总固碳量:不考虑泥沙输移比的情况下,区域固碳量C为1.239t+3.57t=4.809t;在考虑泥沙输移比的情况下,采集侵蚀区和淤积区的土壤样品,测定土壤有机碳的含量分别为G控侵2.3t/a和G淤积1.5t/a,根据公式(9)-(11),计算得到b为0.35,区域内固碳量C为3.4167t。Total carbon sequestration of engineering measures: without considering the sediment transport ratio, the regional carbon sequestration amount C is 1.239t+3.57t=4.809t; in the case of considering the sediment transport ratio, the collection of erosion areas and deposition areas The measured soil organic carbon content is 2.3t/a for G control invasion and 1.5t/a for G deposition respectively. According to the formula (9)-(11), b is calculated as 0.35, and the carbon fixation amount C in the area is 3.4167t.

最后应说明的是,以上仅用以说明本发明的技术方案而非限制,尽管参照较佳布置方案对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围。Finally, it should be noted that the above is only used to illustrate the technical solution of the present invention and not to limit it. Although the present invention has been described in detail with reference to the preferred arrangement, those skilled in the art should understand that the technical solution of the present invention can be modified Or an equivalent replacement without departing from the spirit and scope of the technical solution of the present invention.

Claims (9)

1. The method for measuring and calculating the carbon sink potential of the soil by the water and soil conservation engineering measure is characterized by comprising the following steps of:
step one, classifying adopted water and soil conservation engineering measures: the method is divided into a surface type measure and a point type measure; the surface type measure is an engineering measure for controlling the soil erosion of the surface scale; the point-type measure is a measure for intercepting and storing corroded sediment;
step two, determining the maximum storage capacity Q of the point engineering measure Maximum value And a stratified sedimentation amount of the blocked sediment amount, the sum of the stratified sedimentation amounts being the blocked sediment amount Q of the dot type Cheng Cuoshi Has already been provided with
Calculating sediment storage potential of the part of the point type engineering which is not full of sediment: the sediment blocking amount of the point type tool Cheng Cuoshi is subtracted from the maximum reservoir capacity of the point type tool, namely: q (Q) Maximum value -Q Has already been provided with
Calculating the coefficient alpha of other corrosion reducing effects which can be exerted after the spot-type engineering reaches the maximum reservoir capacity:
measuring the layered organic carbon content of sediment in the fixed-point engineering siltation layer and the organic carbon content of soil in an erosion area;
step six, calculating sediment blocking and carbon fixation amount of the spot type tool Cheng Cuoshi: the calculation is performed by using formulas (1) to (5):
C point(s) =C The dot has been +C Dot is not (1)
Figure QLYQS_1
Figure QLYQS_2
C Dot is not =Q Potential G Erosion of /1000 (4)
Q Potential =(Q Maximum value -Q Has already been provided with )+Q Maximum value α (5)
Wherein: c (C) Point(s) The unit is t for blocking sediment and carbon fixation of the spot type Cheng Cuoshi; c (C) The dot has been The carbon fixation amount of the deposited sediment is t for the point-type engineering measure; c (C) Dot is not The future carbon fixation potential is measured by point-type engineering measures, and the unit is t; q (Q) i The unit of the i layer deposition amount of the deposited sediment is t; g i The unit of the organic carbon content of the ith layer of sediment is kg/t; q (Q) Has already been provided with The sediment quantity is blocked for point-type engineering measures, and the unit is t; q (Q) Potential The total potential is blocked by a dotted type Cheng Nisha, and the unit is t; g Erosion of The organic carbon content of soil in the erosion area is kg/t for the point-type engineering measure; q (Q) Maximum value The maximum storage capacity of the spot-type engineering measure is represented by t; alpha is the coefficient of other corrosion reducing effects which can be exerted after the point type engineering reaches the maximum reservoir capacity, the value range is 0-1, and no unit exists; n is the total number of layers of deposited sediment, and no unit exists;
step seven, calculating the service year Y of engineering measures;
step eight, determining the organic carbon content G of the surface soil in the inner area of the surface type engineering measure In-plane
Step nine, measuring the area A of a control area of the surface engineering measure Flour with a plurality of grooves
Step ten, respectively establishing a non-surface engineering measure and a surface engineering Cheng Cuoshi runoff district, and measuring the annual total sand production of the district; no-surface type Cheng Cuoshi diameterDividing the annual total sand production of the flowing cell by the cell area to obtain the soil erosion modulus X of the non-surface engineering measure 0 The method comprises the steps of carrying out a first treatment on the surface of the Dividing the annual total sand production of the area-shaped Cheng Cuoshi runoff district by the area of the district to obtain a soil erosion modulus X of the area-shaped shaping measure;
step eleven, measuring the sand conveying amount and calculating the sediment conveying ratio SDR;
step twelve, calculating the sediment storage and carbon fixation amount of the surface-shaped tool Cheng Cuoshi, and calculating by adopting formulas (6) - (8):
C flour with a plurality of grooves =Q Flour with a plurality of grooves G In-plane Y/1000 (6)
Q Flour with a plurality of grooves =A Flour with a plurality of grooves (X 0 -X) (7)
X=X 0 E (8)
Wherein: c (C) Flour with a plurality of grooves The unit of the sediment accumulation carbon fixation quantity is t for the surface type Cheng Cuoshi; q (Q) Flour with a plurality of grooves The unit of the soil erosion sand reduction/fixation amount is t/a; g In-plane The organic carbon content of the soil on the surface layer of the inner area of the surface type engineering measure is kg/t; y is the service year of engineering measures, and the unit is a; a is that Flour with a plurality of grooves The area of the control area is km for the surface engineering measure 2 ;X 0 The soil erosion modulus is the unit of t/(km) for the non-planar engineering measure 2 A); x is the soil erosion modulus of the surface engineering measure, and the unit is t/(km) 2 A); e is an engineering measure factor;
thirteenth step, determining soil organic carbon content G in erosion zone in face type engineering measure control zone Intrusion control And sediment organic carbon content G in the sedimentation zone Fouling of
Fourteen, calculating the total carbon fixation amount of engineering measures in the area:
and when the sediment transport ratio is not considered, adopting a formula (9) for calculation:
C=C flour with a plurality of grooves +C Point(s) (9) When the sediment transport ratio is considered, the calculation is carried out by adopting formulas (10) to (11):
C=C flour with a plurality of grooves SDR+C Point(s) +C Flour with a plurality of grooves (1-SDR)b (10)
b=1-G Fouling of /G Intrusion control (11)
Wherein: c is the carbon fixation amount in the region, and the unit is t; SDR is sediment transport ratio of engineering measure management and control area; b is the carbon ratio lost in the carrying-sedimentation process after the soil is eroded, and is between 0 and 1; g Intrusion control The soil organic carbon content of the eroded area in the area is controlled by the surface engineering measure, and the unit is kg/t; g Fouling of The sediment organic carbon content of the silt area in the area is controlled by the surface engineering measures, and the unit is kg/t.
2. The method for measuring and calculating the carbon sequestration potential of soil by soil and water conservation engineering measures according to claim 1, wherein in the first step, the surface type measure comprises: terrace, horizontal step, horizontal ditch, bamboo joint ditch, fish scale hole, large-scale fruit tree hole, sand barrier sand fixation, engineering protection road, the punctiform measure includes: slope small-sized storage and discharge engineering, ditch head protection, check dam and silt dam.
3. The method for measuring and calculating the carbon sequestration potential of soil in soil and water conservation engineering measures according to claim 1, wherein in the second step, the stratified sediment amount of the blocked sediment amount is selected to be 20cm as one layer, and the ith sediment amount of the sediments is represented as Q i
Q i =A i ·h i ·BD i /1000
Wherein: q (Q) i The unit of the i layer deposition amount of the deposited sediment is t; a is that i The bottom area of the ith layer is the unit of square meter for sediment deposition; h is a i The thickness of the ith layer of sediment is m; BD (BD) i The unit of the volume weight of the sediment is kg/m for the ith layer of the sediment 3
4. The method for measuring and calculating the carbon sequestration potential of soil and water conservation engineering measures according to claim 1, wherein in the fourth step, α=1- (S) Full of /S Not yet ) Wherein: s is S Full of The sand production amount is exported for the point-type engineering measure after the point-type engineering measure is full; s is S Not yet And outputting the sand production amount when the point-type engineering measures are not built.
5. The method for measuring and calculating the carbon sequestration potential of soil in soil and water conservation engineering measures according to claim 1, wherein in the fifth step, the specific operation of measuring the layered organic carbon content of the silt in the site-specific engineering siltation layer and the organic carbon content of the soil in the erosion area is as follows: the point-type engineering sample of the sediment of the whole siltation layer is collected, and the organic carbon content G of the ith layer of the siltation sediment is measured by using an organic carbon analyzer in a layered manner i Taking 20cm as a layer; the organic carbon content of the soil in the erosion zone is determined by using an organic carbon analyzer according to a collected 5cm soil sample on the soil surface layer of the erosion zone, and the organic carbon content of the soil in the erosion zone is expressed as G by using point-type engineering measures Erosion of
6. The method for measuring and calculating the carbon sequestration potential of soil in soil and water conservation engineering measures according to claim 1, wherein in the eighth step, the organic carbon content G of the surface soil in the inner area of the surface engineering measures is measured In-plane The specific operation is as follows: collecting a soil sample with the surface layer of 5cm in the inner area of the surface engineering measure, and measuring the organic carbon content G of the soil by using an organic carbon analyzer In-plane
7. The method for measuring and calculating the carbon sequestration potential of soil in soil and water conservation engineering measures according to claim 1, wherein in the step ten, the specific operations of establishing a non-planar engineering measure and a planar engineering Cheng Cuoshi runoff plot and measuring the annual total sand production of the plot are as follows: the non-surface engineering measures are the same as those of the surface engineering Cheng Cuoshi runoff plot except the engineering measures, sediment collecting and accumulating equipment is placed below the plot, sediment in the runoff is collected, the sediment is collected and weighed after each rainfall, each sediment production amount is obtained, and the sum of each sediment production amount in one year is the total sediment production amount in one year.
8. The method for measuring and calculating the carbon sequestration potential of soil in soil and water conservation engineering measures according to claim 1, wherein in the eleventh step, the concrete operation of measuring the sediment transport ratio SDR is as follows: and a bayonet station is arranged at an outlet of the area, sediment collection and storage equipment is arranged, sediment in runoff is collected, the sediment is collected and weighed after each rainfall, each time of sediment transport is obtained, the sum of each time of sediment transport in one year is the total sediment transport in one year, and the ratio of the total sediment transport to the total sediment production is the sediment transport ratio SDR.
9. The method for measuring and calculating the carbon sequestration potential of soil in soil and water conservation engineering measures according to claim 1, wherein in step thirteen, the soil organic carbon content G of the eroded area is measured in the area of control of engineering measures of surface type Corrosion control And sediment organic carbon content G in the sedimentation zone Fouling of The specific operation is as follows: soil samples of 5cm of the soil surface layer of the erosion zone are collected in the erosion zone of the surface engineering measure control zone, sediment samples of the whole sedimentation layer are collected in the sedimentation zone, and the organic carbon content G of the soil of the erosion zone is measured by an organic carbon analyzer respectively Corrosion control And sediment organic carbon content G in the sedimentation zone Fouling of
CN202211070781.2A 2022-09-02 2022-09-02 Method for measuring and calculating carbon sink potential of soil in water and soil conservation engineering measures Active CN115407049B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211070781.2A CN115407049B (en) 2022-09-02 2022-09-02 Method for measuring and calculating carbon sink potential of soil in water and soil conservation engineering measures

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211070781.2A CN115407049B (en) 2022-09-02 2022-09-02 Method for measuring and calculating carbon sink potential of soil in water and soil conservation engineering measures

Publications (2)

Publication Number Publication Date
CN115407049A CN115407049A (en) 2022-11-29
CN115407049B true CN115407049B (en) 2023-05-05

Family

ID=84163615

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211070781.2A Active CN115407049B (en) 2022-09-02 2022-09-02 Method for measuring and calculating carbon sink potential of soil in water and soil conservation engineering measures

Country Status (1)

Country Link
CN (1) CN115407049B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116429723B (en) * 2023-04-04 2023-09-05 中国水利水电科学研究院 Method for evaluating carbon sequestration of soil by soil and water conservation cultivation measures
CN116644261B (en) * 2023-05-30 2024-01-26 中国水利水电科学研究院 Urban bluish green space and construction land storage and discharge balance calculation method
CN118446403B (en) * 2024-04-07 2025-02-14 中国科学院地理科学与资源研究所 Black soil carbon sink potential assessment method and system based on grid precise zoning

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777688A (en) * 2016-12-15 2017-05-31 中国水利水电科学研究院 A kind of erosion and torrent control works runoff and sedimentation reduction method for quantitatively evaluating
CN111260250A (en) * 2020-02-13 2020-06-09 中国科学院地理科学与资源研究所 Method for evaluating ecological service function of water and soil conservation measure
CN113610334A (en) * 2021-05-14 2021-11-05 上海同渠工程咨询有限公司 Intelligent planning method, system, terminal and storage medium for carbon sink capacity measurement monitoring point
CN114169161A (en) * 2021-12-02 2022-03-11 辽宁省农业发展服务中心 Method and system for estimating space-time variation and carbon sequestration potential of soil organic carbon

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110174506B (en) * 2019-06-28 2020-07-07 中国科学院地球化学研究所 Method for estimating organic carbon in soil in karst region

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777688A (en) * 2016-12-15 2017-05-31 中国水利水电科学研究院 A kind of erosion and torrent control works runoff and sedimentation reduction method for quantitatively evaluating
CN111260250A (en) * 2020-02-13 2020-06-09 中国科学院地理科学与资源研究所 Method for evaluating ecological service function of water and soil conservation measure
CN113610334A (en) * 2021-05-14 2021-11-05 上海同渠工程咨询有限公司 Intelligent planning method, system, terminal and storage medium for carbon sink capacity measurement monitoring point
CN114169161A (en) * 2021-12-02 2022-03-11 辽宁省农业发展服务中心 Method and system for estimating space-time variation and carbon sequestration potential of soil organic carbon

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李智广.碳达峰与碳中和目标下水土保持碳汇的机理、途径及特征.《水土保持通报》.2022,第42卷(第3期),312-317. *

Also Published As

Publication number Publication date
CN115407049A (en) 2022-11-29

Similar Documents

Publication Publication Date Title
CN115407049B (en) Method for measuring and calculating carbon sink potential of soil in water and soil conservation engineering measures
Olley et al. Changes in the flux of sediment in the Upper Murrumbidgee catchment, Southeastern Australia, since European settlement
CN102201034B (en) Method for estimating nonpoint source pollution load of village adjacent to river
CN109614655B (en) Method for researching river runoff
CN114580910B (en) Watershed sediment connectivity evaluation method considering road influence
Li et al. Fingerprinting sediment sources in a typical karst catchment of southwest China
CN110838079A (en) Intelligent city flood control and drainage remote monitoring system and method
CN113408870A (en) Comprehensive safety control method for long-term environmental safety in slag disposal engineering
Al-Ansari et al. Flow of River Tigris and its effect on the bed sediment within Baghdad, Iraq
CN116561529A (en) Dam break-up state prediction method based on case statistics and break-up process
CN109736259B (en) Hydraulic method for determining the use range of compound flatlands and flood control safety positions in mountainous rivers
Yan et al. Characteristics and dynamic mechanism of rill erosion driven by extreme rainfall on karst plateau slopes, SW China
CN114818339B (en) A decision-making method for reservoir ecological desilting scheme based on multi-objective constraints
Wei et al. Evolution of historical sediment yield using check-dam systems as carriers: A case study in a restored agricultural catchment on the Loess Plateau, China
CN114239912A (en) A Prediction and Simulation Method of Compound Erosion Process
CN116737989B (en) Urban surface hydrologic monitoring network construction method based on video image
CN110889185A (en) Peak flow analysis method for small watershed and application thereof
Berg et al. Low flow variability in groundwater-fed streams
CN112464474B (en) Low-frequency debris flow early warning method based on risk evaluation
CN111999228B (en) A Method for Measuring and Calculating the Infiltration of New Urban Areas
Zhang et al. Effects of check-dams on sediment storage-release in Chabagou Watershed
CN115015504A (en) Surface water quality monitoring method based on surface land type
CN110188404B (en) Method for determining forest and grass vegetation coverage rate threshold capable of restraining sand production in river basin
Borji Sedimentation and sustainability of hydropower reservoirs: cases of Grand Ethiopian Renaissance Dam on the Blue Nile River in Ethiopia
CN112884279B (en) A rainstorm flood safety risk assessment method for a trench step land remediation project

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