当前位置: 首页 > 专利查询>宁夏大学专利>正文

一种同步反演人为制造技术

技术编号:39521135 阅读:10 留言:0更新日期:2023-11-25 19:00
本发明专利技术涉及二氧化碳反演技术领域,公开了一种同步反演人为

【技术实现步骤摘要】
一种同步反演人为CO2排放与自然CO2通量的同化方法


[0001]本专利技术涉及二氧化碳反演
,具体是一种同步反演人为
CO2
排放与自然
CO2
通量的同化方法


技术介绍

[0002]二氧化碳
(CO2)
是大气中重要的温室气体,准确掌握大气
CO2
的含量及变化,可为气候变化预测及环境决策提供支持,为满足气候研究的需求,需要对
CO2
的排放量进行反演

[0003]中国专利公开了一种二氧化碳通量反演方法

系统
(
授权公告号
CN111487216B)
,该专利技术适用于我国高分五号卫星搭载的
GMI
载荷的二氧化碳通量反演方法,在
CMAQ
模式中引入并耦合了陆地生态系统光合作用和呼吸作用的二氧化碳模块,使该模式成为污染气体和碳循环为一个大气的处理系统,从而得到的二氧化碳通量更准确

提高了反演精度;但是其存在一些问题:其只能得出总的
CO2
排放量,无法甄别出自然
CO2
通量和人为
CO2
排放量,不方便研究人分别以这两种数据进行有效地碳排放管理

因此,本领域技术人员提供了一种同步反演人为
CO2
排放与自然
CO2
通量的同化方法,以解决上述
技术介绍
中提出的问题

专利技术内容
[0004]本专利技术的目的在于提供一种同步反演人为
CO2
排放与自然
CO2
通量的同化方法,以解决上述
技术介绍
中提出的问题

[0005]为实现上述目的,本专利技术提供如下技术方案:
[0006]一种同步反演人为
CO2
排放与自然
CO2
通量的同化方法,包括以下步骤:
[0007]S1、
筛选数据:通过碳含量检测装置观测并选择合适的太阳天,选择合适的气溶胶光学厚度;
[0008]S2、
光谱信息分析:通过碳含量检测装置观测光谱各通道对应的大气
CO2、
温度

湿度和地表压力;
[0009]S3、
光谱数据处理:对反演中的实测或模拟的绝对辐亮度光谱进行重构,获得的辐亮度比值谱;
[0010]S4、
反演计算:采用估计反演算法,对碳含量检测装置观测光谱和模拟光谱开展迭代计算,获得
CO2
的浓度值;
[0011]S5、
总通量计算:通过
CO2
的浓度值计算出总
CO2
的排放总量;
[0012]S6、
分量计算:根据
CO2
的排放总量,通过分量计算法计算出人为
CO2
排放量与自然
CO2
排放量

[0013]作为本专利技术再进一步的方案:所述
S1
中的太阳天顶角小于
70
°
,所述气溶胶光学厚度小于
0.3
,气溶胶光学厚度来自于碳含量检测装置平台上的大气多角度偏振探测仪气溶胶产品

[0014]作为本专利技术再进一步的方案:所述
S2
中大气
CO2、
温度

湿度和地表压力的香农信息含量计算公式为:
[0015][0016]公式
(1)
中,
K
代表权重函数
,
不同参数的权重函数表示对应参数单位扰动导致仿真光谱的变化量;
S
a
代表先验协方差矩阵,该矩阵为对角阵,对角线上元素为各参数波动范围的平方;
S
e
是测量误差协方差矩阵,该矩阵同样为对角阵,对角线上元素为光谱各通道对应仪器噪声的平方;
[0017]利用大气
CO2、
温度

湿度和地表压力的香农信息含量,计算出综合香农信息含量,计算公式为:
[0018]H
tot

w1H
CO2
+w2H
T
+w3H
R
+w4H
P

ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0019]公式
(2)
中,
w
代表各参数信息含量的权重,计算用于
CO2
含量反演的信息时,
w1取3,
w2~
w4均取
‑1;计算用于背景扣除的信息时,
w1取7,
w2~
w4均取
1。
[0020]作为本专利技术再进一步的方案:所述
S3
中辐亮度比值谱计算公式为:
[0021][0022]公式
(3)
中,
i
代表各光谱通道,共
18
个;
Y
on
代表用于
CO2
含量反演的大气
CO2
信息谱,
Y
off
代表用于背景扣除的参考谱

[0023]作为本专利技术再进一步的方案:所述
S4
中估计反演算法的计算公式为:
[0024]y

F(x)+
ε

ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0025]公式
(4)
中,
y
代表碳含量检测装置观测光谱,
F
代表前向模型,状态向量
x
代表大气
CO2
廓线,
ε
代表碳含量检测装置观测误差;
[0026][0027]公式
(5)
中,
R
F(x)
代表由前向模型计算的模拟谱经公式
(3)
转换获得的模拟辐亮度比值谱,
R
y
代表
y
根据公式
(3)
转换获得的实测辐亮度比值谱;观测误差协方差
S
e
代表模型和测量的不确定性;先验误差协方差
S
a
代表先验估计
x
a
的不确定性

[0028]作为本专利技术再进一步的方案:采用
Levenberg

Marquardt
迭代法求解公式
(5)
,迭代形式为:
[0029][0030]公式
(6)
中,下标
i
代表第
i
次迭代;为雅克比矩阵,是
R
F(x)
的偏导数,
...

【技术保护点】

【技术特征摘要】
1.
一种同步反演人为
CO2
排放与自然
CO2
通量的同化方法,其特征在于,包括以下步骤:
S1、
筛选数据:通过碳含量检测装置观测并选择合适的太阳天,选择合适的气溶胶光学厚度;
S2、
光谱信息分析:通过碳含量检测装置观测光谱各通道对应的大气
CO2、
温度

湿度和地表压力;
S3、
光谱数据处理:对反演中的实测或模拟的绝对辐亮度光谱进行重构,获得的辐亮度比值谱;
S4、
反演计算:采用估计反演算法,对碳含量检测装置观测光谱和模拟光谱开展迭代计算,获得
CO2
的浓度值;
S5、
总通量计算:通过
CO2
的浓度值计算出总
CO2
的排放总量;
S6、
分量计算:根据
CO2
的排放总量,通过分量计算法计算出人为
CO2
排放量与自然
CO2
排放量
。2.
根据权利要求1所述的一种同步反演人为
CO2
排放与自然
CO2
通量的同化方法,其特征在于,所述
S1
中的太阳天顶角小于
70
°
,所述气溶胶光学厚度小于
0.3
,气溶胶光学厚度来自于碳含量检测装置平台上的大气多角度偏振探测仪气溶胶产品
。3.
根据权利要求1所述的一种同步反演人为
CO2
排放与自然
CO2
通量的同化方法,其特征在于,所述
S2
中大气
CO2、
温度

湿度和地表压力的香农信息含量计算公式为:公式
(1)
中,
K
代表权重函数
,
不同参数的权重函数表示对应参数单位扰动导致仿真光谱的变化量;
S
a
代表先验协方差矩阵,该矩阵为对角阵,对角线上元素为各参数波动范围的平方;
S
e
是测量误差协方差矩阵,该矩阵同样为对角阵,对角线上元素为光谱各通道对应仪器噪声的平方;利用大气
CO2、
温度

湿度和地表压力的香农信息含量,计算出综合香农信息含量,计算公式为:
H
tot

w1H
CO2
+w2H
T
+w3H
R
+w4H
P

ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
公式
(2)
中,
w
代表各参数信息含量的权重,计算用于
CO2
含量反演的信息时,
w1取3,
w2~
w4均取
‑1;计算用于背景扣除的信息时,
w1取7,
w2~
w4均取
1。4.
根据权利要求3所述的一种同步反演人为
CO2
排放与自然
CO2
通量的同化方法,其特征在于,所述
S3
中辐亮度比值谱计算公式为:公式
(3)
中,
i
代表各光谱通道,共
18
个;
Y
on
代表用于
CO2
含量反演的大气
CO2
信息谱,
Y
off
代表用于背景扣除的参考谱
。5.
根据权利要求4所述的一种同步反演人为
CO2
排放与自然
CO2
通量的同化方法,其特征在于,所述
S4
中估计反演算法的计算公式为:公式
(4)
中,
y
代表碳含量检测装置观测光谱,
F
代表前向模型,状态向量
x
...

【专利技术属性】
技术研发人员:卜晓燕王佳蕊王潇敏董锁成李小伟
申请(专利权)人:宁夏大学
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1