CN105115859B - 用于监测泄漏气体扩散浓度的方法及装置 - Google Patents
用于监测泄漏气体扩散浓度的方法及装置 Download PDFInfo
- Publication number
- CN105115859B CN105115859B CN201510456988.7A CN201510456988A CN105115859B CN 105115859 B CN105115859 B CN 105115859B CN 201510456988 A CN201510456988 A CN 201510456988A CN 105115859 B CN105115859 B CN 105115859B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- msup
- mtr
- 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
- 238000012544 monitoring process Methods 0.000 title claims abstract description 242
- 238000009792 diffusion process Methods 0.000 title claims abstract description 142
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims description 33
- 238000006243 chemical reaction Methods 0.000 claims description 22
- 238000012937 correction Methods 0.000 claims description 17
- 230000005855 radiation Effects 0.000 claims description 11
- 238000012986 modification Methods 0.000 claims description 4
- 230000004048 modification Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 abstract description 4
- 239000006185 dispersion Substances 0.000 abstract 2
- 239000007789 gas Substances 0.000 description 81
- 238000010586 diagram Methods 0.000 description 9
- 238000013468 resource allocation Methods 0.000 description 5
- 238000000605 extraction Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 239000002341 toxic gas Substances 0.000 description 3
- 238000013507 mapping Methods 0.000 description 2
- 231100000572 poisoning Toxicity 0.000 description 2
- 230000000607 poisoning effect Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- 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
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Examining Or Testing Airtightness (AREA)
Abstract
本发明公开了一种用于监测泄漏气体扩散浓度的方法及装置。该方法包括:采集泄漏源和地理区域上与监测点对应的位置点在地理坐标系中的坐标,以及泄漏源的泄漏强度;根据泄漏源和与监测点对应的位置点在地理坐标系中的坐标、以及地理区域的平均风向角,确定泄漏源和监测点在笛卡尔坐标系中的坐标;根据泄漏源在地理坐标系中的坐标、地理区域的气象参数、监测时间、以及监测点在笛卡尔坐标系中的坐标,确定监测点处的实际扩散系数;根据监测点处的实际扩散系数、泄漏源在笛卡尔坐标系中的坐标、监测点在笛卡尔坐标系中的坐标和泄漏强度,确定监测点处的泄漏气体扩散浓度。由此,能够在地形起伏变化较大的条件下,准确监测出泄漏气体的扩散浓度。
Description
技术领域
本发明涉及泄漏气体浓度监测领域,具体地,涉及一种用于监测泄漏气体扩散浓度的方法及装置。
背景技术
随着工业技术的不断进步,人类将面临越来越多的安全生产问题、污染排放问题等,气体(特别是有毒气体)泄漏问题则是其中之一。泄漏的有毒气体会在风力的作用下快速传播,导致大面积人员中毒、伤亡,甚至导致重大、特大事故。因此,一旦有毒气体泄漏事故发生,需要在迅速控制住泄漏的同时,迅速疏散泄漏源周围一定范围内的居民,并进行中毒人员救援,而这些都需要首先确定气体扩散范围、趋势,从而确定警戒与疏散范围、撤离路线以及救援资源分配等。上述这些工作均需要建立在对气体扩散浓度的准确监测的基础上。因此如何准确、有效地监测泄漏气体扩散浓度是社会安全保障的基础,也是提升应急救援的效率的前提之一。
发明内容
本发明的目的是提供一种用于监测泄漏气体扩散浓度的方法及装置,该方法及装置能够在地形起伏变化较大的条件下,准确监测出泄漏气体的扩散浓度。
为了实现上述目的,本发明提供一种用于监测泄漏气体扩散浓度的方法,该方法包括:采集泄漏源和地理区域上与监测点对应的位置点在地理坐标系中的坐标,以及所述泄漏源的泄漏强度;根据所述泄漏源和所述与监测点对应的位置点在所述地理坐标系中的坐标、以及所述地理区域的平均风向角,确定所述泄漏源和所述监测点在笛卡尔坐标系中的坐标,其中,所述笛卡尔坐标系为以所述泄漏源向地面的铅直投影点为原点,以所述地理区域的平均风向为x轴正向,以垂直所述泄漏源所在地的海平面、指向地球外部为z轴正向的右手笛卡尔坐标系;根据所述泄漏源在所述地理坐标系中的坐标、所述地理区域的气象参数、监测时间、以及所述监测点在所述笛卡尔坐标系中的坐标,确定所述监测点处的实际扩散系数;根据所述监测点处的实际扩散系数、所述泄漏源在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标和所述泄漏源的泄漏强度,确定所述监测点处的泄漏气体扩散浓度。
本发明还一种用于监测泄漏气体扩散浓度的装置,该装置包括:采集模块,用于采集泄漏源和地理区域上与监测点对应的位置点在地理坐标系中的坐标,以及所述泄漏源的泄漏强度;笛卡尔坐标确定模块,用于根据所述泄漏源和所述与监测点对应的位置点在所述地理坐标系中的坐标、以及所述地理区域的平均风向角,确定所述泄漏源和所述监测点在笛卡尔坐标系中的坐标,其中,所述笛卡尔坐标系为以所述泄漏源向地面的铅直投影点为原点,以所述地理区域的平均风向为x轴正向,以垂直所述泄漏源所在地的海平面、指向地球外部为z轴正向的右手笛卡尔坐标系;实际扩散系数确定模块,用于根据所述泄漏源在所述地理坐标系中的坐标、所述地理区域的气象参数、监测时间、以及所述监测点在所述笛卡尔坐标系中的坐标,确定所述监测点处的实际扩散系数;泄漏气体扩散浓度确定模块,用于根据所述监测点处的实际扩散系数、所述泄漏源在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标和所述泄漏源的泄漏强度,确定所述监测点处的泄漏气体扩散浓度。
通过上述技术方案,可进行大范围泄漏气体扩散浓度的监测,适用性强,能够与地理信息***相结合。通过获取泄漏源相关参数以及气象参数,能够快速确定出监测点处的泄漏气体扩散浓度,并快速确定出泄漏气体的污染区域,特别是在高起伏地形条件下,也能够准确地确定泄漏气体扩散浓度分布,从而为确定警戒与疏散范围、撤离路线以及救援资源分配等提供准确、可靠的数据支持。
本发明的其他特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明的进一步理解,并且构成说明书的一部分,与下面的具体实施方式一起用于解释本发明,但并不构成对本发明的限制。在附图中:
图1是根据本发明的实施方式的用于监测泄漏气体扩散浓度的方法的流程图;
图2是泄漏源、监测点与一地理区域的示意图;
图3是根据本发明的另一实施方式的用于监测泄漏气体扩散浓度的方法的流程图;
图4是根据本发明的另一实施方式的用于监测泄漏气体扩散浓度的方法的流程图;
图5是根据本发明的另一实施方式的用于监测泄漏气体扩散浓度的方法的流程图;
图6是泄漏源、监测点、泄漏源的镜像位置与一地理区域的示意图;
图7是根据本发明的实施方式的用于监测泄漏气体扩散浓度的装置的框图;
图8是根据本发明的另一实施方式的用于监测泄漏气体扩散浓度的装置的框图;
图9是根据本发明的另一实施方式的用于监测泄漏气体扩散浓度的装置的框图;
图10是根据本发明的另一实施方式的用于监测泄漏气体扩散浓度的装置的框图。
具体实施方式
以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本发明,并不用于限制本发明。
图1示出了根据本发明的实施方式的用于监测泄漏气体扩散浓度的方法的流程图。如图1所示,该方法可以包括:
在步骤S101中,采集泄漏源和地理区域上与监测点对应的位置点在地理坐标系中的坐标,以及所述泄漏源的泄漏强度。其中,所述地理区域的范围可以被预先设定,例如,可以根据实际的监测需求进行设定。
图2示出了泄漏源、监测点与一地理区域的示意图。如图2所示,泄漏源所在位置可以由点L表示。监测点为点P所指示的位置,在本发明中,监测点可以是位于地理区域上方的监测点。在监测点位于地理区域上方的情况下,监测点P在地理区域上存在一个对应的位置点W,例如,监测点P垂直投射到地理区域上的点。
泄漏源L在地理坐标系中的坐标可以表示为(Lo0,La0,H0),其中,Lo0表示泄漏源L的经度(单位为度,°),La0表示泄漏源L的纬度(单位为度,°),H0表示泄漏源L的有效海拔高度(单位为m)。Lo0、La0可以通过GPS(全球定位***)接收机进行测量,或者可以通过手机定位和地图进行数值提取。H0可以通过H00及泄漏源L距地面的有效高度进行叠加得到,其中,H00表示泄漏源L所在位置的地面海拔高度(单位为m),并且H00可以通过GPS接收机进行测量,或通过手机定位和地图进行数值提取。
与监测点P对应的位置点W在地理坐标系中的坐标可以表示为(Loij,Laij,Hij),其中,Loij表示经度(单位为度,°),Laij表示纬度(单位为度,°),Hij表示海拔(单位为m),Loij,Laij,Hij可以通过GPS接收机进行测量,或者可以通过手机定位和地图进行数值提取。
在步骤S102中,根据所述泄漏源和所述与监测点对应的位置点在所述地理坐标系中的坐标、以及所述地理区域的平均风向角,确定所述泄漏源和所述监测点在笛卡尔坐标系中的坐标,其中,所述笛卡尔坐标系为以所述泄漏源向地面的铅直投影点为原点,以所述地理区域的平均风向为x轴正向,以垂直所述泄漏源所在地的海平面、指向地球外部为z轴正向的右手笛卡尔坐标系。
通过这一步骤,可以将泄漏源在地理坐标系中的坐标(Lo0,La0,H0)转换为在笛卡尔坐标系中的坐标(X0,Y0,Z0),以及将与监测点对应的位置点在地理坐标系中的坐标(Loij,Laij,Hij)转换为在笛卡尔坐标系中的坐标(Xij,Yij,Zij),并且,可以根据与监测点对应的位置点在笛卡尔坐标系中的坐标(Xij,Yij,Zij),得出监测点在笛卡尔坐标系中的坐标(Xij,Yij,Zpij)。
在步骤S103中,根据所述泄漏源在所述地理坐标系中的坐标、所述地理区域的气象参数、监测时间、以及所述监测点在所述笛卡尔坐标系中的坐标,确定所述监测点处的实际扩散系数。其中,所述气象参数可以例如包括总云量、低云量、风速等参数。
在步骤S104中,根据所述监测点处的实际扩散系数、所述泄漏源在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标和所述泄漏源的泄漏强度,确定所述监测点处的泄漏气体扩散浓度。
通过上述技术方案,可进行大范围泄漏气体扩散浓度的监测,适用性强,能够与地理信息***相结合。通过获取泄漏源相关参数以及气象参数,能够快速确定出监测点处的泄漏气体扩散浓度,并快速确定出泄漏气体的污染区域,特别是在高起伏地形条件下,也能够准确地确定泄漏气体扩散浓度分布,从而为确定警戒与疏散范围、撤离路线以及救援资源分配等提供准确、可靠的数据支持。
图3示出了根据本发明的另一实施方式的用于监测泄漏气体扩散浓度的方法的流程图。如图3所示,所述步骤S102可以包括:
在步骤S301中,根据所述与监测点对应的位置点在所述地理坐标系中的坐标、所述地理区域的平均海拔、以及第一地理坐标与全球坐标的转换函数,确定所述与监测点对应的位置点在全球坐标系中的坐标。
所述地理区域中可以包括至少一个位置点,其中一个位置点为与所述监测点对应的位置点,如图2所示的位置点W。可以通过GPS接收机测量出(或者可以通过手机定位和地图进行数值提取出)地理区域中每一个位置点在地理坐标系中的坐标,即,可以得到每一个位置点的经度、纬度和海拔。所述每一个位置点的经度可以组成一个M×N的二阶经度矩阵Lo,即,
其中,i=1,2,……,M;j=1,2,……,N,M×N可以表示所述地理区域上的位置点的总数。同理,所述每一个位置点的纬度可以组成一个M×N的二阶纬度矩阵La,所述每一个位置点的海拔可以组成一个M×N的二阶海拔矩阵H,即,
这样,就可以根据地理区域中每一个位置点的海拔确定出该地理区域的平均海拔Havg,即:
其中,Hij可以表示地理区域中任一个位置点的海拔。
之后,就可以根据所述与监测点对应的位置点在所述地理坐标系中的坐标Loij,Laij,Hij、所述地理区域的平均海拔Havg、以及第一地理坐标与全球坐标的转换函数,确定所述与监测点对应的位置点在全球坐标系中的坐标。其中,所述全球坐标系为以所述地理区域内海拔高度为0的地球椭球表面这个近似平面为xoy平面,以与监测点对应的位置点指向地心方向的反方向为z轴正方向的O1xyz坐标系。
在一个实施方式中,所述第一地理坐标与全球坐标的转换函数可以包括:
其中,LXij、LYij、LZij表示所述与监测点对应的位置点在所述全球坐标系中的坐标;Havg表示所述地理区域的平均海拔;Loij、Laij、Hij表示所述与监测点对应的位置点在所述地理坐标系中的坐标;Ngl表示地球的卯酉圈曲率半径,并且agl表示地球椭球长半轴(例如,agl=6378.137×103m),Wgl为中间变量,并且e表示地球椭球的第一偏心率,bgl表示地球椭球的短半轴(例如,bgl=6356.7523141×103m)。
这样,就可以将与监测点对应的位置点在所述地理坐标系中的坐标转换为在全球坐标系O1xyz中的坐标。
在步骤S302中,根据所述泄漏源在所述地理坐标系中的坐标、所述地理区域的平均海拔、以及第二地理坐标与全球坐标的转换函数,确定所述泄漏源在全球坐标系中的坐标。
在一个实施方式中,所述第二地理坐标与全球坐标的转换函数可以包括:
其中,LX0、LY0、LZ0表示所述泄漏源在所述全球坐标系中的坐标;Lo0、La0、H0表示所述泄漏源在所述地理坐标系中的坐标。
这样,就可以将泄漏源在所述地理坐标系中的坐标转换为在全球坐标系O1xyz中的坐标。
在步骤S303中,根据所述与监测点对应的位置点在所述全球坐标系中的坐标、所述泄漏源在所述全球坐标系中的坐标、所述地理区域的平均风向角、以及第一全球坐标与笛卡尔坐标的转换函数,确定所述与监测点对应的位置点在所述笛卡尔坐标系中的坐标。
在一个实施方式中,所述第一全球坐标与笛卡尔坐标的转换函数可以包括:
其中,Xij、Yij、Zij表示所述与监测点对应的位置点在所述笛卡尔坐标系中的坐标;θ0表示所述全球坐标系O1xyz的x轴正向与正东方向的夹角;H00表示所述泄漏源所在位置的地面海拔高度,α表示所述地理区域的平均风向角(即,风向与正东方向的夹角),其中,α∈[0,2π),单位为弧度(rad)。
在一个实施方式中,可以通过风向测定仪测得所述地理区域的平均风向角。或者,在另一个实施方式中,可以通过气象资料或气象预报得到的所述地理区域的风向,来确定所述地理区域的平均风向角,其中,风向所对应的平均风向角α如下:
在一个实施方式中,可以通过以下方式来确定θ0:
首先,采集所述地理区域内的任意两个位置点在地理坐标系中的坐标,例如,第一位置点(Lo11,La11,H11)和第二位置点(Lo12,La12,H12)。
接下来,根据等式(2),可以得出这两个位置点在全球坐标系中的坐标(LX11,LY11,H11)和(LX12,LY12,H12)。
接下来,利用以下等式(5)和等式(6)构建向量v1和v2,其中,v1表示在地理坐标系中第二个位置点到第一个位置点的向量,v2表示在全球坐标系中第二个位置点到第一个位置点的向量:
v1=(Lo12-Lo11,La12-La11)(|v1|≠0) (5)
v2=(LX12-LX11,LY12-LY11)(|v2|≠0) (6)
之后,可以根据向量v1和v2,利用等式(7),确定出全球坐标系O1xyz的x轴正向与正东方向的夹角θ0(单位为度,°):
将O1xyz坐标系绕z轴旋转至x轴正向指向与所述地理区域的平均风向一致的方向,且原点平移至泄漏源向地面的铅直投影点处,得到新坐标系O2xyz,即,所述笛卡尔坐标系。
之后,再利用(4),可以将与监测点对应的位置点在所述全球坐标系O1xyz中的坐标转换为在所述笛卡尔坐标系O2xyz中的坐标。
在步骤S304中,根据所述与监测点对应的位置点在所述笛卡尔坐标系中的坐标,确定所述监测点在所述笛卡尔坐标系中的坐标。
如上所述,监测点在地理区域上存在一个对应的位置点,例如,监测点垂直投射到地理区域上的点。因此,监测点及其相对应的位置点可以具有相同的经、纬度,并且仅在海拔上有所不同。可以根据与监测点对应的位置点在笛卡尔坐标系中的z轴坐标、利用等式(8)来确定监测点在笛卡尔坐标系中的z轴坐标:
Zpij=Zij+hpij (8)
其中,Zij表示与监测点对应的位置点在笛卡尔坐标系中的z轴坐标;hpij表示监测点距所述与监测点对应的位置点的高度(单位为m),并且,hpij∈(0,100],该hpij的值可以根据实际的监测需求进行设定;Zpij表示监测点在笛卡尔坐标系中的z轴坐标。
这样,就可以确定出监测点在笛卡尔坐标系中的坐标为(Xij,Yij,Zpij)。
在步骤S305中,根据所述泄漏源在所述全球坐标系中的坐标、以及第二全球坐标与笛卡尔坐标的转换函数,确定所述泄漏源在所述笛卡尔坐标系中的坐标。
在一个实施方式中,所述第二全球坐标与笛卡尔坐标的转换函数可以包括:
其中,X0、Y0、Z0表示所述泄漏源在所述笛卡尔坐标系中的坐标。
图4示出了根据本发明的另一实施方式的用于监测泄漏气体扩散浓度的方法的流程图。如图4所示,所述步骤S103可以包括:
在步骤S401中,根据所述泄漏源在所述地理坐标系中的坐标、泄漏发生的时刻、监测日期序数、以及太阳高度角计算函数,确定太阳高度角。
在一个实施方式中,太阳高度角计算函数可以包括:
θ=360dn/365
hsun=arcsin[sinLo0sinσ+cosLo0cosσcos(15T0)+La0-300]
其中,hsun表示所述太阳高度角(单位为度,°);dn表示所述监测日期序数(dn=0~364,其中,针对元月1日dn=0,以此类推);θ为中间变量;σ表示泄漏源所在地的太阳倾角;Lo0、La0表示所述泄漏源在所述地理坐标系中的x轴坐标和y轴坐标,即,泄漏源的经度和纬度;T0表示泄漏发生的时刻(单位为小时),例如,针对时刻15时30分,T0=15.5小时。
在步骤S402中,根据所述太阳高度角、所述总云量、所述低云量和所述监测时间,确定太阳辐射等级。
首先,可以根据地理区域的气象条件以及表1,确定出地理区域的天空的总云量和低云量。
表1
之后,可以根据在步骤S401中确定的太阳高度角hsun、根据表1确定的总云量、低云量、以及监测时间,按照表2,确定太阳辐射等级SL。
表2
表2中的“夜间”指日落后1小时至日出前1小时的时间范围。
例如,在总云量≤4并且低云量≤4的情况下,如果监测时间属于夜间,那么确定太阳辐射等级SL=-2。如果监测时间不属于夜间,那么如果太阳高度角hsun≤15°,则确定太阳辐射等级SL=-1。
接下来,在步骤S403中,根据所述地理区域的风速和所述太阳辐射等级,确定大气稳定度。
地理区域的风速U(单位为m/s)可以例如通过风速测定仪获得。在获得地理区域的风速U后,可以结合太阳辐射等级SL、按照表3确定出大气稳定度,其中,表3中的A、A~B、C、C~D、D、E、F表示大气稳定度。
表3
在步骤S404中,根据所述大气稳定度和所述监测点在所述笛卡尔坐标系中的坐标,确定所述监测点处的实际扩散系数。
在一个实施方式中,所述步骤S404可以包括以下步骤:
首先,根据所述大气稳定度、所述监测点在所述笛卡尔坐标系中的坐标、以及基础扩散系数计算函数,确定所述监测点处的基础扩散系数(单位为m,x/y/z分别代表x/y/z方向)。
在一个实施方式中,所述基础扩散系数计算函数可以包括:
其中,表示所述监测点处的基础扩散系数;A、A~B、B、B~C、C、C~D、D、D~E、E、E~F、F表示所述大气稳定度;Xij表示所述监测点在所述笛卡尔坐标系中的x轴坐标。
这样,在确定出大气稳定度后,可以首先根据等式(11)、(12)和(13)确定出监测点处的基础扩散系数
接下来,根据所述大气稳定度、所述地理区域上所有位置点在所述笛卡尔坐标系中的坐标、以及修正系数计算函数,确定修正系数。
在一个实施方式中,所述修正系数计算函数可以包括:
fy=1+a0difmax (14)
其中,fy表示所述修正系数;difmax表示所述地理区域上的所有位置点在所述笛卡尔坐标系中的z轴坐标与所述所有位置点在所述笛卡尔坐标系中的z轴坐标的平均值之间的差值中的最大值,其中,difmax>0;a0表示预定系数,并且a0与所述大气稳定度相对应。
在一个实施方式中,可以通过以下方式确定difmax:
其中,Zij表示地理区域上的任一位置点在笛卡尔坐标系中的z轴坐标;Z表示地理区域上的所有位置点在所述笛卡尔坐标系中的z轴坐标的平均值。
a0可以与所述大气稳定度相对应。例如,可以根据大气稳定度、按照以下等式(17)确定a0的值:
最后,可以根据所述修正系数fy、所述监测点处的基础扩散系数 以及基础扩散系数修正函数,得出所述监测点处的实际扩散系数
在一个实施方式中,所述基础扩散系数修正函数可以包括:
由此,就可以确定出监测点处的实际扩散系数(单位为m,x/y/z分别代表x/y/z方向)。
图5示出了根据本发明的另一实施方式的用于监测泄漏气体扩散浓度的方法的流程图。如图5所示,所述步骤S104可以包括:
在步骤S501中,根据所述监测点处的实际扩散系数、所述泄漏源在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标、所述泄漏源的泄漏强度、以及第一泄漏气体扩散浓度计算函数,确定所述泄漏源在所述监测点处产生的第一泄漏气体扩散子浓度。
在一个实施方式中,所述第一泄漏气体扩散浓度计算函数可以包括:
其中,crij表示所述第一泄漏气体扩散子浓度(单位为mg/m3);Xij、Yij、Zpij表示所述监测点在所述笛卡尔坐标系中的坐标;表示所述监测点处的实际扩散系数;Z0表示所述泄漏源在所述笛卡尔坐标系中的z轴坐标;Q表示所述泄漏源的泄漏强度(单位为kg/min)。
在步骤S502中,以与所述监测点对应的位置点所在的平原为镜面,确定所述泄漏源相对于所述镜面的镜像位置,并得出所述镜像位置在所述笛卡尔坐标系中的坐标。
例如,如图6所示,所述监测点P对应的位置点W所在的平原可以由MF表示。以该MF为镜面,泄漏源L相对于该镜面的镜像位置为点V所示,该点V可以称为泄漏源L的虚源点,该镜像位置在所述笛卡尔坐标系中的坐标可以为(0,0,Zvij),Zvij=2Zij-Z0,其中,Zvij表示泄漏源的镜像位置在笛卡尔坐标系中的z轴坐标,Zij表示与监测点对应的位置点在笛卡尔坐标系中的z轴坐标;Z0表示泄漏源在笛卡尔坐标系中的z轴坐标。
在步骤S503中,根据所述监测点处的实际扩散系数、所述镜像位置在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标、所述泄漏源的泄漏强度、以及第二泄漏气体扩散浓度计算函数,确定所述镜像位置在所述监测点处产生的第二泄漏气体扩散子浓度。
在一个实施方式中,所述第二泄漏气体扩散浓度计算函数可以包括:
其中,cvij表示所述第二泄漏气体扩散子浓度(单位为mg/m3)。
在步骤S504中,根据所述第一泄漏气体扩散子浓度crij、所述第二泄漏气体扩散子浓度cvij和泄漏气体扩散浓度计算函数,确定所述监测点处的泄漏气体扩散浓度。
在一个实施方式中,所述泄漏气体扩散浓度计算函数可以包括:
cij=crij+cvij (21)
其中,cij表示监测点处的泄漏气体扩散浓度(单位为mg/m3)。
这样,就可以监测出监测点处的泄漏气体扩散浓度。可以通过上述方法监测出任意期望的地理区域上方的任意期望的监测点处的泄漏气体扩散浓度,这样,可进行大范围泄漏气体扩散浓度的监测。
此外,通过上述用于监测泄漏气体扩散浓度的方法,能够快速确定出监测点处的泄漏气体扩散浓度,并快速确定出泄漏气体的污染区域,特别是在高起伏地形条件下,也能够准确地确定泄漏气体扩散浓度分布,从而为确定警戒与疏散范围、撤离路线以及救援资源分配等提供准确、可靠的数据支持。
图7是根据本发明的实施方式的用于监测泄漏气体扩散浓度的装置的框图。如图7所示,该装置可以包括:采集模块701,用于采集泄漏源和地理区域上与监测点对应的位置点在地理坐标系中的坐标,以及所述泄漏源的泄漏强度;笛卡尔坐标确定模块702,用于根据所述泄漏源和所述与监测点对应的位置点在所述地理坐标系中的坐标、以及所述地理区域的平均风向角,确定所述泄漏源和所述监测点在笛卡尔坐标系中的坐标,其中,所述笛卡尔坐标系为以所述泄漏源向地面的铅直投影点为原点,以所述地理区域的平均风向为x轴正向,以垂直所述泄漏源所在地的海平面、指向地球外部为z轴正向的右手笛卡尔坐标系;实际扩散系数确定模块703,用于根据所述泄漏源在所述地理坐标系中的坐标、所述地理区域的气象参数、监测时间、以及所述监测点在所述笛卡尔坐标系中的坐标,确定所述监测点处的实际扩散系数;泄漏气体扩散浓度确定模块704,用于根据所述监测点处的实际扩散系数、所述泄漏源在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标和所述泄漏源的泄漏强度,确定所述监测点处的泄漏气体扩散浓度。
图8是根据本发明的另一实施方式的用于监测泄漏气体扩散浓度的装置的框图。如图8所示,所述笛卡尔坐标确定模块702可以包括:第一全球坐标确定子模块801,用于根据所述与监测点对应的位置点在所述地理坐标系中的坐标、所述地理区域的平均海拔、以及第一地理坐标与全球坐标的转换函数,确定所述与监测点对应的位置点在全球坐标系中的坐标;第二全球坐标确定子模块802,用于根据所述泄漏源在所述地理坐标系中的坐标、所述地理区域的平均海拔、以及第二地理坐标与全球坐标的转换函数,确定所述泄漏源在全球坐标系中的坐标;第一笛卡尔坐标确定子模块803,用于根据所述与监测点对应的位置点在所述全球坐标系中的坐标、所述泄漏源在所述全球坐标系中的坐标、所述地理区域的平均风向角、以及第一全球坐标与笛卡尔坐标的转换函数,确定所述与监测点对应的位置点在所述笛卡尔坐标系中的坐标;第二笛卡尔坐标确定子模块804,用于根据所述与监测点对应的位置点在所述笛卡尔坐标系中的坐标,确定所述监测点在所述笛卡尔坐标系中的坐标;以及第三笛卡尔坐标确定子模块805,用于根据所述泄漏源在所述全球坐标系中的坐标、以及第二全球坐标与笛卡尔坐标的转换函数,确定所述泄漏源在所述笛卡尔坐标系中的坐标。
图9是根据本发明的另一实施方式的用于监测泄漏气体扩散浓度的装置的框图。在图9所示的实施方式中,所述地理区域的气象参数包括总云量、低云量和风速;以及,所述实际扩散系数确定模块703可以包括:太阳高度角确定子模块901,用于根据所述泄漏源在所述地理坐标系中的坐标、泄漏发生的时刻、监测日期序数、以及太阳高度角计算函数,确定太阳高度角;太阳辐射等级确定子模块902,用于根据所述太阳高度角、所述总云量、所述低云量和所述监测时间,确定太阳辐射等级;大气稳定度确定子模块903,用于根据所述风速和所述太阳辐射等级,确定大气稳定度;实际扩散系数确定子模块904,用于根据所述大气稳定度和所述监测点在所述笛卡尔坐标系中的坐标,确定所述监测点处的实际扩散系数。
其中,所述实际扩散系数确定子模块904可以包括:基础扩散系数确定单元,用于根据所述大气稳定度、所述监测点在所述笛卡尔坐标系中的坐标、以及基础扩散系数计算函数,确定所述监测点处的基础扩散系数;修正系数确定单元,用于根据所述大气稳定度、所述地理区域上所有位置点在所述笛卡尔坐标系中的坐标、以及修正系数计算函数,确定修正系数;实际扩散系数确定单元,用于根据所述修正系数、所述监测点处的基础扩散系数、以及基础扩散系数修正函数,得出所述监测点处的实际扩散系数。
图10是根据本发明的另一实施方式的用于监测泄漏气体扩散浓度的装置的框图。如图10所示,所述泄漏气体扩散浓度确定模块704可以包括:第一泄漏气体扩散子浓度确定子模块1001,用于根据所述监测点处的实际扩散系数、所述泄漏源在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标、所述泄漏源的泄漏强度、以及第一泄漏气体扩散浓度计算函数,确定所述泄漏源在所述监测点处产生的第一泄漏气体扩散子浓度;镜像位置确定子模块1002,用于以与所述监测点对应的位置点所在的平原为镜面,确定所述泄漏源相对于所述镜面的镜像位置,并得出所述镜像位置在所述笛卡尔坐标系中的坐标;第二泄漏气体扩散子浓度确定子模块1003,用于根据所述监测点处的实际扩散系数、所述镜像位置在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标、所述泄漏源的泄漏强度、以及第二泄漏气体扩散浓度计算函数,确定所述镜像位置在所述监测点处产生的第二泄漏气体扩散子浓度;泄漏气体扩散浓度确定子模块1004,用于根据所述第一泄漏气体扩散子浓度、所述第二泄漏气体扩散子浓度和泄漏气体扩散浓度计算函数,确定所述监测点处的泄漏气体扩散浓度。
本发明的用于监测泄漏气体扩散浓度的装置对应于用于监测泄漏气体扩散浓度的方法,因此相同的内容不再赘述。
综上所述,通过本发明的用于监测泄漏气体扩散浓度的方法及装置,可进行大范围泄漏气体扩散浓度的监测,适用性强,能够与地理信息***相结合。通过获取泄漏源相关参数以及气象参数,能够快速确定出监测点处的泄漏气体扩散浓度,并快速确定出泄漏气体的污染区域,特别是在高起伏地形条件下,也能够准确地确定泄漏气体扩散浓度分布,从而为确定警戒与疏散范围、撤离路线以及救援资源分配等提供准确、可靠的数据支持。
以上结合附图详细描述了本发明的优选实施方式,但是,本发明并不限于上述实施方式中的具体细节,在本发明的技术构思范围内,可以对本发明的技术方案进行多种简单变型,这些简单变型均属于本发明的保护范围。
此外,尽管在附图中以特定顺序描述了本发明方法的操作,但是,这并非要求或者暗示必须按照该特定顺序来执行这些操作,或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地,可以省略某些步骤,将多个步骤合并为一个步骤执行,和/或将一个步骤分解为多个步骤执行。
本领域技术人员知道,本发明的实施方式可以实现为一种***、装置、设备、方法或计算机程序产品。因此,本发明的实施方式可以具体实现为以下形式,即:完全的硬件、完全的软件(包括固件、驻留软件、微代码等),或者硬件和软件结合的形式。
另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合。为了避免不必要的重复,本发明对各种可能的组合方式不再另行说明。
此外,本发明的各种不同的实施方式之间也可以进行任意组合,只要其不违背本发明的思想,其同样应当视为本发明所公开的内容。
Claims (9)
1.一种用于监测泄漏气体扩散浓度的方法,其特征在于,该方法包括:
采集泄漏源和地理区域上与监测点对应的位置点在地理坐标系中的坐标,以及所述泄漏源的泄漏强度;
根据所述泄漏源和所述与监测点对应的位置点在所述地理坐标系中的坐标、以及所述地理区域的平均风向角,确定所述泄漏源和所述监测点在笛卡尔坐标系中的坐标,其中,所述笛卡尔坐标系为以所述泄漏源向地面的铅直投影点为原点,以所述地理区域的平均风向为x轴正向,以垂直所述泄漏源所在地的海平面、指向地球外部为z轴正向的右手笛卡尔坐标系;
根据所述泄漏源在所述地理坐标系中的坐标、所述地理区域的气象参数、监测时间、以及所述监测点在所述笛卡尔坐标系中的坐标,确定所述监测点处的实际扩散系数;
根据所述监测点处的实际扩散系数、所述泄漏源在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标、所述泄漏源的泄漏强度、以及第一泄漏气体扩散浓度计算函数,确定所述泄漏源在所述监测点处产生的第一泄漏气体扩散子浓度;
以与所述监测点对应的位置点所在的平原为镜面,确定所述泄漏源相对于所述镜面的镜像位置,并得出所述镜像位置在所述笛卡尔坐标系中的坐标;
根据所述监测点处的实际扩散系数、所述镜像位置在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标、所述泄漏源的泄漏强度、以及第二泄漏气体扩散浓度计算函数,确定所述镜像位置在所述监测点处产生的第二泄漏气体扩散子浓度;
根据所述第一泄漏气体扩散子浓度、所述第二泄漏气体扩散子浓度和泄漏气体扩散浓度计算函数,确定所述监测点处的泄漏气体扩散浓度。
2.根据权利要求1所述的方法,其特征在于,所述根据所述泄漏源和所述与监测点对应的位置点在所述地理坐标系中的坐标、以及所述地理区域的平均风向角,确定所述泄漏源和所述监测点在笛卡尔坐标系中的坐标包括:
根据所述与监测点对应的位置点在所述地理坐标系中的坐标、所述地理区域的平均海拔、以及第一地理坐标与全球坐标的转换函数,确定所述与监测点对应的位置点在全球坐标系中的坐标;
根据所述泄漏源在所述地理坐标系中的坐标、所述地理区域的平均海拔、以及第二地理坐标与全球坐标的转换函数,确定所述泄漏源在全球坐标系中的坐标;
根据所述与监测点对应的位置点在所述全球坐标系中的坐标、所述泄漏源在所述全球坐标系中的坐标、所述地理区域的平均风向角、以及第一全球坐标与笛卡尔坐标的转换函数,确定所述与监测点对应的位置点在所述笛卡尔坐标系中的坐标;
根据所述与监测点对应的位置点在所述笛卡尔坐标系中的坐标,确定所述监测点在所述笛卡尔坐标系中的坐标;以及
根据所述泄漏源在所述全球坐标系中的坐标、以及第二全球坐标与笛卡尔坐标的转换函数,确定所述泄漏源在所述笛卡尔坐标系中的坐标。
3.根据权利要求2所述的方法,其特征在于,所述第一地理坐标与全球坐标的转换函数包括:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>LX</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>g</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;Lo</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>/</mo>
<mn>180</mn>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;La</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>/</mo>
<mn>180</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>LY</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>g</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;Lo</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>/</mo>
<mn>180</mn>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;La</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>/</mo>
<mn>180</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>LZ</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,LXij、LYij、LZij表示所述与监测点对应的位置点在所述全球坐标系中的坐标;Havg表示所述平均海拔;Loij、Laij、Hij表示所述与监测点对应的位置点在所述地理坐标系中的坐标;Ngl表示地球的卯酉圈曲率半径,并且agl表示地球椭球长半轴,Wgl为中间变量,并且e表示地球椭球的第一偏心率,bgl表示地球椭球的短半轴;
所述第二地理坐标与全球坐标的转换函数包括:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>LX</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>g</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;Lo</mi>
<mn>0</mn>
</msub>
<mo>/</mo>
<mn>180</mn>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;La</mi>
<mn>0</mn>
</msub>
<mo>/</mo>
<mn>180</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>LY</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>g</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;Lo</mi>
<mn>0</mn>
</msub>
<mo>/</mo>
<mn>180</mn>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;La</mi>
<mn>0</mn>
</msub>
<mo>/</mo>
<mn>180</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>LZ</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>H</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,LX0、LY0、LZ0表示所述泄漏源在所述全球坐标系中的坐标;Lo0、La0、H0表示所述泄漏源在所述地理坐标系中的坐标;
所述第一全球坐标与笛卡尔坐标的转换函数包括:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>LX</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>&alpha;</mi>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>LY</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>&alpha;</mi>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>LX</mi>
<mn>0</mn>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>&alpha;</mi>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>LY</mi>
<mn>0</mn>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>&alpha;</mi>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>LY</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>&alpha;</mi>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>LX</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>&alpha;</mi>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>LY</mi>
<mn>0</mn>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>LX</mi>
<mn>0</mn>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>&alpha;</mi>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>H</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mn>00</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,Xij、Yij、Zij表示所述与监测点对应的位置点在所述笛卡尔坐标系中的坐标;θ0表示所述全球坐标系的x轴正向与正东方向的夹角;H00表示所述泄漏源所在位置地面海拔高度,α表示所述地理区域的平均风向角;
所述第二全球坐标与笛卡尔坐标的转换函数包括:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>LZ</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mn>00</mn>
</msub>
<mo>=</mo>
<msub>
<mi>H</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mn>00</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,X0、Y0、Z0表示所述泄漏源在所述笛卡尔坐标系中的坐标。
4.根据权利要求1所述的方法,其特征在于,所述地理区域的气象参数包括总云量、低云量和风速;以及
所述根据所述泄漏源在所述地理坐标系中的坐标、所述地理区域的气象参数、监测时间、以及所述监测点在所述笛卡尔坐标系中的坐标,确定所述监测点处的实际扩散系数包括:
根据所述泄漏源在所述地理坐标系中的坐标、泄漏发生的时刻、监测日期序数、以及太阳高度角计算函数,确定太阳高度角;
根据所述太阳高度角、所述总云量、所述低云量和所述监测时间,确定太阳辐射等级;
根据所述风速和所述太阳辐射等级,确定大气稳定度;
根据所述大气稳定度和所述监测点在所述笛卡尔坐标系中的坐标,确定所述监测点处的实际扩散系数。
5.根据权利要求4所述的方法,其特征在于,所述太阳高度角计算函数包括:
θ=360dn/365
<mrow>
<mi>&sigma;</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0.006918</mn>
<mo>-</mo>
<mn>0.39912</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mn>0.070257</mn>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>-</mo>
<mn>0.006758</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mn>2</mn>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>0.000907</mn>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mn>2</mn>
<mi>&theta;</mi>
<mo>-</mo>
<mn>0.002697</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mn>3</mn>
<mi>&theta;</mi>
<mo>+</mo>
<mn>0.00148</mn>
<mi>sin</mi>
<mn>3</mn>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mn>180</mn>
<mo>/</mo>
<mi>&pi;</mi>
</mrow>
hsun=arcsin[sinLo0sinσ+cosLo0cosσcos(15T0)+La0-300]
其中,hsun表示所述太阳高度角;dn表示所述监测日期序数;θ为中间变量;σ表示泄漏源所在地的太阳倾角;Lo0、La0表示所述泄漏源在所述地理坐标系中的x轴坐标和y轴坐标;T0表示所述泄漏发生的时刻。
6.根据权利要求4所述的方法,其特征在于,所述根据所述大气稳定度和所述监测点在所述笛卡尔坐标系中的坐标,确定所述监测点处的实际扩散系数包括:
根据所述大气稳定度、所述监测点在所述笛卡尔坐标系中的坐标、以及基础扩散系数计算函数,确定所述监测点处的基础扩散系数;
根据所述大气稳定度、所述地理区域上所有位置点在所述笛卡尔坐标系中的坐标、以及修正系数计算函数,确定修正系数;
根据所述修正系数、所述监测点处的基础扩散系数、以及基础扩散系数修正函数,得出所述监测点处的实际扩散系数。
7.根据权利要求6所述的方法,其特征在于,所述基础扩散系数计算函数包括:
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>b</mi>
<mi>x</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>b</mi>
<mi>y</mi>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>b</mi>
<mi>y</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0.22</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mi>A</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mn>0.22</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>0.16</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>A</mi>
<mo>~</mo>
<mi>B</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.16</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mi>B</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mn>0.16</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>0.11</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>B</mi>
<mo>~</mo>
<mi>C</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.11</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mi>C</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mn>0.11</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>0.08</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>C</mi>
<mo>~</mo>
<mi>D</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.08</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mi>D</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mn>0.08</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>0.06</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>D</mi>
<mo>~</mo>
<mi>E</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.06</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mi>E</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mn>0.06</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>0.04</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>E</mi>
<mo>~</mo>
<mi>F</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.04</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.001</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mi>F</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>b</mi>
<mi>z</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0.20</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mi>A</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mn>0.20</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>0.12</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>A</mi>
<mo>~</mo>
<mi>B</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.12</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mi>B</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mn>0.12</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>+</mo>
<mn>0.08</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.002</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>B</mi>
<mo>~</mo>
<mi>C</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.08</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.002</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mi>C</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mn>0.08</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.002</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>0.06</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.0015</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>C</mi>
<mo>~</mo>
<mi>D</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.06</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.0015</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mi>D</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mn>0.06</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.0015</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>0.03</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.0003</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>D</mi>
<mo>~</mo>
<mi>E</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.03</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.0003</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mi>E</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mn>0.03</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.0003</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>+</mo>
<mn>0.016</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.0003</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>E</mi>
<mo>~</mo>
<mi>F</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.016</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>0.0003</mn>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mi>F</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,表示所述监测点处的基础扩散系数;A、A~B、B、B~C、C、C~D、D、D~E、E、E~F、F表示所述大气稳定度;Xij表示所述监测点在所述笛卡尔坐标系中的x轴坐标;
所述修正系数计算函数包括:
fy=1+a0difmax
其中,fy表示所述修正系数;difmax表示所述地理区域上的所有位置点在所述笛卡尔坐标系中的z轴坐标与所述所有位置点在所述笛卡尔坐标系中的z轴坐标的平均值之间的差值中的最大值;a0表示预定系数,并且a0与所述大气稳定度相对应;
所述基础扩散系数修正函数包括:
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>x</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>y</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>y</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>b</mi>
<mi>y</mi>
</mrow>
</msubsup>
<mo>*</mo>
<msub>
<mi>f</mi>
<mi>y</mi>
</msub>
</mrow>
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>z</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mi>z</mi>
<mi>y</mi>
</mrow>
</msubsup>
<mo>*</mo>
<msub>
<mi>f</mi>
<mi>y</mi>
</msub>
</mrow>
其中,表示所述监测点处的实际扩散系数。
8.根据权利要求1所述的方法,其特征在于,所述第一泄漏气体扩散浓度计算函数包括:
<mrow>
<msub>
<mi>cr</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mi>Q</mi>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mn>3</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<mn>2</mn>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>x</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<mn>2</mn>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>y</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>Zp</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>z</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
其中,crij表示所述第一泄漏气体扩散子浓度;Xij、Yij、Zpij表示所述监测点在所述笛卡尔坐标系中的坐标;表示所述监测点处的实际扩散系数;Z0表示所述泄漏源在所述笛卡尔坐标系中的z轴坐标;Q表示所述泄漏源的泄漏强度;
所述第二泄漏气体扩散浓度计算函数包括:
<mrow>
<msub>
<mi>cv</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mi>Q</mi>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mn>3</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<mn>2</mn>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>x</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<mn>2</mn>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>y</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>Zp</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Zv</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>z</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
其中,cvij表示所述第二泄漏气体扩散子浓度;Zvij表示所述镜像位置在所述笛卡尔坐标系中的z轴坐标;
所述泄漏气体扩散浓度计算函数包括:
cij=crij+cvij
其中,cij表示所述泄漏气体扩散浓度。
9.一种用于监测泄漏气体扩散浓度的装置,其特征在于,该装置包括:
采集模块,用于采集泄漏源和地理区域上与监测点对应的位置点在地理坐标系中的坐标,以及所述泄漏源的泄漏强度;
笛卡尔坐标确定模块,用于根据所述泄漏源和所述与监测点对应的位置点在所述地理坐标系中的坐标、以及所述地理区域的平均风向角,确定所述泄漏源和所述监测点在笛卡尔坐标系中的坐标,其中,所述笛卡尔坐标系为以所述泄漏源向地面的铅直投影点为原点,以所述地理区域的平均风向为x轴正向,以垂直所述泄漏源所在地的海平面、指向地球外部为z轴正向的右手笛卡尔坐标系;
实际扩散系数确定模块,用于根据所述泄漏源在所述地理坐标系中的坐标、所述地理区域的气象参数、监测时间、以及所述监测点在所述笛卡尔坐标系中的坐标,确定所述监测点处的实际扩散系数;
泄漏气体扩散浓度确定模块,用于根据所述监测点处的实际扩散系数、所述泄漏源在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标和所述泄漏源的泄漏强度,确定所述监测点处的泄漏气体扩散浓度;其中,所述泄漏气体扩散浓度确定模块包括:
第一泄漏气体扩散子浓度确定子模块,用于根据所述监测点处的实际扩散系数、所述泄漏源在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标、所述泄漏源的泄漏强度、以及第一泄漏气体扩散浓度计算函数,确定所述泄漏源在所述监测点处产生的第一泄漏气体扩散子浓度;
镜像位置确定子模块,用于以与所述监测点对应的位置点所在的平原为镜面,确定所述泄漏源相对于所述镜面的镜像位置,并得出所述镜像位置在所述笛卡尔坐标系中的坐标;
第二泄漏气体扩散子浓度确定子模块,用于根据所述监测点处的实际扩散系数、所述镜像位置在所述笛卡尔坐标系中的坐标、所述监测点在所述笛卡尔坐标系中的坐标、所述泄漏源的泄漏强度、以及第二泄漏气体扩散浓度计算函数,确定所述镜像位置在所述监测点处产生的第二泄漏气体扩散子浓度;
泄漏气体扩散浓度确定子模块,用于根据所述第一泄漏气体扩散子浓度、所述第二泄漏气体扩散子浓度和泄漏气体扩散浓度计算函数,确定所述监测点处的泄漏气体扩散浓度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510456988.7A CN105115859B (zh) | 2015-07-29 | 2015-07-29 | 用于监测泄漏气体扩散浓度的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510456988.7A CN105115859B (zh) | 2015-07-29 | 2015-07-29 | 用于监测泄漏气体扩散浓度的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105115859A CN105115859A (zh) | 2015-12-02 |
CN105115859B true CN105115859B (zh) | 2018-01-12 |
Family
ID=54663899
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510456988.7A Active CN105115859B (zh) | 2015-07-29 | 2015-07-29 | 用于监测泄漏气体扩散浓度的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105115859B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106595757B (zh) * | 2016-11-29 | 2019-04-19 | 西南石油大学 | 一种环境监测方法及*** |
CN107589223B (zh) * | 2017-09-15 | 2020-04-07 | 深圳市卡普瑞环境科技有限公司 | 一种定位的方法、装置以及存储介质 |
CN107526905A (zh) * | 2017-10-16 | 2017-12-29 | 哈尔滨理工大学 | 基于fluent的氯气泄漏绕障碍物扩散的研究方法 |
CN109030289B (zh) * | 2018-05-18 | 2020-06-16 | 北京理工大学 | 一种燃气泄漏扩散范围预测方法 |
CN108760580A (zh) * | 2018-05-21 | 2018-11-06 | 众安仕(北京)科技有限公司 | 一种动态环境监测的气体扩散推测***及方法 |
CN109342658A (zh) * | 2018-10-30 | 2019-02-15 | 河南省环境保护科学研究院 | 基于wrf模型的计算机程序化计算大气环境容量的方法 |
CN110332465A (zh) * | 2019-06-27 | 2019-10-15 | 中石化川气东送天然气管道有限公司 | 一种长输天然气管道泄漏监测决策方法和*** |
CN112016738B (zh) * | 2020-08-13 | 2023-08-18 | 中国人民解放军空军军医大学 | 一种吸入性有毒物质的危害预测方法 |
CN113758891B (zh) * | 2021-08-18 | 2022-06-03 | 清华大学 | 混合气体组分浓度计算方法、装置、设备及存储介质 |
CN113562363B (zh) * | 2021-08-20 | 2022-07-08 | 北京云迹科技股份有限公司 | 一种清理垃圾桶的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101882184A (zh) * | 2010-05-25 | 2010-11-10 | 中冶赛迪工程技术股份有限公司 | 基于gis技术和aermode模型的大气环评***与环评方法 |
US7983844B2 (en) * | 2009-12-18 | 2011-07-19 | General Electric Company | Device and method for determining wind conditions using multiple wind resource grids |
JP5189897B2 (ja) * | 2008-05-29 | 2013-04-24 | 株式会社日立製作所 | パイプライン危険地域監視システム |
CN103413062A (zh) * | 2013-08-29 | 2013-11-27 | 中国测绘科学研究院 | 一种放射性核素扩散的计算方法 |
CN104657573A (zh) * | 2014-01-06 | 2015-05-27 | 北京千安哲信息技术有限公司 | 用于三维空间的泄漏气体扩散预测方法 |
-
2015
- 2015-07-29 CN CN201510456988.7A patent/CN105115859B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5189897B2 (ja) * | 2008-05-29 | 2013-04-24 | 株式会社日立製作所 | パイプライン危険地域監視システム |
US7983844B2 (en) * | 2009-12-18 | 2011-07-19 | General Electric Company | Device and method for determining wind conditions using multiple wind resource grids |
CN101882184A (zh) * | 2010-05-25 | 2010-11-10 | 中冶赛迪工程技术股份有限公司 | 基于gis技术和aermode模型的大气环评***与环评方法 |
CN103413062A (zh) * | 2013-08-29 | 2013-11-27 | 中国测绘科学研究院 | 一种放射性核素扩散的计算方法 |
CN104657573A (zh) * | 2014-01-06 | 2015-05-27 | 北京千安哲信息技术有限公司 | 用于三维空间的泄漏气体扩散预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105115859A (zh) | 2015-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105115859B (zh) | 用于监测泄漏气体扩散浓度的方法及装置 | |
CN105372388B (zh) | 用于有限时间稳态泄漏的气体浓度监测方法及装置 | |
CN107316012B (zh) | 小型无人直升机的火灾检测与跟踪方法 | |
CN108932394A (zh) | 有毒气体泄漏危害范围的确定方法 | |
KR102443718B1 (ko) | 환경시설의 대기오염가스 확산을 모니터링하기 위한 가상물리-기반 3차원 확산 예측지도 작성 시스템 및 그 방법 | |
CN104834017A (zh) | 一种有毒有害气体泄漏扩散事故源定位方法 | |
CN108760580A (zh) | 一种动态环境监测的气体扩散推测***及方法 | |
CN108061572B (zh) | 一种海洋核动力平台综合态势显控***及方法 | |
CN103984310A (zh) | 基于多源遥感数据的化工园区环境污染检测方法和装置 | |
CN109660290A (zh) | 一种基于多卫星多频段的多架构应急保障*** | |
CN107367773A (zh) | 一种基于北斗定位的单点气象要素确定方法 | |
CN108680936A (zh) | 北斗卫星导航***完好性监测评估方法 | |
CN110044338A (zh) | 一种溃堤溃坝场景的无人机监测方法及*** | |
CN114091756A (zh) | 一种基于泰森多边形的乡镇级海啸风险评估方法 | |
CN108181618A (zh) | 一种雷达标定方法 | |
CN113625731A (zh) | 一种基于dem数据的无人机地形匹配仿地飞行方法 | |
CN106918317A (zh) | 一种扫描式卫星对地侧摆观测覆盖带确定方法及装置 | |
CN114374716A (zh) | 地质灾害远程监测***及其监测方法 | |
JP5181379B1 (ja) | 発電ポテンシャル評価装置、及び発電ポテンシャル評価プログラム | |
Xiu et al. | Plane accuracy of aerophotogrammetry without image control points of DJI PHANTOM 4 RTK UAV | |
Xiaochun | Development of Beidou navigation satellite system | |
CN201623731U (zh) | 集成化林火动态监测与精确定位*** | |
CN114413850A (zh) | 一种基于SRTM 90m DEM的单点火点定位方法 | |
Khalaf et al. | Accuracy Assessment of World View-2 Satellite Imagery for Planimetric Maps Production | |
CN113611159B (zh) | 一种障碍物告警方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |