资料介绍
ICS 65.020.20 CCS B 61
21
辽宁 省地 方标 准
DB21/T 4469—2026
主要针叶树种立木含碳量相容性模型及碳
计量
2026-07-07 发布 2026-08-07 实施
辽宁省市场监督管理局发 布
前言
本文件按照GB/T 1.1—2020《标准化工作导则第1部分：标准化文件的结构和起草规则》的规定起草。
请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。
本文件由辽宁省林业和草原局提出并归口。
本文件起草单位：沈阳农业大学、辽宁省林业科学研究院、北票市林业和草原局、岫岩满族自治县自然资源事物服务中心、本溪市林业和草原局、彰武县林业和草原局、宽甸满族自治县林业和草原局、国有清原满族自治县大孤家林场、本溪满族自治县清河城实验林场。
本文件主要起草人：高慧淋、殷有、孙云霞、吕海燕、郑雪敏、张晏恺、单晓宇、白国华、隋炎佐、史秀举、张亮、丛剑锋、刘冰、王小星、连永刚、赵庆喜、于翔、胡延辰、刘家腾、佟艺玟、倪铭岐、邹怡茜、刘芷岐、马雪诺、闫慧鹏、程海嘉、程淇、邰浩宇、马恩普。
本文件发布实施后，任何单位和个人如有问题和意见建议，均可以通过来电和来函等方式进行反馈，我们将及时答复并认真处理，根据实际情况依法进行评估及复审。
归口管理部门通讯地址：辽宁省林业和草原局(辽宁省沈阳市和平区太原北街2号)，联系电话： 024-23448927。
文件起草单位通讯地址： 沈阳农业大学 (辽宁省沈阳市沈河区东陵路120号)，联系电话： 024-88487150。
IV
主要针叶树种立木含碳量相容性模型及碳计量
1 范围
本文件规定了辽宁省落叶松( Larix gmelin ii)、红松( Pinus koraiensis)、樟子松( Pinus sylvestris)等8种针叶树种立木树干、树枝、树叶、树根含碳量相容性模型数据调查及处理方法、相容性模型研建方法、参数估计方法、不同针叶树种相容性模型的参数估计值。
本文件适用于主要针叶树种树干、树枝、树叶、树根含碳量的估算。
2 规范性引用文件
下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件;不注日期的引用文件，其最新版本(包括所有的修改单)适用于本文件。
GB/T 26424 森林资源规划设计调查技术规程
GB/T 43648 主要树种立木生物量模型与碳计量参数
LY/T 2188.2 森林资源数据采集技术规范第2部分：森林资源规划设计调查
LY/T 2258 立木生物量建模方法技术规程
3 术语和定义
下列术语和定义适用于本文件。
3.1
含水率 water content
树木各器官中所含水分的重量和干重之比。
[来源：LY/T 2259—2014，3.2，有修改]
3.2
生物量 biomass
某一时刻单位面积内实际存活的有机物质(干重)总量。
[来源：GB/T 26424—2010，3.17，有修改]
3.3
含碳量相容性模型 compatibility carbon stock model
以立木总含碳量为基础，总含碳量与树干各组成部分的含碳量之和相等的非线性相容性模型。
3.4
碳计量 carbon parameter
1
相容性模型参数的估计结果。
4 外业调查
4.1 标准地设置
选取具有代表性的林分设置标准地，标准地的形状可为矩形或圆形，面积≥0.06hm2。
4.2 解析木的选取
依据标准地调查的结果，所有树木按照胸径由大到小的顺序排序，按照等断面积划分为5个组，分别在第1组、第3组和第5组选取1株平均木作为样木进行解析。
4.3 树干鲜重测定及样品选择
将选取的优势木、平均木、劣势木伐倒去掉树枝后，从树根至树梢按照1m的长度将树干进行分段，上端不足1m的部分为梢头。称量所有区分段鲜重(kg)，计算得到全树干的鲜重。将树木伐倒去掉树枝后，从树干基部到树梢按照1m的长度将树干进行切割，不足1m的部分确定为梢头。为计算树干含水率，在树干每个区分段上端及0m、1.3m处截取厚度约为2cm的圆盘，在每个圆盘上截取楔形木块作为样品，记录鲜重(g)。
4.4 枝条鲜重测定及样品选择
将树干上所有活枝条取下，称量每一个带叶枝条的鲜重(g)。将树冠上、中、下层三层的标准枝条进行枝叶分离，分别记录树枝和树叶的鲜重。为计算树干含水率，将树冠每层的树枝取100g左右作为样品，记录鲜重(g)。
4.5 树叶鲜重测定及样品选择
将树冠等分为长度相等的上、中、下三层，每层选取1个大小中等、健康的枝条作为标准枝条。将标准枝条进行枝叶分离，分别记录树枝和树叶的鲜重(g)，将树冠每层的树叶取100g左右作为样品，记录鲜重(g)。
4.6 树根鲜重测定及样品选择
将树根全部挖出，将树根按照大于5cm、2cm～5cm和小于2cm三个等级分组，分别测定每个组树根的鲜重，将三个组的重量累计得到树根鲜重(kg)。将树根分为大于5cm、2cm～5cm和小于2cm三个等级，为计算树干含水率，每个等级的树根分别选取200g左右作为样品，对其进行保鲜处理后，带回实验室，记录鲜重(g)。
5 内业测定
5.1 生物量计算
用干燥箱将选取的各树木器官的所有样品在85℃的条件下烘干至恒重，记录各样品干重，由各器官样品的干重和鲜重的比例计算含水率，进而计算单木各器官的生物量。
5.2 含碳率测定
将烘干后的树木器官的样品进行研磨，将研磨后的样品重新在85℃的恒温条件下进行烘干处理。所有样品均分别取50mg左右的样品，进行样品含碳率的测定。
2
5.3 含碳量计算
利用各树木器官的生物量乘以各器官的含碳率，计算树木各器官的含碳量。
6 相容性模型构建
6.1 变量的选取
以树木胸径为单一变量构建树木器官一元相容性含碳量预估模型，以胸径、树高为变量构建二元相容性含碳量预估模型。
6.2 相容性模型构建
一元相容性模型：
式中：
Ct ——针叶树全树含碳量(kg) ;
Cs ——针叶树树干含碳量(kg) ;
Cb ——针叶树树枝含碳量(kg) ;
Cl ——针叶树树叶含碳量(kg) ;
Cr ——针叶树树根含碳量(kg) ;
D ——树木胸径(cm) ;
bc0、cc0、rc1、rc2、rc3、kc1、kc2、kc3 ——待估参数。
二元相容性模型：
3
H ——树高(m) ;
b0、c0、r1、r2、r3、k1、k2、k3 ——待估参数。
6.3 参数估计
所有针叶树种的一元、二元立木含碳量相容性模型的参数均采用R语言联立方程组的方法进行拟合。
6.4 模型的评价指标
以调整后的决定系数(Ra(2)dj)作为模型拟合优度的指标，公式如下：
yyi(i)————均(())(;;kg);
n ——样本个数;
p ——参数个数。
7 碳计量方程
基于6.3部分估计得出的相容性模型的参数结果，编制辽宁省主要针叶树种一元立木含碳量相容性模型计量参数用表，按照附录A执行;编制辽宁省主要针叶树种二元立木含碳量相容性模型计量参数用表，按照附录B执行。
4
附录 A
(规范性)
主要针叶树一元立木含碳量相容性模型计量参数
辽宁省主要针叶树种一元立木含碳量相容性模型计量参数按照表 A.1 执行。
表 A.1 主要针叶树一元立木含碳量相容性模型计量参数表
树种
bc0
cc0
rc1
rc2
rc3
kc1
kc2
kc3
落叶松
Larix gmelin ii
2.7255
0.0245
2.6580
3.0004
0.1766
-1.0755
-1.5083
0.1655
0.0245D2.7255
Cs= 1+2.6580D-1.0755+3.0004D-1.5083+0.1766D0.1655
0.0245×2.6580D(-1.0755+2.7255)
Cb= 1+2.6580D-1.0755+3.0004D-1.5083+0.1766D0.1655
0.0245×3.0004D(-1.5083+2.7255)
Cl= 1+2.6580D-1.0755+3.0004D-1.5083+0.1766D0.1655
0.0245×0.1766D(0.1655+2.7255)
1+2.6580D-1.0755+3.0004D-1.5083+0.1766D0.1655
Cr=
红松
Pinus koraiensis
2.5196
0.0384
0.0315
0.1146
0.0540
0.5499
-0.0785
0.5807
0.0384D2.5196
Cs= 1+0.0315D0.5499+0.1146D-0.0785+0.0540D0.5807
0.0384×0.0315D(0.5499+2.5196)
Cb= 1+0.0315D0.5499+0.1146D-0.0785+0.0540D0.5807
0.0384×0.1146D(-0.0785+2.5196)
Cl= 1+0.0315D0.5499+0.1146D-0.0785+0.0540D0.5807
0.0384×0.0540D(0.5807+2.5196)
C =
r 1+0.0315D0.5499+0.1146D-0.0785+0.0540D0.5807
樟子松
Pinus sylvestris
2.1529
0.1035
0.9827
2.7908
0.9955
-0.4280
-1.1112
-0.7936
0.1035D2.1529
Cs= 1+0.9827D-0.4280+2.7908D-1.1112+0.9955D-0.7936
0.1035×0.9827D(-0.4280+2.1529)
Cb= 1+0.9827D-0.4280+2.7908D-1.1112+0.9955D-0.7936
0.1035×2.7908D(-1.1112+2.1529)
Cl= 1+0.9827D-0.4280+2.7908D-1.1112+0.9955D-0.7936
0.1035×0.9955D(-0.7936+2.1529)
Cr=
1+0.9827D-0.4280+2.7908D-1.1112+0.9955D-0.7936
华山松
Pinus armandii
1.8452
0.3062
11.2183
0.7925
2.4683
-1.0215
-0.9827
-0.5069
0.3062D1.8452
Cs= 1+11.2183D-1.0215+0.7925D-0.9827+2.4683D-0.5069
0.3062 × 11.2183D(-1.0215+1.8452)
Cb= 1+11.2183D-1.0215+0.7925D-0.9827+2.4683D-0.5069
0.3062×0.7925D(-0.9827+1.8452)
Cl= 1+11.2183D-1.0215+0.7925D-0.9827+2.4683D-0.5069
0.3062×2.4683D(-0.5069+1.8452)
Cr=
1+11.2183D-1.0215+0.7925D-0.9827+2.4683D-0.5069
5
表 A.1 主要针叶树一元立木含碳量相容性模型计量参数表(续)
侧柏
Platy cladus
orientalis
1.4091
0.5809
2.7069
1.0892
1.5426
-0.8084
-0.8144
-0.3152
Cs=
Cb=
Cl=
Cr=
0.5809D1.4091
1+2.7069D-0.8084+1.0892D-0.8144+1.5426D-0.3152
0.5809×2.7069D(-0.8084+1.4091)
0.5809×1.0892D(-0.8144+1.4091)
0.5809×1.5426D(-0.3152+1.4091)
油松
Pinus tabuliformis
2.1607
0.0865
1.7199
0.3249
2.5066
-0.6478
-0.3377
-0.7331
0.0865D2.1607
1+1.7199D-0.6478+0.3249D-0.3377+2.5066D-0.7331
0.0865 × 1.7199D(-0.6478+2.1607)
0.0865×0.3249D(-0.3377+2.1607)
0.0865×2.5066D(-0.7331+2.1607)
冷杉
Abies fabri
1.6535
0.6209
10.5846
12.0344
3.7853
-1.1325
-1.2039
-1.0227
0.6209D1.6535
Cs= 1+10.5846D-1.1325+12.0344D-1.2039+3.7853D-1.0227
0.6209×10.5846D(-1.1325+1.6535)
Cb= 1+10.5846D-1.1325+12.0344D-1.2039+3.7853D-1.0227
0.6209×12.0344D(-1.2039+1.6535)
Cl= 1+10.5846D-1.1325+12.0344D-1.2039+3.7853D-1.0227
0.6209×3.7853D(-1.0227+1.6535)
1+10.5846D-1.1325+12.0344D-1.2039+3.7853D-1.0227
Cr=
彰武松
Pinus den siflora
2.0262
0.1205
0.1680
0.5852
0.4553
0.2188
-0.5224
-0.3939
0.1205D2.0262
Cs= 1+0.1680D0.2188+0.5852D-0.5224+0.4553D-0.3939
0.1205×0.1680D(0.2188+2.0262)
Cb= 1+0.1680D0.2188+0.5852D-0.5224+0.4553D-0.3939
0.1205×0.5852D(-0.5224+2.0262)
Cl= 1+0.1680D0.2188+0.5852D-0.5224+0.4553D-0.3939
0.1205×0.4553D(-0.3939+2.0262)
Cr= 1+0.1680D0.2188+0.5852D-0.5224+0.4553D-0.3939
6
附录 B
主要针叶树二元立木含碳量相容性模型计量参数
辽宁省主要针叶树种二元立木含碳量相容性模型计量参数参照表 B.1 执行。
表 B.1 主要针叶树二元立木含碳量相容性模型计量参数表
b0
c0
r1
r2
r3
k1
k2
k3
0.9665
0.0161
6.7484
5.9725
0.2144
-0.4651
-0.5845
-0.0352
0.0161×D2×0.9665H0.9665
Cs= 1+6.7484D-2×0.4651H-0.4651+5.9725D-2×0.5845H-0.5845+0.2144D-2×0.0352H-0.0352
0.0161×6.7484D(-2×0.4651+2×0.9665)H(-0.4651+0.9665)
Cb= 1+6.7484D-2×0.4651H-0.4651+5.9725D-2×0.5845H-0.5845+0.2144D-2×0.0352H-0.0352
0.0161×5.9725D(-2×0.5845+2×0.9665)H(-0.5845+0.9665)
Cl= 1+6.7484D-2×0.4651H-0.4651+5.9725D-2×0.5845H-0.5845+0.2144D-2×0.0352H-0.0352
0.0161×0.2144D(-2×0.0352+2×0.9665)H(-0.0352+0.9665)
Cr=
1+6.7484D-2×0.4651H-0.4651+5.9725D-2×0.5845H-0.5845+0.2144D-2×0.0352H-0.0352
1.0285
0.0096
0.0191
0.1151
0.2462
-0.0278
0.2400
0.0096×D2×1.0285H1.0285
Cs= 1+0.0191D-2×0.1365H0.2462+0.1151D-2×0.4468H-0.0278+0.0384D-2×0.2095H0.2400
0.0096×0.0191D(-2×0.1365+2×1.0285)H(0.2462+1.0285)
Cb= 1+0.0191D-2×0.1365H0.2462+0.1151D-2×0.4468H-0.0278+0.0384D-2×0.2095H0.2400
0.0096×0.1151D(-2×0.4468+2×1.0285)H(-0.0278+1.0285)
Cl= 1+0.0191D-2×0.1365H0.2462+0.1151D-2×0.4468H-0.0278+0.0384D-2×0.2095H0.2400
0.0096×0.0384D(-2×0.2095+2×1.0285)H(0.2400+1.0285)
C =
r 1+0.0191D-2×0.1365H0.2462+0.1151D-2×0.4468H-0.0278+0.0384D-2×0.2095H0.2400
0.8138
0.0783
1.0100
3.7782
1.2225
-0.1573
-0.4415
-0.3139
0.0783×D2×0.8138H0.8138
Cs= 1+1.0100D-2×0.1573H-0.1573+3.7782D-2×0.4415H-0.4415+1.2225D-2×0.3139H-0.3139
0.0783×1.0100D(-2×0.1573+2×0.8138)H(-0.1573+0.8138)
Cb= 1+1.0100D-2×0.1573H-0.1573+3.7782D-2×0.4415H-0.4415+1.2225D-2×0.3139H-0.3139
0.0783×3.7782D(-2×0.4415+2×0.8138)H(-0.4415+0.8138)
Cl= 1+1.0100D-2×0.1573H-0.1573+3.7782D-2×0.4415H-0.4415+1.2225D-2×0.3139H-0.3139
0.0783×1.2225D(-2×0.3139+2×0.8138)H(-0.3139+0.8138)
C =
r 1+1.0100D-2×0.1573H-0.1573+3.7782D-2×0.4415H-0.4415+1.2225D-2×0.3139H-0.3139
0.7120
0.1843
31.0080
1.8873
3.7082
-0.4831
-0.4503
-0.2272
0.1843×D2×0.7120H0.7120
Cs= 1+31.0080D-2×0.4831H-0.4831+1.8873D-2×0.4503H-0.4503+3.7083D-2×0.2272H-0.2272
0.1843×31.0080D(-2×0.4831+2×0.7120)H(-0.4831+0.7120)
Cb= 1+31.0080D-2×0.4831H-0.4831+1.8873D-2×0.4503H-0.4503+3.7083D-2×0.2272H-0.2272
0.1843×1.8873D(-2×0.4503+2×0.7120)H(-0.4503+0.7120)
Cl= 1+31.0080D-2×0.4831H-0.4831+1.8873D-2×0.4503H-0.4503+3.7083D-2×0.2272H-0.2272
0.1843×3.7083D(-2×0.2272+2×0.7120)H(-0.2272+0.7120)
C =
r 1+31.0080D-2×0.4831H-0.4831+1.8873D-2×0.4503H-0.4503+3.7083D-2×0.2272H-0.2272
表 B.1 主要针叶树二元立木含碳量相容性模型计量参数表(续)
0.5830
0.4129
3.6153
1.4116
1.9590
-0.3511
-0.3479
-0.1591
7
0.4129×D2×0.5830H0.5830
Cs= 1+3.6153D-2×0.3511H-0.3511+1.4116D-2×0.3479H-0.3479+1.9590D-2×0.1591H-0.1591
0.4129×3.6153D(-2×0.3511+2×0.5830)H(-0.3511+0.5830)
Cb= 1+3.6153D-2×0.3511H-0.3511+1.4116D-2×0.3479H-0.3479+1.9590D-2×0.1591H-0.1591
0.4129×1.4116D(-2×0.3479+2×0.5830)H(-0.3479+0.5830)
Cl= 1+3.6153D-2×0.3511H-0.3511+1.4116D-2×0.3479H-0.3479+1.9590D-2×0.1591H-0.1591
0.4129×1.9590D(-2×0.1591+2×0.5830)H(-0.1591+0.5830)
Cr= 1+3.6153D-2×0.3511H-0.3511+1.4116D-2×0.3479H-0.3479+1.9590D-2×0.1591H-0.1591
0.7702
0.1057
1.9121
0.3844
3.7514
-0.2518
-0.1455
-0.3244
0.1057×D2×0.7702H0.7702
Cs= 1+1.9121D-2×0.2518H-0.2518+0.3844D-2×0.1455H-0.1455+3.7514D-2×0.3244H-0.3244
0.1057×1.9121D(-2×0.2518+2×0.7702)H(-0.2518+0.7702)
Cb= 1+1.9121D-2×0.2518H-0.2518+0.3844D-2×0.1455H-0.1455+3.7514D-2×0.3244H-0.3244
0.1057×0.3844D(-2×0.1455+2×0.7702)H(-0.1455+0.7702)
Cl= 1+1.9121D-2×0.2518H-0.2518+0.3844D-2×0.1455H-0.1455+3.7514D-2×0.3244H-0.3244
0.1057×3.7514D(-2×0.3244+2×0.7702)H(-0.3244+0.7702)
1+1.9121D-2×0.2518H-0.2518+0.3844D-2×0.1455H-0.1455+3.7514D-2×0.3244H-0.3244
0.6469
0.3337
34.4468
33.9541
18.2061
-0.5265
-0.5341
-0.5364
0.3337×D2×0.6469H0.6469
Cs= 1+34.4468D-2×0.5265H-0.5265+33.9541D-2×0.5341H-0.5341+18.2061D-2×0.5364H-0.5364
0.3337×34.4468D(-2×0.5265+2×0.6469)H(-0.5265+0.6469)
Cb= 1+34.4468D-2×0.5265H-0.5265+33.9541D-2×0.5341H-0.5341+18.2061D-2×0.5364H-0.5364
0.3337×33.9541D(-2×0.5341+2×0.6469)H(-0.5341+0.6469)
Cl= 1+34.4468D-2×0.5265H-0.5265+33.9541D-2×0.5341H-0.5341+18.2061D-2×0.5364H-0.5364
0.3337×18.2061D(-2×0.5364+2×0.6469)H(-0.5364+0.6469)
Cr=
1+34.4468D-2×0.5265H-0.5265+33.9541D-2×0.5341H-0.5341+18.2061D-2×0.5364H-0.5364
0.8317
0.0573
0.1102
1.0938
0.9050
0.1294
-0.2686
-0.2282
0.0573×D2×0.8317H0.8317
Cs= 1+0.1102D2×0.1294H0.1294+1.0938D-2×0.2686H-0.2686+0.9050D-2×0.2282H-0.2282
0.0573×0.1102D(2×0.1294+2×0.8317)H(0.1294+0.8317)
Cb= 1+0.1102D2×0.1294H0.1294+1.0938D-2×0.2686H-0.2686+0.9050D-2×0.2282H-0.2282
0.0573×1.0938D(-2×0.2686+2×0.8317)H(-0.2686+0.8317)
Cl= 1+0.1102D2×0.1294H0.1294+1.0938D-2×0.2686H-0.2686+0.9050D-2×0.2282H-0.2282
0.0573×0.9050D(-2×0.2282+2×0.8317)H(-0.2282+0.8317)
r 1+0.1102D2×0.1294H0.1294+1.0938D-2×0.2686H-0.2686+0.9050D-2×0.2282H-0.2282
8
参考 文献
[1] 林学名词审定委员会.林学名词(第二版).科学出版社.2016.
[2] DB11/T 953—2024 林地碳汇计量监测技术规程
[3] LY/T 2259—2014 立木生物量建模样本采集技术规程
[4] GB/T 26424—2010 森林资源规划设计调查技术规程
9
