Aquaculture wastewater disposal through dual microbial fuel cell constructed by coupling anaerobic ammonia oxidation sludge and chlorella
-
摘要:目的
利用厌氧氨氧化污泥和小球藻耦合构建双生物微生物燃料电池(DMFC)，以期在产电的同时有效去除养殖废水中的高浓度氮磷和有机物，为养殖废水处理和生物能源制备提供新途径。
方法将厌氧氨氧化污泥和模拟废水以体积比1∶4混合后投至阳极室，小球藻和BG-11培养基以体积比1∶4混合后加入阴极室，运行24 h后阳极开始出水，并将阳极出水加入阴极室中，保持恒温并设置光暗周期模式，监测电压及废水中氨氮、NO2−−N、NO3−−N、总磷质量浓度，可溶性化学需氧量(SCOD)和pH。
结果厌氧氨氧化污泥和小球藻之间存在明显的电子交换，系统产电性能良好，DMFC稳定运行阶段正负峰值电压分别为45、−125 mV。氨氮质量浓度由1 588.97 mg/L降至602.75 mg/L，去除率达62.07%；NO2−−N、总磷和SCOD去除率分别为88.62%、54.45%和63.72%。系统在pH为9.5的条件下能稳定运行。
结论本研究成功构建了一个厌氧氨氧化污泥耦合小球藻的DMFC系统，该系统不仅能有效脱氮、除磷、降低有机物浓度，还能持续产电，是协同解决环境与能源问题的有效途径。
Abstract:ObjectiveTo construct the dual microbial fuel cell (DMFC) by coupling anaerobic ammonia oxidation sludge and chlorella, hoping to remove the high concentration of nitrogen, phosphorus and organic matter in aquaculture wastewater while generating electricity, and provide a new way for aquaculture wastewater disposal and bioenergy preparation.
MethodThe anaerobic ammonia oxidation sludge was mixed with simulated wastewater at the volume ratio of 1∶4 and added to the anode chamber, while chlorella was mixed with BG-11 medium at the volume ratio of 1∶4 and added to the cathode chamber. After running for 24 h, the anode chamber began to discharge water, which was then added to the cathode chamber. The system was maintained at a constant temperature and a light-dark cycle mode. The voltage, mass concentrations of ammonia nitrogen/NO2−-N/NO3−-N/total phosphorus, soluble chemical oxygen demand (SCOD) and pH of wastewater were monitored during the system running.
ResultThere was an obvious electron exchange between the anaerobic ammonia oxidation sludge and chlorella, and the system had good electricity generation performance, the positive and negative peak voltages of DMFC in the stable operation stage were 45 and −125 mV respectively. The ammonia nitrogen mass concentration decreased from 1 588.97 mg/L to 602.75 mg/L, with the removal rate of 62.07%; The removal rates of NO2−-N, total phosphorus and SCOD were 88.62%, 54.45% and 63.72% respectively. The system operated stably at pH 9.5.
ConclusionThis study successfully establishes a DMFC system constructed by coupling anaerobic ammonia oxidation sludge and chlorella. The system not only effectively removes nitrogen, phosphorus and reduces organic matter concentration, but also continuously generates electricity, providing an effective way to address environmental and energy issues simultaneously.
-
森林生物量指的是森林范围内生物所含有机物的量，常用干质量表示。作为重要的生态学指标，森林生物量直接反映森林的生产力与碳汇能力，对于评价森林生态系统的结构和功能具有重要作用[1]。传统森林生物量的获取方法通常采用直接收获法，主要有皆伐法[2]、标准木法[3]和相对生长法[4]3种方式，通过对样地内树木伐倒、烘干以及称质量等步骤得到单株树木各部分的生物量，进而得到单株树木生物量和样地生物量。这些方法获取数据存在工作量大、费时费力等不足，且对树木造成不可逆转的损伤。因此如何通过有效的途径快速、无损地获取森林生物量一直是生态学领域的研究热点[5]。
无人机激光雷达(Unmanned aerial vehicle laser scanning，UAV-LS)是一种新兴的遥感技术，能够通过高密度的点云数据获取水平以及垂直结构的森林参数[6]，因其高效无损、便捷灵活等优点广泛应用于林业中。目前，基于UAV-LS数据进行森林生物量估测的方法大多依赖于异速生长模型[7]，主要是将激光雷达数据中提取的变量与生长模型结合，从而完成生物量估测。其中，将UAV-LS数据衍生的变量(高度百分位数、强度百分位数等)与实测树高、胸径等参数相结合，利用机器学习方法进行地上生物量估算已取得良好效果[8]，并且，这一方法进一步拓展到基于机载激光雷达冠层高度与信号强度等变量的结合，构建生物量估测模型[9]。相较于传统回归模型，机器学习方法因其不需要描述预测关系和响应变量而被越来越多地运用在森林清查中[10]，特别是在预测碳储量、生物量等方面[11]。除此之外，也有利用UAV-LS提取的单木参数(树高、胸径、冠幅等)在非线性CAR(Constant allometric ratio)模型的基础上进行生物量的估测[12]。
UAV-LS所提取的参数信息已被证实在生物量估测上具有较高的可信度[13]，但是，利用树木定量结构模型(Quantitative structure model，QSM)从UAV-LS数据中提取参数估测生物量则鲜有报道。QSM是利用三维模型定量描述和理解树木结构与生长规律的工具，已广泛应用在地基激光雷达(Terrestrial laser scanning，TLS)中，通过TLS点云数据建模并提取树木三维结构参数，进而估算材积[14]、蓄积量[15]、生物量[16]等森林量化指标。同样，QSM利用无人机激光点云进行生物量估测的应用潜力已得到证实[17]，特别是在考虑到树木结构和生物量构成组分时，无人机QSM可以成为森林管理活动和决策过程中有效的替代方案[18]。但是，如果数据集的质量太差、点云密度不高时，建模精度便无法保证，会出现生物量被低估的现象[19]。
因此，本研究在获取高密度UAV-LS数据的基础上，基于落叶阔叶林树种进行QSM的构建。通过三维模型提取树木参数，并进一步利用随机森林(Random forest，RF)、多层感知机(Multi-layer perceptron，MLP)和支持向量机(Support vector machine，SVM) 3种机器学习方法进行单木生物量估测，筛选最优单木生物量估测模型，为复杂森林环境下生物量估算提供科学参考。
1. 材料与方法
1.1 研究区概况
研究区位于安徽省六安市金寨县天马国家级自然保护区。保护区地处大别山腹地，生物资源丰富，是北亚热带向暖温带的过渡地域，地理坐标介于115°20′~115°50′E，31°10′~31°20′N，海拔800~
1729 m。保护区属于北亚热带湿润季风气候类型，四季分明，年平均气温13.8 ℃，年平均降雨量1400 mm。保护区总面积约28913.7 hm2，其中核心区5745.3 hm2，缓冲区8118.6 hm2，实验区15049.8 hm2。研究区植被类型属暖温带落叶阔叶林向亚热带常绿阔叶林过渡型，植被垂直分布明显，落叶阔叶林属天然次生林，主要有茅栗Castanea seguinii、短柄枹栎Quercus glandulifera、锐齿槲栎Quercus aliena var. acutiserrata等栎类以及化香Platycarya strobilacea、大别山山核桃Carya dabieshanensis、领春木Euptelea pleiosperma等其他硬阔类树种。样地内最高海拔1103 m，最低海拔1055 m，高度差达48 m，坡度主要集中在22°~36°。1.2 地面样地数据
地面数据的获取来源于2023年3月对落叶阔叶林样地的调查。首先在研究区布设75 m×75 m的大样地，划分为9个25 m×25 m的小样地，利用全站仪获取样地4个至角点坐标，并确定界外木。对样地内胸径(DBH)≥5 cm的乔木进行每木检尺，并记录树种信息。使用TruePulse200激光测距仪测量树高，用胸径卷尺测量胸径，用徕卡DistoD510激光测距仪测量树冠垂直投影距离。同时，将树的基部位置作为定位点，利用华测i86RTK(精度15 mm)记录每株树准确的位置信息，定位点的空间参照采用WGS84坐标系和UTM投影。样地基本概况见表1，除枯死木外，共统计样地内563株阔叶树种。根据相关研究并结合样地实际情况[20]，按胸径大小将样地树木划分为4个生活史阶段，分别为小树(5 cm≤DBH<7.5 cm，176株)、中树(7.5 cm≤DBH<22.5 cm，234株)、大树(22.5 cm≤DBH<47.5 cm，140株)以及老树(DBH≥47.5 cm，13株)。在忽略小树后，样地内主要优势树种为茅栗(100株)、短柄枹栎(86株)与化香(63株)。
表 1 样地基本情况Table 1. Basic information of sample plots统计值
Statistical value胸径/cm
DBH树高/m
Tree height平均冠幅/m
Average crown
breadth最大值 Max. 65.50 21.70 12.00 最小值 Min. 5.00 2.00 0.50 平均值 Average 17.03 10.16 3.17 标准差 SD 12.12 4.82 1.67 1.3 UAV-LS数据及预处理
1.3.1 Lidar数据采集
利用大疆M300RTK无人机搭载禅思L1传感器(水平精度：150 m测量距离下误差为5 cm，垂直精度：150 m测量距离下误差为4 cm)，同步获取样地近地面的激光雷达数据。在仿地飞行60 m的高度下进行作业，飞行速度为5 m/s，回波模式为三回波并重复扫描，激光旁向重叠率为60%，点云密度达到631 m−2。点云采集过程中，飞行高度越低、速度越慢，点云密度就越大，高密度点云可以更好地描述树木的结构特征。本试验使用的无人机具备载波相位差分技术(Real-time kinematic，RTK)，能提供精确的位置信息，其坐标系和投影设置与华测i86RTK保持一致，航线设置如图1所示。
1.3.2 归一化处理
通过开源软件CloudCompare (https://www.cloudcompare.org)对Las格式点云进行去噪处理，再利用CSF(Cloth simulation filter)算法[21]对去噪后的点云进行滤波，分离出地面点与非地面点。最后，将地面点转换为分辨率0.1 m大小的栅格数据，使用拉普拉斯算子对栅格平滑后计算非地面点到该栅格的距离，并将计算的距离作为非地面点的高度，从而完成归一化(图2)。
1.3.3 点云分割
使用LiDAR360软件中集成的相对最短路径算法(Comparative shortest-path algorithm，CSP)对点云进行分割，该算法借鉴了生态学理论，采用自下而上的方式，从点云主干进行搜索，并通过分支到主干的最近距离来判断所属关系，从而完成点云分割。这种方法在落叶期的阔叶树中表现出良好的分割效果[22]。本文研究数据采集时期正是树木的落叶期，枝干点云清晰，符合该算法的应用场景。点云分割精度采用召回率(t)、正确率(p)以及t和p的调和平均值(F)进行评价，三者取值范围均为0到1，数值越接近1表示分割效果越好。计算公式如下：
$$ t = \dfrac{{{\mathrm{TP}}}}{{{\mathrm{TP}} + {\mathrm{FN}}}}{\text{，}} $$ (1) $$ p = \dfrac{{{\mathrm{TP}}}}{{{\mathrm{TP}} + {\mathrm{FP}}}}{\text{，}} $$ (2) $$ F = \dfrac{{2t p}}{{t + p}}{\text{，}} $$ (3) 式中，
${\mathrm{TP}}$ 表示正确分割，指一对一的关系，将树木分割为单独的一株树，即有效单木分割；${\text{FP}}$ 表示过分割，指一对多的关系，把原本属于1个整体的单木分成了多株单木；${\text{FN}}$ 表示欠分割，指1株或几株树被划分到其临近的单木中，且未被检测出来。1.4 树木定量结构模型
树木定量结构模型是重建树木三维模型的方法。目前较为常用的模型有PypeTree[23]、TreeQSM[24]、SimpleTree[25]等，本文将TreeQSM (https://github.com/InverseTampere/TreeQSM)作为试验模型，该模型由芬兰坦佩雷理工大学的Pasi Raumonen等开发，在不同森林类型中均具有较好的普适性，其构建单木模型的核心方法是对点云进行体元化分割，实现对树木几何特征的定量表达，如体积、冠幅、分支角度等。其中，PatchDiam1、PatchDiam2Min和PatchDiam2Max是TreeQSM建模所必需的3个重要参数。建模过程主要分为2个阶段，第一阶段是按照PatchDiam1定义的体元大小进行简单分割；第二阶段是按照设置的最大、最小体元实现精细分割，也可同时添加其他参数实现更具体的分割。其中，在所有样本中随机选取20株树对参数进行敏感性分析，并将最优参数作为其他样本的输入参数。具体参数设置可以参照开发者文档(https://github.com/InverseTampere/TreeQSM/blob/master/Manual)。在完成三维模型后，提取胸径、主干体积、枝干体积、分支数、冠层基部高度、冠层面积、冠层体积以及冠幅等参数，并且为保证建模效果具有稳定性，每个样本获取20组树木定量结构模型，各个参数取其平均值作为最终结果。
1.5 单木生物量模型构建及精度评价
1.5.1 地上生物量计算
按照《主要树种立木生物量模型与碳计量参数》标准(GB/T 43648—2024)[26]计算单木生物量，研究区内栎树采用南方其他省份栎树分布区对应的计算公式，而其他阔叶树则采用全国通用的其他硬阔树地上生物量公式，具体如下：
$$ M = 0.213\;6{D^{2.304\;16}}{\text{，}} $$ (4) $$ M = 0.187\;5{D^{2.333\;00}}{\text{，}} $$ (5) 式中，M为地上生物量估计值，D为林木胸径。式(4)为栎树公式，式(5)为其他阔叶树公式。
1.5.2 建模变量筛选
利用Pearson相关系数(r)和方差膨胀因子(Variance inflation factor，VIF)进行特征筛选，移除共线性变量，具体计算公式如下：
$$ r=\dfrac{{\displaystyle\sum }_{i=1}^{n}({X}_{i}-\overline {X})({Y}_{i}-\overline {Y})}{\sqrt{{\displaystyle\sum }_{i=1}^{n}{({X}_{i}-\overline {X})}^{2}}\sqrt{{\displaystyle\sum }_{i=1}^{n}{({Y}_{i}-\overline {Y})}^{2}}}， $$ (6) $$ {{\mathrm{VIF}}}_{i}=\dfrac{1}{1-{r}_{i}^{2}}， $$ (7) 式中，
$ {X}_{i} $ 和$ {Y}_{i} $ 分别是2个变量的第i个观察值，$ \overline {X} $ 和$ \overline {Y} $ 是2个变量的均值，ri表示第i个变量与其他变量的相关系数，VIFi是第i个变量的方差膨胀因子，当VIF小于5时可以认为该变量与其他变量之间不存在明显的共线性。1.5.3 机器学习算法
选取RF、SVM和MLP 3种方法进行生物量建模。RF是集成算法的延伸，其核心是对样本进行随机抽样，通过构造多个决策树来提升算法的性能，对模型具有很好的解释性。SVM通过划分超几何空间并结合核函数，将数据隐式地映射到一个更高维的空间，使其在新的空间中变得线性可分，进而完成模型构建。MLP是前馈神经网络，由输入层、1个或多个隐藏层和输出层组成，数据传入输入层，通过激活函数完成在隐藏层的传输，最终由输出层对数据进行预测。机器学习可以使用各种具有更高精度的大输入数据对复杂的空间模式进行建模，这种方法越来越多地用于森林生物量的估测研究中[27-28]。
1.5.4 模型评价
本文将试验样本数据集按7∶3随机划分成训练集和测试集。在训练集上，通过网格搜索与十折交叉验证的方式寻找模型的最优参数。具体而言，将训练集进一步分为10等份，每次选取其中1份作为验证集，剩余的9份作为训练集，循环10次，确保每一份数据都被用作验证集，并利用交叉验证确定的最优参数，在整个训练集上重新训练模型。模型训练完成后，使用测试集对模型的性能进行评估，以检验模型的泛化能力。其中，采用决定系数(R2)、均方根误差(RMSE)和相对均方根误差(rRMSE)进行模型训练集及测试集的精度评价。R2用于评估模型的整体拟合效果，RMSE和rRMSE用于衡量模型的预测准确性，具体公式如下：
$$ {R}^{2}=1-{ \displaystyle\sum\nolimits _{i=1}^{n}{({y}_{i}-{\overset{\frown}{y}_{i}})}^{2}/}{ \displaystyle\sum\nolimits _{i=1}^{n}{({y}_{i}-\overline{y})}^{2}}， $$ (8) $$ {\mathrm{RMSE}}=\sqrt{{ \displaystyle\sum\nolimits _{i=1}^{n}{({y}_{i}-{\overset{\frown}{y}_{i}})}^{2}/n}} ，$$ (9) $$ {\mathrm{rRMSE}} = \dfrac{{{\mathrm{RMSE}}}}{{\overline y}} \times 100{\text{%}}， $$ (10) 式中，
$n$ 为样本中单木数量，${{{y}}_{{i}}}$ 为单木生物量的估计值，$ {\overset{\frown}{y}_{i}} $ 为模型预测值，$\overline y$ 为样本平均值。2. 结果与分析
2.1 基于机载激光雷达的多变量精度评价
2.1.1 空间位置精度
将非地面点按照TIN内插法生成分辨率为0.1 m的栅格，即数字高程模型(Digital elevation model，DEM)。在相同位置处对比DEM与单木RTK的海拔，以此验证激光雷达的位置精度。结果(图3)显示两者海拔几乎一致(R2=
0.9999 )，说明UAV-LS所采集的数据位置偏差极小，点云数据中显示的单木位置就是其本身在真实空间的位置，不需要再进行空间位置校正。2.1.2 点云分割精度
本文采用边界角点所确定的范围对整个大样地进行划分，共分为9块小样地，当靠近边界的树木提取不完整时，按照树基部位置判断所属的样地，并对缺失部分采取提取、合并的方式进行补全，得到完整的样地点云。由于阔叶林林下植被复杂且UAV-LS采用自上而下的方式获取数据，因此所获取的点云均为中上层枝干较为清晰的点云数据，较难获取样地内中下层小树点云。根据样地内树种分布的复杂性，在进行点云分割时，需针对不同样地采用不同参数进行设置，以期达到最优分割。最后对各个样地的分割结果进行精度评价，具体结果如表2所示。
表 2 样地分割精度统计结果Table 2. Plot segmentation accuracy statistical results样地
Plot实测株数
Actual tree count分割株数
Segmented tree countTP FN FP t P F值
F-score1 20 19 18 2 1 0.90 0.95 0.92 2 23 22 19 4 3 0.83 0.86 0.84 3 27 24 21 6 3 0.78 0.88 0.82 4 23 23 18 5 5 0.78 0.78 0.78 5 36 37 30 6 7 0.83 0.81 0.82 6 27 26 20 6 5 0.77 0.80 0.78 7 19 19 18 1 1 0.95 0.95 0.95 8 17 19 16 3 1 0.84 0.94 0.89 9 28 28 25 3 3 0.89 0.89 0.89 由表2可知，在所有样地的分割中，4号和6号样地调和平均值相对偏低，即F值为0.78。影响其分割精度的主要原因有2点：一是样地中存在丛生的树木点云，二是林分中间层的乔木点云不够清晰。整体来看，共有185株树木被准确地识别并分割出来，占样地内点云树木总量的84%。在阔叶林树种中，分割的综合效果较好，对于一些分割效果较差的树木，通过添加种子点的位置进行再次分割，完成单株树木的精准提取。最终得到220株中大型的乔木点云数据，其胸径范围为15.0~67.2 cm，生活史阶段均为中树及以上树木。
2.1.3 树高提取精度
根据单木分割结果，提取归一化后高度的最大值与最小值，其差值即是树高。与野外实际测量值对比结果(图4)表明，提取的单木树高与测量结果非常接近(R2=0.870 0)，说明激光雷达数据能很好地反演树木的实际高度。
2.2 TreeQSM构建及参数提取
在进行参数敏感性分析后，设置PatchDiam1、PatchDiam2Min、PatchDiam2Max参数值分别为[0.05，0.08]、[0.01，0.015，0.02]与[0.06，0.10]。敏感性分析后的参数能较好地构建出大胸径树木，但针对胸径较小的树木需要在此基础上稍作修改，其中PatchDiam1参数较为影响小胸径树木的模型效果。在剔除树冠不完整以及主干缺失的点云数据后，共构建176个TreeQSM模型。所有模型中，茅栗(65株)数量最多，其次是短柄枹栎(43株)和化香(42株)。槲栎Quercus aliena Blume和锐齿槲栎也略高，分别为8和9株。其他树种如灯台树Cornus controversa、金缕梅Hamamelis mollis、鹅耳枥Carpinus turczaninowii、黄檀Dalbergia hupeana、山合欢Albizia kalkora以及尾叶樱Cerasus dielsiana等数量较少，各自仅有1或2株。总体来看，树种分布不均，其中茅栗、短柄枹栎和化香占据了较大比例，为优势树种。树木具体建模效果以1株分支较多的点云数据(图5)进行展示，图中所构建的TreeQSM模型能较为清晰地模拟出树木各级枝干的特征。
![]()
图 5 TreeQSM建模前(a)、后(b)对比图Figure 5. Comparison images before (a) and after (b) TreeQSM modeling通过构建三维模型从而量化树木参数，并在其中选取7个树木结构参数进行分析，分别统计各个参数的平均值、标准差、最小值、上四分位数、中位数、下四分位数以及最大值，结果见表3。并通过实测胸径与冠幅检验提取精度，其中提取的胸径与实测胸径R2达到0.930 0，冠幅拟合R2达到0.440 0(图6)，说明构建的三维模型能较为准确地反映出树木的三维结构。这与相关研究结果[29-30]相符，当胸径大于30 cm时，TreeQSM能更准确地实现重建[31]。本文参与构建模型的点云树木多为胸径较大的乔木树种，因提取的胸径数据为计算生物量的强关联因子，故不参与生物量的模型构建。
表 3 参数统计信息Table 3. Parameter statistical information统计量
Statistic主干体积/m3
Trunk volume枝干体积/m3
Branch volume分支数
Branch number冠层基部高度/m
Crown base height冠层面积/m2
Crown area冠层体积/m3
Crown volume冠幅/m
Crown width平均值 Average 0.82 0.76 146.97 7.89 16.79 39.72 3.78 标准差 SD 0.58 1.40 166.26 3.17 19.34 103.76 2.31 最小值 Min. 0.12 0.01 8.00 0.44 1.38 0.01 0.67 下四分位 P25 0.41 0.11 41.00 5.47 5.09 0.71 2.10 中位数 P50 0.67 0.25 85.00 8.34 10.49 3.81 3.25 上四分位 P75 1.04 0.68 173.00 10.43 19.27 26.33 4.85 最大值 Max. 3.56 10.35 903.00 14.57 138.73 929.58 13.17 ![]()
图 6 胸径(DBH)及冠幅拟合结果Figure 6. Diameter at breast height (DBH) and crown width fitting results2.3 生物量估算模型构建
2.3.1 变量筛选
经过共线性分析后，共保留5个特征因子用于模型构建，分别是主干体积、分支数、冠层体积、冠层基部高度以及冠幅。其VIF均小于5，主干体积VIF为2.95，分支数3.82，冠层体积3.11，冠层基部高度1.30，冠幅4.62，说明筛选后的变量之间没有较为严重的共线性。进一步对各变量之间进行显著性检验，结果如表4所示。各变量之间均在P<0.001水平下显著相关，除变量冠层基部高度与其他变量负相关外，其他变量之间均为正相关，其中主干体积与生物量相关系数最高，达到0.91。通过相关性分析结果可以看出，冠层结构参数对于生物量研究具有重要意义，如冠层体积(0.72)、冠幅(0.77)等。在陆地生态系统中，植物冠层是介于大气和土壤基质间进行有机物合成代谢的重要部位，对于研究树木生物量起到直接或间接的作用。如单木生物量随冠层基部高度的上升会出现一定程度的下降，这与冠层大小有着密切的关系，当冠层基部位置低时，属于树木冠层的部分就会增多，从而使生物量增加。
表 4 各建模因子相关性统计1)Table 4. Correlation analysis of each modeling factor因素
Factor生物量
Biomass主干体积
Trunk volume冠层体积
Crown volume冠幅
Crown width冠层基部高度
Crown base height分支数
Branch number生物量 Biomass 1.00 主干体积 Trunk volume 0.91*** 1.00 冠层体积 Crown volume 0.72*** 0.71*** 1.00 冠幅 Crown width 0.77*** 0.77*** 0.73*** 1.00 冠层基部高度 Crown base height −0.46*** −0.40*** −0.38*** −0.46*** 1.00 分支数 Branch number 0.64*** 0.61*** 0.77*** 0.81*** −0.44*** 1.00 1) ***表示在P<0.001水平显著相关(Pearson法)。
1) *** indicates significant correlation at P<0.001 level (Pearson method).2.3.2 模型参数及建模结果
利用网格搜索遍历所有可能的参数组合，并在优化过程中通过交叉验证来评估每个参数组合的性能，从而选择最优的参数集。RF中，max_depth代表决策树的最大深度，限制树的深度可以防止过拟合；min_samples_split是每个节点可以分裂的最小样本数量，该参数影响树的分枝过程，较大的值可以防止过拟合；n_estimators控制树的数量，更多的树通常能提高模型的稳定性。MLP中，alpha代表正则化系数，用于防止模型过拟合，较大的alpha会使模型更平滑，但可能降低模型拟合能力；hidden_layer_sizes是隐藏层的大小和数量，影响模型的学习能力和复杂性；learning_rate控制学习率模型参数更新的速度，其中自适应学习率可根据模型收敛情况自动调整；max_iter为模型训练过程中允许的最大迭代步数。SVM中，C代表惩罚系数，较大的C值会强制模型在训练集上更好地拟合，但可能导致过拟合；gamma是核函数中的参数，决定样本的影响范围，较大的gamma会使模型更倾向于对训练数据的过拟合；Kernel是核函数类型(如线性核等)。
在经过参数筛选后，RF的最优参数为max_depth=3，min_samples_split=5，n_estimators=100；MLP的最优参数是alpha=0.1，hidden_layer_sizes=(50, 50)，learning_rate=‘adaptive’，max_iter=
1000 ；SVM的最优参数是C=1，gamma=0.1，kernel=‘linear’。根据最优参数进行生物量估测的结果(图7)表明，RF的测试集R2达到0.880 0，单木生物量的RMSE为192.81 kg，rRMSE为29.88%；MLP的测试集R2达到0.820 0，单木生物量的RMSE为236.48 kg，rRMSE为36.65%；SVM的测试集R2达到0.810 0，单木生物量的RMSE为243.67 kg，rRMSE为37.77%。3种机器学习模型中，RF的效果最佳，MLP与SVM结果较为接近，精度稍低于RF，但也都表现出较高的准确性，证实了树木定量结构模型所提取的参数能够较好地构建出单木生物量模型，为森林资源清查提供了新方法，对后续研究具有重要意义。![]()
图 7 3种机器学习算法的生物量估测模型散点拟合对比Figure 7. Scatter plot fitting comparison of biomass estimation models based on the three machine learning algorithms3. 讨论
本研究使用的RF算法可以训练出最优的单木生物量估测模型，在生物量估测模型的研究中，可以分为参数和非参数模型两大类[10,32]，参数模型分为多元线性模型和非线性模型，非参数模型则包括人工神经网络、支持向量机、深度学习等机器学习方法。针对不同研究区域和研究对象，两类模型估测精度各有所长[33]。本文采用机器学习方法进行模型构建，旨在提高目标变量的预测精度，而参数模型可以通过参数组合评估和解释变量对目标变量的影响程度。后续研究如需揭示目标变量与解释变量之间的特定关系，可以选择多元线性和非线性模型探究各参数对单木生物量的响应关系。如Abd Rahman等[34]根据树高、胸径、树冠大小和树冠基部高度生成新的异速生长方程进行生物量计算。也可以将生物量模型细化，分为主干生物量与枝干生物量，这可能需要加入更多的树木参数进行模型完善，如唐依人等[35]利用枝条因子进行分支生物量模型构建。同样，可以根据主干体积与单木生物量之间的强相关性，探索基于体积与木材密度的生物量计算方式是否优于异速生长模型[36]，值得注意的是，不同树种的木材密度存在差异，因此需要先区分树种的木材密度，才能得到更准确的结果。
点云密度也会对TreeQSM模型效果产生影响，并在一定程度上影响树木分割精度。虽然本研究获取的数据为高密度的点云数据，但是也存在点云冠层以及枝干缺失的情况，UAV-LS采用自上而下的方式获取点云数据，当上层树木出现遮挡时，就无法完整获取中下层的树木点云[37]。本研究的阔叶林为天马国家级自然保护区内的天然阔叶混交林，群落发育比较成熟，具有地带性稳定群落特征。因其冠层、林下环境较为复杂，所获取的树木点云主要为主林层大胸径乔木树种数据，胸径均在15.0 cm以上，而对于林下的小乔木则较难获取。但随着地基激光雷达(TLS)和移动激光扫描仪(Mobile laser scanning，MLS)等便携设备的深入应用，也为多源点云数据融合获取高精度三维数据提供了新的技术手段。通过与地面点云数据相结合，可以很好地弥补无人机点云数据在林下缺失的情况，无人机点云数据也能补充地面点云数据缺失的上层树冠部分，多源点云的融合可以促进对树木参数更可靠的估计，从而有助于对复杂森林场景进行互补解释[38]，起到扩大应用范围与提高应用精度的效果。Qi等[39]研究发现，在冠层覆盖率增加的情况下，无人机激光点云与MLS相融合的数据在树木定量结构模型中表现最好。除多源点云数据的融合外，探究无人机点云最佳的建模密度也是一种新思路。本文研究情境下表现为仿地飞行高度60 m，后续的研究可以尝试更低的飞行高度以获取更丰富的树木结构，进而深入探究哪种高度下无人机点云数据的三维建模效果最好。
树木定量结构模型不仅用来估测森林蓄积量、生物量等参数，在树种识别方面也取得了不错的结果[40-41]。Hui等[42]利用几何分形特征和QSM特征进行树种分类，探索特征向量的分类效果，获得了较好的分类结果。QSM可以定量化树木结构，统计树木几何和拓扑属性，如本研究中的主干体积等参数。树木结构可以作为树种分类的重要依据，本文主要基于树木参数进行回归训练，未涉及分类任务。针对植被结构复杂的落叶阔叶林，如何精准识别树种显得尤为重要。特别在树木落叶期，无法根据树叶等信息判断种属，如果能通过树冠、分支等树木结构信息直接进行树种识别，将大大降低野外工作中树种识别的难度。
4. 结论
本文基于UAV-LS和QSM开展天马自然保护区阔叶树种单木生物量研究。其中，空间位置精度R2达到
0.9999 ，树高提取精度R2为0.870 0，表明激光雷达数据精度很高，并在点云总体分割精度较高的情况下，进行单木点云的精确提取。在此基础上，构建单木的三维结构模型，结果显示，胸径、冠幅拟合精度R2分别为0.930 0和0.440 0，并进一步选择3种机器学习算法对树木参数进行生物量估测模型构建。3种模型中，RF算法结果最优(R2=0.880 0, RMSE=192.81 kg, rRMSE=29.88%)，证明QSM所提取的树木参数可以构建出精度较高的落叶期阔叶树种生物量估测模型，为以后相关研究提供启发。本文使用的方法在一定程度上大大减轻了野外工作量和难度，点云数据的获取比地基激光雷达更方便，在森林资源调查中具有巨大的应用价值。 -
![]()
图 1 试验装置模拟图(A)和实物图(B)
1：阳极室(厌氧氨氧化污泥)，2：阴极室(小球藻)，3：阳离子交换膜，4：石墨棒电极，5：数据采集卡；废水从进水口a处进入装置，停留24 h后从出水口a处排出并收集，循环利用至进水口b，停留24 h后从出水口b排出。
Figure 1. Simulation image (A) and actual object image (B) of experiment apparatus
1: Anode chamber (anaerobic ammonia oxidation sludge), 2: Cathode chamber (chlorella), 3: Cation exchange membrane, 4: Graphite rod electrode, 5: Data acquisition card; The wastewater enters the device from inlet a, and is then discharged and collected from outlet a after 24 h, it is recycled and used again at inlet b, where it stays for another 24 h before being discharged from outlet b.
![]()
图 3 出水中氨氮(A)和NO2−−N(B)质量浓度的变化
Figure 3. Mass concentration changes of ammonium nitrogen (A) and NO2−-N (B) in discharged wastewater
![]()
图 4 出水中可溶性化学需氧量的变化
Figure 4. Changes of soluble chemical oxygen demand in discharged wastewater
![]()
图 5 出水中总磷质量浓度的变化
Figure 5. Mass concentration changes of total phosphorus in discharged wastewater
-
[1] ZHENG H, WU X, ZOU G, et al. Cultivation of Chlorella vulgaris in manure-free piggery wastewater with high-strength ammonium for nutrients removal and biomass production: Effect of ammonium concentration, carbon/nitrogen ratio and pH[J]. Bioresource Technology, 2019, 273: 203-211. doi: 10.1016/j.biortech.2018.11.019
[2] AHMAD A L, CHIN J Y, HARUN M H Z M, et al. Environmental impacts and imperative technologies towards sustainable treatment of aquaculture wastewater: A review[J]. Journal of Water Process Engineering, 2022, 46: 102553. doi: 10.1016/j.jwpe.2021.102553
[3] LI X, WU S, YANG C, et al. Microalgal and duckweed based constructed wetlands for swine wastewater treatment: A review[J]. Bioresource Technology, 2020, 318: 123858. doi: 10.1016/j.biortech.2020.123858
[4] ZENG W S, WANG D H, WU Z Y, et al. Recovery of nitrogen and phosphorus fertilizer from pig farm biogas slurry and incinerated chicken manure fly ash[J]. Science of the Total Environment, 2021, 782: 146856. doi: 10.1016/j.scitotenv.2021.146856
[5] KABUBA J, LEPHALLO J, RUTTO H. Comparison of various technologies used to eliminate nitrogen from wastewater: A review[J]. Journal of Water Process Engineering, 2022, 48: 102885. doi: 10.1016/j.jwpe.2022.102885
[6] ZHANG W, CHU H, YANG L, et al. Technologies for pollutant removal and resource recovery from blackwater: A review[J]. Frontiers of Environmental Science & Engineering, 2023, 17(7): 83. doi: 10.1007/s11783-023-1683-3
[7] CHEN X, LIU L, BI Y, et al. A review of anammox metabolic response to environmental factors: Characteristics and mechanisms[J]. Environmental Research, 2023, 223: 115464. doi: 10.1016/j.envres.2023.115464
[8] ABMA W, SCHULTZ C, MULDER J W, et al. The advance of anammox[J]. Water, 2007: 36-37.
[9] 李彩林, 刘扬, 李月. 菌藻微生物燃料电池处理模拟养殖废水的研究[J]. 青海大学学报(自然科学版), 2023, 41(6): 1-8. [10] 苟珍琼, 郑道会, 罗发文. 光催化耦合微生物燃料电池协同处理废水的研究进展[J]. 现代化工, 2023, 43(S1): 101-104. [11] 王子义, 张嫄, 刘根深, 等. 微生物燃料电池污废水处理及能源化研究进展[J]. 净水技术, 2024, 43(11): 29-38. [12] WANG Y M, LIN Z Y, SU X S, et al. Cost-effective domestic wastewater treatment and bioenergy recovery in an immobilized microalgal-based photoautotrophic microbial fuel cell (PMFC)[J]. Chemical Engineering Journal, 2019, 372: 956-965. doi: 10.1016/j.cej.2019.05.004
[13] SONU K, SOGANI M, SYED Z, et al. Improved decolorization of dye wastewater and enhanced power output in the electrically stacked microbial fuel cells with H2O2 modified corncob anodes[J]. Environmental Progress & Sustainable Energy, 2021, 40(5): ep13638. doi: 10.1002/ep.13638
[14] MIN B, KIM J R, OH S E, et al. Electricity generation from swine wastewater using microbial fuel cells[J]. Water Research, 2005, 39(20): 4961-4968. doi: 10.1016/j.watres.2005.09.039
[15] 王佳璇, 段嘉琪, 刘喆, 等. 藻类微生物燃料电池的构型发展及应用现状[J]. 环境科学与技术, 2023, 46(4): 61-71. [16] HASSAN M, WEI H, QIU H, et al. Power generation and pollutants removal from landfill leachate in microbial fuel cell: Variation and influence of anodic microbiomes[J]. Bioresource Technology, 2018, 247: 434-442. doi: 10.1016/j.biortech.2017.09.124
[17] YADAV G, SHARMA I, GHANGREKAR M, et al. A live bio-cathode to enhance power output steered by bacteria-microalgae synergistic metabolism in microbial fuel cell[J]. Journal of Power Sources, 2020, 449: 227560. doi: 10.1016/j.jpowsour.2019.227560
[18] 严茜. 藻膜阴极微生物电化学系统处理垃圾渗滤液效能及机理研究[D]. 烟台: 烟台大学, 2023. [19] DING A, ZHAO D, DING F, et al. Effect of inocula on performance of bio-cathode denitrification and its microbial mechanism[J]. Chemical Engineering Journal, 2018, 343: 399-407. doi: 10.1016/j.cej.2018.02.119
[20] ZHANG Y, ZHAO Y, ZHOU M. A photosynthetic algal microbial fuel cell for treating swine wastewater[J]. Environmental Science and Pollution Research, 2019, 26(6): 6182-6190. doi: 10.1007/s11356-018-3960-4
[21] ZHANG J, ZHANG Z, RONG K, et al. Simultaneous anaerobic ammonium oxidation and electricity generation in microbial fuel cell: Performance and electrochemical characteristics[J]. Processes, 2022, 10(11): 2379. doi: 10.3390/pr10112379
[22] SAM T, LE ROES-HILL M, HOOSAIN N, et al. Strategies for controlling filamentous bulking in activated sludge wastewater treatment plants: The old and the new[J]. Water, 2022, 14(20): 3223. doi: 10.3390/w14203223
[23] GUO J, WANG S, WANG Z, et al. Effects of feeding pattern and dissolved oxygen concentration on microbial morphology and community structure: The competition between floc-forming bacteria and filamentous bacteria[J]. Journal of Water Process Engineering, 2014, 1: 108-114. doi: 10.1016/j.jwpe.2014.03.011
[24] DALIRY S, HALLAJISANI A, MOHAMMADI ROSHANDEH J, et al. Investigation of optimal condition for Chlorella vulgaris microalgae growth[J]. Global Journal of Environmental Science and Management, 2017, 3(2): 217-230.
[25] ZHANG C, LI S, HO S. Converting nitrogen and phosphorus wastewater into bioenergy using microalgae-bacteria consortia: A critical review[J]. Bioresource Technology, 2021, 342: 126056. doi: 10.1016/j.biortech.2021.126056
[26] 黄山, 杨莹莹, 仇志峰, 等. 蛋白核小球藻-硝化细菌共培养条件优化与氮磷去除效果[J]. 青岛理工大学学报, 2023, 44(3): 15-25. doi: 10.3969/j.issn.1673-4602.2023.03.003 [27] 中国环境监测总站. 水质 化学需氧量的测定 重铬酸盐法: HJ 828—2017[S]. 1989. [28] 北京市环保监测中心, 上海市环境监测中心. 水质 总磷的测定 钼酸铵分光光度法: GB/T 11893—1989[S]. 1989. [29] 江苏省环境监测站. 水质 铵的测定 纳氏试剂比色法: GB/T 7479—1987[S]. 1987. [30] 中国科学院南京土壤研究所, 中国科学院亚热带农业生态研究所, 中国科学院西双版纳热带植物园, 等. 土壤硝态氮的测定 紫外分光光度法: GB/T 32737—2016[S]. 北京: 中国标准出版社, 2016. [31] 中国环境监测总站. 大气降水中亚硝酸盐的测定 N-(1-萘基)-乙二胺光度法: GB/T 13580.7—1992[S]. 1992. [32] ZHANG H, GE C, YU M, et al. Performance of cathodic nitrate reduction driven by electricity generated from ANAMMOX sludge in anode[J]. Process Biochemistry, 2020, 90: 148-155. doi: 10.1016/j.procbio.2019.11.013
[33] ZHANG G, ZHAO Q, JIAO Y, et al. Efficient electricity generation from sewage sludge using biocathode microbial fuel cell[J]. Water Research, 2012, 46(1): 43-52. doi: 10.1016/j.watres.2011.10.036
[34] WANG J T, SONG A I, HUANG Y, et al. Domesticating Chlorella vulgaris with gradually increased the concentration of digested piggery wastewater to bio-remove ammonia nitrogen[J]. Algal Research-Biomass Biofuels and Bioproducts, 2021, 60: 102526. doi: 10.1016/j.algal.2021.102526
[35] COLLOS Y, HARRISON P J. Acclimation and toxicity of high ammonium concentrations to unicellular algae[J]. Marine Pollution Bulletin, 2014, 80(1/2): 8-23.
[36] 李宝珍, 范晓荣, 徐国华. 植物吸收利用铵态氮和硝态氮的分子调控[J]. 植物生理学通讯, 2009, 45(1): 80-88. [37] MAESTRINI S Y, ROBERT, J M, LEFTLEY J W, et al. Ammonium thresholds for simutaneous uptake of ammonium and nitrate by oyster-pond algae[J]. Journal of Experimental Marine Biology and Ecology, 1986, 102(1): 75-98. doi: 10.1016/0022-0981(86)90127-9
[38] ZHANG L, TIAN Z, QIAN Y, et al. Long-term effects of phosphorus deficiency on one-stage partial nitrification-anammox system and recovery strategies[J]. Journal of Cleaner Production, 2023, 402: 136820. doi: 10.1016/j.jclepro.2023.136820
[39] HUANG W, ZHOU J, HE X, et al. Simultaneous nitrogen and phosphorus removal from simulated digested piggery wastewater in a single-stage biofilm process coupling anammox and intracellular carbon metabolism[J]. Bioresource Technology, 2021, 333: 125152. doi: 10.1016/j.biortech.2021.125152
[40] ZHA X, MA J, LU X. Use of a low-cost and energy-efficient device for treating low-strength wastewater at low temperatures focusing on nitrogen removal and microbial community[J]. Science of the Total Environment, 2020, 722: 137916. doi: 10.1016/j.scitotenv.2020.137916
[41] 谭顺. 驯化对小球藻在沼液废水处理中的强化效果研究[D]. 南宁: 广西大学, 2022. [42] PRAVEEN P, GUO Y C, KANG H, et al. Enhancing microalgae cultivation in anaerobic digestate through nitrification[J]. Chemical Engineering Journal, 2018, 354: 905-912. doi: 10.1016/j.cej.2018.08.099
[43] LU S M, LIU X G, LIU C, et al. A review of ammonia-oxidizing archaea and anaerobic ammonia-oxidizing bacteria in the aquaculture pond environment in China[J]. Frontiers in Microbiology, 2021, 12: 775794. doi: 10.3389/fmicb.2021.775794
[44] HU P, STROM P F. Effect of pH on fungal growth and bulking in laboratory-activated sludges[J]. Research Journal of the Water Pollution Control Federation, 1991, 63(3): 276-277.
下载:
计量
- 文章访问数: 67
- HTML全文浏览量: 14
- PDF下载量: 29
