CN112540052A - 测算待测储罐的污染物排放量的测算方法及系统 - Google Patents
测算待测储罐的污染物排放量的测算方法及系统 Download PDFInfo
- Publication number
- CN112540052A CN112540052A CN201910894314.3A CN201910894314A CN112540052A CN 112540052 A CN112540052 A CN 112540052A CN 201910894314 A CN201910894314 A CN 201910894314A CN 112540052 A CN112540052 A CN 112540052A
- Authority
- CN
- China
- Prior art keywords
- storage tank
- pollutant
- detected
- light path
- tested
- 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.)
- Granted
Links
- 239000003344 environmental pollutant Substances 0.000 title claims abstract description 206
- 231100000719 pollutant Toxicity 0.000 title claims abstract description 206
- 238000003860 storage Methods 0.000 title claims abstract description 203
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012544 monitoring process Methods 0.000 claims abstract description 124
- 238000004364 calculation method Methods 0.000 claims abstract description 35
- 238000009792 diffusion process Methods 0.000 claims description 54
- 238000005033 Fourier transform infrared spectroscopy Methods 0.000 claims description 38
- 230000003287 optical effect Effects 0.000 claims description 33
- 238000012360 testing method Methods 0.000 claims description 12
- 238000012806 monitoring device Methods 0.000 claims description 9
- 238000004836 empirical method Methods 0.000 claims description 3
- 239000000356 contaminant Substances 0.000 description 10
- 238000005259 measurement Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000007716 flux method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000002329 infrared spectrum Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000004642 transportation engineering Methods 0.000 description 1
- 239000012855 volatile organic compound Substances 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3504—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
-
- 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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N2021/3595—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using FTIR
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Environmental & Geological Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及污染物监测领域，公开了一种测算待测储罐的污染物排放量的测算方法及系统，所述测算方法包括：获取所述待测储罐所在位置的气象参数；获取所述待测储罐的罐体参数；获取所述待测储罐的下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度；以及将所述气象参数、所述罐体参数、所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度输入预先建立的待测储罐的污染物排放量计算模型，以测算出所述待测储罐的污染物排放量，通过上述技术方案，排除待测储罐位置处环境中原有污染物对待测储罐污染物排放量测算过程的影响，并实现了对待测储罐污染物排放量的实时监测和计算，并提高测算精确度。
Description
技术领域
本发明涉及污染物监测领域，具体地涉及测算待测储罐的污染物排放量的测算方法及系统。
背景技术
储罐是石化企业中最常用的生产装置和最主要的无组织排放源之一，由于温度、压力以及现场的操作工艺等的影响，在油品的储存过程中会产生大量的呼吸损耗，一方面造成了资源的大量浪费，另一方面排放的大气污染物严重危害了人们的生命健康和生态环境，石化企业罐区无组织排放VOCs的源强计算、定量研究，对于保护环境不受污染、营造良好的大气生活环境具有关键的指导作用，已经成为国内外石油储运工程、化学工程、环境工程的重点研究课题。
现有技术中常用的核算方法主要有实测法、公式法、物料衡算法、经验估算或类比法、模型法、通量法、扩散模式反推法等，虽然这些方法得到了广泛的应用，但就其本身而言，仍然有许多不足之处：核算方法适用范围小，没有综合考虑对污染物排放量影响的因素，不能实时对污染物浓度进行测算。
发明内容
本发明的目的是为了克服现有技术存在的无法实现实时在线自动监测的问题，提供一种测算待测储罐的污染物排放量的测算方法，该测算待测储罐的污染物排放量的测算方法综合考虑了气象参数以及储罐参数，基于污染物积分浓度，实时对污染物排放量进行监测和计算。
为了实现上述目的，本发明一方面提供一种测算待测储罐的污染物排放量的测算方法，所述测算待测储罐的污染物排放量的测算方法包括：获取所述待测储罐所在位置的气象参数；获取所述待测储罐的罐体参数；获取所述待测储罐的下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度；以及将所述气象参数、所述罐体参数、所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度输入预先建立的待测储罐的污染物排放量计算模型，以测算出所述待测储罐的污染物排放量。
优选的，在将所述气象参数、所述罐体参数、所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度输入预先建立的待测储罐的污染物排放量计算模型之前，所述测算待测储罐的污染物排放量的测算方法还包括：剔除所述气象参数、所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度中的无效数据。
优选的，所述获取所述待测储罐的下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度包括：在所述待测储罐的上风向距离所述待测储罐的距离为第一设定距离处布置第一遥感傅里叶变换红外光谱仪，在所述待测储罐的下风向距离所述待测储罐的距离为第二设定距离处布置第二遥感傅里叶变换红外光谱仪；以及利用所述第一遥感傅里叶变换红外光谱仪和所述第二遥感傅里叶变换红外光谱仪进行监测，分别得到所述待测储罐的下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度。
优选的，所述罐体参数包括所述待测储罐的罐体高度和储罐直径，且所述气象参数包括所述罐体高度对应的平均风速、总云量及低云量；所述预先建立的待测储罐的污染物排放量计算模型为：
其中，QP表示所述待测储罐的污染物排放量；H表示所述待测储罐的罐体高度；y表示测量所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度的光程长度；u表示H高度处平均风速；CL表示所述下风向监测光路积分浓度；Cb表示所述上风向背景监测光路积分浓度；σy表示所述待测储罐的污染物的水平方向扩散参数；σz表示所述待测储罐的污染物的垂直方向扩散参数；
表示所述待测储罐的污染物的水平方向初始扩散参数；以及
表示所述待测储罐的污染物的垂直方向初始扩散参数，其中，根据所述气象参数和所述储罐直径确定所述待测储罐的污染物扩散参数σy、σz、
优选的，所述根据所述气象参数和所述储罐参数中的储罐直径确定所述待测储罐的污染物扩散参数σy、σz、
包括：以所述待测储罐的中心为原点，以主导风向为x轴，以横风向为y轴，建立平面直角坐标系；根据经验法确定所述待测储罐的污染物的水平方向初始扩散参数
和所述待测储罐的污染物的垂直方向初始扩散参数
以及根据所述待测储罐所在位置处的气象参数中的风向、风速、总云量、低云量、监测所述下风向监测光路污染物积分浓度的监测距离、所述待测储罐的污染物的水平方向初始扩散参数
和所述待测储罐的污染物的垂直方向初始扩散参数
计算得到所述水平方向扩散参数σy和所述垂直方向扩散参数σz。
根据本发明的另一方面，还提供一种测算待测储罐的污染物排放量的测算系统，所述测算待测储罐的污染物排放量的测算系统包括：气象站，用于获取所述待测储罐所在位置的气象参数；测量装置，用于获取所述待测储罐的罐体参数；监测装置，用于获取所述待测储罐的下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度；以及计算装置，用于将所述气象参数、所述罐体参数、所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度、输入预先建立的待测储罐的污染物排放量计算模型，以测算出所述待测储罐的污染物排放量。
优选的，所述监测装置包括：第一遥感傅里叶变换红外光谱仪，布置在所述待测储罐的上风向距离所述待测储罐的距离为第一设定距离处，用于监测所述上风向背景监测光路污染物积分浓度；以及第二遥感傅里叶变换红外光谱仪，布置在所述待测储罐的下风向距离所述待测储罐的距离为第二设定距离处，用于监测所述待测储罐的下风向监测光路污染物积分浓度。
优选的，所述计算装置包括：建模模块，用于预先建立所述待测储罐的污染物排放量计算模型；计算模块，用于将所述气象参数、所述罐体参数、所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度输入预先建立的待测储罐的污染物排放量计算模型，以测算出所述待测储罐的污染物排放量。
优选的，所述罐体参数包括所述待测储罐的罐体高度和储罐直径，且所述气象参数包括所述罐体高度对应的平均风速、总云量及低云量，并且所述建模模块预先建立的所述待测储罐的污染物排放量计算模型为：
其中，QP表示所述待测储罐的污染物排放量；H表示所述待测储罐的罐体高度；y表示测量所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度的光程长度；u表示H高度处平均风速；CL表示所述下风向监测光路积分浓度；Cb表示所述上风向背景监测光路积分浓度；σy表示所述待测储罐的污染物的水平方向扩散参数；σz表示所述待测储罐的污染物的垂直方向扩散参数；
表示所述待测储罐的污染物的水平方向初始扩散参数；以及
表示所述待测储罐的污染物的垂直方向初始扩散参数，其中，根据所述气象参数和所述储罐直径确定所述待测储罐的污染物扩散参数σy、σz、
此外，本发明实施例还提供一种机器可读存储介质，该机器可读存储介质上存储有指令，该指令用于使得机器执行上述测算待测储罐的污染物排放量的测算方法。
通过上述技术方案，综合考虑了气象参数以及储罐参数，基于污染物积分浓度，同时获取待测储罐的上风向背景监测光路污染物积分浓度对下风向监测光路污染物积分浓度，排除待测储罐位置处环境中原有污染物对待测储罐污染物排放量测算过程的影响，实现对污染物排放量的实时监测和计算，并提高测算精确度。
本发明实施例的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明实施例的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明实施例，但并不构成对本发明实施例的限制。
图1是本发明实施例一提供的测算待测储罐的污染物排放量的测算方法的流程图；
图2是本发明实施例二提供的测算待测储罐的污染物排放量的测算系统的框图；
图3是本发明实施例三提供的测算待测储罐的污染物排放量的测算系统的监测装置的框图；以及
图4是本发明实施例四提供的测算待测储罐的污染物排放量的测算系统的计算装置的框图。
附图标记说明
1、气象站 2、监测装置 3、计算装置
21、第一遥感傅里叶变换红外光谱仪
22、第二遥感傅里叶变换红外光谱仪
31、建模模块 32、计算模块
具体实施方式
以下结合附图对本发明实施例的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明实施例，并不用于限制本发明实施例。
实施例一
图1是本发明实施例一提供的测算待测储罐的污染物排放量的测算方法的流程图；如图1所示，可以包括：
S100、获取待测储罐所在位置的气象参数。
优选的，在待测储罐区域内进行一定时间的气象条件观测并收集储罐相关数据。
具体来讲，在进行储罐现场污染物测算的数据采集之前，应首先进行至少一周的气象观测和实地考察，以获取待测储罐所在位置的气象参数。
其中，气象参数主要通过在待测储罐位置处所架设的气象站得到，气象站的架设高度应为10米且不受储罐及周围建筑物的遮挡。获取气象参数主要包括确定待测储罐所在位置处的污染物排放量测算过程中的主导风向以及罐体高度对应的平均风速。
此外，需要对待测储罐所在位置及周围地形进行实地考察，确定待测储罐的主导风向上风向处和下风向处均有足够空间，以保证待测储罐附近能够架设监测仪器，例如遥感傅里叶变换红外光谱仪和气象站。以及，在本发明后续计算大气扩散参数时，需要考虑到地形的复杂程度。
S200、获取待测储罐的罐体参数。
优选的，罐体参数包括待测储罐的罐体高度和储罐直径，可通过储罐生产厂家提供或实际测量获得。
S300、获取待测储罐的下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度。
优选的，在待测储罐的上风向距离待测储罐的距离为第一设定距离处布置第一遥感傅里叶变换红外光谱仪，在待测储罐的下风向距离待测储罐的距离为第二设定距离处布置第二遥感傅里叶变换红外光谱仪；以及利用第一遥感傅里叶变换红外光谱仪和第二遥感傅里叶变换红外光谱仪进行监测，分别得到待测储罐的下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度。
具体来讲，布置第一遥感傅里叶变换红外光谱仪和第二遥感傅里叶变换红外光谱仪时，首先要尽量保证遥感傅里叶变换红外光谱仪发射光路与主导风向垂直，并近似接近于地面。由于在监测过程中，风向是不断变化的，所以要根据前期的气象观测，确定主导风向，因此本发明中并不要求遥感傅里叶变换红外光谱仪发射光路与主导风向一定要垂直，而只要近似垂直即可。
其中，设定第一设定距离和第二设定距离的设置，是为了保证第一遥感傅里叶变换红外光谱仪和第二遥感傅里叶变换红外光谱仪所采集的数据免受待测储罐自身扩散的影响。设定第一设定距离时要保证待测储罐的中心至第一遥感傅里叶变换红外光谱仪发射光路的距离足够长。设定第二设定距离时要保证待测储罐中心至第二遥感傅里叶变换红外光谱仪发射光路的距离足够长，例如，监测距离应至少>储罐直径。其次，对第一设定距离和第二设定距离进行设定时，还需要分别充分考虑第一遥感傅里叶变换红外光谱仪和第二遥感傅里叶变换红外光谱仪的监测距离，在保证前述条件基础上，尽量缩短第一设定距离和第二设定距离，以保证遥感傅里叶变换红外光谱仪的监测精度。如此，能够保证遥感傅里叶变换红外光谱仪的自身灵敏度，保证其所采集的数据的精度，进而使得对待测储罐污染物排放量的监测和计算更为精确。。
需要说明的是，在该步骤中在待测储罐的上风向和下风处分别设置遥感傅里叶变换红外光谱仪，主要是因为某些物质在待测储罐所处的环境中本身可能就有一定浓度，以下风向所监测到的污染物浓度减去上风向处所监测到的污染物浓度的差值来计算该待测储罐的污染物排放量更为准确，具有实际参考价值。
S400、测算出待测储罐的污染物排放量。
优选的，将气象参数、罐体参数、下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度输入预先建立的待测储罐的污染物排放量计算模型，以测算出待测储罐的污染物排放量。
在本发明的优选实施中，在将气象参数、罐体参数、下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度输入预先建立的待测储罐的污染物排放量计算模型之前，进行如下操作：剔除气象参数、下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度中的无效数据。
具体来讲，气象参数、罐体参数、下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度形成了多组监测数据，其中一个数据不符合监测和计算要求时会对待测储罐的污染物排放量测算造成较大影响，因此，当其中一个或多个参数不符合要求时，需要剔除该参数对应的一组数据。
举例说明，当平均风速<1m/s；和/或风向波动角度>30°；和/或大气稳定度为A/B状态；和/或遥感傅里叶变换红外光谱仪的信号强度过低时所监测到的数据均会对待测储罐污染物的测算带来较大影响，导致测算结果不准确，因此，需要剔除相应的一组数据。其中，风向波动角度是指风向与光路的垂线之间的夹角，大气稳定度可以根据P.S(Pasquill-Turner stability cl-assification，帕斯奎尔—特纳大气稳定度分类法)大气稳定度分类法确定。
优选的，预先建立的待测储罐的污染物排放量计算模型为：
其中，QP表示待测储罐的污染物排放量；H表示待测储罐的罐体高度；y表示测量下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度的光程长度；u表示H高度处平均风速；CL表示下风向监测光路积分浓度；Cb表示上风向背景监测光路积分浓度；σy表示待测储罐的污染物的水平方向扩散参数；σz表示待测储罐的污染物的垂直方向扩散参数；
表示待测储罐的污染物的水平方向初始扩散参数；以及
表示待测储罐的污染物的垂直方向初始扩散参数，其中，根据气象参数和储罐直径确定待测储罐的污染物扩散参数σy、σz、
上风向背景监测光路积分浓度Cb和下风向监测光路积分浓度CL可分别通过第一遥感傅里叶变换红外光谱仪和第二遥感傅里叶变换红外光谱仪进行监测得到。
在本发明的优选实施例中，通过以下方式来确定待测储罐的污染物扩散参数σy、σz、
以待测储罐的中心为原点，以主导风向为x轴，以横风向为y轴，建立平面直角坐标系；根据经验法确定待测储罐的污染物的水平方向初始扩散参数
和待测储罐的污染物的垂直方向初始扩散参数
以及根据待测储罐所在位置处的气象参数中的风向、风速、待测储罐的污染物的水平方向初始扩散参数
和待测储罐的污染物的垂直方向初始扩散参数
下风向监测距离(下风向遥感傅里叶变换红外光谱仪发射光束与污染源中心的垂直距离)，计算得到水平方向扩散参数σy和垂直方向扩散参数σz。此外，在计算σy、σz时，需要根据地形复杂状况进行修正。
实施例二
图2是本发明实施例二提供的测算待测储罐的污染物排放量的测算系统的框图，如图2所示，可以包括：气象站1，用于获取待测储罐所在位置的气象参数；测量装置，用于获取待测储罐的罐体参数；监测装置2，用于获取待测储罐的下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度；以及计算装置3，用于将气象参数、罐体参数、下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度、输入预先建立的待测储罐的污染物排放量计算模型，以测算出待测储罐的污染物排放量。
实施例三
图3是本发明实施例三提供的测算待测储罐的污染物排放量的测算系统的监测装置的框图，如图3所示，监测装置2包括：第一遥感傅里叶变换红外光谱仪21，布置在待测储罐的上风向距离待测储罐的距离为第一设定距离处，用于监测上风向背景监测光路污染物积分浓度；以及第二遥感傅里叶变换红外光谱仪22，布置在待测储罐的下风向距离待测储罐的距离为第二设定距离处，用于监测待测储罐的下风向监测光路污染物积分浓度。
实施例四
图4是本发明实施例四提供的测算待测储罐的污染物排放量的测算系统的计算装置的框图，如图4所示，计算装置3包括：建模模块31，用于预先建立待测储罐的污染物排放量计算模型；计算模块32，用于将气象参数、罐体参数、下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度输入预先建立的待测储罐的污染物排放量计算模型，以测算出待测储罐的污染物排放量。
优选的，罐体参数包括待测储罐的罐体高度和储罐直径，且气象参数包括罐体高度对应的平均风速，并且建模模块31预先建立的待测储罐的污染物排放量计算模型为：
其中，QP表示待测储罐的污染物排放量；H表示待测储罐的罐体高度；y表示测量下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度的光程长度；u表示H高度处平均风速；CL表示下风向监测光路积分浓度；Cb表示上风向背景监测光路积分浓度；σy表示待测储罐的污染物的水平方向扩散参数；σz表示待测储罐的污染物的垂直方向扩散参数；
表示待测储罐的污染物的水平方向初始扩散参数；以及
表示待测储罐的污染物的垂直方向初始扩散参数，其中，根据气象参数和储罐直径确定待测储罐的污染物扩散参数σy、σz、
测算待测储罐的污染物排放量的测算系统的其他具体实施细节和有益效果同上述测算待测储罐的污染物排放量的测算方法，此处不再赘述。
通过上述技术方案，综合考虑了气象和地形参数以及储罐参数，基于污染物积分浓度，同时获取待测储罐的上风向背景监测光路污染物积分浓度对下风向监测光路污染物积分浓度，排除待测储罐位置处环境中原有污染物对待测储罐污染物排放量测算过程的影响，实现对污染物排放量的实时监测和计算，并提高测算精确度。
以上结合附图详细描述了本发明实施例的可选实施方式，但是，本发明实施例并不限于上述实施方式中的具体细节，在本发明实施例的技术构思范围内，可以对本发明实施例的技术方案进行多种简单变型，这些简单变型均属于本发明实施例的保护范围。
另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明实施例对各种可能的组合方式不再另行说明。
本领域技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序存储在一个存储介质中，包括若干指令用以使得单片机、芯片或处理器(processor)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
此外，本发明实施例的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明实施例的思想，其同样应当视为本发明实施例所公开的内容。
Claims (10)
1.一种测算待测储罐的污染物排放量的测算方法，其特征在于，所述测算待测储罐的污染物排放量的测算方法包括：
获取所述待测储罐所在位置的气象参数；
获取所述待测储罐的罐体参数；
获取所述待测储罐的下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度；以及
将所述气象参数、所述罐体参数、所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度输入预先建立的待测储罐的污染物排放量计算模型，以测算出所述待测储罐的污染物排放量。
2.根据权利要求1所述的测算待测储罐的污染物排放量的测算方法，其特征在于，在将所述气象参数、所述罐体参数、所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度输入预先建立的待测储罐的污染物排放量计算模型之前，所述测算待测储罐的污染物排放量的测算方法还包括：
剔除所述气象参数、所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度中的无效数据。
3.根据权利要求1所述的测算待测储罐的污染物排放量的测算方法，其特征在于，所述获取所述待测储罐的下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度包括：
在所述待测储罐的上风向距离所述待测储罐的距离为第一设定距离处布置第一遥感傅里叶变换红外光谱仪，在所述待测储罐的下风向距离所述待测储罐的距离为第二设定距离处布置第二遥感傅里叶变换红外光谱仪；以及
利用所述第一遥感傅里叶变换红外光谱仪和所述第二遥感傅里叶变换红外光谱仪进行监测，分别得到所述待测储罐的下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度。
4.根据权利要求1所述的测算待测储罐的污染物排放量的测算方法，其特征在于，所述罐体参数包括所述待测储罐的罐体高度和储罐直径，且所述气象参数包括所述罐体高度对应的平均风速、总云量及低云量；
所述预先建立的待测储罐的污染物排放量计算模型为：
其中，QP表示所述待测储罐的污染物排放量；
H表示所述待测储罐的罐体高度；
y表示测量所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度的光程长度；
u表示H高度处平均风速；
CL表示所述下风向监测光路积分浓度；
Cb表示所述上风向背景监测光路积分浓度；
σy表示所述待测储罐的污染物的水平方向扩散参数；
σz表示所述待测储罐的污染物的垂直方向扩散参数；
其中，根据所述气象参数和所述储罐直径确定所述待测储罐的污染物扩散参数σy、σz、σy0、
5.根据权利要求4所述的测算待测储罐的污染物排放量的测算方法，其特征在于，所述根据所述气象参数和所述储罐参数中的储罐直径确定所述待测储罐的污染物扩散参数σy、σz、σy0、
包括：
以所述待测储罐的中心为原点，以主导风向为x轴，以横风向为y轴，建立平面直角坐标系；
根据经验法确定所述待测储罐的污染物的水平方向初始扩散参数
和所述待测储罐的污染物的垂直方向初始扩散参数
以及
根据所述待测储罐所在位置处的气象参数中的风向、风速、总云量、低云量、监测所述下风向监测光路污染物积分浓度的监测距离、所述待测储罐的污染物的水平方向初始扩散参数
和所述待测储罐的污染物的垂直方向初始扩散参数
计算得到所述水平方向扩散参数σy和所述垂直方向扩散参数σz。
6.一种测算待测储罐的污染物排放量的测算系统，其特征在于，所述测算待测储罐的污染物排放量的测算系统包括：
气象站，用于获取所述待测储罐所在位置的气象参数；
测量装置，用于获取所述待测储罐的罐体参数；
监测装置，用于获取所述待测储罐的下风向监测光路污染物积分浓度和上风向背景监测光路污染物积分浓度；以及
计算装置，用于将所述气象参数、所述罐体参数、所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度、输入预先建立的待测储罐的污染物排放量计算模型，以测算出所述待测储罐的污染物排放量。
7.根据权利要求6所述的测算待测储罐的污染物排放量的测算系统，其特征在于，所述监测装置包括：
第一遥感傅里叶变换红外光谱仪，布置在所述待测储罐的上风向距离所述待测储罐的距离为第一设定距离处，用于监测所述上风向背景监测光路污染物积分浓度；以及
第二遥感傅里叶变换红外光谱仪，布置在所述待测储罐的下风向距离所述待测储罐的距离为第二设定距离处，用于监测所述待测储罐的下风向监测光路污染物积分浓度。
8.根据权利要求6所述的测算待测储罐的污染物排放量的测算系统，其特征在于，所述计算装置包括：
建模模块，用于预先建立所述待测储罐的污染物排放量计算模型；
计算模块，用于将所述气象参数、所述罐体参数、所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度输入预先建立的待测储罐的污染物排放量计算模型，以测算出所述待测储罐的污染物排放量。
9.根据权利要求8所述的测算待测储罐的污染物排放量的测算系统，其特征在于，所述罐体参数包括所述待测储罐的罐体高度和储罐直径，且所述气象参数包括所述罐体高度对应的平均风速、总云量及低云量，并且所述建模模块预先建立的所述待测储罐的污染物排放量计算模型为：
其中，QP表示所述待测储罐的污染物排放量；
H表示所述待测储罐的罐体高度；
y表示测量所述下风向监测光路污染物积分浓度和所述上风向背景监测光路污染物积分浓度的光程长度；
u表示H高度处平均风速；
CL表示所述下风向监测光路积分浓度；
Cb表示所述上风向背景监测光路积分浓度；
σy表示所述待测储罐的污染物的水平方向扩散参数；
σz表示所述待测储罐的污染物的垂直方向扩散参数；
其中，根据所述气象参数和所述储罐直径确定所述待测储罐的污染物扩散参数σy、σz、
10.一种机器可读存储介质，该机器可读存储介质上存储有指令，该指令用于使得机器执行本申请权利要求1-5中任一项所述的测算待测储罐的污染物排放量的测算方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910894314.3A CN112540052B (zh) | 2019-09-20 | 2019-09-20 | 测算待测储罐的污染物排放量的测算方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910894314.3A CN112540052B (zh) | 2019-09-20 | 2019-09-20 | 测算待测储罐的污染物排放量的测算方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112540052A true CN112540052A (zh) | 2021-03-23 |
| CN112540052B CN112540052B (zh) | 2024-07-30 |
Family
ID=75012526
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910894314.3A Active CN112540052B (zh) | 2019-09-20 | 2019-09-20 | 测算待测储罐的污染物排放量的测算方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112540052B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116541654A (zh) * | 2023-07-06 | 2023-08-04 | 中汽研汽车检验中心(昆明)有限公司 | 农用收获机械污染物排放量测算方法及装置 |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101694461A (zh) * | 2009-10-16 | 2010-04-14 | 中国科学院安徽光学精密机械研究所 | 一种污染源气体排放通量红外多组份监测方法与系统 |
| CN201503545U (zh) * | 2009-06-25 | 2010-06-09 | 中国农业大学 | 污染源实时监控系统 |
| RU2012149045A (ru) * | 2012-11-19 | 2014-05-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Государственный университет по землеустройству" | Способ оперативного контроля атмосферных загрязнений локальных территорий |
| KR20150031577A (ko) * | 2013-09-16 | 2015-03-25 | 주식회사 에니텍 | 대기오염물질 배출량 역산출 방법 |
| CN106841087A (zh) * | 2017-04-17 | 2017-06-13 | 中国石油化工股份有限公司 | 石油化工企业无组织排放源vocs排放量测算方法 |
| CN107064046A (zh) * | 2017-04-17 | 2017-08-18 | 中国石油化工股份有限公司 | 挥发性有机物排放总量的计算方法 |
| CN107331125A (zh) * | 2017-08-03 | 2017-11-07 | 深圳市盛路物联通讯技术有限公司 | 一种交通污染预警方法及系统 |
| CN107436343A (zh) * | 2017-07-31 | 2017-12-05 | 南京南瑞集团公司 | 一种模拟计算敏感区域污染物浓度的方法 |
| CN109115966A (zh) * | 2018-11-11 | 2019-01-01 | 西藏神州瑞霖环保科技股份有限公司 | 一种可实时监测污染物浓度及扩散范围的监测设备 |
| US20190113445A1 (en) * | 2017-10-16 | 2019-04-18 | Weimin Zhang | Air pollution monitoring system and air pollution monitoring method |
| CN109855898A (zh) * | 2019-02-27 | 2019-06-07 | 中国石油大学(华东) | 单拉井井场储罐、敞开液面污油池VOCs排放模拟装置 |
-
2019
- 2019-09-20 CN CN201910894314.3A patent/CN112540052B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201503545U (zh) * | 2009-06-25 | 2010-06-09 | 中国农业大学 | 污染源实时监控系统 |
| CN101694461A (zh) * | 2009-10-16 | 2010-04-14 | 中国科学院安徽光学精密机械研究所 | 一种污染源气体排放通量红外多组份监测方法与系统 |
| RU2012149045A (ru) * | 2012-11-19 | 2014-05-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Государственный университет по землеустройству" | Способ оперативного контроля атмосферных загрязнений локальных территорий |
| KR20150031577A (ko) * | 2013-09-16 | 2015-03-25 | 주식회사 에니텍 | 대기오염물질 배출량 역산출 방법 |
| CN106841087A (zh) * | 2017-04-17 | 2017-06-13 | 中国石油化工股份有限公司 | 石油化工企业无组织排放源vocs排放量测算方法 |
| CN107064046A (zh) * | 2017-04-17 | 2017-08-18 | 中国石油化工股份有限公司 | 挥发性有机物排放总量的计算方法 |
| CN107436343A (zh) * | 2017-07-31 | 2017-12-05 | 南京南瑞集团公司 | 一种模拟计算敏感区域污染物浓度的方法 |
| CN107331125A (zh) * | 2017-08-03 | 2017-11-07 | 深圳市盛路物联通讯技术有限公司 | 一种交通污染预警方法及系统 |
| US20190113445A1 (en) * | 2017-10-16 | 2019-04-18 | Weimin Zhang | Air pollution monitoring system and air pollution monitoring method |
| CN109115966A (zh) * | 2018-11-11 | 2019-01-01 | 西藏神州瑞霖环保科技股份有限公司 | 一种可实时监测污染物浓度及扩散范围的监测设备 |
| CN109855898A (zh) * | 2019-02-27 | 2019-06-07 | 中国石油大学(华东) | 单拉井井场储罐、敞开液面污油池VOCs排放模拟装置 |
Non-Patent Citations (7)
| Title |
|---|
| CHAMBERS, ALLAN K.等: "Direct measurement of fugitive emissions of hydrocarbons from a refinery", JOURNAL OF THE AIR&WASTE MANAGEMENT ASSOCIATION, 31 October 2008 (2008-10-31), pages 1047 - 1056 * |
| 朱亮;贾润中;高少华;邹兵;: "基于OP-FTIR技术的炼油装置VOCs排放量监测", 安全、健康和环境, no. 05, 31 May 2017 (2017-05-31), pages 24 - 28 * |
| 杨旭: "工业废气无组织排放污染的监测和评价方法", 中国卫生工程学, vol. 9, no. 1, pages 30 - 38 * |
| 蒋维楣主编: "《空气污染气象学》", vol. 2003, 30 April 2004, 南京大学出版社, pages: 93 - 98 * |
| 薛诚;毕俊;: "某电解铝厂氟化物无组织排放量的估算", 中国卫生工程学, vol. 8, no. 5, pages 257 * |
| 邹兵;王国龙;赵东风;: "石化罐区挥发性有机物源强反演技术的研究", 环境工程, no. 10, 31 October 2017 (2017-10-31), pages 87 - 90 * |
| 霍玉侠;李发荣;仝纪龙;袁九毅;: "石化企业储罐区无组织排放大气环境影响及对策研究", 环境科学与技术, no. 07, 31 July 2011 (2011-07-31), pages 195 - 199 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116541654A (zh) * | 2023-07-06 | 2023-08-04 | 中汽研汽车检验中心(昆明)有限公司 | 农用收获机械污染物排放量测算方法及装置 |
| CN116541654B (zh) * | 2023-07-06 | 2023-09-12 | 中汽研汽车检验中心(昆明)有限公司 | 农用收获机械污染物排放量测算方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112540052B (zh) | 2024-07-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20200110019A1 (en) | Atmospheric pollution source mapping and tracking of pollutants by using air quality monitoring networks having high space-time resolution | |
| KR101751642B1 (ko) | 대기 관측용 라이다의 소산계수 보정방법 | |
| CN111798928B (zh) | 大气颗粒物污染源解析方法和装置 | |
| CN105334147A (zh) | 基于β射线法和光散射法的颗粒物在线监测系统及方法 | |
| CA2843753C (en) | Methods for calibrating a fugitive emission rate measurement | |
| CN107612999A (zh) | 大气网格化精准监控系统 | |
| Brown et al. | Standardisation of a European measurement method for organic carbon and elemental carbon in ambient air: results of the field trial campaign and the determination of a measurement uncertainty and working range | |
| CN115060850B (zh) | 空地双场耦合的大气污染源追踪和通量测量装置及方法 | |
| CN106092841A (zh) | 以pm2.5质量浓度为约束条件反演气溶胶消光系数吸湿增长因子与相对湿度函数关系的方法 | |
| CN112364940A (zh) | 基于多源数据的大气污染物源解析方法、存储介质及设备 | |
| Margaritis et al. | Calibration of low-cost gas sensors for air quality monitoring | |
| CN112540052A (zh) | 测算待测储罐的污染物排放量的测算方法及系统 | |
| CN113218912A (zh) | 一种基于远程质控监控平台的烟气在线监测方法及系统 | |
| JP2740718B2 (ja) | ガス、蒸気等の漏洩地点および漏洩量推定システム | |
| CN116258101B (zh) | 一种基于气体浓度监测数据快速计算排放源强的方法 | |
| CN117455110A (zh) | 一种考虑碳排放点源扩散的碳排监测及评估方法 | |
| CN110907318B (zh) | 一种近地面大气总悬浮颗粒物质量浓度遥感物理估算方法 | |
| CN113029989B (zh) | 基于激光传感器的气体检测方法、装置、设备及存储介质 | |
| KR20010097331A (ko) | 레이저 레이더 및 확산모델을 이용한 대기오염 정보시스템 | |
| CN112259172B (zh) | 一种存在空气污染物化学反应情况的污染源逆向辨识方法 | |
| CN115705510A (zh) | 厂区气体污染分布预测方法、系统、电子设备及存储介质 | |
| CN207410366U (zh) | 大气网格化精准监控系统 | |
| CN117849279B (zh) | 一种温室气体排放源位置及排放量的反演方法 | |
| KR102538345B1 (ko) | 원격분광측정차량을 이용한 미세먼지 원인 물질 배출원 추적방법 | |
| CN115950797B (zh) | 污染物溯源方法和系统 |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20240701 Address after: 100728 No. 22 North Main Street, Chaoyang District, Beijing, Chaoyangmen Applicant after: CHINA PETROLEUM & CHEMICAL Corp. Country or region after: China Applicant after: Sinopec Safety Engineering Research Institute Co.,Ltd. Address before: 100728 No. 22 North Main Street, Chaoyang District, Beijing, Chaoyangmen Applicant before: CHINA PETROLEUM & CHEMICAL Corp. Country or region before: China Applicant before: SINOPEC Research Institute OF SAFETY ENGINEERING |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |