CN116881604A - 一种碳汇造林项目碳汇动态计量方法 - Google Patents
一种碳汇造林项目碳汇动态计量方法 Download PDFInfo
- Publication number
- CN116881604A CN116881604A CN202310203806.XA CN202310203806A CN116881604A CN 116881604 A CN116881604 A CN 116881604A CN 202310203806 A CN202310203806 A CN 202310203806A CN 116881604 A CN116881604 A CN 116881604A
- Authority
- CN
- China
- Prior art keywords
- carbon
- monitoring
- reserves
- area
- project
- 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.)
- Granted
Links
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 292
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 290
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000012544 monitoring process Methods 0.000 claims abstract description 190
- 238000005070 sampling Methods 0.000 claims abstract description 34
- 238000011835 investigation Methods 0.000 claims abstract description 15
- 238000011156 evaluation Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims description 48
- 239000002028 Biomass Substances 0.000 claims description 45
- 239000002023 wood Substances 0.000 claims description 44
- 238000005259 measurement Methods 0.000 claims description 27
- 238000013461 design Methods 0.000 claims description 18
- 230000008859 change Effects 0.000 claims description 14
- 238000012937 correction Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 7
- 230000009919 sequestration Effects 0.000 claims description 6
- 238000007689 inspection Methods 0.000 claims description 5
- 230000000737 periodic effect Effects 0.000 claims description 3
- 101000872083 Danio rerio Delta-like protein C Proteins 0.000 claims description 2
- 238000004891 communication Methods 0.000 description 44
- 238000010295 mobile communication Methods 0.000 description 13
- 241000894007 species Species 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 6
- 230000012010 growth Effects 0.000 description 6
- QVFWZNCVPCJQOP-UHFFFAOYSA-N chloralodol Chemical compound CC(O)(C)CC(C)OC(O)C(Cl)(Cl)Cl QVFWZNCVPCJQOP-UHFFFAOYSA-N 0.000 description 5
- 230000001419 dependent effect Effects 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 244000166124 Eucalyptus globulus Species 0.000 description 2
- 235000011611 Pinus yunnanensis Nutrition 0.000 description 2
- 241000018652 Pinus yunnanensis Species 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 2
- 210000000481 breast Anatomy 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000005431 greenhouse gas Substances 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 238000003825 pressing Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 238000004804 winding Methods 0.000 description 2
- 229910000838 Al alloy Inorganic materials 0.000 description 1
- 241000218691 Cupressaceae Species 0.000 description 1
- 240000007263 Pinus koraiensis Species 0.000 description 1
- 235000011615 Pinus koraiensis Nutrition 0.000 description 1
- QUWBSOKSBWAQER-UHFFFAOYSA-N [C].O=C=O Chemical compound [C].O=C=O QUWBSOKSBWAQER-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000013479 data entry Methods 0.000 description 1
- -1 dead falling matters Substances 0.000 description 1
- 230000034373 developmental growth involved in morphogenesis Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Y—INFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
- G16Y10/00—Economic sectors
- G16Y10/10—Forestry
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Business, Economics & Management (AREA)
- Economics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Human Resources & Organizations (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Health & Medical Sciences (AREA)
- Computational Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Algebra (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Forests & Forestry (AREA)
- Accounting & Taxation (AREA)
- Development Economics (AREA)
- Computing Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种碳汇造林项目碳汇动态计量方法，包括如下步骤：(1)物联网监测样地设计布置；(2)物联网样地自动数据采集；设定物联网样地数据采集频率，对采集数据完成解析和存储任务；(3)碳汇造林森林碳储量监测：依次进行监测区碳储量计算、监测区碳储量精度计算与检验以及监测区碳储量精度扣减校正；(4)碳汇价值评估及监测出数。本发明提供了一种碳汇造林项目碳汇动态计量方法，通过布设物联网监测样地，按固定周期性自动获取监测数据计算森林碳储量，对比基线情况计量项目碳汇价值，可大幅降低监测成本，减少调查错误，提升碳汇连续监测精度和监测时效性，对项目区碳汇实时监测。
Description
技术领域
本发明涉及物联网及地理空间数据采集、计算、处理领域，具体涉及一种碳汇造林项目碳汇动态计量方法。
背景技术
碳汇造林是以增加森林碳汇为主要目标，对造林和林木生长全过程实施碳汇计量和监测而进行的有特殊要求的项目活动。准确监测森林碳储量和计量碳汇价值是碳汇造林项目最基础和必不可少的环节，目前对碳汇计量方法分为两类，一类是调查生物量计算碳储量，另一类是利用微气象原理测定森林二氧化碳碳通量。现有技术中的碳汇造林项目中碳汇计量和森林碳储量监测存在以下缺点：
测量成本高，容易产生错测。固定样地每木调查需完成现地测量、拍照、挂牌、登记、计算流程才能获得碳汇量，需要专业调查人员实地测量，容易由于测量错误、数据录入和操作错误，导致数据不准确。
监测时效性差，数据为静态数据。由于现地调查工作量大、费用高，开展一次成本较高，监测时效性较差，无法准确掌握树木生长动态，而林木碳汇同林木生长密切相关并且动态变化。
碳储量变化量时间分辨率低。限于实测工作量，监测周期较长，对于碳汇量变化量在一段时间内被认定为是线性的，无法体现出生物量随季节、林分结构、气候条件的动态变化情况。
综上所述，亟需提供一种可大幅降低监测成本，减少调查错误，提升碳汇连续监测精度和监测时效性，对项目区碳汇实时监测的碳汇造林项目碳汇动态计量方法。
发明内容
本发明的目的是提供一种可大幅降低监测成本，减少调查错误，提升碳汇连续监测精度和监测时效性，对项目区碳汇实时监测的碳汇造林项目碳汇动态计量方法。
上述目的是通过如下技术方案实现：一种碳汇造林项目碳汇动态计量方法，包括如下步骤：
(1)物联网监测样地设计布置：依次进行项目边界监测、抽样设计与布点和物联网监测样地布设，用于完成碳汇造林项目区内监测样地地上生物量调查和周期性监测；
(2)物联网样地自动数据采集；设定物联网样地数据采集频率，对采集数据完成解析和存储任务；
(3)碳汇造林森林碳储量监测：
(3.1)监测区碳储量计算：根据步骤(2)中采集的物联网样地数据计算碳储量；
(3.2)监测区碳储量精度计算与检验：计算监测区项目总体样地平均碳储量及其方差，计算监测区平均单位面积碳储量的不确定性，依此对监测结果开展精度检验，对不符合精度要求数据查找原因并执行步骤(3.3)进行校正；
(3.3)监测区碳储量精度扣减校正；
(4)碳汇价值评估及监测出数：依次进行计量基线碳汇量、项目碳汇量计算以及项目碳汇价值计算，输出成果指标、统计数表以及本监测期成果数据库。
进一步的技术方案是，所述步骤(3.1)中包括如下步骤：
(3.1.1)根据采集的物联网样地数据计算样地所有活立木样木直径和高度，进而计算出样地的生物量和林木生物量碳含量；
(3.1.2)计算各样地碳储量和监测区内平均单位面积碳储量，样地碳储量为样地内活立木样木林木碳储量合计，样地平均碳储量为样地碳储量除以样地面积，计算公式如下：
C＝∑c
式中，C为样地碳储量；c为活立木样木的碳储量；为监测区平均单位面积碳储量；s为样地面积；
(3.1.3)计算监测区内总碳储量，为监测区内平均单位面积碳储量与监测区总面积的乘积，计算公式如下：
式中，Ctotal为监测区总碳储量；为监测区内平均单位面积碳储量；S为监测区总面积；
(3.1.4)计算监测区内的碳汇量，为监测末期碳储量同监测初期碳储量的差值，计算公式如下：
式中，ΔCtotal,t2,t1为监测区t1日至t2日的碳汇量；为第t2日的监测区碳储量；/>为第t1日的监测区碳储量。
进一步的技术方案是，所述步骤(3.2)中包括如下步骤：
(3.2.1)计算监测区项目总体样地平均碳储量及其方差，计算公式如下：
式中，为样地平均碳储量；C为样地碳储量；ci为第i个样地碳储量；/>为样地平均碳储量；X为项目总体平均单位面积碳储量估计值的方差；
(3.2.2)计算监测区平均单位面积碳储量的不确定性，计算公式如下：
式中，u为样地评价单位面积碳储量的不确定性；t为可靠性指数；X为项目总体平均单位面积碳储量估计值的方差；为样地平均碳储量；
(3.2.3)判断u是否满足抽样设定精度要求，若满足，当前计算监测区碳储量即为本期监测结果，若不满足，则执行步骤(3.3)。
进一步的技术方案是，所述步骤(3.3)中包括如下步骤：
(3.3.1)通过步骤(3.2.3)计算的u值和抽样设定精度，判断是否超过设定精度最大阈值，若是，则不能扣减校正，进行修正或重新开展样地抽样设计，增设或重新部署物联网监测样地，当前计算监测区碳储量不可以作为本期监测结果；若否，则执行步骤(3.3.2)；
(3.3.2)采取扣减系数，在当前计算监测区碳储量结果上进行系数扣减，扣减结果为本期监测结果。
进一步的技术方案是，所述步骤(4)的具体步骤如下：
(4.1)计量基线碳汇量；
(4.2)项目碳汇量计算：计算项目边界内总林木生物质碳储量及按监测频率计算的变化量；
(4.3)项目碳汇价值计算：监测区内碳储量变化量乘以固碳价格即为碳汇价值，计算公式如下：
式中，Utotal,t2,t1为监测区t1日至t2日的监测区碳汇价值；ΔCtotal,t2,t1为监测区t1日至t2日的碳汇量；为第t2日的监测区碳储量；/>为第t1日的监测区碳储量；Pc为固碳价格；
(4.4)输出成果指标及统计数表：
(4.5)输出本监测期成果数据库，并加时间戳进行归档存储。
进一步的技术方案是，所述步骤(1)的具体步骤如下：
(1.1)项目边界监测：使用GPS、遥感、地形图为辅助，以项目实际边界为准动态确定项目边界；
(1.2)抽样设计与布点：确定抽样设计的可靠性水平和抽样精度，计算项目总体样地数量并确定样地大小，其中，样地数量计算公式如下：
式中：n为项目边界内项目总体样地数量；t为可靠性指数；Y为预估生物质碳储量变动系数；E为抽样允许相对误差；
(1.3)物联网监测样地布设：对监测样地完成首次测量和测量设备布设，完成设备测量频率设定。
相比于现有技术，本发明具备如下优势：
(1)降低调查错误。每木调查时，对于树径的测量由树径测量传感器自动完成，无需测量人员读数，避免了错测，并由于对已经测量的立木绑定有树径测量传感器，减少了重测的发生，而且树径测量传感器标记了已测立木和未测立木，大幅减少漏测的可能性。
(2)减少数据记录和计算出错概率。对于连续监测，由于采用自动数据采集，不用再进行数据录入，计算通过软件复用执行，减少数据记录和计算出错的概率。
(3)直径连续监测准确性高。由于树径测量传感器一直绑定与固定位置，故每次测量时均是同一位置，不会产生测量位置不一致及围尺绕测误差。
(4)大幅降低连续监测的工作量。本发明对于连续监测不需要到达现地，即可稳定的获得活立木样木碳储量以及碳储量变化量数据，由服务器端直接进行计算和展示，大大降低了连续监测的工作量。
(5)监测时效性大幅提高。本发明中监测的频率根据任务需要和电池容量共同决定，由于采用低功耗设计，采用每天测量1次，每周回传一次工作可达10年，使用寿命显著延长，可以准确掌握树木生长动态和碳汇变化量。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
图1为本发明一种实施方式所涉及的碳汇造林项目碳汇动态计量方法的流程示意图；
图2为本发明一种实施方式所涉及的物联网监测样地设计布置的流程示意图；
图3为本发明一种实施方式所涉及的物联网监测样地布设的流程示意图；
图4为本发明一种实施方式所涉及的物联网样地自动数据采集的流程示意图；
图5为本发明一种实施方式所涉及的碳汇造林森林碳储量监测的流程示意图；
图6为本发明一种实施方式所涉及的监测区碳储量计算的流程示意图；
图7为本发明一种实施方式所涉及的监测区碳储量精度计算与检验的流程示意图；
图8为本发明一种实施方式所涉及的监测区碳储量精度校正的流程示意图；
图9为本发明一种实施方式所涉及的碳汇价值评估及监测出数的流程示意图。
具体实施方式
下面结合附图对本发明进行详细描述，本部分的描述仅是示范性和解释性，不应对本发明的保护范围有任何的限制作用。此外，本领域技术人员根据本文件的描述，可以对本文件中实施例中以及不同实施例中的特征进行相应组合。
本发明实施例如下，参照图1，一种碳汇造林项目碳汇动态计量方法，包括如下步骤：
(1)物联网监测样地设计布置：依次进行项目边界监测、抽样设计与布点和物联网监测样地布设，用于完成碳汇造林项目区内监测样地地上生物量调查和周期性监测；
物联网监测样地由树径测量传感器、通信中继装置、数据采集终端、北斗卫星及北斗指挥机、移动通信基站及通信服务器、计算服务器及网络服务器组成。
设置数据采集频率：数据采集终端同通信中继装置通信，数据采集终端通过使用WiFi、蓝牙、LORA等无线传输协议同通信中继装置连接，设置数据采集频率，确定通信中继装置和树径测量传感器的启动时间和启动时长。
录入树种、立木类型：树径测量传感器同数据采集终端通信，数据采集终端通过使用WiFi、蓝牙、LORA等无线传输协议同树径测量传感器连接，树径测量传感器获得并存储数据采集终端录入树种、立木类型。
时钟同步和获取唤醒时间：树径测量传感器同通信中继装置通信，树径测量传感器使用WiFi、蓝牙、LORA等无线传输协议同通信中继装置通信连接，同通信中继装置进行时钟同步，获得树径测量传感器的下一次启动时间和启动时长。
数据汇入：树径测量传感器同通信中继装置通信，在树径测量传感器进入睡眠状态前，将测量信息、录入每木调查信息、电池容量信息，使用WiFi、蓝牙、LORA等无线传输协议，传输汇入通信中继装置存储。
数据下载和检查：树径测量传感器同通信中继装置通信，结束每木调查后，再次通过使用WiFi、蓝牙、LORA等无线传输协议同通信中继装置连接，下载数据到数据采集终端进行查看和检查，查看是否存在缺项数据。
数据回传：通信中继装置数据同卫星/移动通信基站通信，通信中继装置的信号发射单元，如使用北斗短报文通信方式，通过卫星接收北斗短报文通信中继转置信息后，转发至北斗短报文通信接收机。如使用2/3/4/5G通信方式，通过移动通信基站设备接收通信中继转置信息后，使用移动运行商通信网络回传至通信服务器。
数据接收和存储：通信接收装置和服务器端通信，通信服务器或是北斗短报文通信接收机通过光纤网络将接收数据转换处理后，通过光纤网络存储进入服务器端网络服务器中。
如图2，具体步骤如下：
(1.1)项目边界监测：使用GPS、遥感、地形图为辅助，以项目实际边界为准动态确定项目边界；
(1.2)抽样设计与布点：确定抽样设计的可靠性水平和抽样精度，计算项目总体样地数量并确定样地大小，其中，样地数量计算公式如下：
式中：n为项目边界内项目总体样地数量；t为可靠性指数；Y为预估生物质碳储量变动系数；E为抽样允许相对误差；
《碳汇造林项目方法学》要求为90％可靠性水平下90％的抽样精度，样地面积为0.04-0.06公顷，本发明中样地大小建议为方形样地，面积1亩，即0.0667公顷。开展样地布设，采用随机起点，机械布点的系统抽样方案。其中机械布点是指间隔固定距离抽取设定数量的点位，点位的坐标作为样地的西南角点。
(1.3)物联网监测样地布设：对监测样地完成首次测量和测量设备布设，完成设备测量频率设定，测试设备运行正常，数据采集和传输稳定可靠。使用树径测量传感器代替人工围尺测量，使用数据采集终端数据采集代替人工记录，使移动通信中继装置或北斗短报文网关定时采集、传输树径测量传感器数据，实现样地数据自动、连续监测。如图3，物联网样地布设流程如下：
(1.3.1)到达布设样地，完成样地测设，使用数据采集终端测试样地信号类型，选择网关类型。如果有2/3/4/5G等移动通信信号，选择移动通信中继装置，如果没有2/3/4/5G等移动通信信号，选择北斗短报文网关。
(1.3.2)在样地中央位置选择一株样木布设铝合金支架，固定通信中继装置并进行信号测试，成功后保持通信中继装置开启状态。
(1.3.3)通过数据采集终端连接通信中继装置，设置数据采集频率，确定通信中继装置和树径测量传感器的下次自动启动时间和启动时长。
(1.3.4)选择样木直径测量位置，固定树径测量传感器于被测样木上。
(1.3.5)按下树径测量传感器启动按钮，指示灯亮起，使用数据采集终端连接树径测量传感器，保证显示连接树径测量传感器编码同树径测量传感器外壳标签编码一致。
(1.3.6)从树径测量传感器的拉线出口将拉绳拉出，绕样木一圈后扣入防拆拉线固定口。
(1.3.7)数据采集终端中查看样木直径是否有测量数值或是否明显偏离实际值，是，则按下树径测量传感器启动按钮，在数据采集终端清零后，重新拉动拉绳安装。
(1.3.8)点击连接树径测量传感器，录入树种、检尺类型。其中树种和检尺类型均按技术规定编码要求填写代码，例如树种为柏木，代码填写601，检尺类型为活立木，代码填写1。
(1.3.9)树径测量传感器同通信中继装置通信，传输测量数值至通信中继装置并进行时钟同步，以及获得树径测量传感器的下次启动时间和启动时长。
(1.3.10)树径测量传感器在无连接情况下，间隔固定时间后自动进入休眠状态，指示灯熄灭。
(1.3.11)重复步骤1.3.4-1.3.10，完成样地内全部样木的测量和树径测量传感器安装。
(1.3.12)通信中继装置通过卫星/移动通信基站将存储信息传回通信接收装置，并汇总进入平台层的网络服务器。对通过卫星回传的信息，传回北斗指挥机汇总进入网络服务器，对通过移动通信基站回传的信息，传回通信服务器汇总进入网络服务器。
(1.3.13)断开数据采集终端和通信中继装置连接，间隔固定时间后通信中继装置自动进入休眠状态，指示灯熄灭。
(2)物联网样地自动数据采集；设定物联网样地数据采集频率，对采集数据完成解析和存储任务；如图4，具体步骤如下：
(2.1)判断是否到达通信中继装置唤醒时间，否，等待到达唤醒时间。
(2.2)到达唤醒时间，北斗指挥机广播北斗短报文信息，北斗短报文网关接收短报文信息，获取北斗短报文网关和树径测量传感器唤醒时间、唤醒时长和唤醒频率。
(2.3)移动通信中继装置通过移动通信基站访问服务器端的网络服务器，获取移动通信中继装置和树径测量传感器唤醒时间、唤醒时长和唤醒频率。
(2.4)判断是否到达树径测量传感器唤醒时间，否，等待到达唤醒时间。
(2.5)到达唤醒时间，树径测量传感器测量样木直径，传输至通信中继装置并同通信中继装置进行时钟同步，获取唤醒时间、唤醒时长和唤醒频率，成功后进入睡眠状态。
(2.6)通信中继装置通过卫星/移动通信基站将通信中继装置存储信息传回通信接收装置，并汇总进入平台层的网络服务器。
(2.7)通信中继装置间隔固定时间后自动进入休眠状态，指示灯熄灭，结束本次流程。
(3)碳汇造林森林碳储量监测，联网监测样地布设流程包括监测区碳储量计算、监测区碳储量精度计算与检验、监测区碳储量精度扣减校正三个部分，如图5，具体流程如下：
(3.1)监测区碳储量计算：根据步骤(2)中采集的物联网样地的数据计算碳储量；
依据样地监测数据回传的实测样地内样木树种、活立木的胸径数据，计算各树种林木生物量，样地生物量、样地碳储量以及项目区生物量、碳储量。如图6，具体流程如下：
(3.1.1)依据物联网样地采集的样地内样木树种、胸径和立木类型数据，使用本区域当前树种的树高曲线，将直径作为自变量，树高作为因变量，计算样地所有活立木样木的高度。云南松可采用树高曲线经验方程如下：
h＝1.9293+0.753375d-0.017924d2
h—活立木样木的树高；d—活立木样木的直径；
(3.1.2)计算样地内活立木林木生物量，如果使用异速生长方程法，执行步骤(3.1.5)，使用生物量扩展因子法，执行步骤(3.1.3)。
(3.1.3)计算活立木样木材积，使用本区域二元材积表计算，将直径和树高作为自变量，材积作为因变量进行计算。黑龙江红松可采用二元材积方程如下：
v＝0.000063446d1.9760h0.8843
v—活立木样木的蓄积；h—活立木样木的树高；d—活立木样木的直径；
(3.1.4)使用生物量扩展因子法计算生物量，材积作为自变量，生物量作为因变量，查表获取木材密度、根冠比、生物量扩展因子三个缺省数据进行计算，计算公式如下：
w＝v×WD×BEF×(1+R)
w—活立木样木的生物量；v—活立木样木的材积；WD—活立木样木的木材密度；BEF—活立木样木的树干生物量转换到地上生物量的生物量扩展因子，无量纲；R—活立木样木的根冠比(即地下生物量与地上生物量的比值，无量纲)；
如广东桉树采用生物量扩展因子法的计算方程如下：
w＝v×0.578×1.263×(1+0.221)
WD＝0.578，BEF＝1.263，R＝0.221；
(3.1.5)使用异速生长方程计算生物量，直径和树高作为自变量，生物量作为因变量进行计算。云南松可采用生物量异速生长方程如下：
w＝0.211(d2h)1.0371
w—活立木样木的生物量；h—活立木样木的树高；d—活立木样木的直径；
(3.1.6)利用林木生物量碳含量将林木生物量换算为碳储量计算碳储量，再利用CO2与C的分子量(44/12)比将碳储量转换为二氧化碳当量，计算公式如下：
c—活立木样木的碳储量；w—活立木样木的生物量；CF—活立木样木的碳含量；
如广东桉树使用生物量扩展因子法的计算方程为如下：
CF＝0.5144；
(3.1.7)计算各样地碳储量和监测区内平均单位面积碳储量，样地碳储量是样地内活立木样木林木碳储量合计，样地平均碳储量为样地碳储量除以样地面积，计算公式如下：
C＝∑c
C—样地碳储量；c—样地活立木样木的碳储量；—监测区平均单位面积碳储量；s—样地面积；
(3.1.8)计算监测区内总碳储量，为监测区内平均单位面积碳储量与监测区总面积的乘积，计算公式如下：
Ctotal—监测区总碳储量；—监测区内平均单位面积碳储量；S—监测区总面积；
(3.1.9)计算监测区内碳储量变化量，即碳汇量，为监测末期碳储量同监测初期碳储量的差值，计算公式如下：
ΔCtotal,t2,t1—监测区t1日至t2日碳储量变化量，本发明中监测频率由人为设定，一般为每日测量，7天传输一次，碳储量变化量为日变化量；—第t2日的监测区碳储量；/>—第t1日的监测区碳储量。
(3.2)监测区碳储量精度计算与检验：计算监测区项目总体样地平均碳储量及其方差，计算监测区平均单位面积碳储量的不确定性，依此对监测结果开展精度检验，对不符合精度要求数据查找原因并执行步骤(3.3)进行校正；
如图7，具体步骤如下：
(3.2.1)计算监测区项目总体样地平均碳储量及其方差，计算公式如下：
—样地平均碳储量；C—样地碳储量；ci—第i个样地碳储量；/>—样地平均碳储量；X—项目总体平均单位面积碳储量估计值的方差；
(3.2.2)计算监测区平均单位面积碳储量的不确定性。
u—样地评价单位面积碳储量的不确定性；t—可靠性指数；X—项目总体平均单位面积碳储量估计值的方差；—样地平均碳储量；
(3.2.3)判断u是否满足抽样设定精度要求，满足，当前计算监测区碳储量即为本期监测结果，不满足，执行监测区碳储量精度扣减校正，校正结果作为本期监测结果。例如抽样精度设定为90％，则u小于或等于10％即满足精度要求，u大于10％即不满足精度要求，需进行精度扣减校正流程。
(3.3)监测区碳储量精度扣减校正，如图8，具体步骤如下：
(3.3.1)通过u值和抽样设定精度，判读是否超过设定精度最大阈值，是，不能扣减校正，需修正或是重新开展样地抽样设计，增设或是重新部署物联网监测样地，当前计算监测区碳储量不可以作为本期监测结果。
(3.3.2)否，在最大阈值内，采取扣减系数，在当前计算监测区碳储量结果上进行系数扣减，扣减结果为本期监测结果。
例如抽样精度设定为90％，30％即为精度最大阈值，u如果大于等于30％，不能扣减校正，需修正或是重新开展样地抽样设计，增设或是重新部署物联网监测样地，当前计算监测区碳储量不可以作为本期监测结果。u大于10小于20％且前后两期变化量大于0，扣减6％，u大于10小于20％且前后两期变化量小于0，扣减-6％。u大于20小于30％且前后两期变化量大于0，扣减11％，u大于20小于30％且前后两期变化量小于0，扣减-11％。
(4)碳汇价值评估及监测出数：依次进行计量基线碳汇量、项目碳汇量计算以及项目碳汇价值计算，输出成果指标、统计数表以及本监测期成果数据库。如图9，具体流程如下：
(1)计量基线碳汇量。碳汇造林项目为无林地上造林，基线情况下地上生物量、地下生物量、枯死木、枯落物、土壤有机质、和木产品碳库的变化率统一视为0，并在项目开始后不再进行监测。
(2)项目碳汇量计算，计算项目边界内总林木生物质碳储量及按监测频率计算的变化量。此时林木生物量的变化可以精确到周或是天，并且可以选取任意监测期间计算项目碳储量的变化率，由于本发明中基线情况碳库视为0，不考虑目边界内温室气体排放增加量和泄露，故碳储量的变化量即为项目碳汇量。
(3)项目碳汇价值计算，监测区内碳储量变化量乘以固碳价格即为碳汇价值。
Utotal,t2,t1—监测区t1日至t2日的监测区碳汇价值
ΔCtotal,t2,t1—监测区t1日至t2日碳储量变化量，本发明中监测频率由人为设定，一般为每日测量，7天传输一次，碳储量变化量为日变化量。
—第t2日的监测区碳储量
—第t1日的监测区碳储量
Pc—固碳价格
(4)按需打印输出各成果主要指标及统计数表。主要指标包括但不限于总体名称、可靠性，以及按监测周期更新的抽样精度、样地生物量、样地碳储量、监测区生物量、监测区碳储量以及按年/月/日计算的碳汇量和碳汇价值。统计数表包括但不限于按监测周期更新统计的各类土地面积统计表、各类森林碳储量统计表、各类土地面积动态表、总体特征数计算表，以及按年/月/日计算的各类森林按碳汇量表、各类森林碳汇价值表。
(5)输出本监测期成果数据库，并加时间戳进行归档存储。本监期成果数据库包括但不限于现地调查数据库、样地样木监测数据库、样地样木变化数据库、碳储量及碳汇价值数据库。
本发明针对碳汇造林项目碳汇计量和森林碳储量监测过程中普遍存在的测量成本高，监测时效性差，碳储量变化量被线性化处理的缺点，本发明提供了一种碳汇造林项目碳汇动态计量方法，通过布设物联网监测样地，按固定周期性自动获取监测数据计算森林碳储量，对比基线情况计量项目碳汇价值，可大幅降低监测成本，减少调查错误，提升碳汇连续监测精度和监测时效性，对项目区碳汇实时监测。
需说明的本发明中碳汇计量中只考虑林木生物质碳储量及变化，灌木、枯死木、枯落物、土壤有机碳和木产品碳库采用缺省值。本发明中不考虑项目边界内温室气体排放增加量的监测。本发明中碳汇造林项目活动无潜在泄露，识为0。
对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (6)
1.一种碳汇造林项目碳汇动态计量方法，其特征在于，包括如下步骤：
(1)物联网监测样地设计布置：依次进行项目边界监测、抽样设计与布点和物联网监测样地布设，用于完成碳汇造林项目区内监测样地地上生物量调查和周期性监测；
(2)物联网样地自动数据采集；设定物联网样地数据采集频率，对采集数据完成解析和存储任务；
(3)碳汇造林森林碳储量监测：
(3.1)监测区碳储量计算：根据步骤(2)中采集的物联网样地数据计算碳储量；
(3.2)监测区碳储量精度计算与检验：计算监测区项目总体样地平均碳储量及其方差，计算监测区平均单位面积碳储量的不确定性，依此对监测结果开展精度检验，对不符合精度要求数据查找原因并执行步骤(3.3)进行校正；
(3.3)监测区碳储量精度扣减校正；
(4)碳汇价值评估及监测出数：依次进行计量基线碳汇量、项目碳汇量计算以及项目碳汇价值计算，输出成果指标、统计数表以及本监测期成果数据库。
2.根据权利要求1所述的碳汇造林项目碳汇动态计量方法，其特征在于，所述步骤(3.1)中包括如下步骤：
(3.1.1)根据采集的物联网样地数据计算样地所有活立木样木直径和高度，进而计算出样地的生物量和林木生物量碳含量；
(3.1.2)计算各样地碳储量和监测区内平均单位面积碳储量，样地碳储量为样地内活立木样木林木碳储量合计，样地平均碳储量为样地碳储量除以样地面积，计算公式如下：
C＝∑c
式中，C为样地碳储量；c为活立木样木的碳储量；为监测区平均单位面积碳储量；s为样地面积；
(3.1.3)计算监测区内总碳储量，为监测区内平均单位面积碳储量与监测区总面积的乘积，计算公式如下：
式中，Ctotal为监测区总碳储量；为监测区内平均单位面积碳储量；S为监测区总面积；
(3.1.4)计算监测区内的碳汇量，为监测末期碳储量同监测初期碳储量的差值，计算公式如下：
式中，ΔCtotal,t2,t1为监测区t1日至t2日的碳汇量；为第t2日的监测区碳储量；为第t1日的监测区碳储量。
3.根据权利要求2所述的碳汇造林项目碳汇动态计量方法，其特征在于，所述步骤(3.2)中包括如下步骤：
(3.2.1)计算监测区项目总体样地平均碳储量及其方差，计算公式如下：
式中，为样地平均碳储量；C为样地碳储量；ci为第i个样地碳储量；/>为样地平均碳储量；X为项目总体平均单位面积碳储量估计值的方差；
(3.2.2)计算监测区平均单位面积碳储量的不确定性，计算公式如下：
式中，u为样地评价单位面积碳储量的不确定性；t为可靠性指数；X为项目总体平均单位面积碳储量估计值的方差；为样地平均碳储量；
(3.2.3)判断u是否满足抽样设定精度要求，若满足，当前计算监测区碳储量即为本期监测结果，若不满足，则执行步骤(3.3)。
4.根据权利要求3所述的碳汇造林项目碳汇动态计量方法，其特征在于，所述步骤(3.3)中包括如下步骤：
(3.3.1)通过步骤(3.2.3)计算的u值和抽样设定精度，判断是否超过设定精度最大阈值，若是，则不能扣减校正，进行修正或重新开展样地抽样设计，增设或重新部署物联网监测样地，当前计算监测区碳储量不可以作为本期监测结果；若否，则执行步骤(3.3.2)；
(3.3.2)采取扣减系数，在当前计算监测区碳储量结果上进行系数扣减，扣减结果为本期监测结果。
5.根据权利要求1～4任意一项所述的碳汇造林项目碳汇动态计量方法，其特征在于，所述步骤(4)的具体步骤如下：
(4.1)计量基线碳汇量；
(4.2)项目碳汇量计算：计算项目边界内总林木生物质碳储量及按监测频率计算的变化量；
(4.3)项目碳汇价值计算：监测区内碳储量变化量乘以固碳价格即为碳汇价值，计算公式如下：
式中，Utotal,t2,t1为监测区t1日至t2日的监测区碳汇价值；ΔCtotal,t2,t1为监测区t1日至t2日的碳汇量；为第t2日的监测区碳储量；/>为第t1日的监测区碳储量；Pc为固碳价格；
(4.4)输出成果指标及统计数表：
(4.5)输出本监测期成果数据库，并加时间戳进行归档存储。
6.根据权利要求5所述的碳汇造林项目碳汇动态计量方法，其特征在于，所述步骤(1)的具体步骤如下：
(1.1)项目边界监测：使用GPS、遥感、地形图为辅助，以项目实际边界为准动态确定项目边界；
(1.2)抽样设计与布点：确定抽样设计的可靠性水平和抽样精度，计算项目总体样地数量并确定样地大小，其中，样地数量计算公式如下：
式中：n为项目边界内项目总体样地数量；t为可靠性指数；Y为预估生物质碳储量变动系数；E为抽样允许相对误差；
(1.3)物联网监测样地布设：对监测样地完成首次测量和测量设备布设，完成设备测量频率设定。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310203806.XA CN116881604B (zh) | 2023-03-06 | 2023-03-06 | 一种碳汇造林项目碳汇动态计量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310203806.XA CN116881604B (zh) | 2023-03-06 | 2023-03-06 | 一种碳汇造林项目碳汇动态计量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116881604A true CN116881604A (zh) | 2023-10-13 |
| CN116881604B CN116881604B (zh) | 2024-03-08 |
Family
ID=88259205
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310203806.XA Active CN116881604B (zh) | 2023-03-06 | 2023-03-06 | 一种碳汇造林项目碳汇动态计量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116881604B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117789045A (zh) * | 2023-11-07 | 2024-03-29 | 国网宁夏电力有限公司电力科学研究院 | 一种构建科学的碳汇模型的碳汇监测方法、系统及装置 |
| CN117871792A (zh) * | 2024-03-13 | 2024-04-12 | 河北省建筑科学研究院有限公司 | 用于对园区绿植碳汇的动态监测方法及系统 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105225160A (zh) * | 2015-10-12 | 2016-01-06 | 温州科技职业学院 | 大型藻类渔业碳汇计量方法 |
| IN201911004661A (zh) * | 2019-02-06 | 2019-04-05 | ||
| CN110032611A (zh) * | 2019-04-03 | 2019-07-19 | 海南省林业科学研究所 | 一种碳汇计量监测方法 |
| CN112287287A (zh) * | 2020-11-06 | 2021-01-29 | 东北林业大学 | 一种林木碳汇量的测定方法、系统及装置 |
| CN114894254A (zh) * | 2022-07-13 | 2022-08-12 | 四川省林业和草原调查规划院(四川省林业和草原生态环境监测中心) | 一种单株木碳汇动态计量方法 |
| US20220374912A1 (en) * | 2021-05-16 | 2022-11-24 | Ning Zeng | Method and system for wood harvest and storage, carbon sequestration and carbon management |
| CN115700370A (zh) * | 2021-07-23 | 2023-02-07 | 久瓴(上海)智能科技有限公司 | 碳储量计算方法、装置、电子设备和存储介质 |
-
2023
- 2023-03-06 CN CN202310203806.XA patent/CN116881604B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105225160A (zh) * | 2015-10-12 | 2016-01-06 | 温州科技职业学院 | 大型藻类渔业碳汇计量方法 |
| IN201911004661A (zh) * | 2019-02-06 | 2019-04-05 | ||
| CN110032611A (zh) * | 2019-04-03 | 2019-07-19 | 海南省林业科学研究所 | 一种碳汇计量监测方法 |
| CN112287287A (zh) * | 2020-11-06 | 2021-01-29 | 东北林业大学 | 一种林木碳汇量的测定方法、系统及装置 |
| US20220374912A1 (en) * | 2021-05-16 | 2022-11-24 | Ning Zeng | Method and system for wood harvest and storage, carbon sequestration and carbon management |
| CN115700370A (zh) * | 2021-07-23 | 2023-02-07 | 久瓴(上海)智能科技有限公司 | 碳储量计算方法、装置、电子设备和存储介质 |
| CN114894254A (zh) * | 2022-07-13 | 2022-08-12 | 四川省林业和草原调查规划院(四川省林业和草原生态环境监测中心) | 一种单株木碳汇动态计量方法 |
Non-Patent Citations (1)
| Title |
|---|
| 刘欢;武曙红;于天飞;: "森林保护碳汇项目方法学研究", 世界林业研究, no. 05 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117789045A (zh) * | 2023-11-07 | 2024-03-29 | 国网宁夏电力有限公司电力科学研究院 | 一种构建科学的碳汇模型的碳汇监测方法、系统及装置 |
| CN117871792A (zh) * | 2024-03-13 | 2024-04-12 | 河北省建筑科学研究院有限公司 | 用于对园区绿植碳汇的动态监测方法及系统 |
| CN117871792B (zh) * | 2024-03-13 | 2024-05-14 | 河北省建筑科学研究院有限公司 | 用于对园区绿植碳汇的动态监测方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116881604B (zh) | 2024-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116881604B (zh) | 一种碳汇造林项目碳汇动态计量方法 | |
| CN116070080B (zh) | 一种森林碳储量与碳汇价值监测系统及动态评估方法 | |
| CN114894254B (zh) | 一种单株木碳汇动态计量方法 | |
| US11961150B2 (en) | National forest resources continuous inventory cloud platform and sample plot monitoring method | |
| US11869023B1 (en) | Continuous monitoring method and system for forest stock and execution method therefor | |
| König-Langlo et al. | The baseline surface radiation network and its World radiation monitoring centre at the Alfred Wegener Institute | |
| Shah et al. | Testing the near-field Gaussian plume inversion flux quantification technique using unmanned aerial vehicle sampling | |
| US11221323B2 (en) | System and method for instantaneously determining uniform distribution of water, salinity, conductivity, temperature and other conditions in soil | |
| CN116071419B (zh) | 一种碳汇造林和森林经营碳汇项目碳汇动态计量方法 | |
| Conil et al. | Continuous atmospheric CO 2, CH 4 and CO measurements at the Observatoire Pérenne de l'Environnement (OPE) station in France from 2011 to 2018 | |
| US11913928B2 (en) | System for estimating amount of absorbed carbon dioxide and method for estimating amount of absorbed carbon dioxide | |
| CN113051653B (zh) | 一种基于多维度数据分析的城市规划道路修建施工评估管理系统 | |
| CN101158576A (zh) | 内置程序式电子经纬仪角规自动立体测树技术 | |
| CN104849219A (zh) | 一种作物氮素营养诊断仪及方法 | |
| CN107680148A (zh) | 一种基于监测数据实时生成等值线图的方法 | |
| JP5826120B2 (ja) | 歩容計測装置、方法及びプログラム | |
| CN115879835B (zh) | 一种国家储备林建设成效监测方法、系统及云平台 | |
| CN111766613B (zh) | Gps功能测试方法及其装置、存储介质 | |
| CN116307868B (zh) | 一种森林生态系统服务功能动态评估方法、应用系统 | |
| CN105717517B (zh) | 一种车载北斗多模gnss高精度道路基础数据采集方法 | |
| JP7511721B2 (ja) | 単木炭素吸収源の動的計量方法 | |
| CN214200311U (zh) | 一种加油机检定计量罐 | |
| CN110824529A (zh) | 一种用于环境γ剂量率测量的便携式智能装置及测量方法 | |
| CN117114174A (zh) | 一种联合物联网概念土木工程测量的优化方法 | |
| KR102616969B1 (ko) | 인공지능 기반 총질소 및 총인 측정용 스마트 센서 시스템 |
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 |