本发明属于污染物监测，尤其涉及一种无人机双半球采样装置与大气污染源排放定量方法。
背景技术：
1、获取工业企业排气筒大气污染物排放情况是大气环境监管的重要依据，采用无人机进行污染物浓度采样监测是生态环境领域常用方法。已有技术大多采用无人机进行排气筒周边烟羽采样分析，直接以烟羽浓度表征排放情况，或通过高斯扩散模型等模拟手段反演污染物浓度或计算排放量，受环境空气小微尺度平流、湍流及排气筒烟羽排放的复杂三维气象场，以及实测浓度数据代表性和空气质量模型不确定性影响，上述方法难以准确表征排气筒烟羽排放量。
技术实现思路
1、发明目的：本发明所要解决的技术问题是针对现有排气筒大气污染物排放定量技术的不足，提供一种无人机双半球采样装置，包括搭载在无人机上的在线监测装置、双旋臂采样装置、板载计算机数据传输控制系统；
2、所述双旋臂采样装置用于采集环境样品，实现在一个采样位置同时完成双球面点位大气环境样品采集；
3、所述板载计算机数据传输控制系统包括基于单片机系统设计的控制板，所述板载控制板中集成有双半球飞行路径建模系统，所述双半球飞行路径建模系统用于产生飞行轨迹与监测点位数据，所述板载计算机数据传输控制系统读取双半球飞行路径建模系统产生的飞行轨迹与监测点位数据，并发送控制飞行所需参数，同时与在线监测装置、双旋臂采样装置联动，将在线监测数据和采样数据传输回地面站。
4、所述在线监测装置包括高分辨率大气环境六参数仪器、辅助无人机定位的基于全球卫星导航系统的实时载波相位差分定位rtk模块；
5、所述高分辨率大气环境六参数仪器用于同时测量no2、co、o3、so2浓度气体浓度和pm2.5和pm10颗粒物质量浓度，no2、co、o3、so2浓度气体浓度的测量分辨率为1ppb，pm2.5和pm10颗粒物质量浓度的测量分辨率为0.01ug/m3；
6、所述实时载波相位差分定位rtk模块用于实时获取无人机的gps坐标和飞行高度（rtk模块采用大疆d-rtk，包括天空端和地面端），测量分辨率为1厘米。
7、所述双旋臂采样装置包括采样管旋转控制装置和恒流恒温的样品采集控制系统；
8、所述采样管旋转控制装置包括第一采样管和第二采样管、金属连接件、行星齿轮减速机构；
9、所述采样管包括石英进样头、石英过滤棉、塑料连接部件、聚四氟采样管、碳纤维管；所述塑料连接部件用于连接聚四氟采样管和石英进样头；所述石英进样头入口处填充石英过滤棉；
10、所述行星齿轮减速机构包括两组行星减速齿轮组、带旋转编码器的双出轴减速驱动电机、固定座、固定块；所述固定块与采样管固定连接；每组行星减速齿轮组包括一个太阳齿轮、三个行星齿轮和一个外圈齿轮；所述两组行星减速齿轮组、带旋转编码器的双出轴减速驱动电机、固定块和固定座相连接；
11、所述固定座与无人机的碳纤维板连接，通过旋转编码器对电机实时测量，保证了采样管的精确旋转，对预定空间位置的环境样品实现精确采集；
12、所述行星减速齿轮组的动力来源为双出轴减速驱动电机，双出轴减速驱动电机与行星减速齿轮组中心的太阳齿轮同轴固定安装；分布在行星减速齿轮组的中间的三个行星齿轮等间距分布，通过齿轮啮合方式，实现将太阳齿轮的动力输入，转换到外圈齿轮输出；
13、所述旋转编码器提供高精度位置信号，从而实现多级减速下的采样管转动角度的控制；
14、所述双旋臂采样装置通过行星齿轮减速机构实现采样管精确旋转，所述行星齿轮减速机构通过太阳齿轮输入，实现前向0~90°正转和后向0~负90°的旋转；所述双旋臂采样分析装置通过旋转采样管实现前后两个不同半径球面的取样。
15、所述双出轴减速驱动电机采样直流为12v，减速箱减速比为12：1；所述行星齿轮的内齿圈齿数为30，外齿圈齿数为60，行星齿轮传动比为3；所述双旋臂采样装置中双出轴减速驱动电机到采样管旋转的减速比为36：1。
16、所述恒流恒温的样品采集控制系统包括恒温控制系统和恒流控制系统；
17、所述恒温控制系统实现两个温控通道，分别对每根采样管进行温度控制；每个温控通道包括独立的温度传感器、聚酰亚胺薄膜加热套和保温材料；
18、所述恒流控制系统实现恒流样品采集，包括样品电磁阀、质量流量控制器、直流采样泵；
19、所述恒流恒温的样品采集控制系统用于采集排气筒周边的环境气体样本，降低采样管壁低温对样品气体的吸附和冷凝，对采集流量进行恒流控制；
20、所述采样管采用碳纤维管作为外套管，外径10mm，内径8mm，采样管长度可根据监测要求设计，所述外套管内嵌1/4英寸聚四氟采样管；
21、所述聚四氟采样管穿过金属连接件，通过样品电磁阀，最后再与在线监测装置连接；
22、所述采样管的外壁缠绕聚酰亚胺薄膜加热套和酚醛树脂保温套；所述聚酰亚胺薄膜加热套采用直流12vdc加热，电流小于2a；可实现飞行采样期间采样管恒温50摄氏度。
23、无人机飞控实时返回目标点位的rtk坐标数据给板载计算机数据传输控制系统进行核对坐标，到达坐标位置后通知在线监测装置进行采样分析。
24、设定监测排气筒污染物排放时，监测点位周边大气达到平衡状态，则污染物的输送主要由边界层内湍流造成。无人机在两个半球面中间位置飞行，通过监测排气筒外围同一球心的两个半球面上大气污染物浓度阵列，计算每个网格内的大气污染物排放通量，得到固定污染源大气污染物排放强度；
25、所述双半球飞行路径建模系统基于无人机飞行路径建模算法，生成双半球监测所需的点位与路径，基于监测球面半径、采样管长度、水平与垂直方向监测点位间隔，以飞行时间与距离最短为约束条件，通过双半球飞行路径建模系统进行航迹规划，生成航迹点和停留时间，完成无人机飞行与采样监测所有参数设置，所述参数包括飞行参数和采样参数，所述飞行参数包括空间点位坐标序列、飞行高度、飞行速度，航向角；所述采样参数包括采样温度、采样流量、采样时间、两个采样管角度；每个点位均对应飞行参数和采样参数。
26、所述恒流恒温的样品采集控制系统通过板载质量流量控制器对采样流量进行稳定控制和精确测量，对采样管外壁的聚酰亚胺薄膜加热套进行加热，由靠近保温层的pt100（铂电阻传感器，温度测量精度1摄氏度）传感器反馈温度信号，实现恒温控制；
27、所述恒流恒温的样品采集控制系统在采样开始前手动设置参数，或者在飞行中进行远程控制。
28、所述装置将获取的实时流量、实时温度、采样时间、rtk坐标通过数据发送模块实时传输至地面站，方便地面人员监控。
29、所述板载计算机数据传输控制系统接收基于全球卫星导航系统的实时载波相位差分定位模块提供的高精度gps信号，并执行导航计算以实现准确定位和航迹规划。
30、本发明还提供了一种大气污染源排放定量方法，包括如下步骤：
31、步骤1，所述双半球飞行路径建模系统通过三维建模计算排气筒周边烟羽监测的双半球坐标；
32、通过双半球监测路径建模系统，得到飞行任务规划文件，起飞前将飞行任务规划文件上传至板载计算机数据传输控制系统，所述飞行任务规划文件包括飞行参数和序列的采样参数；
33、步骤2，所述实时载波相位差分定位rtk模块通过rtk高精度双定位方式，确定飞行坐标与飞行轨迹，读取离线规划好的监测数据，参照无人机飞控格式要求，将飞行轨迹数据实时传输至板载计算机数据传输控制系统；无人机按照规定飞行路径与监测点位完成双球面大气污染物浓度监测；
34、步骤3，根据球面参数计算无人机飞行球面面积，基于无人机双球面大气污染物监测结果，计算每个网格的大气污染物排放通量并求和，结合飞行时间，计算得到排气筒整体大气污染物排放量。
35、步骤1包括：
36、步骤1-1，确定无人机飞行与采样监测点位：
37、以监测的排气筒口中心为球心，确定无人机飞行球面半径radius，采样监测的双半球内圈半径和外圈半径，无人机在内外半球中间位置飞行，radius计算公式如公式（1）所示，第一采样管和第二采样管的采样头之间距离d即为两个半球的半径差：
38、（1）
39、确定垂直方向监测层数m，内外半球的层数需一致且均匀覆盖半球面；每层以球面横截面圆心为中心，确定每个垂直监测层水平截面监测间隔度数（以横截面圆心为中心的放射状间隔），水平监测点位阵列需要形成闭合圆形，每层采样点数量为n；最终形成封闭的双半球面监测网格；
40、考虑到通量积分的有效性和达标性，每个球面的垂直层数量不小于5，每层的水平截面监测点位数量不少于10；
41、步骤1-2，建立采样点位三维坐标并转换为无人机飞行的经纬度坐标。
42、使用已知的排气筒口中心的gps坐标（wgs84）和高度信息作为球心，将球心的经纬度转换成弧度，根据给定的内外球半径、和监测层数layer（包括水平和垂直层的监测层数），在球面上生成采样点，根据球心经纬度和高度计算得到每个采样点的球面坐标，再将球面坐标转换成wgs84坐标，最后输出所有坐标的经纬度和高度；
43、步骤1-3，将半球面采样点经纬度与高度参数输入板载计算机数据传输控制系统，通过双半球飞行路径建模系统规划航迹，生成航迹点和停留时间；完成无人机飞行与采样监测所有参数设置，生成无人机飞行与采样配置文件，上传至无人机的板载计算机数据传输控制系统；
44、步骤2包括：
45、步骤2-1，启动高分辨率大气环境六参数仪器、双旋臂采样装置、实时载波相位差分定位rtk模块和板载计算机数据传输控制系统；
46、步骤2-2，使用实时载波相位差分定位rtk模块提供的高精度gps信号执行球形网格点导航飞行，无人机飞至每个采样点位调整转向航向角后，由板载计算机数据传输控制系统发送控制命令给双旋臂采样装置，完成样品采集与分析监测，每个悬停点位停留时间为t1（一般为10s），实时获得污染物排放浓度，遍历步骤1中所有采样点位，完成最终任务；
47、步骤2-3，通过板载计算机数据传输控制系统，将高分辨率大气环境六参数仪器、恒流恒温的样品采集控制系统采集的数据实时传输至地面数据接收系统，同时通过地面控制站对无人机的电池电压、续航时间、飞行状态（飞行速度、航向角、高度、信号强度等等）进行实时监控；
48、步骤3包括：
49、步骤3-1，为计算排气筒污染物排放通量，将原经纬度坐标系转换为球面坐标系，连线原点和飞行点位，构建新坐标系；
50、步骤3-2，进行无人机监测并记录飞行所使用的时间t，单位为s；
51、步骤3-3，仅考虑湍流通量，则排气筒污染物排放强度为双采样管监测的差值：首先求得无人机飞行坐标系内第i,j个网格对应面积：
52、（2）
53、其中，单位为m2；、、、、、、和为网格的四角坐标；
54、步骤3-4，根据每个网格监测所得污染物浓度，µg/m3，对半圆进行求和，并求两个球面上的监测浓度差，µg/m：
55、（3）
56、其中，n和m分别为半球水平方向网格数和垂直方向网格数，代表外半球第i，j网格的污染物浓度，µg/m3；代表内半球第i，j网格的污染物浓度，µg/m3，则最终排气筒排放强度i (单位为：µg s-1)为：
57、。
58、有益效果：本发明的无人机双旋臂采样监测装置，能同时对双点位进行监测采样，配合无人机监测路径规划系统，能快速获得排气筒周边烟羽规定点位阵列的大气污染物（no2、co、o3、so2、pm2.5和pm10）浓度数据，大大提升了无人机大气环境采样监测的时效性。基于双半球污染物实测浓度的排放通量计算方法，具有较好的流场闭合性和稳定性，规避了环境空气小微尺度平流、湍流及排气筒烟羽排放的复杂三维气象场影响以及空气质量模拟方法的不确定性影响。本发明的采样装置与排放定量方法能有效提升排气筒烟羽排放定量的有效性、准确性和稳定性。
1.一种无人机双半球采样装置，其特征在于，包括搭载在无人机上的在线监测装置（4）、双旋臂采样装置、板载计算机数据传输控制系统；
2.根据权利要求1所述的一种无人机双半球采样装置，其特征在于，所述在线监测装置包括高分辨率大气环境六参数仪器、辅助无人机定位的基于全球卫星导航系统的实时载波相位差分定位rtk模块；
3.根据权利要求2所述的一种无人机双半球采样装置，其特征在于，所述双旋臂采样装置包括采样管旋转控制装置和恒流恒温的样品采集控制系统；
4.根据权利要求3所述的一种无人机双半球采样装置，其特征在于，所述恒流恒温的样品采集控制系统包括恒温控制系统和恒流控制系统；
5.根据权利要求4所述的装置，其特征在于，所述采样管采用碳纤维管（21）作为外套管，所述外套管内嵌聚四氟采样管（20）；
6.根据权利要求5所述的装置，其特征在于，设定监测排气筒污染物排放时，监测点位周边大气达到平衡状态，无人机在两个半球面中间位置飞行，通过监测排气筒外围同一球心的两个半球面上大气污染物浓度阵列，计算每个网格内的大气污染物排放通量，得到固定污染源大气污染物排放强度；
7.根据权利要求6所述的装置，其特征在于，所述恒流恒温的样品采集控制系统通过板载质量流量控制器对采样流量进行稳定控制和精确测量，对采样管外壁的聚酰亚胺薄膜加热套进行加热，由传感器反馈温度信号，实现恒温控制；
8.根据权利要求7所述的装置，其特征在于，所述板载计算机数据传输控制系统接收基于全球卫星导航系统的实时载波相位差分定位模块提供的高精度gps信号，并执行导航计算以实现准确定位和航迹规划。
9.一种大气污染源排放定量方法，采用如权利要求2~8任一项所述的一种无人机双半球采样装置实现，其特征在于，包括如下步骤：
10.根据权利要求9所述的方法，其特征在于，步骤1包括：
