本发明涉及介质阻挡放电(dielectricbarrierdischarge:dbd)反应器供电能量检测计算技术领域，特别是一种在功率密度调制(powerdensitymodulation:pdm)工作模式的高压电源激励dbd反应器供电能量检测系统和方法。
背景技术：
近几年来，大气压等离子体被广泛应用于环境修复，污染控制，生物医学，流动控制和材料处理等方面。dbd是有绝缘介质插入放电空间的一种非平衡态气体放电，具有放电稳定和放电面积大等特点，是在常压下产生低温等离子体的有效方法之一，适合于大规模工业生产，具有广阔的应用前景。在实际应用中，dbd的供电电压、供电电流、供电能量等电学参量，是低温等离子体研究和应用中极其重要的参量，在放电过程中它们相互作用，共同影响多重微放电的产生、放电的效果和等离子体的特性。准确测量并计算介质阻挡放电系统的有效能量将成为提高介质阻挡放电性能的关键所在，对于提高介质阻挡放电在应用中的放电效率具有重要意义。
由于dbd通过大量短暂的微放电来消耗能量，具有时空不确定性，故供电功率和能量很难测定。对于这类放电的能量计算，通常采用电压-电荷李萨如(lissajous)图形法。
专利201310132909.8公开了一种介质阻挡放电电路参数检测装置和检测方法。根据李萨如图形法设计了一种检测和计算dbd反应器等效参数的检测装置和检测方法。但是，该专利主要根据李萨如图形对反应器的等效参数进行计算，未对非连续、供电电压峰值不固定的放电状态进行分析，没有提出相应情况下供电能量计算的算法思想和放电效果评估方法。
目前，使用电压-电荷李萨如图形法能够计算在连续的高压交流电源激励下的dbd反应器供电能量。但是当放电不连续或者供电电压电流峰值不固定时，尚未有方法对反应器的供电能量进行准确计算。若将高压交流电源激励下dbd供电能量的计算方法直接运用到功率密度调制(powerdensitymodulation:pdm)电源激励下的反应器进行供电能量计算，存在以下问题：
(1)根据pdm高压电源的放电特性，在一段工作时间内，电压电流波形是不连续的；
(2)每个供电周期的放电强度和供电电压电流峰值不同；
(3)在一个调功周期内，放电时间是可调节的，因此在工作时间内的放电时间不固定。
因此，传统的李萨如图形法计算dbd供电能量不适合直接应用到pdm高压电源驱动dbd系统进行能量消耗分析，若能够提供一种针对使用pdm高压电源激励dbd的供电能量计算方法，将对dbd的研究、应用具有指导意义。
技术实现要素：
本发明要解决的技术问题为：针对使用pdm高压电源激励的dbd在放电过程中的电气和光学参数，获取相关参数，进而借助李萨如图形进行能量消耗分析，以更加准确地获得经过反应器的能量转换效率。
本发明采用的技术方案为：一种pdm工作模式的高压电源激励dbd反应器供电能量的检测系统，包括可编程电源、pdm高压电源、被测dbd阵列、电气参数检测模块、光纤光谱仪和上位机；
可编程电源的输出端连接pdm高压电源的输入端，pdm高压电源的输出端连接dbd阵列的激励电压输入端；上位机与可编程电源之间连接通信，以通过可编程电源控制pdm高压电源的工作，进而使得dbd阵列在pdm高压电源的激励下放电；
dbd阵列的放电回路包括一端连接地线的积分电容；
电气参数检测模块包括供电电压检测单元、积分电压检测单元和放电电流检测单元；供电电压检测单元的电压输入端采集dbd阵列激励电压输入端的供电电压；积分电压检测单元的电压输入端采集所述积分电容上的积分电压；放电电流检测单元的电流输入端采集dbd阵列的放电电流；供电电压检测单元、积分电压检测单元和放电电流检测单元的输出端分别连接上位机；
光纤光谱仪采集dbd阵列中放电区域中活性物种发射光谱的相对强度，传输至上位机中；
上位机根据获取到的供电电压、积分电容电压和放电电流，计算dbd的供电能量；上位机根据获取到的放电区域中活性物种发射光谱的相对强度，计算dbd的相对光量子产率能效比。
本发明所述可编程电源对pdm高压电源的控制包括输入电压控制、供电的功率密度控制和放电强度控制，分别为现有技术。
优选的，上位机利用李萨如图形算法进行dbd供电能量的计算：
上位机根据dbd阵列在pdm高压电源激励下产生的微放电数据，进行李萨如图形重构：选择供电电压上升到零的两个相邻点之间的时间段作为一个供电周期；定义单个供电周期中产生微放电的时段有ab时段和cd时段，分别位于单个供电周期的负半周期和正半周期中，则对应的李萨如图形的四个顶点坐标分别为a(ux1,uy1)、b(ux2,uy2)、c(ux3,uy3)、d(ux4,uy4)，其中各顶点的横坐标和纵坐标分别为相应时间点上的积分电容电压和供电电压；
计算结果包括：单个供电周期的能量为：
即单个供电周期的能量表示为积分电容与该供电周期对应的李萨如图形面积的乘积；其中，ed,i是第i个供电周期的供电能量，下标d,i表示供电周期的次数，td,i表示第i个供电周期的时间，u(t)表示供电电压，i(t)表示放电电流，um(t)指积分电容上的累积电荷电压，cm为积分电容的电容值，sd,i表示第i个供电周期的李萨如图形面积；
整个运行时间的dbd系统消耗总能量为：
式中，tm,j为pdm高压电源一个调功周期的时间；下标m,j表示功率密度调制周期的次数；ton,j为一个调功周期内的dbd放电时间，由多个供电周期td,j组成；占空比
dbd放电区域中活性物种的相对光量子产率能效比eer为：
其中i为光纤光谱仪采集的dbd放电区域中活性物种发射光谱的相对强度；em,j为一个调功周期的供电能量，即
优选的，所述可编程电源模块采用可编程ac/dc电源模块。可编程ac/dc电源用于设定pdm高压电源的输入电压及整个系统的功率。
优选的，所述电气参数检测模块采用数字示波器，数字示波器包括至少3个检测通道，3个检测通道分别设置高压探头、电压探头和电流互感器，并分别用于检测dbd阵列的高压电源输入端电压、积分电容上的电压和放电回路的电流。
优选的，所述电流互感器采用罗氏线圈。罗氏线圈以感应方式接在放电回路的地线上，获得放电电流，为现有技术。
本发明还公开基于上述pdm工作模式的高压电源激励dbd反应器供电能量的检测系统的供电能量检测方法，包括：
s1，系统初始化：将dbd的激励电源调节为pdm供电模式，设定pdm供电模式下的放电条件，使得dbd被激励放电；
s2，利用电压电流检测模块获取dbd阵列的放电参数，包括供电电压、放电电流，和积分电容上的积分电压；利用光纤光谱仪获取dbd阵列放电区域中活性物种发射光谱的相对强度；
s3，上位机从电压电流检测模块和光纤光谱仪获取已采集到的数据，基于获取到的数据进行dbd供电能量、等效参数、有效放电时间以及相对光量子产率能效比的计算和输出。
优选的，s1中，所述放电条件包括pdm高压电源供电电压、调功周期长度、供电频率和功率密度。
优选的，步骤s3包括：
s31，计算dbd的供电能量：
定义tm,j为pdm高压电源的一个调功周期，ton,j为dbd在pdm高压电源一个调功周期内的放电时间，由多个dbd供电周期组成；则dbd的单个供电周期能量ed,i表示为积分电容和李萨如图形面积的乘积，为：
式中，td,i表示单个供电周期的时间，u(t)是供电电压，i(t)是放电电流，d表示放电，i表示放电次数；
一个调功周期tm,j中供电时间的占空比为：
一个调功周期内的供电周期个数为：
一个调功周期的供电能量em,j即ndm,j个单个供电能量之和：
定义ed,a是调功周期内单个供电周期的平均供电能量，则：
当dbd系统工作在连续状态时，运行时间内总供电周期个数non,t可以表示为：
式中，tt是放电系统的运行总时间；
则系统运行时间内的dbd总供电能量et表示为：
当tt＝tm,j时，
根据式(1)、(6)、(7)和(8)、可得到et：
根据dbd阵列在pdm高压电源激励下产生的微放电数据，绘制激励dbd的供电能量的电压电荷李萨如图形，则单个供电周期的能量与相应李萨如图形的面积成正比；
由于：
其中cm是回路中外接积分电容的容值，um(t)为积分电容上的累积电荷电压；
则将(10)式带入(1)式可得单个供电周期的供电能量ed,i，表示为积分电容和相应李萨如图形面积sd,i的乘积，即：
将(11)代入(9)，利用李萨如图形面积来表示供电能量，可以得到运行时间内dbd总供电能量et为：
s32，计算dbd的等效电容：
定义单个供电周期中产生微放电的时段有ab时段和cd时段，分别位于单个供电周期的负半周期和正半周期中，对应的李萨如图形的四个顶点坐标分别为a(ux1,uy1)、b(ux2,uy2)、c(ux3,uy3)、d(ux4,uy4)，各顶点的横坐标和纵坐标分别为相应时间点上的积分电容电压和供电电压；
根据单个供电周期的李萨如图形，可以计算对应供电周期中dbd阵列的总等效电容c、介质等效电容cd和放电间隙等效电容cg分别为：
在一个调功周期中，dbd阵列的等效电容看做该调功周期内所有供电周期对应的等效电容的平均值。
s33计算量子产率能效比eer：
其中i为光纤光谱仪采集的dbd放电区域中活性物种发射光谱的相对强度。
优选的，s31在进行dbd的供电能量的计算时，首先获取检测到的供电电压数据和积分电容电压数据，然后根据获取到的数据，进行单个供电周期的李萨如图形重构，图形重构时，包括：
以电容积分电压为横坐标数据，以“供电电压”为纵坐标数据；得到所有供电周期堆叠的李萨如图形；
进行单个周期李萨如图形的分离，得到单个供电周期李萨如图形，分离规则为：选择供电电压上升到零的两个相邻点之间的时间段作为一个供电周期；
对分离得到的单个李萨如图形分别计算得到李萨如图形的面积，进而根据公式(11)得到相应供电周期的供电能量ed,i。
进一步的，s3还包括放电参数处理步骤：
根据获取到的供电电压参数，得到供电电压峰-峰值；
根据获取到的回路放电电流，得到dbd的平均微放电强度和有效放电时间，包括：
对获取的回路放电电流进行波形重构；
对重构形成的放电电流波形进行平滑处理；优选为采用savitzky-golay平滑算法；
将回路放电电流数据与平滑处理后的相应电流数据相减，得到放电过程中的微放电电流的波形数据；
根据微放电的电流波形的产生，得到每个供电周期的放电起始时间和放电截止时间，从而计算单个供电周期的有效放电时间；
将所有供电周期的有效放电时间累加，得到有效放电总时间；
根据每个供电周期的微放电脉冲峰值，计算放电过程中的平均微放电强度。
有益效果
与现有技术相比，本发明具有以下优点和进步：
(1)既能直接测反应器每个供电周期能量，也可以对单个供电周期的能量进行累加统计，得到任意时间段内消耗的能量。因此，本方法既能应用在连续放电状态中，又能应用在非连续放电状态中，此外，还可应用在供电峰-峰值不稳定的放电状态。
(2)定义了在获得的pdm供电波形后如何进行李萨如图形的重构，例如如何定义每个周期起始点与终止点。
(3)传统的系统运行效率评估是使用高压电源的输出端的能量和高压电源输入端的能量的比值作为能量转换效率，本文中的方法可以获得反应器上消耗的能量，因此更加准确地获得经过反应器的能量转换效率。
附图说明
图1所示为本发明在pdm高压电源激励下的dbd供电能量参数检测示意图；
图2所示为本发明在pdm高压电源激励下的dbd供电能量检测系统示意图；
图3所示为dbd阵列结构示意图；
图4所示为功率密度调制电源典型供电波形图；
图5所示为高压交流电源激励下的供电波形；
图6(a)所示为dbd典型李萨如图；图6(b)是调功周期内堆叠的李萨如图；
图7所示为检测步骤主流程示意图；
图8所示为参数处理主程序流程示意图；
图9所示为放电电流处理子程序流程示意图；
图10所示为供电能量计算子程序流程示意图；
图11所示为等效参数计算子程序流程示意图；
图12所示为光量子产率能效比计算子程序流程示意图；
图13所示为牛顿爬山算法示意图。
具体实施方式
以下结合附图和具体实施例进一步描述。
实施例1
参考图1和图2所示，pdm工作模式的高压电源激励的dbd反应器供电能量的检测系统，包括可编程电源、pdm高压电源、被测dbd阵列、电气参数检测模块和上位机(电脑)；
可编程电源的输出端连接pdm高压电源的输入端，pdm高压电源的输出端连接dbd阵列的激励电压输入端；上位机与可编程电源之间连接通信，以通过可编程电源控制pdm高压电源的工作，进而使得dbd阵列在pdm高压电源的激励下放电；
dbd阵列的放电回路包括一端连接地线的积分电容；
电气参数检测模块包括供电电压检测单元、积分电压检测单元和放电电流检测单元；供电电压检测单元的电压输入端采集dbd阵列激励电压输入端的供电电压；积分电压检测单元的电压输入端采集所述积分电容上的积分电压；放电电流检测单元感应dbd阵列的放电电流；供电电压检测单元、积分电压检测单元和放电电流检测单元的输出端分别连接上位机；
上位机根据获取到的供电电压、积分电容电压和放电电流，计算dbd的供电能量。
上位机利用李萨如图形算法进行dbd供电能量的计算，具体为，根据dbd阵列在pdm高压电源激励下产生的供电电压和积分电压数据，绘制电压电荷李萨如图形：选择供电电压上升到零的两个相邻点之间的时间段作为一个供电周期；定义单个供电周期中产生微放电的时段有ab时段和cd时段，分别位于单个供电周期的负半周期和正半周期中，对应的李萨如图形的四个顶点坐标分别为a(ux1,uy1)、b(ux2,uy2)、c(ux3,uy3)、d(ux4,uy4)，其中各顶点的横坐标和纵坐标分别为相应时间点上的积分电容电压和供电电压；
计算结果包括：单个供电周期的能量为：
即单个供电周期的能量表示为积分电容与该供电周期对应的李萨如图形面积的乘积。其中，ed,i是第i个供电周期的供电能量，td,i表示第i个供电周期的时间，u(t)表示供电电压，i(t)表示放电电流，um(t)指积分电容上的累积电荷电压，cm为积分电容的电容值，sd,i表示第i个供电周期的李萨如图形面积；
整个运行时间的dbd系统消耗总能量为：
式中，tm,j为pdm高压电源一个调功周期的时间；ton,j为一个调功周期内的dbd放电时间，由多个供电周期td,j组成；占空比
dbd放电区域内活性物种的量子产率能效比eer为：
其中i为光纤光谱仪采集的dbd放电区域中活性物种发射光谱的相对强度；em,j为一个调功周期的供电能量，即
当李萨如图形的各顶点坐标已定，其图形面积的计算为现有技术。
本发明所述可编程电源对pdm高压电源的控制包括输入电压控制、供电的功率密度控制和放电强度控制，分别为现有技术。
如图2所示，本实施例中，可编程电源模块采用可编程ac/dc电源模块。可编程ac/dc电源用于设定pdm高压电源的输入电压及整个系统的功率。
电气参数检测模块采用数字示波器，数字示波器包括至少3个检测通道，3个检测通道分别设置高压探头、电压探头和电流互感器，并分别用于检测dbd阵列的高压电源输入端电压、积分电容上的电压和放电回路的电流。电流互感器采用罗氏线圈。罗氏线圈以感应方式接在放电回路的地线上，获得放电电流，为现有技术。
实施例2
参考图7，基于实施例1所示pdm工作模式的高压电源激励dbd反应器供电能量的检测系统的供电能量检测方法，包括：
s1，系统初始化：将dbd的激励电源调节为pdm供电模式，设定pdm供电模式下的放电条件，使得dbd被激励放电；所述放电条件包括pdm高压电源供电电压、调功周期长度、供电频率和功率密度；
s2，利用电压电流检测模块获取dbd阵列的放电参数，包括供电电压、放电电流，和积分电容上的积分电压；
s3，上位机从电压电流检测模块获取已采集到的放电参数，基于获取到的放电参数进行dbd供电能量、等效参数和有效放电时间的计算和输出。
步骤s3包括：
s31，计算dbd的供电能量：
定义tm,j为pdm高压电源的一个调功周期，ton,j为dbd在pdm高压电源一个调功周期内的放电时间，由多个dbd供电周期组成；则dbd的单个供电周期能量ed,i表示为积分电容和李萨如图形面积的乘积，为：
式中，td,i表示单个供电周期的时间，u(t)是供电电压，i(t)是放电电流，d表示放电，i表示放电次数；
一个调功周期tm,j中供电时间的占空比为：
一个调功周期内的供电周期个数为：
一个调功周期的供电能量em,j即ndm,j个单个供电能量之和：
定义ed,a是调功周期内单个供电周期的平均供电能量，则：
当dbd系统工作在连续状态时，运行时间内总供电周期个数non,t可以表示为：
式中，tt是放电系统的运行总时间；
则系统运行时间内的dbd总供电能量et表示为：
当tt＝tm,j时，
根据式(1)、(6)、(7)和(8)、可得到et：
根据dbd阵列在pdm高压电源激励下产生的微放电数据，绘制激励dbd的供电能量的电压电荷李萨如图形，则单个供电周期的能量与相应李萨如图形的面积成正比；
由于：
其中cm是回路中外接积分电容的容值，um(t)为积分电容上的累积电荷电压；
则将(10)式带入(1)式可得单个供电周期的供电能量ed,i，表示为积分电容和相应李萨如图形面积sd,i的乘积，即：
将(11)代入(9)，利用李萨如图形面积来表示供电能量，可以得到运行时间内dbd总供电能量et为：
s32，计算dbd反应器的等效电容：
定义单个供电周期中产生微放电的时段有ab时段和cd时段，分别位于单个供电周期的负半周期和正半周期中，对应的李萨如图形的四个顶点坐标分别为a(ux1,uy1)、b(ux2,uy2)、c(ux3,uy3)、d(ux4,uy4)，各顶点的横坐标和纵坐标分别为相应时间点上的积分电容电压和供电电压；
则单个供电周期td,i对应的dbd反应器的总等效电容ci、介质等效电容cd,i和放电间隙等效电容cg,i分别为：
在一个调功周期tm.j中，形成如图6(b)所示的堆叠的李萨如图，根据李萨如图形，dbd阵列的等效电容看做该调功周期内所有供电周期对应的等效电容的平均值。
式(16)-(18)分别是根据式(13)-(15)得到的一个调功周期内对应的等效电容平均值。
s33，计算量子产率能效比eer：
其中i为光纤光谱仪采集的dbd放电区域中活性物种发射光谱的相对强度。
s31在进行dbd的供电能量的计算时，首先获取检测到的供电电压数据和积分电容电压数据，然后根据获取到的数据，进行单个供电周期的李萨如图形重构，图形重构时，包括：
以电容积分电压为横坐标数据，以“供电电压”为纵坐标数据；得到所有供电周期堆叠的李萨如图形；
进行单个周期李萨如图形的分离，得到单个供电周期李萨如图形，分离规则为：选择供电电压上升到零的两个相邻点之间的时间段作为一个供电周期；
对分离得到的单个李萨如图形分别计算得到李萨如图形的面积，进而根据公式(11)得到相应供电周期的供电能量ed,i。
更进一步的，s3还包括放电参数处理步骤：
根据获取到的供电电压参数，得到供电电压峰-峰值；
根据获取到的回路放电电流，得到dbd的平均微放电强度和有效放电时间，包括：
对获取的回路放电电流进行波形重构；
对重构形成的放电电流波形进行平滑处理；优选为采用savitzky-golay平滑算法；
将回路放电电流数据与平滑处理后的相应电流数据相减，得到放电过程中的微放电电流的波形数据；
根据微放电的电流波形的产生，得到每个供电周期的放电起始时间和放电截止时间，从而计算单个供电周期的有效放电时间；
将所有供电周期的有效放电时间累加，得到有效放电总时间；
根据每个供电周期的微放电脉冲峰值，计算放电过程中的平均微放电强度。
实施例3
本实施例综合本发明的所有研究内容进行描述，检测系统和检测方法主要解决以下各个方面问题：
1)dbd放电检测系统结构
参考图1，测试系统是采用功率密度调制(powerdensitymodulation:pdm)电源驱动体积dbd反应器产生等离子体，同时通过高压探头、电压探头、电流互感器、示波器等测试仪器获得放电参数的原始数据，设计一种pdm高压电源激励dbd反应器供电能量的计算方法，在放电过程中，将测试仪器获取到的信号数据通过特定的数据采集口(如usb接口、rs232接口等)传输到上位机，在上位机中利用文中设计的计算方法进行处理，对放电状态进行显示、监测和评估。
dbd阵列结构如图3所示，单个电极由两块石英板中间夹着铜箔电极组成。正负电极交错排列。d为正负电极之间的放电间隙，约为2到6mm；l1为石英板的长度，约170到180mm；l2为石英板的宽度，约为105到115mm；d为石英板的厚度，约为1到2mm。
如图2是pdm高压电源激励的dbd测试系统图。由dbd阵列、可编程ac/dc电源、pdm高压激励电源、高压探头、电流互感器、电压探头、光纤光谱仪、数字示波器、以及pc机(上位机)组成。
可编程ac/dc电源用于设定pdm高压电源的输入电压及整个系统的功率，pdm高压电源作为dbd反应器的激励电源，其供电的功率密度、电压、放电强度都可进行调节，因此反应器的供电条件随着pdm高压电源的工作模式而改变。高压探头连接在pdm高压电源的高压端(即dbd阵列的高压供电端)，以测量dbd反应器的供电电压；电流互感器接在dbd供电回路的地线，以测量流过整个回路的放电电流变化情况；在dbd反应器与地线之间连接一个电容作为微放电积分电容，通过一个电压探头获取积分电容两端的电压信号，将检测到的三路信号同时接入示波器的三个通道。使用交互软件将示波器与pc机相连，可直接在pc机上观察并获取到放电参数的原始数据。另外，光纤光谱仪的输入端朝向dbd放电区域，输出端与pc机相连用来检测放电空间内活性物种的相对浓度。pc机上的数据处理模块直接将原始数据导入，使用参数计算程序进行处理，对放电参数进行实时检测及输出。
2)dbd反应器供电能量检测
如图4为pdm高压电源的典型供电波形，图中的tm,j为一个功率密度调制周期，简称调功周期，对于一个固定频率的pdm高压电源，在调节电源的输出功率和电压时，tm,j是固定的。定义图中的ton,j为一个调功周期内的放电时间，由多个单个供电周期组成，在tm,j固定时，通过改变ton,j改变供电持续时间，即可改变供电能量。因此，对于pdm高压电源，供电电压、单个供电周期的频率、一个供电周期中的放电次数都可以进行调节。由于pdm高压电源具有放电不连续，放电频率可调，供电周期个数可变、单个供电峰-峰值不固定等特点，需要一种针对于pdm高压电源激励下dbd反应器供电能量的计算方法。本发明根据这一状况，设计了一种可对单个供电周期能量进行累加计算，得到总能量的方法，此方法既能计算连续工作状态下的供电能量，又能计算非连续、供电峰-峰值不固定工作状态下的供电能量。
结合图4的典型供电波形图，本发明通过改变ton,j来调整功率密度，因此采用以下计算方法获得供电能量：
公式(1)是单个供电周期能量的计算公式，表示为积分电容和李萨如图形面积的乘积。其中，ed,i是单个供电周期的供电能量，td,i表示单个供电周期的时间，u(t)是供电电压，i(t)是供电电流，此外，d和i分别表示了放电和放电次数。
一个调功周期tm,j中供电时间的占空比是供电持续时间与调功周期的比值，通过下式表示：
在一个调功周期内的供电周期个数可以通过以下方法计算：
其中ndm,j是一个调功周期内的供电周期个数。
一个调功周期内单个供电能量之和即为该调功周期的供电能量，即：
em,j是一个调功周期消耗的总能量，m表示调功周期，j表示调功周期的次数。
一个调功周期的供电能量还可通过单个供电周期的平均能量得到：
em,j＝ndm,jed,a(5)
其中，ed,a是调功周期内单个供电周期的平均供电能量，可以表示如下：
当dbd系统工作在连续状态时，运行时间内总供电周期个数可以表示为：
non,t是运行时间内总供电周期个数，tt是放电系统的运行总时间。
系统运行时间内的总供电能量et表示如下：
根据(1)，(6)，(7)和(8)式可得到et：
图5显示了激励dbd单元的典型电压电流波形，是供电周期的波形(为图4的供电部分的放大图)。在供电电流波形中可以观察到微放电的产生，对应图6(a)中标记的ab和cd时段，并且在标记为bc和da的时段中具有位移电流。图6(a)是激励dbd单元的供电能量的电压电荷李萨如图形。一个供电周期的能量与李萨如图形的面积成正比。通过李萨如图形可以得到部分放电特性。
由于：
其中cm是回路中外接积分电容的容值。
其中供电能量表示如下：
将(10)式带入(1)式可得：
公式(11)是一个供电周期的能量，表示为积分电容和李萨如图形面积的乘积。其中，ed,i是一个供电周期的供电能量，td,i表示一个供电周期的时间，u(t)是供电电压，i(t)是供电电流，um(t)指累积电荷电压，cm为辅助电容的电容值，sd,i是相应供电周期李萨如图形面积，
将(11)代入(9)，利用李萨如图形面积来表示供电能量，可以得到运行时间内dbd系统消耗的总能量为：
如图4所示，单个供电周期中产生微放电的时段有ab时段和cd时段，分别位于单个供电周期的负半周期和正半周期中，结合图5，对应的李萨如图形的四个顶点坐标分别为a(ux1,uy1)、b(ux2,uy2)、c(ux3,uy3)、d(ux4,uy4)，各顶点的横坐标和纵坐标分别为相应时间点上的积分电容电压ux和供电电压uy；
则单次供电周期td,i对应的dbd反应器的总等效电容ci、介质等效电容cd,i和放电间隙等效电容cg,i分别为：
在一个调功周期tm.j中，形成如图6(b)所示的堆叠的李萨如图，根据李萨如图形，dbd阵列的等效电容看做该调功周期内所有供电周期对应的等效电容的平均值。
式(16)-(18)分别是根据式(13)-(15)得到的一个调功周期内对应的等效电容平均值。
其中c，cd和cg分别是总等效电容，介质等效电容和放电等效间隙电容。这些等效电容值随着加载在dbd反应器上的供电电压和能量改变。
3)计算相对光量子产率能效比
本发明根据光谱仪测试获得的相对光谱强度对dbd装置的放电效果进行评估。用放电区域内羟基自由基(·oh)的相对光谱强度除以一个调功周期消耗能量表示相对光量子产率能效比(eer)。因此，来自放电区域的活性物种的相对光量子产率能效比可以通过下式求出：
其中i是放电区域中活性物种的发射光谱相对强度，由图2的光纤光谱仪采集得到。等离子体区域中活性物种的相对浓度与其发射光谱的相对强度呈正相关。因此，可以通过使用相对光量子产率eer来评价活性物种的产率。
4)放电参数检测与处理流程
如图7为检测步骤主流程图。首先，dbd的激励电源调节到pdm供电模式，在pdm模式下设定放电条件，包括pdm高压电源供电电压、调功周期的频率，供电频率，功率密度。使用第3节介绍的检测仪器进行放电参数采集，包括高压探头接在dbd反应器的高压供电端获得供电电压；电压探头接在积分电容两端获得积分电压；电流互感器接在回路的地线上获得放电电流。将这三组信号接到示波器进行原始数据采集。除此之外，使用光谱仪测量放电区域内活性物种的相对光谱强度。将示波器和光谱仪都与pc机相连，使用交互软件将示波器与pc机相连，可直接在pc机上观察并获取到放电参数的原始数据，在pc机上保存采集到的放电参数的原始数据。之后，将pc机采集到的原始数据自动导入到数据处理模块。经过一系列的数据处理程序得到供电能量、等效参数、有效放电时间等放电参数，并将结果输出显示。
如图8，参数处理程序是将pc机读取并保存的数据进行计算处理，通过供电电压处理子程序、放电电流处理子程序、供电能量计算子程序、等效参数计算子程序、量子产率能效比子程序，分别得到供电电压峰-峰值、有效放电时间、微放电平均强度、系统运行时间内总供电能量、单个供电周期平均能量、反应器等效电容、量子产率能效比。
第一个进行的是供电电压处理子程序，在这一部分，只需自动读取pc机上存储的“供电电压”采集通道的数据，之后得到供电电压的峰-峰值进行显示输出。
图9是放电电流处理子程序，该部分得到装置运行过程中的平均微放电强度和有效放电总时间。自动读取pc机上存储的“放电电流”采集通道的数据。然后进行对放电电流波形的重构。将得到的放电电流波形进行平滑处理，进行平滑处理时，本文采用的方法是savitzky-golay平滑算法，这是一种在时域内基于最小二乘法原理拟合的滤波方法，在对放电电流进行平滑的同时可以确保信号的形状和宽度不变。在matlab中可直接调用sgolayfilt函数进行处理：
data＝xlsread('放电电流.xlsx')；％读入数据
b＝data(:,1)；
y＝sgolayfilt(b,3,7)；％平滑处理
plot(y,'r')；
holdon
将平滑处理后的电流波形数据进行存储。之后，原始的放电电流数据与平滑处理过的电流数据相减，得到放电过程中的微放电电流的波形数据。对微放电电流数据进行检测，得到每一供电周期的放电起始时间和放电截止时间，计算出单个供电周期的有效放电时间，有效放电时间即电流波形上发生微放电的时间。将所有供电周期的有效放电时间累加，得到有效放电总时间。另外，对得到的微放电电流数据进行检测，得到每一供电周期的微放电脉冲峰值，计算放电过程中的平均微放电强度。最后对放电过程中有效放电总时间和平均微放电强度进行输出显示。
图10是供电能量计算子程序，该部分得到装置运行过程中的单个供电周期的能量平均值和总供电能量。自动读取pc机上存储的“供电电压”与“积分电压”采集通道的数据。根据得到的数据，进行单个周期的李萨如图形重构。由于李萨如图形是两个正交矢量在做周期性振荡时合成的曲线，在本方案中，两个矢量的振动频率相同，能够合成封闭的图形。但是在振动过程中，矢量的模是不固定，因此每个周期合成的图形大小会有差异。在进行供电能量计算时，需要对每个周期李萨如图形的面积进行计算，因此要进行图形的重构并对每个周期进行分离，确立图形分离的规则。重构过程如下：
以“积分电容电压”为横坐标数据，“供电电压”为纵坐标数据，进行图形重构。此时得到的是所有供电周期堆叠的李萨如图形。由于单个供电周期对应的是单个李萨如图形，因此需进行单个周期图形分离。确立单个供电周期李萨如图形重构的数据选取范围规则是：在绘制李萨如图时，选择施加的电压上升到零的两个相邻点的时间作为一个周期。
在对单个供电周期的李萨如图形进行分离重构后，对重构的李萨如图形进行供电电压与积分电压的积分计算，得到李萨如图形的面积sd,i。根据李萨如的面积，结合公式(11)得到单个供电周期的能量ed,i。根据公式(3)和(7)得到一个调功周期内的供电周期个数ndm,j，之后根据公式(4)得到该调功周期的供电能量em,j。单个供电周期的平均能量通过公式(6)得到。之后根据运行时间内总供电周期个数non,t和公式(12)对每个供电周期的能量进行累加，得到系统运行时间内的总供电能量。对ed,a和et进行输出显示。
图11是等效参数计算子程序，该部分得到反应器的等效电容参数。自动读取pc机上存储的“供电电压”与“积分电压”采集通道的数据。根据得到的数据，进行单个周期的李萨如图形重构。根据公式(13)-(15)得到装置单个供电周期的等效电容数据c，cd和cg。然后对所有周期的电容值进行平均，得到等效电容平均值。对计算结果进行输出显示。
图12是相对光量子产率能效比(eer)计算子程序，自动读取pc机上存储的“相对光谱强度”数据，读取“供电能量处理子程序”中经过处理存储的“调功周期供电能量em,j”数据，通过公式(16)计算得到相对光量子产率能效比，并对结果进行显示输出。
5)相对光量子产率能效比评估
通过前面的步骤得到放电参数后，本发明与牛顿爬山算法相结合，设计了最佳放电效果评估方法。根据光量子产率的变化规律，得到在最佳放电效果时对应的放电参数。根据牛顿爬山算法找出最佳eer(能效比)对应的放电条件，确定相应的参数范围。
如图13所示是牛顿爬山法的xx-eer图，横轴表示一个调功周期的供电能量em,j，纵轴表示相对光量子产率eer。牛顿爬山法又称扰动观察法，本文中是通过不断调节放电反应系统的放电条件来比较调整前后dbd产生eer的变化情况，再根据变化情况来调整放电条件，包括可编程ac/dc电源的输出电压、pdm高压电源的供电电压、供电能量的参数,使放电反应器工作在最佳能效比附近。根据上图可以分析牛顿爬山法具体工作情况如下：
(1)在a点加一个扰动变量，如改变供电能量，使反应器的eer达到b点；
(2)检测到之前提高供电能量使得反应器的相对光量子产率能效比增加，继续原来的方向增加扰动变量，使反应器工作在c点；
(3)继续原来的方向加扰动变量，使反应器工作在m点；
(4)继续原来的方向加扰动变量，使反应器工作在d点；
(5)此时检测到之前的扰动变量使得反应器的相对光量子产率能效比减小，改变原来的方向加扰动变量，使反应器的eer重新达到m点；
(6)继续原来的方向加扰动变量，使反应器工作在c点；
(7)最后，反应器在c点，m点，d点三个工作点间波动。
算法中的扰动变量可取气体流速，供电电压，一个调功周期的供电能量。使用这种方法可以确定最佳的相对光量子产率能效比，而且可以得到所对应的放电条件，因此可以确定最佳放电参数调节范围。
在牛顿爬山算法中，根据扰动变量的步长所确定的m点不一定是eer的最高点，在确定了最佳参数调节范围c到d后，用优选法对扰动变量的步长进行重新设定，根据下列步骤找出最高点：
(1)在(c,m)区间内取中点p1，(m,d)区间内取中点p2；
(2)当p1对应的函数值大于p2对应的函数值时，相对光量子产率能效比的极大值在(p1,m)区间内，(m,p2)的区间舍去；
(3)反之，极大值在(m,p2)区间内，舍去(p1,m)区间；
(4)当p1和p2对应的函数值相等时，极大值在(p1,p2)范围内，舍去(c,p1)和(p2,d)的区间；
(5)在剩余的区间内重新取中点，找到p3，p4，以步骤(1)到(4)的方式继续进行迭代计算；直到剩余区间范围小于设定值
当自变量为气体流速时，
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。