一种基于CO2-DIAL模拟测量的点源CO2排放的新型估算方法

一种基于co2‑
dial模拟测量的点源co2排放的新型估算方法
技术领域
1.本发明涉及大气遥感领域，尤其涉及一种基于co2‑
dial模拟测量的点源co2排放的新型估算方法。


背景技术：

2.人类活动排放的co2是大气中co2增加的主要原因。过去几个世纪大气中的co2急剧增加，造成了全球变暖，气候变化剧烈等各种影响。然而，大功率的发电厂的排放量超过人为排放量的30％，大气中co2含量的上升主要是由化石燃料燃烧产生的co2排放增加引起的。《联合国气候变化公约》(unfccc)以控制人为二氧化碳的排放。监测co2排放的传统方法取决于排放清单。尽管排放清单的准确性和可靠性很高，但它们的质量高度依赖于许多发展中国家所缺乏的完善的法律体系和监管体系。此外，快速的经济发展导致二氧化碳排放量的巨大变化，这很难通过排放量清单报告。通过卫星获取co2的浓度，例如通过卫星oco
‑
2的xco2(co2的柱平均干空气摩尔分数)，虽然可以获取大面积的co2浓度，但是限制了强点源co2浓度的估算。因此，有必要开发一种新颖的基于测量的强co2点源监测方法，以不依赖于排放清单的编制，以补充当前的碳排放验证系统。
3.随着dial(差分吸收激光雷达)的快速发展，地面co2‑
dial可以部署在任何可以到达的位置，以测量大气中co2的浓度场，这种设备依靠气溶胶反向散射的红外激光信号来获取至少3km内的大气co2浓度。
4.将co2浓度场转为co2排放量的反演方法的基本思想是，从点源排放的co2分子的空间分布服从与气象和地形因素有关的扩散模型。高斯扩散模型被广泛的应用于co2的排放量估计。因此我们以高斯扩散模型为反演核，以co2‑
dial模拟观测的距离分辨的大气co2浓度为关键输入，计算点源的co2排放。


技术实现要素：

5.本发明要解决的技术问题在于针对现有技术中卫星观测忽略强点源排放的缺陷，提供一种基于co2‑
dial模拟测量的点源co2排放的新型估算方法。
6.本发明解决其技术问题所采用的技术方案是：
7.本发明提供一种基于co2‑
dial模拟测量的点源co2排放的新型估算方法，该方法包括以下步骤：
8.步骤s1、获取co2‑
dial对烟囱排放co2的模拟观测数据，设定观测的初始条件；
9.步骤s2、建立坐标系，计算每个co2‑
dial测量co2浓度的测量点位置；
10.步骤s3、首先选取co2‑
dia的一个观测方向作为估算方向，将co2‑
dial中每一个测量点的co2浓度，作为估算点；利用线性化的高斯扩散模型建立线性化方程；
11.步骤s4、路径检验，在选取一个观测方向作为验证方向，该方向得到co2的浓度，作为验证点；
12.步骤s5、根据有效排放高度的不同选取，将验证点与估算点的co2浓度计算值进行
比较，误差最小时，得到烟囱的co2排放量。
13.进一步地，本发明的步骤s1中设定的初始条件包括：假设风的平均流场稳定，风速均匀，风向笔直；假设co2的排放服从正态分布，污染物浓度在y轴和z轴方向上均符合正态分布；污染物在运输和扩散过程中被大量保存；污染源统一，连续。
14.进一步地，本发明的步骤s2中的具体方法为：
15.计算测量点位置的公式为：
[0016][0017]
其中，每个测量点的位置为(x，y，z)；σy，σz的弥散系数由位置与顺风方向上的坐标原点的距离确定；为测量垂直于x轴的浓度，保证每个测量点的σy和σz相同；顺风方向被视为x轴，烟囱高度方向被视为z轴，y轴垂直于xoz平面，得到dial和x轴之间的角度；s是co2‑
dial位置与烟囱之间的距离，α是x轴与s方向之间的角度，β是xoy平面与l之间的角度，α和β由激光测角仪测量，l是被测点与co2‑
dial之间的距离；
[0018]
计算不同高度的风速，平均风速与不同高度之间的关系采用迪肯的幂定律描述：
[0019]
u＝u1(z/z1)
n
[0020]
其中，u以及u1表示在高度z以及已知高度z1处的平均风速；n是与大气稳定度有关的指数，n＝1/7，本步骤得到每个co2‑
dial测量co2浓度的测量点位置(x,y,z)。
[0021]
进一步地，本发明的步骤s3中的具体方法为：
[0022]
计算co2‑
dial中每一个测量点的co2浓度c的公式为：
[0023][0024]
其中，q是烟囱中co2的排放速率，us是烟囱高度处的平均风速，h是烟囱有效高度，h＝h+δh；h是烟囱高度，δh是羽流上升高度，σy是水平方向的弥散系数，σz是垂直方向的弥散系数；当co2分子向下到达表面时，它可能会反射到大气中，α表示反射指数，范围为0到1，当分子完全被表面吸收时将为0，如果co2分子完全被表面完全反射，则它将为1；
[0025]
为了计算方便，将上式简化为线性方程；则有：
[0026][0027][0028]
本步骤得到co2‑
dial中每一个测量点的co2浓度q，作为估算点。
[0029]
进一步地，本发明的步骤s4中的具体方法为：
[0030]
使用jacobi迭代算法求解以下方程，求解方程的雅各比矩阵b，然后对x1,x2,x3进行求解；
[0031][0032]
其中，b＝(d
‑
l)
‑1u,f＝(d
‑
l)
‑1b
[0033]
求解过程中引入以下定理：定理1：简单迭代法收敛的充分必要条件是迭代矩阵b的谱半径ρ(b)<1；定理2：若迭代矩阵b的某种范数||b||<1则确定的迭代法对任意初值x
(0)
均收敛于方程组x＝bx+f的唯一解x0；对任意初值x
(0)
均收敛于方程组x＝bx+f的唯一解x0，得到对应的x1,x2,x3，参数可用于计算验证方向的co2浓度。
[0034]
进一步地，本发明的步骤s5中的具体方法为：
[0035]
对于h，假设其范围为1m至100m，间隔为1m；n分别取值1,2,3
…
n，hn的假定值将对反演模型进行迭代，得到参数(x1n，x2n，x3n)，对应也将得到qn，σny，σnz；之后，检验路径将更改为另一个方向，l1—>l2，两个检验路径的x保持一致；将测量此路径中的co2分布，因为σy，σz仅与x的距离有关，并且风速q，σy，σz均与先前的测量路径相同；qn，σny，σnz将在新方向上迭代到已验证的测量点，将其选择为将参数(x1n，x2n，x3n)的co2浓度计算值与实际检测结果进行比较；当获得比较小的平均误差时，认为hn将是最终值；更进一步地，根据确定的x1,x2,x3值获得烟囱排放的真实二氧化碳强度。
[0036]
本发明产生的有益效果是：本发明的基于co2‑
dial模拟测量的点源co2排放的新型估算方法，基于co2‑
dial的模拟观测值和高斯扩散模型能够准确的反演强点源的co2排放，弥补卫星观测忽略的强点源的碳排放和补充排放清单。
附图说明
[0037]
下面将结合附图及实施例对本发明作进一步说明，附图中：
[0038]
图1是本发明实施例的总流程图；
[0039]
图2是本发明实施例的co2‑
idal测量示意图。
具体实施方式
[0040]
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
[0041]
本发明实施例的基于co2‑
dial模拟测量的点源co2排放的新型估算方法，包括以下步骤：
[0042]
步骤s1、如图1所示，首先假设风的平均流场稳定，风速均匀，风向笔直；污染物浓度在y轴和z轴方向上均符合正态分布；污染物在运输和扩散过程中被大量保存；污染源统一，连续。
[0043]
步骤s2、建立坐标系，如图1所示。在坐标系中，通过公式(1)计算每个测量点的位置(x，y，z)。σy，σz的弥散系数由位置与顺风方向上的坐标原点的距离确定。为测量垂直于x轴的浓度，可以保证每个测量点的σy和σz相同。顺风方向被视为x轴，烟囱高度方向被视为z轴，y轴垂直于xoz平面。我们可以得到dial和x轴之间的角度。s是co2‑
dial位置与烟囱之间的距离，α是x轴与s方向之间的角度，β是xoy平面与l之间的角度，α和β由激光测角仪测量，l是被测点与co2‑
dial之间的距离。
[0044][0045]
计算不同高度的风速，平均风速与不同高度之间的关系一般可以用迪肯(deacon)的幂定律描述：
[0046]
u＝u1(z/z1)
n
ꢀꢀꢀ
(2)
[0047]
式中u以及u1表示在高度z以及已知高度z1处的平均风速，m/s。n是与大气稳定度有关的指数，n＝1/7。本步骤得到每个co2‑
dial测量co2浓度的测量点位置(x,y,z)。
[0048]
步骤s3、通过公式(3)计算每一点的浓度c(mg/m3)。其中q(mg/s)是烟囱中co2的排放速率，us是烟囱高度处的平均风速(m/s)，h是烟囱有效高度(m)＝h+δh；h是烟囱高度(m)，δh是羽流上升高度(m)，σy是水平方向(m)的弥散系数，σz是垂直方向(m)的弥散系数。当co2分子向下到达表面时，它可能会反射到大气中，α表示反射指数，范围为0到1，当分子完全被表面吸收时将为0，如果co2分子完全被表面完全反射，则它将为1。本步骤得到co2‑
dial中每一个测量点的co2浓度q，作为估算点。
[0049]
[0050][0051][0052][0053]
ax
‑
l＝v
ꢀꢀꢀ
(7)
[0054]
步骤s4、使用jacobi迭代算法求解方程(5)，求解方程的雅各比矩阵b，然后对x1,x2,x3进行求解。
[0055]
x
(k+1)
＝bx
(k)
+f,(k＝0,1,2,
…
)
ꢀꢀꢀ
(8)
[0056]
其中，
[0057]
b＝(d
‑
l)
‑1u,f＝(d
‑
l)
‑1b
ꢀꢀꢀ
(9)
[0058]
下面为一些定理：
[0059]
定理1：简单迭代法收敛的充分必要条件是迭代矩阵b的谱半径ρ(b)<1。
[0060]
定理2：若迭代矩阵b的某种范数||b||<1则确定的迭代法对任意初值x
(0)
均收敛于方程组x＝bx+f的唯一解x0。对任意初值x
(0)
均收敛于方程组x＝bx+f的唯一解x0，得到对应的x1,x2,x3，参数可用于计算验证方向的co2浓度。
[0061]
步骤s5、对于h，假设其范围为1m至100m，间隔为1m，(1,2,3
…
n)，hn的假定值将对反演模型进行迭代，可以得到参数(x1n，x2n，x3n)，对应也将得到qn，σny，σnz。之后，检验路径将更改为方向(l1—>l2)，两个检验路径的x保持一致。将测量此路径中的co2分布，因为σy，σz仅与x的距离有关，并且风速q，σy，σz均与先前的测量路径相同。qn，σny，σnz将在新方向上迭代到已验证的测量点，将其选择为将参数(x1n，x2n，x3n)的co2浓度计算值与实际检测结果进行比较。当获得比较小的平均误差时，我们认为hn将是最终值,并根据hn对应的x1n，x2n，x3n求解强点源的排放强度。
[0062]
应当理解的是，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，而所有这些改进和变换都应属于本发明所附权利要求的保护范围。