CN115032332B - 一种基于车载***测量强点源碳排放的方法 - Google Patents
一种基于车载***测量强点源碳排放的方法 Download PDFInfo
- Publication number
- CN115032332B CN115032332B CN202210469695.2A CN202210469695A CN115032332B CN 115032332 B CN115032332 B CN 115032332B CN 202210469695 A CN202210469695 A CN 202210469695A CN 115032332 B CN115032332 B CN 115032332B
- Authority
- CN
- China
- Prior art keywords
- value
- point source
- vehicle
- feasible region
- strong point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 37
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 20
- 239000003016 pheromone Substances 0.000 claims abstract description 38
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 18
- 238000001914 filtration Methods 0.000 claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 13
- 238000011065 in-situ storage Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 37
- 230000007704 transition Effects 0.000 claims description 15
- 238000009792 diffusion process Methods 0.000 claims description 12
- 210000003746 feather Anatomy 0.000 claims description 6
- 230000008901 benefit Effects 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 230000007774 longterm Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000006386 neutralization reaction Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000010792 warming Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/0004—Gaseous mixtures, e.g. polluted air
- G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
- G01N33/0062—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method or the display, e.g. intermittent measurement or digital display
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/02—Instruments for indicating weather conditions by measuring two or more variables, e.g. humidity, pressure, temperature, cloud cover or wind speed
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/84—Greenhouse gas [GHG] management systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Pure & Applied Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Health & Medical Sciences (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Environmental & Geological Engineering (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computational Linguistics (AREA)
- Combustion & Propulsion (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Computer Networks & Wireless Communication (AREA)
Abstract
本发明提供一种基于车载***测量强点源碳排放的方法,首先利用车载***(包括原位测量仪、GNSS接收机、车载气象站)获取带有地理坐标的CO2浓度以及风速、风向信息,根据先验强点源排放信息,建立坐标系,转换CO2浓度坐标至强点源坐标系,根据测量信息分别以相关系数和均方根误差作为信息素,借助蚂蚁算法分别计算两种信息素标准下的可行域与初值,对可行域与初始值进行均值滤波,将滤波结果作为内点罚函数的初值与均值进行求解,获取使罚函数最小的最优解,得到兼具高精确度、高时效性性和高可用特征的强点源碳排放产品。
Description
技术领域
本发明涉及环境监测领域,根据车载测量***在强点源附近测量所获得的实测数据,结合高斯烟羽模型定量模拟CO2排放,实现对强点源碳排放进行高准确度的实时监测。
背景技术
全球温室效应已成为世界各国关心的重要议题,巴黎协定的签署意味着世界所有主要经济体在应对气候变化、实现碳中和这一全人类面临的重大问题上达成了共识。CO2作为主要的温室气体,其造成的全球变暖效应被越来越多人关注和研究,控制CO2排放成为有效改善气候变化的有力手段。作为命运共同体的提出者,我国十分重视环境治理问题。有关资料证明,强点源碳排放占人为总排放的31%,如何准确实时检测这一类碳排放成为碳中和首要面临的问题之一。
目前现有手段主要包括通过统计年鉴对其进行粗略性的年际评估方式和基于卫星遥感监测手段测算碳排放方式。其中,通过统计年鉴对其进行粗略性的年际评估方式其时效性与精度都难以满足碳排放监测目的,而基于卫星遥感监测手段测算碳排放方式虽然具有高时效性和高分辨率的特点,但由于卫星受限于轨道和大气情况影响,只能在特定飞行时间对同一地点在大气环境良好的情况下进行碳排放测算,并且在无法提供准确大气信息情况时,现有卫星反演手段会出现难以收敛现象,同时在云层覆盖情况下也难以进行碳排放反演工作,这些因素在很大程度上限制了通过卫星监测方式实现碳排放的检测。
发明内容
本发明针对现有碳排放监测手段不足,提供一种基于车载***测量强点源碳排放的方法,充分利用车载***高灵活性、气象数据的高时效性以及现有高斯烟羽模型,得到兼具高精确度、高时效性性和高可用特征的强点源碳排放产品。
本专利提出一种基于车载***测量强点源碳排放的方法,可以精确且实时地测算强点源碳排放。得益于车载***的灵活简便性,在实际测量过程中,可以通过合理的测量轨迹规划,从而获取大量CO2浓度信息。通过高精度的CO2浓度信息、气象参数等辅助数据以及高斯烟羽模型,本专利所提方法可以快速精确计算强点源碳排放数据,该方法的推广将为我国碳排放监测事业提供一种新的检测手段。
为了达到上述目的,本发明提供的技术方案是一种基于车载***测量强点源碳排放的方法,包括以下步骤:
步骤1,通过车载***测量必要观测信息,包括车载***实时GNSS定位信息、CO2浓度信息、大气观测信息;
步骤2,根据先验强点源排放信息,获取高斯烟羽模型所需信息,建立新坐标系,将步骤1输出的CO2浓度的坐标转换到新坐标系下;
步骤3,根据步骤2输出的新坐标系下的CO2浓度测量值及风速,分别以相关系数以及均方根误差作为适应函数,计算高斯烟羽模型未知参数对应的可行域1、可行域2以及初始值1和初始值2,所述可行域1表示通过相关性系数计算方式计算得到的未知参数的最小值与最大值,所述可行域2表示通过均方根误差计算方式计算得到的未知参数的最小值与最大值,初始值1和初始值2分别表示通过相关性系数和均方根误差计算方式计算得到的未知参数的均值;
步骤4,将步骤3解算的可行域1、可行域2和初始值1、初始值2分别进行均值滤波,获取高斯烟羽模型中未知参数的可行域以及初始值,通过内点罚函数求解未知参数,输出所测强点源排放处的碳排放量。
进一步的,所述车载***包括原位测量仪、GNSS接收机、车载气象站。
进一步的,步骤1中测量必要观测信息的具体实现方式如下;
首先根据强排放点源长久的累积气象资料,充分利用车载***灵活性高的优势,规划车载***测量路径以及车载***行驶速度,车载***在规划路径上按照规定速度行驶,通过车载GNSS***获取实时定位信息,通过车载原位测量仪获取CO2浓度信息,通过车载气象仪获取风速、风向信息,最终输出带有地理信息的CO2测量浓度、风速以及风向。
进一步的,步骤2中,根据先验强点源排放信息,以强点源排放点地理坐标作为坐标轴原点,以测量风向为X轴,水平方向上垂直于X轴的方向为Y轴,以垂直于XOY轴方向为Z轴方向,输出新的XYZ坐标;根据以先验强点源排放信息建立的XYZ坐标系与原始GNSS坐标系关系转换GNSS地理坐标到新的 XYZ坐标系下,输出新的XYZ坐标系下的CO2浓度测量值和风速信息。
进一步的,步骤3中所述高斯烟羽模型由公式1给出:
其中C(x,y,z)为(x,y,z)坐标处的CO2浓度,(x,y,z)为步骤2所建立的新坐标系下的坐标值,C(x,y,z)为该坐标系下测量的CO2浓度,q为强点源的排放强度,u为风速,H为强点源处的碳排放的有效排放高度,σy和σz分别为水平扩散参数和垂直扩散参数,其中σy=a·xb、σz=c·xd,a和b为水平扩散系数,c和d为垂直扩散系数,α为地面反射系数,B为背景CO2浓度,在本式中, C(x,y,z)和u为已知值,待求解量为(q,H,a,b,c,d,B,α,z)。
进一步的,步骤3中通过蚁群算法计算高斯烟羽模型未知参数的最优估计值,蚁群算法是一种在给定信息素条件下存在求解全局最优解(q,H,a,b,c,d,B,α,z)的智能算法,其原理为:
首先随机产生若干个起始点初始化未知数参量[qi,Hi,ai,bi,ci,di,Bi,αi,zi],根据每个未知数参量计算该参量的信息素τ;
根据相关性计算的第i个参量[qi,Hi,ai,bi,ci,di,Bi,αi,zi]的信息素τ1i由公式2给出;
其中Ci(xj,yj,zj)′为第i个参量在高斯烟羽模型下所计算出在(xj,yj,zj)处的CO2浓度,n为测量值的总数量,Var为方差;
根据均方根误差所计算出的第i个参量[qi,Hi,ai,bi,ci,di,Bi,αi,zi]的信息素τ2i由公式3给出;
之后,根据两种不同计算方式的信息素分别计算两种不同计算方式下的状态转移概率;
其中Pi 1为以相关性作为信息素计算方式下第i个参量的转移概率,Pi 2为以均方根误差作为信息素计算方式下第i个参量的转移概率,τ1为τ1i的集合,τ2为τ2i的集合;
根据转移概率与信息素的比较,可以计算出第i个参量在iter+1迭代数下的数值,当第i个参量的信息素小于其转移概率时,第i个参量在iter+1迭代数下的数值由公式6给出,否则以公式7给出;
[qi,Hi,ai,bi,ci,di,Bi,αi,zi]iter+1=[qi,Hi,ai,bi,ci,di,Bi,αi,zi]iter+r1·step·λ (6)
[qi,Hi,ai,bi,ci,di,Bi,αi,zi]iter+1=[qi,Hi,ai,bi,ci,di,Bi,αi,zi]iter+r2·range (7)
其中,r1为介于[-1,1]的随机数,step为局部搜索步长,λ为当前迭代次数iter 的倒数,r2为[-0.5,0.5]的随机数,range为自变量[q,H,a,b,c,d,B,α,z]的区间大小;
计算更新后参量的信息素τ(iter+1)′,若其小于τiter,则参量不进行更改;其中 Rho为挥发因子,可根据先验知识设置;
τ(iter+1)′=(1-Rho)·τiter+τiter+1 (8)
等待蚁群算法在设定迭代次数中迭代完毕之后,可获得蚁群算法在本次迭代中的最优估计值[q,H,a,b,c,d,B,α,z];分别获取两种不同信息素计算方式下的各组最优估计值,分别对其求均值、最小值和最大值,便可输出可行域1的下界 Xl1=[ql1,Hl1,al1,bl1,cl1,dl1,Bl1,αl1,zl1],即通过相关性计算方式计算得到的最优估计值的最小值,上界Xu1=[qu1,Hu1,au1,bu1,cu1,du1,Bu1,αu1,zu1],即通过相关性计算方式计算得到的最优估计值的最大值,可行域2的下界 Xl2=[ql2,Hl2,al2,bl2,cl2,dl2,Bl2,αl2,zl2],即通过均方根误差计算方式计算得到的最优估计值的最小值,上界Xu2=[qu2,Hu2,au2,bu2,cu2,du2,Bu2,αu2,zu2],即通过均方根误差计算方式计算得到的最优估计值的最大值,初始解1
X10=[q10,H10,a10,b10,c10,d10,B10,α10,z10],即通过相关性计算方式计算得到的最优估计值的均值,初始解2X20=[q20,H20,a20,b20,c20,d20,B20,α20,z20],即通过均方根误差计算方式计算得到的最优估计值的均值。
进一步的,步骤4中是将步骤3输出的可行域1的下界与可行域2的下界作均值滤波,可行域1的上界与可行域2的上届作均值滤波,初始值1与初始值2 作均值滤波后,分别获得可行域的下界xl=[ql,Hl,al,bl,cl,dl,Bl,αl,zl],可行域上界 xu=[qu,Hu,au,bu,cu,du,Bu,αu,zu]以及初始值x0=[q0,H0,a0,b0,c0,d0,B0,α0,z0],根据可行域可设置惩罚函数c;
c=f(X)+σmp(X) (9)
通过最小化惩罚函数c,可获得第m+1次迭代中的X(m+1);
其中,X=[q,H,a,b,c,d,B,α,z],X1=q,X2=H,X3=a,X4=b,X5=c, X6=d,X7=B,X8=α,X9=z,将公式10求解出的X(m+1)带入σmp(X(m+1)),若σmp(X(m+1))<ε,则停止迭代输出最优解否则更新σm+1=βσm继续迭代,其中ε和β为预设常数。
与现有技术相比,本发明具有如下优点:1)充分利用车载***高灵活度特点,合理规划探测路线,实现智能探测CO2排放;2)创新提出一种基于车载***测量强点源碳排放的方法,结合蚁群算法全局最优搜索和内点罚函数精确搜索特性,实现具有更高精度潜力的碳排放检测。
附图说明
图1为本发明实施例的流程图。
图2为本发明实施例坐标系建立图。
具体实施方式
本发明提供一种基于车载***测量强点源碳排放的方法,首先利用车载***(包括原位测量仪、GNSS接收机、车载气象站)获取带有地理坐标的CO2浓度以及风速、风向信息,根据先验强点源排放信息,建立坐标系,转换CO2浓度坐标至强点源坐标系,根据测量信息分别以相关系数和均方根误差作为信息素,借助蚂蚁算法分别计算两种信息素标准下的可行域与初值,对可行域与初始值进行均值滤波,将滤波结果作为内点罚函数的初值与均值进行求解,获取使罚函数最小的最优解,得到兼具高精确度、高时效性性和高可用特征的强点源碳排放产品。
下面结合附图和实施例对本发明的技术方案作进一步说明。
如图1所示,本发明实施例的流程包括以下步骤:
步骤1,通过车载***测量必要观测信息,包括(车载***实时GNSS定位信息、CO2浓度信息、大气观测信息)。
首先根据强排放点源长久的累积气象资料,充分利用车载***灵活性高的优势,规划车载***测量路径以及车载***行驶速度。车载***在规划路径上按照规定速度行驶,通过车载GNSS***获取实时定位信息,通过车载原位测量仪获取CO2浓度信息,通过车载气象仪获取风速、风向信息。最终输出带有地理信息的CO2测量浓度、风速以及风向。
步骤2,根据先验强点源排放信息,获取高斯烟羽模型所需信息,建立坐标系。
根据先验强点源排放信息,以强点源排放点地理坐标作为坐标轴原点,以测量风向为X轴,水平方向上垂直于X轴的方向为Y轴,以垂直于XOY轴方向为Z轴方向,输出新的XYZ坐标。根据以先验强点源排放信息建立的XYZ坐标系与原始GNSS坐标系关系转换GNSS地理坐标到新的XYZ坐标系下,输出新的XYZ坐标系下的CO2浓度测量值和风速信息。
步骤3,根据步骤2输出的现有坐标系下的CO2浓度测量值及风速,分别以相关系数以及均方根误差作为适应函数,计算高斯烟羽模型未知参数的可行域 1、可行域2以及初始值1和初始值2。
步骤3中所述高斯烟羽模型可由公式1给出:
其中C(x,y,z)为(x,y,z)坐标处的CO2浓度,在本专利中,(x,y,z)为步骤2 所建立的新坐标系下的坐标值,C(x,y,z)为该坐标系下测量的CO2浓度,q为强点源的排放强度,u为风速,H为强点源处的碳排放的有效排放高度,σy和σz分别为水平扩散参数和垂直扩散参数,其中σy=a·xb、σz=c·xd,a和b为水平扩散系数,c和d为垂直扩散系数,α为地面反射系数,B为背景CO2浓度。在本式中,C(x,y,z)和u为已知值,待求解量为[q,H,a,b,c,d,B,α,z]。
蚁群算法是一种在给定信息素条件下存在求解全局最优解(q,H,a,b,c,d,B,α,z)的智能算法,其原理为:
首先随机产生300个起始点初始化未知数参量[qi,Hi,ai,bi,ci,di,Bi,αi,zi](i∈(1,300)),根据每个未知数参量计算该参量的信息素τ。
根据相关性计算的第i个参量[qi,Hi,ai,bi,ci,di,Bi,αi,zi](i∈(1,300))的信息素τ1i由公式2给出。
其中Ci(xj,yj,zj)'为第i个参量在高斯烟羽模型下所计算出在(xj,yj,zj)处的CO2浓度,n为测量值的总数量,Var为方差
根据均方根误差所计算出的第i个参量[qi,Hi,ai,bi,ci,di,Bi,αi,zi](i∈(1,300))的信息素τ2i由公式3给出。
之后,根据两种不同计算方式的信息素分别计算两种不同计算方式下的状态转移概率。
其中Pi 1为以相关性作为信息素计算方式下第i个参量的转移概率,Pi 2为以均方根误差作为信息素计算方式下第i个参量的转移概率,τ1=[τ11,τ12,τ13......τ1300],τ2=[τ21,τ22,τ23......τ2300]。
根据转移概率与信息素的比较,可以计算出第i个参量在iter+1迭代数下的数值。当第i个参量的信息素小于其转移概率时,第i个参量在iter+1迭代数下的数值由公式6给出,否则以公式7给出。
[qi,Hi,ai,bi,ci,di,Bi,αi,zi]iter+1=[qi,Hi,ai,bi,ci,di,Bi,αi,zi]iter+r1·step·λ (6)
[qi,Hi,ai,bi,ci,di,Bi,αi,zi]iter+1=[qi,Hi,ai,bi,ci,di,Bi,αi,zi]iter+r2·range (7)
其中,r1为介于[-1,1]的随机数,step为局部搜索步长,λ为当前迭代次数iter 的倒数,r2为[-0.5,0.5]的随机数,range为自变量[q,H,a,b,c,d,B,α,z]的区间大小。
计算更新后参量的信息素τ(iter+1)',若其小于τiter,则参量不进行更改。其中 Rho为挥发因子,可根据先验知识设置。
τ(iter+1)'=(1-Rho)·τiter+τiter+1 (8)
等待蚁群算法在设定迭代次数中迭代完毕之后,可获得蚁群算法在本次迭代中的最优估计值[q,H,a,b,c,d,B,α,z],两种不同信息素计算方式下可分别获得两种不同计算方式的最优估计值。在本专利中,进行10000次蚁群算法的迭代,可分别获取两种不同信息素计算方式下的各10000组最优估计值,分别对其求均值、最小值和最大值,便可输出可行域1的下界 Xl1=[ql1,Hl1,al1,bl1,cl1,dl1,Bl1,αl1,zl1],即通过相关性计算方式计算得到的最优估计值的最小值,上界Xu1=[qu1,Hu1,au1,bu1,cu1,du1,Bu1,αu1,zu1],即通过相关性计算方式计算得到的最优估计值的最大值,可行域2的下界 Xl2=[ql2,Hl2,al2,bl2,cl2,dl2,Bl2,αl2,zl2],即通过均方根误差计算方式计算得到的最优估计值的最小值,上界Xu2=[qu2,Hu2,au2,bu2,cu2,du2,Bu2,αu2,zu2],即通过均方根误差计算方式计算得到的最优估计值的最大值,初始解1 X10=[q10,H10,a10,b10,c10,d10,B10,α10,z10],即通过相关性计算方式计算得到的最优估计值的均值,初始解2X20=[q20,H20,a20,b20,c20,d20,B20,α20,z20],即通过均方根误差计算方式计算得到的最优估计值的均值。
步骤4,将步骤3解算的可行域1、可行域2和初始值1、初始值2分别进行均值滤波,获取高斯烟羽模型中未知参数的可行域以及初始值,通过内点罚函数求解未知参数,输出所测强点源排放处的碳排放量。
步骤4中所述通过均值滤波获取可行域与初始值是将步骤3输出的可行域1 的下界与可行域2的下界作均值滤波,可行域1的上界与可行域2的上界作均值滤波,初始值1与初始值2作均值滤波后,分别可以获得可行域的下界 Xl=[ql,Hl,al,bl,cl,dl,Bl,αl,zl],可行域上界Xu=[qu,Hu,au,bu,cu,du,Bu,αu,zu]以及初始值 X0=[q0,H0,a0,b0,c0,d0,B0,α0,z0],根据可行域可设置惩罚函数c。
c=f(X)+σmp(X) (9)
其中, p(X)=-ln(X-Xl)-ln(Xu-X),CX(xj,yj,zj)为高斯烟羽模型在参数X情况下模拟位于(xj,yj,zj)处的CO2浓度,σm为第m次迭代中的罚参数,一般设为大于0的数。
通过最小化惩罚函数c,可获得第m+1次迭代中的X(m+1)。
其中,X=[q,H,a,b,c,d,B,α,z],X1=q,X2=H,X3=a,X4=b,X5=c, X6=d,X7=B,X8=α,X9=z,将公式10求解出的X(m+1)带入σmp(X(m+1)),若σmp(X(m+1))<ε,则停止迭代输出最优解否则更新σm+1=βσm继续迭代。其中ε和β为预设常数。
具体实施时,以上流程可采用计算机软件技术实现自动运行流程。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (5)
1.一种基于车载***测量强点源碳排放的方法,其特征在于,包括如下步骤:
步骤1,通过车载***测量必要观测信息,包括车载***实时GNSS定位信息、CO2浓度信息、大气观测信息;
步骤2,根据先验强点源排放信息,获取高斯烟羽模型所需信息,建立新坐标系,将步骤1输出的CO2浓度的坐标转换到新坐标系下;
步骤2中所述高斯烟羽模型由公式1给出:
其中C(x,y,z)为(x,y,z)坐标处的CO2浓度,(x,y,z)为步骤2所建立的新坐标系下的坐标值,C(x,y,z)为该坐标系下测量的CO2浓度,q为强点源的排放强度,u为风速,H为强点源处的碳排放的有效排放高度,σy和σz分别为水平扩散参数和垂直扩散参数,其中σy=a·xb、σz=c·xd,a和b为水平扩散系数,c和d为垂直扩散系数,α为地面反射系数,B为背景CO2浓度,在本式中,C(x,y,z)和u为已知值,待求解量为(q,H,a,b,c,d,B,α,z);
步骤3,根据步骤2输出的新坐标系下的CO2浓度测量值及风速,分别以相关系数以及均方根误差作为适应函数,计算高斯烟羽模型未知参数对应的可行域1、可行域2以及初始值1和初始值2,所述可行域1表示通过相关性系数计算方式计算得到的未知参数的最小值与最大值,所述可行域2表示通过均方根误差计算方式计算得到的未知参数的最小值与最大值,初始值1和初始值2分别表示通过相关性系数和均方根误差计算方式计算得到的未知参数的均值;
步骤3中通过蚁群算法计算高斯烟羽模型未知参数的最优估计值,蚁群算法是一种在给定信息素条件下存在求解全局最优解(q,H,a,b,c,d,B,α,z)的智能算法,其原理为:
首先随机产生若干个起始点初始化未知数参量[qi,Hi,ai,bi,ci,di,Bi,αi,zi],根据每个未知数参量计算该参量的信息素τ;
根据相关性计算的第i个参量[qi,Hi,ai,bi,ci,di,Bi,αi,zi]的信息素τ1i由公式2给出;
其中Ci(xj,yj,zj)'为第i个参量在高斯烟羽模型下所计算出在(xj,yj,zj)处的CO2浓度,n为测量值的总数量,Var为方差;
根据均方根误差所计算出的第i个参量[qi,Hi,ai,bi,ci,di,Bi,αi,zi]的信息素τ2i由公式3给出;
之后,根据两种不同计算方式的信息素分别计算两种不同计算方式下的状态转移概率;
其中Pi 1为以相关性作为信息素计算方式下第i个参量的转移概率,Pi 2为以均方根误差作为信息素计算方式下第i个参量的转移概率,τ1为τ1i的集合,τ2为τ2i的集合;
根据转移概率与信息素的比较,可以计算出第i个参量在iter+1迭代数下的数值,当第i个参量的信息素小于其转移概率时,第i个参量在iter+1迭代数下的数值由公式6给出,否则以公式7给出;
[qi,Hi,ai,bi,ci,di,Bi,αi,zi]iter+1=[qi,Hi,ai,bi,ci,di,Bi,αi,zi]iter+r1·step·λ (6)
[qi,Hi,ai,bi,ci,di,Bi,αi,zi]iter+1=[qi,Hi,ai,bi,ci,di,Bi,αi,zi]iter+r2·range (7)
其中,r1为介于[-1,1]的随机数,step为局部搜索步长,λ为当前迭代次数iter的倒数,r2为[-0.5,0.5]的随机数,range为自变量[q,H,a,b,c,d,B,α,z]的区间大小;
计算更新后参量的信息素τ(iter+1)',若其小于τiter,则参量不进行更改;其中Rho为挥发因子,可根据先验知识设置;
τ(iter+1)'=(1-Rho)·τiter+τiter+1 (8)
等待蚁群算法在设定迭代次数中迭代完毕之后,可获得蚁群算法在本次迭代中的最优估计值[q,H,a,b,c,d,B,α,z];分别获取两种不同信息素计算方式下的各组最优估计值,分别对其求均值、最小值和最大值,便可输出可行域1的下界Xl1=[ql1,Hl1,al1,bl1,cl1,dl1,Bl1,αl1,zl1],即通过相关性计算方式计算得到的最优估计值的最小值,上界Xu1=[qu1,Hu1,au1,bu1,cu1,du1,Bu1,αu1,zu1],即通过相关性计算方式计算得到的最优估计值的最大值,可行域2的下界Xl2=[ql2,Hl2,al2,bl2,cl2,dl2,Bl2,αl2,zl2],即通过均方根误差计算方式计算得到的最优估计值的最小值,上界Xu2=[qu2,Hu2,au2,bu2,cu2,du2,Bu2,αu2,zu2],即通过均方根误差计算方式计算得到的最优估计值的最大值,初始解1X10=[q10,H10,a10,b10,c10,d10,B10,α10,z10],即通过相关性计算方式计算得到的最优估计值的均值,初始解2X20=[q20,H20,a20,b20,c20,d20,B20,α20,z20],即通过均方根误差计算方式计算得到的最优估计值的均值;
步骤4,将步骤3解算的可行域1、可行域2和初始值1、初始值2分别进行均值滤波,获取高斯烟羽模型中未知参数的可行域以及初始值,通过内点罚函数求解未知参数,输出所测强点源排放处的碳排放量。
2.如权利要求1所述的一种基于车载***测量强点源碳排放的方法,其特征在于:所述车载***包括原位测量仪、GNSS接收机、车载气象站。
3.如权利要求1所述的一种基于车载***测量强点源碳排放的方法,其特征在于:步骤1中测量必要观测信息的具体实现方式如下;
首先根据强排放点源长久的累积气象资料,充分利用车载***灵活性高的优势,规划车载***测量路径以及车载***行驶速度,车载***在规划路径上按照规定速度行驶,通过车载GNSS***获取实时定位信息,通过车载原位测量仪获取CO2浓度信息,通过车载气象仪获取风速、风向信息,最终输出带有地理信息的CO2测量浓度、风速以及风向。
4.如权利要求1所述的一种基于车载***测量强点源碳排放的方法,其特征在于:步骤2中,根据先验强点源排放信息,以强点源排放点地理坐标作为坐标轴原点,以测量风向为X轴,水平方向上垂直于X轴的方向为Y轴,以垂直于XOY轴方向为Z轴方向,输出新的XYZ坐标;根据以先验强点源排放信息建立的XYZ坐标系与原始GNSS坐标系关系转换GNSS地理坐标到新的XYZ坐标系下,输出新的XYZ坐标系下的CO2浓度测量值和风速信息。
5.如权利要求1所述的一种基于车载***测量强点源碳排放的方法,其特征在于:步骤4中是将步骤3输出的可行域1的下界与可行域2的下界作均值滤波,可行域1的上界与可行域2的上界作均值滤波,初始值1与初始值2作均值滤波后,分别获得可行域的下界xl=[ql,Hl,al,bl,cl,dl,Bl,αl,zl],可行域上界xu=[qu,Hu,au,bu,cu,du,Bu,αu,zu]以及初始值x0=[q0,H0,a0,b0,c0,d0,B0,α0,z0],根据可行域可设置惩罚函数c;
c=f(X)+σmp(X) (9)
通过最小化惩罚函数c,可获得第m+1次迭代中的X(m+1);
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210469695.2A CN115032332B (zh) | 2022-04-28 | 2022-04-28 | 一种基于车载***测量强点源碳排放的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210469695.2A CN115032332B (zh) | 2022-04-28 | 2022-04-28 | 一种基于车载***测量强点源碳排放的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115032332A CN115032332A (zh) | 2022-09-09 |
CN115032332B true CN115032332B (zh) | 2023-04-28 |
Family
ID=83118998
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210469695.2A Active CN115032332B (zh) | 2022-04-28 | 2022-04-28 | 一种基于车载***测量强点源碳排放的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115032332B (zh) |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20050004316A (ko) * | 2003-06-27 | 2005-01-12 | 재단법인 포항산업과학연구원 | 복수 기상관측자료를 이용한 대기오염 확산 예측 방법 |
CN104181276A (zh) * | 2013-05-28 | 2014-12-03 | 东北大学 | 一种基于无人机的企业碳排放量检测方法 |
CN107436343A (zh) * | 2017-07-31 | 2017-12-05 | 南京南瑞集团公司 | 一种模拟计算敏感区域污染物浓度的方法 |
CN112711002B (zh) * | 2020-12-23 | 2023-05-16 | 武汉大学 | 一种基于co2-dial模拟测量的点源co2排放的新型估算方法 |
CN112686531B (zh) * | 2020-12-29 | 2021-11-09 | 生态环境部卫星环境应用中心 | 一种结合卫星遥感与车载观测的大气污染企业识别方法 |
CN113624921B (zh) * | 2021-08-04 | 2024-03-29 | 杭州谱育科技发展有限公司 | 多模式综合污染溯源方法 |
CN114088633B (zh) * | 2021-11-19 | 2022-09-27 | 生态环境部卫星环境应用中心 | 一种星地协同监测的煤矿区甲烷排放异常识别和核算方法 |
-
2022
- 2022-04-28 CN CN202210469695.2A patent/CN115032332B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN115032332A (zh) | 2022-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109580003B (zh) | 一种静止气象卫星热红外数据估算近地面大气温度方法 | |
CN114037911B (zh) | 一种顾及生态分区的大尺度森林高度遥感反演方法 | |
Pan et al. | The potential of CO2 satellite monitoring for climate governance: A review | |
CN112684520B (zh) | 一种气象预报的订正方法、装置、计算机设备及存储介质 | |
Reuter et al. | A joint effort to deliver satellite retrieved atmospheric CO 2 concentrations for surface flux inversions: the ensemble median algorithm EMMA | |
CN111487216B (zh) | 一种二氧化碳通量反演方法、*** | |
Mendoza-Dominguez et al. | Estimation of emission adjustments from the application of four-dimensional data assimilation to photochemical air quality modeling | |
Coman et al. | Assimilation of IASI partial tropospheric columns with an Ensemble Kalman Filter over Europe | |
CN106597416A (zh) | 一种地面GPS辅助的LiDAR数据高程差的误差修正方法 | |
CN110389087B (zh) | 一种污染天气下的pm2.5浓度卫星遥感估算方法 | |
CN112711002B (zh) | 一种基于co2-dial模拟测量的点源co2排放的新型估算方法 | |
CN110595968A (zh) | 一种基于静止轨道卫星的pm2.5浓度估算方法 | |
CN110595960B (zh) | 一种基于机器学习的pm2.5浓度遥感估算方法 | |
CN114417580B (zh) | 一种观测***对全球电离层数据同化性能的影响评估方法 | |
CN114898823B (zh) | 一种高时空分辨率遥感近地面no2浓度预估方法和*** | |
CN115480032A (zh) | 基于地面遥感测量的点源排放强度预测方法 | |
CN114595876A (zh) | 一种区域风场预测模型生成方法和装置、电子设备 | |
Chen et al. | Proactive quality control: Observing system simulation experiments with the Lorenz’96 model | |
CN113468799A (zh) | 一种静止气象卫星观测中获取近地面pm2.5浓度的方法及*** | |
CN115032332B (zh) | 一种基于车载***测量强点源碳排放的方法 | |
CN107656905A (zh) | 一种利用误差传递的空气质量数据实时校准方法 | |
Yang et al. | Toward Establishing a Low-cost UAV Coordinated Carbon Observation Network (LUCCN): First Integrated Campaign in China | |
CN116029883A (zh) | 基于二氧化碳激光雷达中尺度碳排放定位核查方法及*** | |
AU2021100706A4 (en) | An Integrated Assessment Method of Climate Change Based on Multi-Source Carbon Data | |
Izanlou et al. | GNSS-derived precipitable water vapor modeling using machine learning methods |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |