Extraction, antioxidant and antibacterial activities of essential oils from Magnolia wufengensis
-
摘要:目的研究红花玉兰花蕾挥发油的萃取工艺，并对其化学组成、抗氧化及抗菌活性进行鉴定和评价，旨在为红花玉兰花蕾挥发油的萃取建立一套可行的工艺体系，并为其后续的开发应用提供技术支持和理论参考。方法以‘娇红1号’红花玉兰阴干花蕾为原料，在单因素试验基础上，采用响应面分析法优化超临界CO2萃取工艺，结合自动化质谱解卷积和保留指数对最佳工艺条件下得到的挥发油进行成分定性和定量；采用DPPH自由基清除法和还原能力法评价挥发油抗氧化活性，纸片扩散法和微量肉汤稀释法评价挥发油抗菌活性。结果阴干预处理条件下，红花玉兰花蕾挥发油萃取的最佳条件为萃取温度47 ℃，萃取压力40.7 MPa，动态萃取时间102 min，该条件下挥发油的萃取率为（1.217 ± 0.014）%。GC-MS共鉴定出30种化学成分，占检测到的挥发油总量的70.65%。DPPH自由基清除法测定得到挥发油IC50值为12.82 mg/mL AEAC值为198.64 mg，表明挥发油具有一定的抗氧化活性，还原力测定也得到了这一结果。抗菌试验表明，挥发油对4种供试菌抑制作用强弱依次为金黄色葡萄球菌 > 枯草芽孢杆菌 > 大肠杆菌 > 鼠伤寒沙门氏菌。结论首次建立了红花玉兰阴干花蕾挥发油的超临界CO2萃取工艺，确定了该工艺条件下挥发油的主要成分，并验证其具有一定的抗氧化和抗菌活性。Abstract:ObjectiveThe aim of this study is to establish a set of reliable process system for the extraction of essential oils from Magnolia wufengensis buds and provide technical support as well as theoretical reference for its subsequent development and application. The extraction process of essential oils from M. wufengensis buds was studied, and chemical compositions, antioxidant and antibacterial activities were identified and evaluated.MethodBased on single-factor experiment, the supercritical CO2 fluid extraction conditions of essential oils from air-dried buds of M. wufengensis ‘Jiaohong No.1’ were optimized by response surface methodology. Components of essential oils were qualitatively and quantitatively analyzed by automated mass spectral deconvolution and identification system combined with retention index. DPPH free radical scavenging method and reducing power method were used to analyze the antioxidant activity of essential oils. Disk diffusion method and broth microdilution method were used to measure the antibacterial activities.ResultUnder the air-dried pretreatment condition, the optimum processing conditions for extracting essential oils from M. wufengensis were as follows: extraction temperature 47 ℃, extraction pressure 40.7 MPa, dynamic extraction time 102 min. The extraction yield of essential oils under this condition was (1.217 ± 0.014)%. A total of 30 components were identified, accounting for 70.65% of the total components. DPPH free radical scavenging method showed that the IC50 value of essential oils was 12.82 mg/mL, and the AEAC value was 198.64 mg. The essential oils had certain antioxidant activity, and the reducing power method showed the same results. The inhibitory effect of essential oils on the four tested bacteria was in order of Staphylococcus aureus > Bacillus subtilis > Escherichia coli > Salmonella typhimurium.ConclusionA feasible supercritical CO2 extraction process of the essential oils from the air-dried buds of M. wufengensis was established for the first time. Furthermore, the main components of the essential oils were identified and its antioxidant and antibacterial activities were verified.
-
在快速城镇化发展进程中，人类较大程度地改变了自然地形地貌，以不透水地面铺砌代替原有透水土壤和植被，造成下渗量与蒸发量的显著减少，从而引起同强度暴雨形成的地表径流量增大，峰值流量增加等现象[1]。为解决城市内涝等问题，“海绵城市”和低影响开发系统(Low impact development, LID)概念应运而生。目前对低影响开发系统的研究多数集中于市政道路[2]、住宅小区[3-4]、公园绿地[5-6]以及城市绿道[7]中绿地的渗透滞留能力，对其后续储蓄利用的量化研究较少。为解决上述问题，本文以位于北京市的首钢西十为研究对象，结合景观需求建设渗透滞留与储蓄利用相结合的低影响开发系统。
1. 研究区域概况
北京市位于华北平原北部，为典型的北温带半湿润大陆性季风气候，年平均降雨量(1950—2012年期间)为585.0 mm，80%降雨量集中于汛期(6—9月)[8]，对城市内涝造成极大压力。本文选取首钢北区西北角的西十片区作为研究区域，西十内部及周边场地硬质化程度高达65.3%，现状下垫面对自然雨水下渗疏解能力非常有限，具有典型特征。
西十地块原为首钢料场，运输线两侧土壤含铁较高，研究表明土壤污染深度为0.5~2.0 m[15]。改造过程中将表层以下2 m的土壤替换为壤土，渗透系数为3.5×10-6 m/s[9]。
如图 1、图 2所示，西十地块四面环街，占地4.24 hm2，地势整体呈西北高东南低，内部有筒仓、转运站、料仓、泵站等多组建筑体，总占地面积约1.5 hm2。场地周边区域(包括秀池东路、秀池北街两条道路和筒仓办公楼、北七筒建筑遗址、停车场、N3-17办公楼、武警用房、员工宿舍、干法除尘等建筑及其外环境，共5.8 hm2)地势略低于西十地块，东侧和北侧有多组由厂房改建的配套服务建筑，强降雨时存在较为严重的内涝问题亟待解决，因此将其与西十地块一同作为本文的研究对象，称为外源雨水地块。
2. 低影响开发系统空间布局设计
2.1 设计目标与原则
根据《新首钢高端产业综合服务区控制性详细规划》要求，西十地块及外源雨水地块低影响开发系统雨水消纳能力应达到北京市3年一遇24 h降雨量(108 mm[9])的标准。设计在满足功能需求和景观效果的前提下，制定合理的低影响开发系统设计策略，在有限的绿地空间中因地制宜地叠加低影响开发设施层以应对降雨后地表径流产生速度加快、峰值流量加大、峰现时间提前等问题，以达到上位规划对径流消纳的要求，同时考虑将收集的地表径流用于后期灌溉。
2.2 设计策略
根据《海绵城市建设技术指南》低影响开发主要通过“源头削减、中途传输、末端调蓄”等多种手段[10]。西十地块及外源雨水场地提出的设计策略为：根据场地竖向、土地利用等条件优先构建渗透滞留系统，在不具备条件的情况下考虑将径流收集再利用，从源头、中途、末端3个方面进行低影响开发系统设计(见图 3)。
2.2.1 西十地块低影响开发系统
西十地块以建立渗透滞留系统为目标，针对“源头削减”本研究提出在满足功能使用的前提下，尽可能多地利用透水铺装、建设绿色屋顶、增加绿地面积，从源头有效削减地表径流。此次设计的“中途传输”途径为植草沟，沿硬质边缘形成线性景观带，达到收集、传输、滞留及下渗径流的目的。“末端调蓄”设施有下沉式绿地、生物滞留池两种。植草沟内的雨水排入末端调蓄设施内，雨水在此进行滞留下渗，当雨量过大时末端调蓄设施内的雨水通过溢流口排入市政雨水管网。
2.2.2 外源雨水地块低影响开发系统
外源雨水地块绿地面积较小(占总面积的28.3%)，不适宜做渗透滞留设施，且建筑密度较大、人群活动量大，雨水应尽快排出场地，因此将建立储存利用系统作为设计策略。“末端调蓄”设施只有雨水调蓄池，调蓄池内收集的雨水用于西十地块浇灌。
2.3 设计内容
综上所述，西十地块及外源雨水地块低影响开发系统分为渗透滞留和储蓄利用两部部分。渗透滞留系统包含种植屋面、透水铺装、下凹式绿地、植草沟、生物滞留池5类LID设施；储蓄利用系统以雨水调蓄池为储蓄设施收集地表径流。低影响系统设计见图 4，各类LID设施设计要点见表 1，各类用地面积及用地比例见表 2和表 3。
表 1 各类LID设施设计要点Table 1. Design points of LID facilities区域Area 低影响开发设施LID facility 设计要点Design point 西十地块West 10 Plot 种植屋面Planting roof 将具备条件的平屋顶(6 495.1 m2)改建为可上人的屋顶花园，为办公人员提供多样活动空间的同时，在急降雨时可滞存一部分雨水以达到推迟峰现时间的目的The flat roof (6 495.1 m2) which can be converted into a roof garden provides office workers with a variety of activity spaces and retains rainwater in case of emergency rainfall to delay the peak time 透水铺装Permeable pavement 将场地中工业废弃物加工成透水混凝土砖，应用于地块内主园路(6 647.8 m2)Permeable concrete bricks which are processed from industrial wastes in the site applied to the main garden road (6 647.8 m2) in the plot 植草沟Vegetative swale 沿硬质边缘形成线性景观带(425.4 m2)，以接纳铺装及绿地上的地表径流，将径流输送至雨水花园内部A linear landscape belt (425.4 m2) is formed along the hard edge to accept the surface runoff from pavement and green space and transport the runoff to the rain garden 下凹式绿地Rain garden 选择地块内部地势较低的区域集中布置两处下凹绿地(1 090.5 m2)，一处位于建筑下沉庭院内，可就近接纳周边绿地及铺装上的地表径流；另一处利用下洼料场基址设计成90余米长的带状雨水花园，其为地表终端LID设施，接纳场地内全部径流Two concave green spaces (1 090.5 m2) are arranged centrally in the low-lying area inside the plot. One is located in the sunken courtyard of the building, which can receive the runoff of surrounding green space and pavement. Another site is designed as a belt rain garden more than 90 m long, which is the surface terminal LID facility and accepts all runoff in the site 生物滞留池Bio-retention cell 在人群活动量相对较小的地势低洼区域布置了1 035.8 m2的生物滞留池，当强降雨发生后，生物滞留池中一定时间内的雨水积蓄不影响场地的正常使用In low-lying areas with relatively small amount of population activity, 1 035.8 m2 of bio-retention cell was arranged. When heavy rainfall occurred, rainwater accumulation in the bio-retention cell within a certain period of time did not affect the normal use of the site 外源雨水地块Exogenous rainwater area 透水铺装Permeable pavement 将场地中工业废弃物加工成透水混凝土砖，应用于秀池北街、秀池东路人行道(7 422.0 m2)Permeable concrete bricks which are processed from industrial wastes in the site applied to the sidewalks of Xiuchi North Street and Xiuchi East Road (7 422.0 m2) 植草沟Vegetative swale 沿硬质边缘形成线性景观带(443.1 m2)，以接纳铺装及绿地上的地表径流，将径流输送至雨水调蓄池中A linear landscape belt (443.1 m2) is formed along the hard edge to accept the surface runoff from pavement and green space and transport the runoff to the rain retention barrel 雨水调蓄池Rain retention barrel 外源雨水场地低处设置有雨水调蓄池，以接纳建筑屋顶雨水及全部地表径流。以设计目标为根据进行计算，外源雨水地块雨水调蓄池总容量为4 189 m3，方可达到消纳北京市3年重现期24 h降雨量这一目标。雨水调蓄池内雨水用于西十地块植物灌溉，并设有溢水口，当降雨强度超过地块内部LID设施消纳能力时，雨水通过溢水口排入市政雨水管网A rain retention barrel is installed at the lower level of the exogenous rainwater area to receive rainwater from the roof of the building and all surface runoff. According to the calculation based on the design objective, the total capacity of rain retention barrel in the exogenous rainwater area was 4 189 m3, which can reach the goal of absorbing 24-hour rainfall in the 3-year recurrence period of Beijing. Rainwater in the rain retention barrel is used for plant irrigation in the west 10 plot, and there is an overflow outlet. When the rainfall intensity exceeds the capacity of LID facilities inside the block, rainwater was discharged into the municipal rainwater pipe network through the overflow outlet 表 2 西十地块用地面积及用地比例表Table 2. Table of West 10 Plot of land area and land use ratio用地类型Type of land use 项目类型Item type 面积Area/m2 总面积Total area/m2 子项比例Subitem proportion/% 比例Proportion/% 绿化用地Green space 普通绿地Common green space 15 737.9 18 371.8 37.1 43.4 植草沟Vegetative swale 425.4 1.1 生物滞留池Bio-retention cell 1 035.8 2.4 下沉式绿地Rain garden 1 172.7 2.8 建筑占地Building space 屋顶花园绿地Green space of roof garden 4 871.3 14 987.3 11.5 35.7 屋顶花园铺装Impermeable pavement of roof garden 1 623.8 3.8 不上人硬质屋顶Hard roof of no people 8 492.2 20.4 园路及铺装场地用地Path and square 透水铺装Permeable pavement 6 647.8 8 954.4 15.6 20.8 不透水铺装Impermeable pavement 2 306.6 5.2 水体Water space 景观水体(不具备调蓄功能) Landscape water space (no regulating function) 52.3 52.3 0.1 0.1 合计Total 42 365.8 100 表 3 外源雨水地块用地面积及用地比例表Table 3. Table of exogenous site of land area and land use ratio用地类型Type of land use 项目类型Item type 面积Area/m2 总面积Total area/m2 子项比例Subitem proportion/% 比例Proportion/% 绿化用地Green space 普通绿地Common green space 15 967.6 16 410.7 27.5 28.3 植草沟Vegetative swale 443.1 0.8 建筑占地Building space 不上人硬质屋顶Hard roof of no people 16 865.1 16 865.1 29.1 29.1 园路及铺装场地用地Path and square 透水铺装Permeable pavement 7 422.0 24 740.1 12.8 42.6 不透水铺装Impermeable pavement 17 318.1 29.8 雨水调蓄池(地下) Rain retention barrel 雨水调蓄池(地下) Rain retention barrel (underground) 4 189 4 189 合计Total 58 015.9 100 3. 研究方法
按照《新首钢高端产业综合服务区控制性详细规划》要求完成西十地块及外源雨水地块低影响开发系统设计，并对上述设计方案进行SWMM模型构建。构建西十地块与外源雨水地块两套体系，西十地块体系用于对比传统模型和低影响开发模型，以分析低影响开发系统的雨洪调控效果；外源雨水地块体系用于对比雨水调蓄池蓄水量与植物灌溉需水量，以分析两者之间的数量关系。
3.1 渗透滞留系统模型建立
通过SWMM构建低影响开发模型和传统开发模型，对比两者之间数据以分析首钢西十低影响开发系统对区域径流总量、峰值流量和峰现时间的调控效果。
3.1.1 低影响开发模型构建
研究区域为西十地块及外源雨水地块，主要包含3类下垫面——屋顶、铺装、绿地。其中：不透水下垫面有硬质屋顶、不透水铺装；透水下垫面有绿色屋顶、透水铺装、绿地。根据上述低影响开发系统设计策略，遵循概化原则，将西十地块及外源雨水地块共划分为21个子汇水单元，末端出水口为2个。
3.1.2 传统开发模型构建
传统开发模式为通过竖向设计将绿地及屋顶径流排到铺装或道路上，再通过雨排水管线排入市政管网(图 5)。为方便后期对比分析，传统开发模型构建时绿地及铺装面积与低影响开发模型中一致(表 4)。雨排水管线图见图 6，SWMM模型概化图见图 7。
![]()
图 6 首钢西十传统雨水系统平面图Figure 6. Plan of traditional storm-water system in Shougang West 10 Plot![]()
图 7 首钢西十传统雨水系统概化图Figure 7. Generalized map of traditional storm-water system in Shougang West 10 Plot表 4 传统开发模式用地面积及用地比例表Table 4. Table of land area and land use ratio in traditional development model项目类型Item type 面积Area/m2 比例Proportion/% 绿化用地Green space 34 782.5 34.7 建筑占地Building space 31 852.4 31.7 园路及铺装场地用地Path and square 33 694.5 33.5 水体Water space 52.3 0.1 雨水调蓄池(地下) Rain retention barrel (underground) 4 189 合计Total 100 381.7 100 3.1.3 模块参数设置
本研究在推求各子汇水区的产流过程时，选用霍顿(Horton)公式进行下渗计算。根据方案设计对不同下垫面进行设定，各项参数参考SWMM用户手册[12]中的推荐值及相关文献。汇水模型为非线性水库模型，水力模型选用动力波模型[11]。根据场地土壤特性，入渗率为12.6 mm/h，衰弱常数为4 h-1。透水和不透水地表糙率系数分别为0.15和0.011，洼蓄量均为3 mm；植草沟宽度为2 000 mm，深度为500 mm；下沉式绿地平均深度为800 mm；生物滞留池平均深度为1 000 mm。
3.2 储蓄利用系统模型建立
根据北京市2017年全年降雨量日值统计(http://www.meteoblue.com/en/weather/archive/export/beijing-capital-international-airport_china_6301354)，通过SWMM模拟得出西十地块及外源雨水场地每日径流总量及峰值流量；整合数据得到每周雨水调蓄池内收集的径流总量，同计算得出的浇灌需水量进行对比，以周为单位分析外源雨水地块收集的径流与西十地块的浇灌需求之间的关系。
4. 结果与分析
4.1 渗透滞留系统结果与分析
4.1.1 设计降雨选择
SWMM模型可使用的降雨数据为两种类型：(1)实测历时降雨资料；(2)设计降雨[13]。本研究根据北京《雨水控制与利用工程设计规范》(DB11/685—2013)北京地区1年一遇、2年一遇、3年一遇、5年一遇、10年一遇、20年一遇重现期降雨量分别为45、81、108、141、209、270 mm，降雨历时为24 h，模拟时间48 h。
4.1.2 结果与分析
渗透滞留系统在重现期为1、2、3、5、10、20年的24 h降雨情况下，西十地块和外源雨水地块的低影响开发模型与传统开发模型的径流总量、峰值径流量和峰现时间分别见表 5和表 6；西十地块出口径流曲线对比图见图 8~13；外源雨水地块雨水调蓄池内部容积变化曲线对比图见图 14~19。根据上述数据结构可得出以下结果：
表 5 西十地块不同降雨强度下的出水口径流量模拟结果Table 5. Simulated results of West 10 Plot of the runoff under different rainfall intensities模拟项目Simulation item 1年一遇(45 mm)24 h 1-year-once(45 mm)24 hour 2年一遇(81 mm)24 h 2-year-once(81 mm)24 hour 3年一遇(108mm)24 h 3-year-once(108 mm)24 hour 5年一遇(141 mm)24 h 5-year-once(141 mm)24 hour 10年一遇(209 mm)24 h 10-year-once(209 mm)24 hour 20年一遇(270 mm)24 h 20-year-once(270 mm)24 hour 低影响开发模式LID model 传统开发模式Traditional development model 差值Difference 低影响开发模式LID model 传统开发模式Traditional development model 差值Difference 低影响开发模式LID model 传统开发模式Traditional development model 差值Difference 低影响开发模式LID model 传统开发模式Traditional development model 差值Difference 低影响开发模式LID model 传统开发模式Traditional development model 差值Difference 低影响开发模式LID model 传统开发模式Traditional development model 差值Difference 径流总量Total runoff/m3 0 1 189 1 189 0 2 411 2 411 0 3 490 3 490 1 028 4 846 3 818 3 818 7 699 3 881 6 712 10 284 3 572 峰值流量Peak flow/(m3·s-1) 0 0.027 0.027 0 0.065 0.065 0 0.092 0.092 0.037 0.124 0.087 0.128 0.192 0.064 0.203 0.252 0.049 峰现时间Time of peak flow occurrence 消失Disappea- rance 09：01 消失Disappea- rance 09：00 消失Disappea- rance 09：00 15：15 09：00 375 min 14：10 09：00 310 min 13：00 09：00 240 min 表 6 外源雨水地块不同降雨强度下的出水口径流量模拟结果Table 6. Simulated results of exogenous site of the runoff under different rainfall intensities模拟项目Simulation item 1年一遇(45 mm)24 h 1-year-once(45 mm)24 hour 2年一遇(81 mm)24 h 2-year-once(81 mm)24 hour 3年一遇(108 mm)24 h 3-year-once(108 mm)24 hour 5年一遇(141 mm)24 h 5-year-once(141 mm)24 hour 10年一遇(209 mm)24 h 10-year-once(209 mm)24 hour 20年一遇(270 mm)24 h 20-year-once(270 mm)24 hour 低影响开发模式LID model 传统开发模式Traditional development model 差值Difference 低影响开发模式LID model 传统开发模式Traditional development model 差值Difference 低影响开发模式LID model 传统开发模式Traditional development model 差值Difference 低影响开发模式LID model 传统开发模式Traditional development model 差值Difference 低影响开发模式LID model 传统开发模式Traditional development model 差值Difference 低影响开发模式LID model 传统开发模式Traditional development model 差值Difference 径流总量Total runoff/m3 800 1 806 1 006 2 633 3 599 966 4 189 5 156 967 6 145 7 108 963 10 249 11 205 956 13 970 14 912 942 峰值流量Peak flow/(m3·s-1) 0.026 0.040 0.014 0.079 0.093 0.014 0.116 0.132 0.016 0.162 0.179 0.017 0.256 0.275 0.019 0.341 0.362 0.021 峰现时间Time of peak flow occurrence 12：20 09：01 119 min 10：45 09：00 105 min 10：15 09：00 75 min 10：10 09：00 70 min 10：00 09：00 60 min 09：55 08：58 57 min ![]()
图 8 西十地块1年一遇24小时降雨出口径流曲线Figure 8. Runoff curve of 1-year-once 24 hour rainfall event in West 10 Plot![]()
图 9 西十地块2年一遇24小时降雨出口径流曲线Figure 9. Runoff curve of 2-year-once 24 hour rainfall event in West 10 Plot![]()
图 10 西十地块3年一遇24小时降雨出口径流曲线Figure 10. Runoff curve of 3-year-once 24 hour rainfall event in West 10 Plot![]()
图 11 西十地块5年一遇24小时降雨出口径流曲线Figure 11. Runoff curve of 5-year-once 24 hour rainfall event in West 10 Plot![]()
图 12 西十地块10年一遇24小时降雨出口径流曲线Figure 12. Runoff curve of 10-year-once 24 hour rainfall event in West 10 Plot![]()
图 13 西十地块20年一遇24小时降雨出口径流曲线Figure 13. Runoff curve of 20-year-once 24 hour rainfall event in West 10 Plot![]()
图 14 外源雨水地块1年一遇24小时降雨雨水调蓄池瞬时容积变化曲线Figure 14. nstantaneous volume change curve of rain retention barrel of 1-year-once 24 hour rainfall event in exogenous site![]()
图 15 外源雨水地块2年一遇24小时降雨雨水调蓄池瞬时容积变化曲线Figure 15. Instantaneous volume change curve of rain retention barrel of 2-year-once 24 hour rainfall event in exogenous site![]()
图 16 外源雨水地块3年一遇24小时降雨雨水调蓄池瞬时容积变化曲线Figure 16. Instantaneous volume change curve of rain retention barrel of 3-year-once 24 hour rainfall event in exogenous site![]()
图 17 外源雨水地块5年一遇24小时降雨雨水调蓄池瞬时容积变化曲线Figure 17. Instantaneous volume change curve of rain retention barrel of 5-year-once 24 hour rainfall event in exogenous site![]()
图 18 外源雨水地块10年一遇24小时降雨雨水调蓄池瞬时容积变化曲线Figure 18. Instantaneous volume change curve of rain retention barrel of 10-year-once 24 hour rainfall event in exogenous site![]()
图 19 外源雨水地块20年一遇24小时降雨雨水调蓄池瞬时容积变化曲线Figure 19. Instantaneous volume change curve of rain retention barrel of 20-year-once 24 hour rainfall event in exogenous site(1) 根据表 3，通过低影响开发模型与传统开发模型结果数据对比，可知经过西十地块地表低影响开发设施后，在1年一遇、2年一遇、3年一遇、5年一遇、10年一遇及20年一遇降雨强度下出水口的径流总量分别削减100%、100%、100%、78.79%、50.41%、34.73%，峰值流量分别减少100%、100%、100%、70.12%、33.33%、19.44%，峰现时间为1年一遇、2年一遇、3年一遇峰值消失，5年一遇、10年一遇及20年一遇降雨强度下分别延迟375、310、240 min。由此得出结论：低影响开发模式可减少径流总量、峰值径流量，推迟峰现时间。外源雨水场地内LID设施除透水铺装外只有植草沟。由表 4可知，对比传统开发模式，低影响开发模式对场地径流总量的削减量较小，1年一遇、2年一遇、3年一遇、5年一遇、10年一遇、20年一遇重现期径流总量分别为800、2 633、4 189、6 145、10 249、13 970 m3，因此可将这部分径流收集到雨水调蓄池中用于西十地块浇灌。
(2) 根据图 8~19可知：低影响开发模式在降雨前期对地表瞬时径流量的削减能力较强；随着时间的推迟削减能力逐渐减弱，在降雨后期低影响开发设施承载容量趋于饱和，此时低影响开发模式下的出口径流量高于传统开发模式。
4.2 储蓄利用系统结果与分析
4.2.1 储蓄结果与分析
根据SWMM模拟结果，外源雨水地块全年总蓄水量为9 617 m3，径流收集主要集中在第25~35周；单周最大收水总量为2 078 m3，发生在第27周(6月初)；第1~20周、43~52周收水量基本为0。
4.2.2 浇灌利用计算与分析
根据植物种类及北京气候类型，可推算出园林植物需水量：
PWR=ET0×KL (1) 式中：PWR为园林植物需水量，mm/d；ETO为参照作物需水量，mm/d；KL为园林系数。
根据彭曼-蒙蒂斯公式计算参照作物需水量：
ETo=0.408Δ(Rn−G)+γ900T+273u2(es−ea)Δ+γ(1+0.34u2) (2) 式中：ETO为参照作物需水量，mm/d；Rn为输入冠层净辐射量，MJ/m2·d；G为土壤热通量，MJ/m2·d；T为2 m高处日平均温度，℃；u2为2 m高处风速，m/s；es为饱和水汽压，kPa；ea为实际水汽压，kPa；Δ为饱和水汽压与温度关系曲线在某处得斜率，kPa/℃；γ为干湿温度计常数，kPa/℃。
园林系数(KL)的计算公式为：
KL=KsKdKme (3) 式中：Ks为园林植物种类因子；Kd为园林种植密度因子；Kmc为园林小气候因子。
园林植物净需水量：
PWRnet=PWR−Re (4) 式中：PWRnet为园林植物净需水量，mm/d；Re为有效降雨量，mm/d。
园林植物净灌溉需水量(IWRnet)是指存储于植物根区，可添加于土壤水库以满足净植物需水量(PWRnet)的那部分灌溉水量。
IWRnet=PWRnet (5) 为简化运算过程，将已报道的不同类型植物的种类系数、密度因子、小气候因子平均值[14](表 7)带入公式(3)得到冷季型草园林系数(KL)为0.8，其余植被类型园林系数(KL)为0.5。本研究参照北京市房山地区气象站20年(1982—2002年)常规观测气象数据，根据公式(2)计算出每月平均ETO值[16]，估算出园林植物需水量PWR(表 8)。园林植物在不降雨或降雨不足时需要人工灌溉。将上述结果带入公式(1)计算得出不降雨时园林植物需水量速率；结合不同植被类型面积求得2017年每周西十场地园林植物灌溉需水量；与雨水调蓄池周容积量进行对比可知每周雨水利用率及雨水调蓄池空置率如图 20所示。由此，计算得出第25周、第27周、第29~32周、第34~35周、第41周雨水调蓄池中收集水量在满足当月灌溉需求后仍有余量；全年径流收集的利用率为43.65%；在冬季雨水调蓄池空置率较大。
![]()
图 20 2017年每周植物需水量与雨水调蓄池容积量关系图Figure 20. Diagram of weekly plant water requirement and rainwater storage of rain retention barrel in 2017表 7 各类园林植物园林系数相关指标及面积Table 7. Relevant indicators and area of KL for all kinds of garden plants植被类型Vegetation type 种类系数Coefficient of species(Ks) 密度因子Density factor(Kd) 小气候因子Microclimatic factor(Kmc) 面积Area/m2 乔木Tree 0.9 0.5 0.2 1.3 1.0 0.5 1.4 1.0 0.5 灌木Shrub 0.7 0.5 0.2 1.1 1.0 0.5 1.3 1.0 0.5 地被植物Ground plant 0.9 0.5 0.2 1.1 1.0 0.5 1.2 1.0 0.5 12 400.5 树木、灌木、地被植物混合Mixed trees, shrubs and ground plants 0.9 0.5 0.2 1.3 1.0 0.6 1.4 1.0 0.5 冷季型草Cold season type grass 0.8 1.0 1.0 0.6 1.2 1.0 0.8 5 971.3 合计Total 18 371.8 表 8 2017年每月植物需水量Table 8. Monthly plant water requirement in 2017时间Time 1月Jan. 2月Feb. 3月Mar. 4月Apr. 5月May. 6月Jun. 7月Jul. 8月Aug. 9月Sep. 10月Oct. 11月Nov. 12月Dec. ETO/(mm·d-1) 0.72 1.20 2.17 3.38 4.26 4.72 4.21 3.74 3.02 2.00 1.19 0.72 PWR/(mm·d-1) 0.43 0.72 1.30 2.02 2.55 2.82 2.52 2.24 1.81 1.20 0.71 0.43 5. 结论
按照西十地块及外源雨水地块低影响开发系统雨水消纳能力应达到北京市3年一遇24 h降雨量(108 mm[9])的标准，经过上述设计研究及对SWMM模拟结果的定量分析可得出以下结论：
(1) 面对相同的雨洪调控目标，本研究根据基底条件和景观需求的差异制定不同的设计策略。西十地块优先构建渗透滞留系统，结合场地现状基底构建425.4 m2植草沟、1 035.8 m2生物滞留池、1 172.7 m2下沉式绿地、4 871.3 m2绿色屋顶以及6 647.8 m2透水铺装等5类LID设施；外源雨水地块以储蓄利用为目的，末端构建4 189 m3雨水调蓄池，最终均达到设计目标要求。
(2) 低影响开发设施对雨洪问题具有良好的调控效果，可以减少径流总量、峰值流量，推迟峰现时间，但随着时间的推迟削减能力逐渐减弱。由此可知，低影响开发系统的雨洪调控能力具有上限。
(3) 根据北京市2017年全年降雨量日值统计，外源雨水地块雨水调蓄池中径流收集的利用率为43.65%，有30周(210 d)雨水调蓄池空置率为100%，较为浪费。因此在进行低影响开发系统设计时，应综合考虑成本造价等因素。
-
![]()
图 1 3种因素对红花玉兰花蕾挥发油萃取率影响的单因素试验结果
字母不同表示各因素对挥发油萃取率影响差异显著（P < 0.05）。Different letters indicate that the effects of various factors on the extraction rate of essential oils are significantly different (P < 0.05).
Figure 1. Single factor experiment results of three factors on extraction yield of essential oils from the buds of Magnolia wufengens
![]()
图 2 红花玉兰花蕾挥发油GCMS总离子流图
Figure 2. Total ion current map of essential oils from the buds of M. wufengensis
表 1 响应面试验设计方案及结果
Table 1 Design and results of response surface experiments
试验号
Experiment No.萃取温度
Extraction temperature/℃萃取压力
Extraction pressure/MPa动态萃取时间
Dynamic extraction time/min挥发油萃取率
Extraction yield of essential oil/%1 − 1 (40) 0 (40) − 1 (60) 1.034 2 − 1 0 1 (120) 1.091 3 0 (45) 0 0 (90) 1.207 4 0 1 (45) − 1 1.116 5 1 (50) 0 − 1 1.145 6 0 0 0 1.203 7 − 1 1 0 1.066 8 0 − 1 (35) 1 1.177 9 0 1 1 1.172 10 1 0 1 1.173 11 1 1 0 1.178 12 1 − 1 0 1.126 13 0 0 0 1.218 14 0 − 1 − 1 1.103 15 − 1 − 1 0 1.048 注：括号内数字表示各编码值代表的实际参数值。Note：the numbers in parentheses represent the actual parameter values represented by each coded value. 
下载: 导出CSV
表 2 回归模型方差分析及显著性检验
Table 2 Variance analysis on regression models and significance test
方差来源
Source of variance平方和
Sum of squares自由度
DF均方
Mean squareF值
F valueP值
P value模型 Model 0.048 9 5.357 × 10− 3 32.39 0.000 7 A 0.018 1 0.018 110.77 0.000 1 B 7.508 × 10− 4 1 7.508 × 10− 4 4.54 0.086 3 C 5.784 × 10− 3 1 5.784 × 10− 3 34.97 0.002 0 AB 2.907 × 10− 4 1 2.907 × 10− 4 1.76 0.242 2 AC 2.088 × 10− 4 1 2.088 × 10− 4 1.26 0.312 2 BC 7.921 × 10− 5 1 7.921 × 10− 5 0.48 0.519 7 A2 0.017 1 0.017 103.24 0.000 2 B2 5.062 × 10− 3 1 5.062 × 10− 3 30.61 0.002 6 C2 3.407 × 10− 3 1 3.407 × 10− 3 20.60 0.006 2 残差 Residual 8.268 × 10− 4 5 1.654 × 10− 4 失拟项 Lack of fit 6.957 × 10− 4 3 2.319 × 10− 4 3.54 0.228 3 纯误差 Pure error 1.312 × 10− 4 2 6.559 × 10− 5 总变异 Total variation 0.049 14
下载: 导出CSV
表 3 红花玉兰花蕾挥发油化学成分
Table 3 Chemical compositions of essential oils from the buds of Magnolia wufengensis
序号
No.化合物名称
Compound name保留时间
Retention time/minaRI/bRI 相对含量
Relative content/%1 8-羟基芳樟醇 8-Hydroxylinalool 5.733 1 364/1 361 0.08 2 反式橙花叔醇 Nerolidol 8.21 1 570/1 564 0.06 3 3,4-二甲氧基苯甲酸甲酯
Benzoic acid, 3,4-dimethoxy-, methyl ester8.535 1 597/1 592 0.13 4 石竹烯氧化物 Caryophyllene oxide 8.636 1 606/1 581 0.14 5 脱氢风毛菊内脂 Dehydrosaussurea lactone 11.364 1 856/1 838 0.19 6 苯甲酸苯乙酯 Benzoic acid, 2-phenylethyl ester 11.534 1 872/1 856 0.06 7 十九烷 Nonadecane 11.822 1 900/1 900 0.23 8 二十烷 Eicosane 12.806 2 000/2 000 0.35 9 二十一烷 Heneicosane 13.754 2 100/2 100 5.38 10 二十二烷 Docosane 14.647 2 200/2 200 1.61 11 裂叶苣荚莱内酯 Santamarine 14.72 2 209/2 227 0.36 12 二十烷醛 Eicosanal 14.894 2 229/2 221 0.24 13 喘诺木烯内酯 Reynosin 14.994 2 241/2 266 0.22 14 诱虫烯 9-Tricosene,（9Z）- 15.305 2 277/2 272 0.35 15 二十三烷 Tricosane 15.521 2 300/2 300 8.77 16 4H-xanthalongia 15.684 2 322/- 5.41 17 对甲氧基肉桂酸异辛酯
2-Ethylhexyl trans-4-methoxycinnamate15.776 2 333/2 339 0.07 18 二十四烷 Tetracosane 16.336 2 400/2 400 1.34 19 二十二烷醛 Docosanal 16.598 2 434/2 430 0.51 20 二十五碳烯 Pentacos-1-ene 16.949 2 478/2 488 0.44 21 二十五烷 Pentacosane 17.142 2 500/2 500 7.79 22 二十六烷 Hexacosane 17.897 2 600/2 600 0.88 23 二十四烷醛 Tetracosanal 18.176 2 636/2 632 0.24 24 二十七碳烯 Heptacos-1-ene 18.565 2 686/2 688 1.44 25 二十八烷 Octacosane 19.574 2 800/2 800 1.15 26 二十六烷醛 Hexacosanal 19.979 2 839/2 815 0.47 27 二十八烷醛 Octacosanal 22.498 3 042/3 039 0.32 28 维生素E dl-α-Tocopherol 24.232 3 146/3 150 0.21 29 16-三十一酮 16-Hentriacontanone 27.283 3 291/3 304 7.82 30 γ-谷淄醇 γ-Sitosterol 28.398 3 332/3 351 24.39 注：aRI/bRI代表计算所得目标成分的保留指数/谱库中目标成分的保留指数。Note：aRI/bRI represents calculated retention index of target component/retention index of target component in the library.
下载: 导出CSV
表 4 红花玉兰花蕾挥发油抗氧化活性
Table 4 Antioxidant activities of essential oils from Magnolia wufengensis buds
方法
Method样品
Sample线性回归方程
Regression equation相关系数R2 IC50/（mg·mL− 1） AEAC /mg DPPH自由基清除法
DPPH free radical scavenging挥发油 Essential oils y = 3.478 7x + 5.408 1 0.995 4 12.82 ± 0.243 7 198.64 抗坏血酸 Ascorbic acid y = 218 4x − 5.609 5 0.998 8 0.025 ± 0.000 7 — 还原能力法
Reducing power挥发油 Essential oil y = 0.006x + 0.089 2 0.977 9 — — 抗坏血酸 Ascorbic acid y = 16.19x + 0.062 0 0.999 8 — —
下载: 导出CSV
表 5 红花玉兰花蕾挥发油抑菌活性测定结果
Table 5 Antibacterial activities of essential oils from Magnolia wufengensis buds
细菌 Bacteria species 纸片扩散法 Disk diffusion method 微量肉汤稀释法 Broth microdilution method 类别 Category 名称 Name 挥发油
Essential oil/mm氨苄青霉素 Ampicillin/mm MIC/(mg·mL− 1) MBC/(mg·mL− 1) 革兰氏阳性
Gram-positive枯草芽孢杆菌 B. subtilis 10.69 ± 1.06 ab 29.83 ± 0.88 b 12.5 > 25 金黄色葡萄球菌
S. aureus12.22 ± 0.83 a 37.01 ± 1.70 a 12.5 革兰氏阴性
Gram-negative鼠伤寒沙门氏菌
S. typhimurium8.26 ± 0.34 b 19.41 ± 0.92 c 25 大肠杆菌 E. coli 8.99 ± 0.38 b 14.29 ± 0.43 d 25 注：字母不同表示氨苄或挥发油对不同细菌抑制作用差异显著（P < 0.05）。Note: different letters indicate that the inhibitory effects of ampicillin or essential oils on different bacteria are significantly different (P < 0.05).
下载: 导出CSV
-
[1] Bayala B, Bassole I H N, Gnoula C, et al. Chemical composition, antioxidant, anti-Inflammatory and anti-proliferative activities of essential oils of plants from Burkina Faso[J/OL]. PLoS One, 2014, 9(3): e92122[2019−05−25]. https://doi.org/10.1371/journal.pone.0092122.
[2] Ahmad S, Ullah F, Sadiq A, et al. Chemical composition, antioxidant and anticholinesterase potentials of essential oil of Rumex hastatus D. Don collected from the north west of Pakistan[J/OL]. BMC Complementary and Alternative Medicine, 2016, 16: 29[2019−05−25]. https://doi.org/10.1186/s12906-016-0998-z.
[3] Almadiy A A, Nenaah G E, Al Assiuty B A, et al. Chemical composition and antibacterial activity of essential oils and major fractions of four Achillea species and their nanoemulsions against foodborne bacteria[J]. LWT - Food Science and Technology, 2016, 69: 529−537. doi: 10.1016/j.lwt.2016.02.009
[4] Boulanouar B, Abdelaziz G, Aazza S, et al. Antioxidant activities of eight Algerian plant extracts and two essential oils[J]. Industrial Crops and Products, 2013, 46: 85−96. doi: 10.1016/j.indcrop.2013.01.020
[5] Zhang L Y, Yang Z W, Wei J W, et al. Essential oil composition and bioactivity variation in wild-growing populations of Curcuma phaeocaulis Valeton collected from China[J]. Industrial Crops and Products, 2017, 103: 274−282. doi: 10.1016/j.indcrop.2017.04.019
[6] Riahi L, Elferchichi M, Ghazghazi H, et al. Phytochemistry, antioxidant and antimicrobial activities of the essential oils of Mentha rotundifolia L. in Tunisia[J]. Industrial Crops and Products, 2013, 47: 883−889.
[7] Abdelli M, Moghrani H, Aboun A, et al. Algerian Mentha pulegium L. leaves essential oil: chemical composition, antimicrobial, insecticidal and antioxidant activities[J]. Industrial Crops and Products, 2016, 94: 197−205. doi: 10.1016/j.indcrop.2016.08.042
[8] Benito-Román O, Rodríguez-Perrino M, Sanz M T, et al. Supercritical carbon dioxide extraction of quinoa oil: study of the influence of process parameters on the extraction yield and oil quality[J]. The Journal of Supercritical Fluids, 2018, 139: 62−71. doi: 10.1016/j.supflu.2018.05.009
[9] 秦娇, 姚永红, 侯智霞, 等. 响应面分析法优化毛竹叶挥发油提取工艺[J]. 食品科学, 2010, 31(6):1−5. Qin J, Yao Y H, Hou Z X, et al. Response surface methodology as an approach to optimization of extraction of volatile oil from leaves of Phyllostachys heterocycla (Carr.) Mitford cv. Pubescens Mazel ex H. de leh[J]. Food Science, 2010, 31(6): 1−5.
[10] 马履一, 王罗荣, 贺随超, 等. 中国木兰科木兰属一新变种(英文)[J]. 植物研究, 2006(5):516−519. doi: 10.3969/j.issn.1673-5102.2006.05.002 Ma L Y, Wang L R, He S C, et al. A new species of magnolia from China[J]. Bulletin of Botanical Research, 2006(5): 516−519. doi: 10.3969/j.issn.1673-5102.2006.05.002
[11] 桑子阳, 马履一, 陈发菊, 等. 五峰红花玉兰种质资源保护现状与开发利用对策[J]. 湖北农业科学, 2011, 50(8):1564−1567. doi: 10.3969/j.issn.0439-8114.2011.08.017 Sang Z Y, Ma L Y, Chen F J, et al. Protection status and utilization countermeasure of germplasm resources of Magnolia wufengensis in Wufeng County[J]. Hubei Agricultural Sciences, 2011, 50(8): 1564−1567. doi: 10.3969/j.issn.0439-8114.2011.08.017
[12] 马履一. 红花玉兰的选育及在国土绿化中的应用[J]. 国土绿化, 2019(3):54−56. Ma L Y. Breeding and application in land greening of Magnolia wufengensis[J]. Land Greening, 2019(3): 54−56.
[13] 张婷婷, 郭夏丽, 黄学勇, 等. 辛夷挥发油GC-MS分析及其抗氧化、抗菌活性[J]. 食品科学, 2016, 37(10):144−150. doi: 10.7506/spkx1002-6630-201610025 Zhang T T, Guo X L, Huang X Y, et al. GC-MS analysis and antioxidant and antimicrobial properties of volatile oil from flos Magnoliae[J]. Food Science, 2016, 37(10): 144−150. doi: 10.7506/spkx1002-6630-201610025
[14] Morshedloo M R, Quassinti L, Bramucci M, et al. Chemical composition, antioxidant activity and cytotoxicity on tumour cells of the essential oil from flowers of Magnolia grandiflora cultivated in Iran[J]. Natural Product Research, 2017, 31(24): 2857−2864. doi: 10.1080/14786419.2017.1303699
[15] Cherrat L, Espina L, Bakkali M, et al. Chemical composition and antioxidant properties of Laurus nobilis L. and Myrtus communis L. essential oils from Morocco and evaluation of their antimicrobial activity acting alone or in combined processes for food preservation[J]. Journal of the Science of Food and Agriculture, 2014, 94(6): 1197−1204. doi: 10.1002/jsfa.6397
[16] 韦献雅, 殷丽琴, 钟成, 等. DPPH法评价抗氧化活性研究进展[J]. 食品科学, 2014, 35(9):317−322. doi: 10.7506/spkx1002-6630-201409062 Wei X Y, Yin L Q, Zhong C, et al. Advances in the DPPH radical scavenging assay for antioxidant activity evaluation[J]. Food Science, 2014, 35(9): 317−322. doi: 10.7506/spkx1002-6630-201409062
[17] Harkat-Madouri L, Asma B, Madani K, et al. Chemical composition, antibacterial and antioxidant activities of essential oil of Eucalyptus globulus from Algeria[J]. Industrial Crops and Products, 2015, 78: 148−153. doi: 10.1016/j.indcrop.2015.10.015
[18] 崔丽丽, 逄世峰, 李亚丽, 等. 响应面试验优化超临界CO2萃取人参挥发油的工艺[J]. 食品科学, 2016, 37(4):58−61. doi: 10.7506/spkx1002-6630-201604010 Cui L L, Pang S F, Li Y L, et al. Optimization of supercritical CO2 fluid extraction of essential oil from Panax ginseng roots by response surface methodology[J]. Food Science, 2016, 37(4): 58−61. doi: 10.7506/spkx1002-6630-201604010
[19] 葛发欢, 李莹, 谢健鸣, 等. 超临界CO2从柴胡中萃取挥发油及其皂甙的研究[J]. 中国中药杂志, 2000(3):21−25. Ge F H, Li Y, Xie J M, et al. Studies on technology of supercritical CO2 fluid extraction for volatile oils and saikosaponins in Bupleurum chinense DC.[J]. China Journal of Chinese Materia Medica, 2000(3): 21−25.
[20] 李莉, 张赛, 何强, 等. 响应面法在试验设计与优化中的应用[J]. 实验室研究与探索, 2015, 34(8):41−45. doi: 10.3969/j.issn.1006-7167.2015.08.011 Li L, Zhang S, He Q, et al. Application of response surface methodology in experiment design and optimization[J]. Research and Exploration in Laboratory, 2015, 34(8): 41−45. doi: 10.3969/j.issn.1006-7167.2015.08.011
[21] Dhouioui M, Boulila A, Chaabane H, et al. Seasonal changes in essential oil composition of Aristolochia longa L. ssp. paucinervis Batt. (Aristolochiaceae) roots and its antimicrobial activity[J]. Industrial Crops & Products, 2016, 83: 301−306.
[22] 傅大立, 赵东欣, 孙金花, 等. 10种国产玉兰属植物挥发油成分及系统学意义[J]. 林业科学, 2005, 41(3):68−74. doi: 10.3321/j.issn:1001-7488.2005.03.012 Fu D L, Zhao D X, Sun J H, et al. Component of volatile oils in 10 Chinese species of Yulania genus and its phylogenetic implication[J]. Scientia Silvae Sinicae, 2005, 41(3): 68−74. doi: 10.3321/j.issn:1001-7488.2005.03.012
[23] Zheng Y F, Ren F, Liu X M, et al. Comparative analysis of essential oil composition from flower and leaf of Magnolia kwangsiensis Figlar and Noot[J]. Natural Product Letters, 2016, 30(13): 1552−1556. doi: 10.1080/14786419.2015.1113410
[24] Mitsuo M, Yoshimi N, Hiroshi N, et al. Volatile compounds with characteristic odor of essential oil from Magnolia obovata leaves by hydrodistillation and solvent-assisted flavor evaporation[J]. Journal of Oleo Science, 2015, 64(9): 999−1007. doi: 10.5650/jos.ess15114
[25] Pirbalouti A G, Oraie M, Pouriamehr M, et al. Effects of drying methods on qualitative and quantitative of the essential oil of Bakhtiari savory (Satureja bachtiarica Bunge.)[J]. Industrial Crops and Products, 2013, 46: 324−327. doi: 10.1016/j.indcrop.2013.02.014
[26] Rahimmalek M, Goli S A H. Evaluation of six drying treatments with respect to essential oil yield, composition and color characteristics of Thymys daenensis subsp. daenensis. Celak leaves[J]. Industrial Crops and Products, 2013, 42(1): 613−619.
[27] 胡长鹰. 当归挥发油的提取与成分分析[J]. 食品与机械, 2006, 22(2):24−26. doi: 10.3969/j.issn.1003-5788.2006.02.007 Hu C Y. Extraction and chemical components of volatile oil from Angelica sinensis[J]. Food and Machinery, 2006, 22(2): 24−26. doi: 10.3969/j.issn.1003-5788.2006.02.007
[28] 张冬阳, 郭雪松. 响应面法优化薏米中β-谷甾醇的提取工艺[J]. 食品工业科技, 2017, 38(22):206−209, 241. Zhang D Y, Guo X S. Optimization of extraction technology of β-sitosterol from Coix seed by response surface methodology[J]. Science and Technology of Food Industry, 2017, 38(22): 206−209, 241.
[29] 尚红梅, 郭玮, 潘丹, 等. 干燥方式对菊苣根多酚含量和抗氧化活性的影响[J]. 食品科学, 2015, 36(1):84−88. doi: 10.7506/spkx1002-6630-201501016 Shang H M, Guo W, Pan D, et al. Effects of drying methods on polyphenol contents and antioxidant activities of chicory root[J]. Food Science, 2015, 36(1): 84−88. doi: 10.7506/spkx1002-6630-201501016
[30] Teixeira B, Marques A, Ramos C, et al. Chemical composition and antibacterial and antioxidant properties of commercial essential oils[J]. Industrial Crops and Products, 2013, 43: 587−595. doi: 10.1016/j.indcrop.2012.07.069
[31] Vallverdú-Queralt A, Regueiro J, Martínez-Huélamo M, et al. A comprehensive study on the phenolic profile of widely used culinary herbs and spices: rosemary, thyme, oregano, cinnamon, cumin and bay[J]. Food Chemistry, 2014, 154: 299−307. doi: 10.1016/j.foodchem.2013.12.106
-
期刊类型引用(6)
1. 莫崇杏，董明亮，李荣生，余纽，郑显澄，杨锦昌. 米老排杂交子代苗期生长性状遗传变异及选择. 森林与环境学报. 2023(05): 555-560 . 
百度学术
2. Shuchun Li，Jiaqi Li，Yanyan Pan，Xiange Hu，Xuesong Nan，Dan Liu，Yue Li. Variation analyses of controlled pollinated families and parental combining ability of Pinus koraiensis. Journal of Forestry Research. 2021(03): 1005-1011 . 必应学术
3. 潘艳艳，许贵友，董利虎，王成录，梁德洋，赵曦阳. 日本落叶松全同胞家系苗期生长性状遗传变异. 南京林业大学学报(自然科学版). 2019(02): 14-22 . 百度学术
4. 秦光华，宋玉民，乔玉玲，于振旭，彭琳. 旱柳苗高年生长与气象因子的灰色关联度. 东北林业大学学报. 2019(05): 42-45+51 . 百度学术
5. 李峰卿，陈焕伟，周志春，楚秀丽，徐肇友，肖纪军. 红豆树优树种子和幼苗性状的变异分析及优良家系的初选. 植物资源与环境学报. 2018(02): 57-65 . 百度学术
6. 张素芳，张磊，赵佳丽，张莉，张含国. 长白落叶松小RNA测序和其靶基因预测. 北京林业大学学报. 2016(12): 64-72 . 本站查看
其他类型引用(6)
计量
- 文章访问数: 7224
- HTML全文浏览量: 866
- PDF下载量: 58
- 被引次数: 12
