1、ICS65.020.20CCS B 6121辽宁省地方标准DB21/T 44692026主要针叶树种立木含碳量相容性模型及碳计量2026-07-07 发布2026-08-07 实施辽宁省市场监督管理局发 布DB21/T 44692026IV前言本文件按照GB/T 1.12020标准化工作导则 第1部分:标准化文件的结构和起草规则的规定起草。请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由辽宁省林业和草原局提出并归口。本文件起草单位:沈阳农业大学、辽宁省林业科学研究院、北票市林业和草原局、岫岩满族自治县自然资源事物服务中心、本溪市林业和草原局、彰武县林业和草原局
2、宽甸满族自治县林业和草原局、国有清原满族自治县大孤家林场、本溪满族自治县清河城实验林场。本文件主要起草人:高慧淋、殷有、孙云霞、吕海燕、郑雪敏、张晏恺、单晓宇、白国华、隋炎佐、史秀举、张亮、丛剑锋、刘冰、王小星、连永刚、赵庆喜、于翔、胡延辰、刘家腾、佟艺玟、倪铭岐、邹怡茜、刘芷岐、马雪诺、闫慧鹏、程海嘉、程淇、邰浩宇、马恩普。本文件发布实施后,任何单位和个人如有问题和意见建议,均可以通过来电和来函等方式进行反馈,我们将及时答复并认真处理,根据实际情况依法进行评估及复审。归口管理部门通讯地址:辽宁省林业和草原局(辽宁省沈阳市和平区太原北街2号),联系电话:024-23448927。文件起草单位
3、通讯地址:沈阳农业大学(辽宁省沈阳市沈河区东陵路120号),联系电话:024-88487150。DB21/T 446920261主要针叶树种立木含碳量相容性模型及碳计量1范围本文件规定了辽宁省落叶松(Larix gmelinii)、红松(Pinus koraiensis)、樟子松(Pinussylvestris)等8种针叶树种立木树干、树枝、树叶、树根含碳量相容性模型数据调查及处理方法、相容性模型研建方法、参数估计方法、不同针叶树种相容性模型的参数估计值。本文件适用于主要针叶树种树干、树枝、树叶、树根含碳量的估算。2规范性引用文件下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。
4、其中,注日期的引用文件,仅该日期对应的版本适用于本文件;不注日期的引用文件,其最新版本(包括所有的修改单)适用于本文件。GB/T 26424 森林资源规划设计调查技术规程GB/T 43648 主要树种立木生物量模型与碳计量参数LY/T 2188.2 森林资源数据采集技术规范第2部分:森林资源规划设计调查LY/T 2258 立木生物量建模方法技术规程3术语和定义下列术语和定义适用于本文件。3.1含水率 water content树木各器官中所含水分的重量和干重之比。来源:LY/T 22592014,3.2,有修改3.2生物量 biomass某一时刻单位面积内实际存活的有机物质(干重)总量。来源:
5、GB/T 264242010,3.17,有修改3.3含碳量相容性模型 compatibility carbon stock model以立木总含碳量为基础,总含碳量与树干各组成部分的含碳量之和相等的非线性相容性模型。3.4碳计量 carbon parameterDB21/T 446920262相容性模型参数的估计结果。4外业调查4.1标准地设置选取具有代表性的林分设置标准地,标准地的形状可为矩形或圆形,面积0.06hm2。4.2 解析木的选取依据标准地调查的结果,所有树木按照胸径由大到小的顺序排序,按照等断面积划分为5个组,分别在第1组、第3组和第5组选取1株平均木作为样木进行解析。4.3 树
6、干鲜重测定及样品选择将选取的优势木、平均木、劣势木伐倒去掉树枝后,从树根至树梢按照1m的长度将树干进行分段,上端不足1m的部分为梢头。称量所有区分段鲜重(kg),计算得到全树干的鲜重。将树木伐倒去掉树枝后,从树干基部到树梢按照1m的长度将树干进行切割,不足1m的部分确定为梢头。为计算树干含水率,在树干每个区分段上端及0m、1.3m处截取厚度约为2cm的圆盘,在每个圆盘上截取楔形木块作为样品,记录鲜重(g)。4.4 枝条鲜重测定及样品选择将树干上所有活枝条取下,称量每一个带叶枝条的鲜重(g)。将树冠上、中、下层三层的标准枝条进行枝叶分离,分别记录树枝和树叶的鲜重。为计算树干含水率,将树冠每层的树
7、枝取100g左右作为样品,记录鲜重(g)。4.5 树叶鲜重测定及样品选择将树冠等分为长度相等的上、中、下三层,每层选取1个大小中等、健康的枝条作为标准枝条。将标准枝条进行枝叶分离,分别记录树枝和树叶的鲜重(g),将树冠每层的树叶取100g左右作为样品,记录鲜重(g)。4.6 树根鲜重测定及样品选择将树根全部挖出,将树根按照大于5cm、2cm5cm和小于2cm三个等级分组,分别测定每个组树根的鲜重,将三个组的重量累计得到树根鲜重(kg)。将树根分为大于5cm、2cm5cm和小于2cm三个等级,为计算树干含水率,每个等级的树根分别选取200g左右作为样品,对其进行保鲜处理后,带回实验室,记录鲜重(
8、g)。5内业测定5.1 生物量计算用干燥箱将选取的各树木器官的所有样品在85的条件下烘干至恒重,记录各样品干重,由各器官样品的干重和鲜重的比例计算含水率,进而计算单木各器官的生物量。5.2 含碳率测定将烘干后的树木器官的样品进行研磨,将研磨后的样品重新在85的恒温条件下进行烘干处理。所有样品均分别取50mg左右的样品,进行样品含碳率的测定。DB21/T 4469202635.3 含碳量计算利用各树木器官的生物量乘以各器官的含碳率,计算树木各器官的含碳量。6相容性模型构建6.1 变量的选取以树木胸径为单一变量构建树木器官一元相容性含碳量预估模型,以胸径、树高为变量构建二元相容性含碳量预估模型。6
9、2 相容性模型构建一元相容性模型:式中:Ct针叶树全树含碳量(kg);Cs针叶树树干含碳量(kg);Cb针叶树树枝含碳量(kg);Cl针叶树树叶含碳量(kg);Cr针叶树树根含碳量(kg);D树木胸径(cm);bc0、cc0、rc1、rc2、rc3、kc1、kc2、kc3 待估参数。二元相容性模型:式中:Ct针叶树全树含碳量(kg);Cs针叶树树干含碳量(kg);Cb针叶树树枝含碳量(kg);Cl针叶树树叶含碳量(kg);DB21/T 446920264Cr针叶树树根含碳量(kg);D树木胸径(cm);H树高(m);b0、c0、r1、r2、r3、k1、k2、k3 待估参数。6.3 参数估计所
10、有针叶树种的一元、二元立木含碳量相容性模型的参数均采用R语言联立方程组的方法进行拟合。6.4 模型的评价指标以调整后的决定系数(Radj2)作为模型拟合优度的指标,公式如下:Radj2=1-n-1n-p-1 i=1N(yi-yi)i=1N(yi-y)式中:yi 实测值(kg);yi预估值(kg);y yi的平均值(kg);n 样本个数;p 参数个数。7碳计量方程基于6.3部分估计得出的相容性模型的参数结果,编制辽宁省主要针叶树种一元立木含碳量相容性模型计量参数用表,按照附录A执行;编制辽宁省主要针叶树种二元立木含碳量相容性模型计量参数用表,按照附录B执行。DB21/T 446920265附录A
11、规范性)主要针叶树一元立木含碳量相容性模型计量参数辽宁省主要针叶树种一元立木含碳量相容性模型计量参数按照表 A.1 执行。表 A.1 主要针叶树一元立木含碳量相容性模型计量参数表树种bc0cc0rc1rc2rc3kc1kc2kc3落叶松Larix gmelinii2.72550.02452.65803.00040.1766-1.0755-1.50830.1655Cs=0.0245D2.72551+2.6580D-1.0755+3.0004D-1.5083+0.1766D0.1655Cb=0.02452.6580D(-1.0755+2.7255)1+2.6580D-1.0755+3.0004D
12、1.5083+0.1766D0.1655Cl=0.02453.0004D(-1.5083+2.7255)1+2.6580D-1.0755+3.0004D-1.5083+0.1766D0.1655Cr=0.02450.1766D(0.1655+2.7255)1+2.6580D-1.0755+3.0004D-1.5083+0.1766D0.1655红松Pinus koraiensis2.51960.03840.03150.11460.05400.5499-0.07850.5807Cs=0.0384D2.51961+0.0315D0.5499+0.1146D-0.0785+0.0540D0.580
13、7Cb=0.03840.0315D(0.5499+2.5196)1+0.0315D0.5499+0.1146D-0.0785+0.0540D0.5807Cl=0.03840.1146D(-0.0785+2.5196)1+0.0315D0.5499+0.1146D-0.0785+0.0540D0.5807Cr=0.03840.0540D(0.5807+2.5196)1+0.0315D0.5499+0.1146D-0.0785+0.0540D0.5807樟子松Pinus sylvestris2.15290.10350.98272.79080.9955-0.4280-1.1112-0.7936Cs=
14、0.1035D2.15291+0.9827D-0.4280+2.7908D-1.1112+0.9955D-0.7936Cb=0.10350.9827D(-0.4280+2.1529)1+0.9827D-0.4280+2.7908D-1.1112+0.9955D-0.7936Cl=0.10352.7908D(-1.1112+2.1529)1+0.9827D-0.4280+2.7908D-1.1112+0.9955D-0.7936Cr=0.10350.9955D(-0.7936+2.1529)1+0.9827D-0.4280+2.7908D-1.1112+0.9955D-0.7936华山松Pinu
15、s armandii1.84520.306211.21830.79252.4683-1.0215-0.9827-0.5069Cs=0.3062D1.84521+11.2183D-1.0215+0.7925D-0.9827+2.4683D-0.5069Cb=0.306211.2183D(-1.0215+1.8452)1+11.2183D-1.0215+0.7925D-0.9827+2.4683D-0.5069Cl=0.30620.7925D(-0.9827+1.8452)1+11.2183D-1.0215+0.7925D-0.9827+2.4683D-0.5069Cr=0.30622.4683D
16、0.5069+1.8452)1+11.2183D-1.0215+0.7925D-0.9827+2.4683D-0.5069DB21/T 446920266表 A.1 主要针叶树一元立木含碳量相容性模型计量参数表(续)树种bc0cc0rc1rc2rc3kc1kc2kc3侧柏Platycladusorientalis1.40910.58092.70691.08921.5426-0.8084-0.8144-0.3152Cs=0.5809D1.40911+2.7069D-0.8084+1.0892D-0.8144+1.5426D-0.3152Cb=0.58092.7069D(-0.8084+1.4
17、091)1+2.7069D-0.8084+1.0892D-0.8144+1.5426D-0.3152Cl=0.58091.0892D(-0.8144+1.4091)1+2.7069D-0.8084+1.0892D-0.8144+1.5426D-0.3152Cr=0.58091.5426D(-0.3152+1.4091)1+2.7069D-0.8084+1.0892D-0.8144+1.5426D-0.3152油松Pinus tabuliformis2.16070.08651.71990.32492.5066-0.6478-0.3377-0.7331Cs=0.0865D2.16071+1.719
18、9D-0.6478+0.3249D-0.3377+2.5066D-0.7331Cb=0.08651.7199D(-0.6478+2.1607)1+1.7199D-0.6478+0.3249D-0.3377+2.5066D-0.7331Cl=0.08650.3249D(-0.3377+2.1607)1+1.7199D-0.6478+0.3249D-0.3377+2.5066D-0.7331Cr=0.08652.5066D(-0.7331+2.1607)1+1.7199D-0.6478+0.3249D-0.3377+2.5066D-0.7331冷杉Abies fabri1.65350.620910
19、584612.03443.7853-1.1325-1.2039-1.0227Cs=0.6209D1.65351+10.5846D-1.1325+12.0344D-1.2039+3.7853D-1.0227Cb=0.620910.5846D(-1.1325+1.6535)1+10.5846D-1.1325+12.0344D-1.2039+3.7853D-1.0227Cl=0.620912.0344D(-1.2039+1.6535)1+10.5846D-1.1325+12.0344D-1.2039+3.7853D-1.0227Cr=0.62093.7853D(-1.0227+1.6535)1+1
20、0.5846D-1.1325+12.0344D-1.2039+3.7853D-1.0227彰武松Pinus densiflora2.02620.12050.16800.58520.45530.2188-0.5224-0.3939Cs=0.1205D2.02621+0.1680D0.2188+0.5852D-0.5224+0.4553D-0.3939Cb=0.12050.1680D(0.2188+2.0262)1+0.1680D0.2188+0.5852D-0.5224+0.4553D-0.3939Cl=0.12050.5852D(-0.5224+2.0262)1+0.1680D0.2188+0
21、5852D-0.5224+0.4553D-0.3939Cr=0.12050.4553D(-0.3939+2.0262)1+0.1680D0.2188+0.5852D-0.5224+0.4553D-0.3939DB21/T 446920267附录B(规范性)主要针叶树二元立木含碳量相容性模型计量参数辽宁省主要针叶树种二元立木含碳量相容性模型计量参数参照表 B.1 执行。表 B.1 主要针叶树二元立木含碳量相容性模型计量参数表树种b0c0r1r2r3k1k2k3落叶松Larix gmelinii0.96650.01616.74845.97250.2144-0.4651-0.5845-0.0352
22、Cs=0.0161D20.9665H0.96651+6.7484D-20.4651H-0.4651+5.9725D-20.5845H-0.5845+0.2144D-20.0352H-0.0352Cb=0.01616.7484D(-20.4651+20.9665)H(-0.4651+0.9665)1+6.7484D-20.4651H-0.4651+5.9725D-20.5845H-0.5845+0.2144D-20.0352H-0.0352Cl=0.01615.9725D(-20.5845+20.9665)H(-0.5845+0.9665)1+6.7484D-20.4651H-0.4651+5.
23、9725D-20.5845H-0.5845+0.2144D-20.0352H-0.0352Cr=0.01610.2144D(-20.0352+20.9665)H(-0.0352+0.9665)1+6.7484D-20.4651H-0.4651+5.9725D-20.5845H-0.5845+0.2144D-20.0352H-0.0352红松Pinus koraiensis1.02850.00960.01910.11510.03840.2462-0.02780.2400Cs=0.0096D21.0285H1.02851+0.0191D-20.1365H0.2462+0.1151D-20.4468
24、H-0.0278+0.0384D-20.2095H0.2400Cb=0.00960.0191D(-20.1365+21.0285)H(0.2462+1.0285)1+0.0191D-20.1365H0.2462+0.1151D-20.4468H-0.0278+0.0384D-20.2095H0.2400Cl=0.00960.1151D(-20.4468+21.0285)H(-0.0278+1.0285)1+0.0191D-20.1365H0.2462+0.1151D-20.4468H-0.0278+0.0384D-20.2095H0.2400Cr=0.00960.0384D(-20.2095+
25、21.0285)H(0.2400+1.0285)1+0.0191D-20.1365H0.2462+0.1151D-20.4468H-0.0278+0.0384D-20.2095H0.2400樟子松Pinus sylvestris0.81380.07831.01003.77821.2225-0.1573-0.4415-0.3139Cs=0.0783D20.8138H0.81381+1.0100D-20.1573H-0.1573+3.7782D-20.4415H-0.4415+1.2225D-20.3139H-0.3139Cb=0.07831.0100D(-20.1573+20.8138)H(-0
26、1573+0.8138)1+1.0100D-20.1573H-0.1573+3.7782D-20.4415H-0.4415+1.2225D-20.3139H-0.3139Cl=0.07833.7782D(-20.4415+20.8138)H(-0.4415+0.8138)1+1.0100D-20.1573H-0.1573+3.7782D-20.4415H-0.4415+1.2225D-20.3139H-0.3139Cr=0.07831.2225D(-20.3139+20.8138)H(-0.3139+0.8138)1+1.0100D-20.1573H-0.1573+3.7782D-20.44
27、15H-0.4415+1.2225D-20.3139H-0.3139华山松Pinus armandii0.71200.184331.00801.88733.7082-0.4831-0.4503-0.2272Cs=0.1843D20.7120H0.71201+31.0080D-20.4831H-0.4831+1.8873D-20.4503H-0.4503+3.7083D-20.2272H-0.2272Cb=0.184331.0080D(-20.4831+20.7120)H(-0.4831+0.7120)1+31.0080D-20.4831H-0.4831+1.8873D-20.4503H-0.4
28、503+3.7083D-20.2272H-0.2272Cl=0.18431.8873D(-20.4503+20.7120)H(-0.4503+0.7120)1+31.0080D-20.4831H-0.4831+1.8873D-20.4503H-0.4503+3.7083D-20.2272H-0.2272Cr=0.18433.7083D(-20.2272+20.7120)H(-0.2272+0.7120)1+31.0080D-20.4831H-0.4831+1.8873D-20.4503H-0.4503+3.7083D-20.2272H-0.2272表 B.1 主要针叶树二元立木含碳量相容性模型
29、计量参数表(续)树种b0c0r1r2r3k1k2k3侧柏0.58300.41293.61531.41161.9590-0.3511-0.3479-0.1591DB21/T 446920268PlatycladusorientalisCs=0.4129D20.5830H0.58301+3.6153D-20.3511H-0.3511+1.4116D-20.3479H-0.3479+1.9590D-20.1591H-0.1591Cb=0.41293.6153D(-20.3511+20.5830)H(-0.3511+0.5830)1+3.6153D-20.3511H-0.3511+1.4116D-20
30、3479H-0.3479+1.9590D-20.1591H-0.1591Cl=0.41291.4116D(-20.3479+20.5830)H(-0.3479+0.5830)1+3.6153D-20.3511H-0.3511+1.4116D-20.3479H-0.3479+1.9590D-20.1591H-0.1591Cr=0.41291.9590D(-20.1591+20.5830)H(-0.1591+0.5830)1+3.6153D-20.3511H-0.3511+1.4116D-20.3479H-0.3479+1.9590D-20.1591H-0.1591油松Pinus tabulif
31、ormis0.77020.10571.91210.38443.7514-0.2518-0.1455-0.3244Cs=0.1057D20.7702H0.77021+1.9121D-20.2518H-0.2518+0.3844D-20.1455H-0.1455+3.7514D-20.3244H-0.3244Cb=0.10571.9121D(-20.2518+20.7702)H(-0.2518+0.7702)1+1.9121D-20.2518H-0.2518+0.3844D-20.1455H-0.1455+3.7514D-20.3244H-0.3244Cl=0.10570.3844D(-20.14
32、55+20.7702)H(-0.1455+0.7702)1+1.9121D-20.2518H-0.2518+0.3844D-20.1455H-0.1455+3.7514D-20.3244H-0.3244Cr=0.10573.7514D(-20.3244+20.7702)H(-0.3244+0.7702)1+1.9121D-20.2518H-0.2518+0.3844D-20.1455H-0.1455+3.7514D-20.3244H-0.3244冷杉Abies fabri0.64690.333734.446833.954118.2061-0.5265-0.5341-0.5364Cs=0.333
33、7D20.6469H0.64691+34.4468D-20.5265H-0.5265+33.9541D-20.5341H-0.5341+18.2061D-20.5364H-0.5364Cb=0.333734.4468D(-20.5265+20.6469)H(-0.5265+0.6469)1+34.4468D-20.5265H-0.5265+33.9541D-20.5341H-0.5341+18.2061D-20.5364H-0.5364Cl=0.333733.9541D(-20.5341+20.6469)H(-0.5341+0.6469)1+34.4468D-20.5265H-0.5265+3
34、3.9541D-20.5341H-0.5341+18.2061D-20.5364H-0.5364Cr=0.333718.2061D(-20.5364+20.6469)H(-0.5364+0.6469)1+34.4468D-20.5265H-0.5265+33.9541D-20.5341H-0.5341+18.2061D-20.5364H-0.5364彰武松Pinus densiflora0.83170.05730.11021.09380.90500.1294-0.2686-0.2282Cs=0.0573D20.8317H0.83171+0.1102D20.1294H0.1294+1.0938D
35、20.2686H-0.2686+0.9050D-20.2282H-0.2282Cb=0.05730.1102D(20.1294+20.8317)H(0.1294+0.8317)1+0.1102D20.1294H0.1294+1.0938D-20.2686H-0.2686+0.9050D-20.2282H-0.2282Cl=0.05731.0938D(-20.2686+20.8317)H(-0.2686+0.8317)1+0.1102D20.1294H0.1294+1.0938D-20.2686H-0.2686+0.9050D-20.2282H-0.2282Cr=0.05730.9050D(-20.2282+20.8317)H(-0.2282+0.8317)1+0.1102D20.1294H0.1294+1.0938D-20.2686H-0.2686+0.9050D-20.2282H-0.2282DB21/T 446920269参考文献1 林学名词审定委员会.林学名词(第二版).科学出版社.2016.2 DB11/T 9532024 林地碳汇计量监测技术规程3 LY/T 22592014 立木生物量建模样本采集技术规程4 GB/T 264242010 森林资源规划设计调查技术规程_