CN108121892A - 增加甲烷化学动态的中尺度大气光化学污染模拟预测算法 - Google Patents
增加甲烷化学动态的中尺度大气光化学污染模拟预测算法 Download PDFInfo
- Publication number
- CN108121892A CN108121892A CN201711085469.XA CN201711085469A CN108121892A CN 108121892 A CN108121892 A CN 108121892A CN 201711085469 A CN201711085469 A CN 201711085469A CN 108121892 A CN108121892 A CN 108121892A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- chemical
- formula
- methane
- 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.)
- Pending
Links
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 title claims abstract description 63
- 239000000126 substance Substances 0.000 title claims abstract description 49
- 238000004088 simulation Methods 0.000 title claims abstract description 14
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 8
- 238000006243 chemical reaction Methods 0.000 claims abstract description 43
- 238000004062 sedimentation Methods 0.000 claims abstract description 13
- 239000013626 chemical specie Substances 0.000 claims abstract description 7
- 239000006185 dispersion Substances 0.000 claims abstract description 5
- 239000000926 atmospheric chemistry Substances 0.000 claims abstract description 3
- 150000003254 radicals Chemical class 0.000 claims description 25
- 230000008859 change Effects 0.000 claims description 19
- 238000009792 diffusion process Methods 0.000 claims description 19
- 238000000034 method Methods 0.000 claims description 15
- 230000007246 mechanism Effects 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000006866 deterioration Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000010521 absorption reaction Methods 0.000 claims description 2
- 230000004913 activation Effects 0.000 claims description 2
- 239000012080 ambient air Substances 0.000 claims description 2
- 238000010574 gas phase reaction Methods 0.000 claims description 2
- 238000006303 photolysis reaction Methods 0.000 claims description 2
- 238000006862 quantum yield reaction Methods 0.000 claims description 2
- 238000009826 distribution Methods 0.000 abstract description 11
- 230000002123 temporal effect Effects 0.000 abstract description 3
- 239000003570 air Substances 0.000 description 13
- 239000007789 gas Substances 0.000 description 9
- 241000894007 species Species 0.000 description 6
- 238000000469 dry deposition Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 239000005436 troposphere Substances 0.000 description 5
- 239000003344 environmental pollutant Substances 0.000 description 4
- 230000003647 oxidation Effects 0.000 description 4
- 238000007254 oxidation reaction Methods 0.000 description 4
- 231100000719 pollutant Toxicity 0.000 description 4
- WSFSSNUMVMOOMR-UHFFFAOYSA-N Formaldehyde Chemical compound O=C WSFSSNUMVMOOMR-UHFFFAOYSA-N 0.000 description 3
- 238000006552 photochemical reaction Methods 0.000 description 3
- WTFNSXYULBQCQV-UHFFFAOYSA-N $l^{1}-oxidanyloxymethane Chemical compound CO[O] WTFNSXYULBQCQV-UHFFFAOYSA-N 0.000 description 2
- 241000209094 Oryza Species 0.000 description 2
- 235000007164 Oryza sativa Nutrition 0.000 description 2
- 238000001311 chemical methods and process Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- GRVDJDISBSALJP-UHFFFAOYSA-N methyloxidanyl Chemical compound [O]C GRVDJDISBSALJP-UHFFFAOYSA-N 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 239000002243 precursor Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 235000009566 rice Nutrition 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 description 1
- 238000003915 air pollution Methods 0.000 description 1
- 239000012491 analyte Substances 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000029087 digestion Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 150000002894 organic compounds Chemical class 0.000 description 1
- 239000007800 oxidant agent Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 239000013618 particulate matter Substances 0.000 description 1
- 230000005180 public health Effects 0.000 description 1
- 239000000376 reactant Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 235000013311 vegetables Nutrition 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2219/00—Indexing scheme relating to application aspects of data processing equipment or methods
- G06F2219/10—Environmental application, e.g. waste reduction, pollution control, compliance with environmental legislation
Landscapes
- Chemical & Material Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
本发明公开了一种增加甲烷化学变化的中尺度大气光化学污染模拟预测算法，其基于CALGRID化学模式，考虑大气化学反应、大气输送与扩散、沉降、地面面源和高架排放源的影响，并对化学物种原始浓度中物理量的平均量分解为脉动量，在该基础上，将OH自由基和CH4作为变量，引入OH自由基和甲烷(CH4)的化学反应，根据OH自由基和CH4的化学反应列表计算OH自由基和CH4的浓度变化。本发明利用空气质量模式模拟OH自由基和CH4的时空分布，不仅有利于了解模拟区域的大气氧化能力，而且有利于深入了解光化学反应的细节。
Description
技术领域
本发明属于大气环境污染物检测技术领域，尤其涉及中尺度大气光化学污染的改进型预测模型。
背景技术
空气质量模式，是在对污染物排入大气环境后传输、扩散、转化和清除等一系列物理和化学过程的认识基础上，利用气象、环境、物理、化学等学科的研究方法和计算机技术，实现模拟和预报不同空间尺度上空气污染物浓度分布状况及变化趋势的方法，在空气质量预报、大气污染控制、环境规划与管理、城市建设及公共卫生等方面均有重要的实际应用价值，具有广阔的发展前景。
CALGRID是由美国加州的ARB(Air Resources Board)开发的，是欧拉型的中尺度大气光化学模式，针对二次污染源如臭氧的模拟具有较好效果，主要适用于晴空条件下的光化学反应的模拟，包含了大气输送与扩散、气相化学反应、人为排放的点面线源、干沉降等过程。
原有的CALGRID模式还是存在一定的局限性，包括假设甲烷为浓度不变物种，考虑甲烷的物理化学过程不完善。在清洁地区，甲烷排放特别是其自然源排放，对OH自由基和气相化学反应有重要的影响。因此，CALGRID在实际应用方面受到一定的限制，且未考虑自然源排放中的甲烷排放对OH自由基的影响，可能造成最终的模拟失真或偏离。
发明内容
为解决上述技术问题，本发明采用了以下技术方案：一种增加甲烷化学动态的中尺度大气光化学污染模拟预测算法，包括以下步骤：
步骤1：采用CALGRID化学模式，考虑大气化学反应、大气输送与扩散、沉降、底面积高架排放源的影响，并对化学物种原始浓度中物理量的平均量分解为脉动量，得到化学物种浓度变化方程如式(1)，
式中，C是化学物种平均浓度，V是三维风矢量的平均量，K是湍流扩散系数，E是污染源排放，是由于沉降引起的物种浓度变化，PCHEM是化学产生率，LCHEM是化学损失率；式(1)中二阶湍流扩散项通过湍流扩散系数K理论闭合转化得到；式(1)中右侧多项式依次分别为平流项、扩散项、源项、沉降项、化学变化项；
步骤2：由于平流项、扩散项、源项、沉降项、化学变化项各过程的特征时间不一致，通过算子分列时间积分对各过程进行积分，
Cn+1＝AxAyAzAcAzAyAxCn (2)
式中，Ax、Ay是水平输送扩散算子；Az是垂直输送扩散、物质输入和物理损耗算子，Ac是化学反应算子；
步骤3：CALGRID化学模式中的化学变化项以SAPRC-90的气相化学机制为基础，并在该基础上引入OH自由基和CH4的化学反应，根据OH自由基和CH4的化学反应列表计算OH自由基和CH4的浓度变化。
进一步改进，步骤3中所述OH自由基和CH4的化学反应列表如表1所示：
表1改进化学模块涉及的与OH、CH4有关的化学反应
进一步的，通过下式计算甲烷CH4的去除速率，
式中，是甲烷的去除速率；[OH]是OH自由基的浓度，R(27)根据气象条件和化学机制中参数计算得到，[OH]可根据表1其它化学反应方程联立并作拟稳态假设求得；CH4的浓度随时间变化则由下式计算：
式中，为甲烷的浓度。
进一步的，步骤3中所述CALGRID化学模式中的化学变化项以SAPRC-90的气相化学机制的化学物质如表2：
其中，光分解反应的光解速率系数通过以下计算表达式求解：
式中，σ[λ,T(h)]是同波长λ和温度T(h)有关的吸收截面系数；φ[λ,T(h)]是量子产额；I[λ,N(h),X]是同波长、大气状态N(h)和空间位置X有关的入射光强；
对于化学反应速率系数随温度升高而变大的双分子气相反应，则用Arrhenius(阿伦尼乌斯)经验公式计算速率系数：
式中，A是前置系数、T是环境气温(K)、Tr是参考气温(300K)、Ea是活化能、B是无量纲数。对于CALGRID模式中每种化学反应的A、B、Ea、σ[λ,T(h)]、φ[λ,T(h)]均由气相化学机制(如SAPRC-90)给出。
图说明
图1为本发明实施例中获得的中国地区近地面CH4平均浓度分布示意图；
图2为本发明实施例中获得的中国地区近地面OH自由基平均浓度分布示意图。
具体实施方式
下面结合附图并以具体实施例，进一步阐明本发明。应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
原CALGRID模式中化学变化项PCHEM-LCHEM是可以采用的SAPRC-90的气相化学机制，而该气相化学反应机制将CH4项作为定常数处理，并未考虑OH自由基与CH4的化学反应变化以及对其他物质的影响。但是CH4是大气中浓度最高的有机痕量气体，同时也是对流层中浓度最高的化学活性气体，其全球平均浓度约为1.803ppm。有关CH4物理和化学特性的研究一直很受重视，一方面它是重要的温室气体，对红外辐射的吸收能力比CO2强15～30倍，因此对地气系统能量收支以及气候有重要影响；另一方面，CH4在大气中易被OH自由基氧化，产生一系列有机化合物，影响对流层OH、CH2O、CO和O3等的浓度，从而影响许多大气成分的化学转化：CH4+OH→CH3+H2O(R1)CH3+O2+M→CH3O2+M(R2)CH3O2+NO→CH3O+NO2(R3)CH3O+O2→CH2O+HO2(R4)CH2O+hυ→H+CHO(R5)H+O2+M→HO2+M(R6)CHO+O2→HO2+CO(R7)HO2+NO→OH+NO2(R8)
作为大气中O3的前体物，虽然CH4的反应活性低于O3的另两种重要前体物NMHC和NOx(朱彬等，2000)，但是由于在大气中寿命较长，它们影响的范围更广，在清洁地区以及区域和全球尺度光化学反应中有重要作用，影响着对流层的氧化能力(王明星等，1999；秦瑜等，2003)。大气CH4很大一部分来自稻田的厌氧分解，研究表明稻田是中国地区最主要的CH4排放源之一(王明星等，2001)；最新研究表明，有氧环境下陆地生态植被也会排放可观的CH4(Keppler等，2006)。鉴于上述原因，利用空气质量模式模拟甲烷自然源排放对中国地区对流层化学的影响，不仅有利于了解我国大气甲烷的空间分布，而且有助于深刻理解甲烷自然源在光化学反应中的重要性。然而在CALGRID模式中，CH4被设定为定常物种，即认为CH4是充足的，其浓度不随化学反应而改变。这种设定在其它一些模式中也有采用，如WRF－CHEM。但是对于本研究而言，由于希望得到甲烷的空间分布以及由甲烷自然源引起的光化学反应的差异，因此有必要对化学模块进行一定的改进，将定常的CH4改变为时变的。
另外，OH自由基在光化学反应中也占有非常重要的地位，它控制着众多微量气体和自由基的氧化和清除过程，是大气中最重要的氧化剂(贾龙等，2006)。从OH自由基生成和反应过程，可见对流层白昼化学是以OH自由基为中心的，高活性的OH自由基使大气微量成分氧化并经链反应实现物种化学转化，其存在对S、N、奇氢、奇氧的循环有重要意义。因此利用空气质量模式模拟OH自由基的时空分布，不仅有利于了解模拟区域的大气氧化能力，而且有利于深入了解光化学反应的细节。在CALGRID模式中，OH自由基被设定为拟稳态物种，但是它的浓度并没有作为模拟结果输出。鉴于本文对对流层氧化特性研究的需要，因此在模式中增加通道将OH浓度输出。
本发明针对CALGRID模式，增加了甲烷化学动态求解过程。CALGRID模式有化学反应、大气输送与扩散、沉降、地面及高架排放源等。其中，化学物种浓度变化方程如下：
上式中已经将原始浓度方程中的物理量分解为平均量为脉动量，然后做了变化，二阶湍流扩散项通过K理论闭合转化为上式第二项。因此(1)式中C是化学物种平均浓度，V是三维风矢量的平均量，K是湍流扩散系数。E是污染源排放。是由于沉降引起的物种浓度变化。PCHEM是化学产生率，LCHEM是化学损失率。式右边分别为平流项、扩散项、源项、沉降项和化学变化项。
鉴于上述各个过程的特征时间不一致，模式采用了算子分裂时间积分方案对各个过程进行积分：
Cn+1＝AxAyAzAcAzAyAxCn (2)
(2)式中Ax、Ay是水平输送扩散算子；Az是垂直输送扩散、物质输入和物理损耗算子，Ac是化学反应算子。
本发明中，(1)平流项模式，采用高阶chapeau函数数值处理方案，并且采用了Forester(1977)的非线性局地滤波，有效地保证了物质守恒并防止出现负浓度和低数值的扩散。(2)垂直输送和扩散项，由于垂直方向的湍流扩散作用远远大于垂直平流，其数值处理方案也不同于水平输送和扩散，这些过程的具体数值解法可以通过Yamartino(1992)法求解。(3)源项，分为面源、周期变化的固定点源、无周期变化的固定点源、流动点源，模式处理时，人为源以面源和无周期变化的固定点源的形式输入；自然源作为面源输入；(4)沉降项，为了能够比较细致准确的反映干沉降过程，CALGRID模式对于不同类型的污染物给出了不同的处理方案，并在计算中将干沉降速率作为地理参数、天气状况、污染特征等因子的函数。具体操作时，用户可以在控制文件中对不同物种选择沉降方案，也可以自己定义或者不考虑干沉降。
对于气体的干沉降，沉降速率被定义为：
Vd＝(ra+rd+rc)-1
式中，ra、rd、rc分别表示湍流层、粘附层和植被层的阻抗。对于固体颗粒物的干沉降，沉降速率被定义为：
Vd＝(ra+rd+ra·rd·vg)-1+vg
式中，ra和rd分别表示湍流层和粘附层的阻抗、vg是重力沉降速度。
以下是对化学反应项具体说明：本发明的CALGRID化学模式中的化学变化项以SAPRC-90的气相化学机制为基础，并在该基础上引入OH自由基和CH4的化学反应，根据OH自由基和CH4的化学反应列表计算OH自由基和CH4的浓度变化。表1列出了改进CALGRID模式化学模块时涉及到的与OH自由基和CH4有关的化学反应。根据其中27式的化学反应速率系数R(27)，可以得到
式中，是甲烷的去除速率；[OH]是OH自由基的浓度。R(27)可根据气象条件和化学机制中参数计算得到，[OH]可根据表1其它化学反应方程联立并作拟稳态假设求得。根据上式，则CH4浓度随时间的变化可以由下式计算：
根据Hybrid或QSSA方法求解甲烷的浓度变化方程，可以得到甲烷的浓度。
表1改进化学模块涉及的与OH、CH4有关的化学反应
(1)以上反应来自Cater等(1990)；“＝”用以区分反应物和生成物；“#数字”表示化学(当量)计算系数。
(2)表中所列光解速率系数为太阳高度为0时的值，热反应速率系数为温度为300K时的值。光解速率系数由(2-10)式计算，热反应速率系数由(2-11)式计算，相关参数见Cater等(1990)。
如图1和2，利用改进化学模块后的CALGRID模式得到的中国地区近地面CH4和OH自由基的年均浓度的空间分布，可见改进后的模式增强了对大气CH4和OH时空分布的模拟能力。一般区域空气质量模式把CH4作为常数，很少给出其分布。一般模型也不输出OH的分布。
Claims (4)
1.一种增加甲烷化学变化的中尺度大气光化学污染模拟预测算法，包括以下步骤：
步骤1：采用CALGRID化学模式，考虑大气化学反应、大气输送与扩散、沉降、地面面源和高架排放源的影响，并对化学物种原始浓度中物理量的平均量分解为脉动量，得到化学物种浓度变化方程如式(1)，
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>C</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mo>&dtri;</mo>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mi>C</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mo>&dtri;</mo>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mi>R</mi>
<mi>C</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>E</mi>
<mo>+</mo>
<msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>C</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>D</mi>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>C</mi>
<mi>H</mi>
<mi>E</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>C</mi>
<mi>H</mi>
<mi>E</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，C是化学物种平均浓度，V是三维风矢量的平均量，K是湍流扩散系数，E是污染源排放，是由于沉降引起的物种浓度变化，PCHEM是化学产生率，LCHEM是化学损失率；式(1)中二阶湍流扩散项▽·(K▽C)通过湍流扩散系数K理论闭合转化得到；式(1)中右侧多项式依次分别为平流项、扩散项、源项、沉降项、化学变化项；
步骤2：由于平流项、扩散项、源项、沉降项、化学变化项各过程的特征时间不一致，通过算子分列时间积分对各过程进行积分，
Cn+1＝AxAyAzAcAzAyAxCn (2)
式中，Ax、Ay是水平输送扩散算子；Az是垂直输送扩散、物质输入和物理损耗算子，Ac是化学反应算子；
步骤3：CALGRID化学模式中的化学变化项以SAPRC-90的气相化学机制为基础，并在该基础上引入OH自由基和甲烷(CH4)的化学反应，根据OH自由基和CH4的化学反应列表计算OH自由基和CH4的浓度变化。
2.根据权利要求1所述增加甲烷化学动态的中尺度大气光化学污染模拟预测算法，其特征在于：步骤3中所述OH自由基和CH4的化学反应列表如表1所示：
表1改进化学模块涉及的与OH、CH4有关的化学反应
3.根据权利要求2所述增加甲烷化学动态的中尺度大气光化学污染模拟预测算法，其特征在于：
通过下式计算甲烷CH4的去除速率，
<mrow>
<msub>
<mi>L</mi>
<mrow>
<msub>
<mi>CH</mi>
<mn>4</mn>
</msub>
</mrow>
</msub>
<mo>=</mo>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mn>27</mn>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mo>&lsqb;</mo>
<mi>O</mi>
<mi>H</mi>
<mo>&rsqb;</mo>
</mrow>
式中，是甲烷的去除速率；[OH]是OH自由基的浓度，R(27)根据气象条件和化学机制中参数计算得到，[OH]可根据表1其它化学反应方程联立并作拟稳态假设求得；CH4的浓度随时间变化则由下式计算：
<mrow>
<mfrac>
<mrow>
<msub>
<mi>dC</mi>
<mrow>
<msub>
<mi>CH</mi>
<mn>4</mn>
</msub>
</mrow>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>1.803</mn>
<mo>-</mo>
<msub>
<mi>L</mi>
<mrow>
<msub>
<mi>CH</mi>
<mn>4</mn>
</msub>
</mrow>
</msub>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>CH</mi>
<mn>4</mn>
</msub>
</mrow>
</msub>
</mrow>
式中，为甲烷的浓度。
4.根据权利要求2所述增加甲烷化学动态的中尺度大气光化学污染模拟预测算法，其特征在于：步骤3中所述CALGRID化学模式中的化学变化项以SAPRC-90的气相化学机制的化学物质如表2：
其中，光分解反应的光解速率系数通过以下计算表达式求解：
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>A</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mi>&infin;</mi>
<mn>0</mn>
</msubsup>
<mi>&sigma;</mi>
<mo>&lsqb;</mo>
<mi>&lambda;</mi>
<mo>,</mo>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mi>&phi;</mi>
<mo>&lsqb;</mo>
<mi>&lambda;</mi>
<mo>,</mo>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mi>I</mi>
<mo>&lsqb;</mo>
<mi>&lambda;</mi>
<mo>,</mo>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>X</mi>
<mo>&rsqb;</mo>
<mi>d</mi>
<mi>&lambda;</mi>
</mrow>
式中，σ[λ,T(h)]是同波长λ和温度T(h)有关的吸收截面系数；φ[λ,T(h)]是量子产额；I[λ,N(h),X]是同波长、大气状态N(h)和空间位置X有关的入射光强；
对于化学反应速率系数随温度升高而变大的双分子气相反应，则用Arrhenius(阿伦尼乌斯)经验公式计算速率系数：
<mrow>
<mi>k</mi>
<mo>=</mo>
<mi>A</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>B</mi>
</msup>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>E</mi>
<mi>a</mi>
</msub>
<mrow>
<mn>0.0019877</mn>
<mi>T</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
式中，A是前置系数、T是环境气温(K)、Tr是参考气温(300K)、Ea是活化能、B是无量纲数。对于CALGRID模式中每种化学反应的A、B、Ea、σ[λ,T(h)]、φ[λ,T(h)]均由气相化学机制(如SAPRC-90)给出。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711085469.XA CN108121892A (zh) | 2017-11-07 | 2017-11-07 | 增加甲烷化学动态的中尺度大气光化学污染模拟预测算法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711085469.XA CN108121892A (zh) | 2017-11-07 | 2017-11-07 | 增加甲烷化学动态的中尺度大气光化学污染模拟预测算法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108121892A true CN108121892A (zh) | 2018-06-05 |
Family
ID=62227753
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711085469.XA Pending CN108121892A (zh) | 2017-11-07 | 2017-11-07 | 增加甲烷化学动态的中尺度大气光化学污染模拟预测算法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108121892A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113704954A (zh) * | 2021-04-13 | 2021-11-26 | 北京师范大学 | 化学物质环境系统行为模拟中大气行为的简化模拟方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003207579A (ja) * | 2002-01-16 | 2003-07-25 | Mitsubishi Heavy Ind Ltd | 大気汚染物質の濃度予測方法及び装置 |
| CN102819661A (zh) * | 2012-06-19 | 2012-12-12 | 中国科学院大气物理研究所 | 一种运用区域空气质量模式的大气环境容量新算法 |
-
2017
- 2017-11-07 CN CN201711085469.XA patent/CN108121892A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003207579A (ja) * | 2002-01-16 | 2003-07-25 | Mitsubishi Heavy Ind Ltd | 大気汚染物質の濃度予測方法及び装置 |
| CN102819661A (zh) * | 2012-06-19 | 2012-12-12 | 中国科学院大气物理研究所 | 一种运用区域空气质量模式的大气环境容量新算法 |
Non-Patent Citations (4)
| Title |
|---|
| 张羽 等: "亚热带农业小流域水系溶存甲烷浓度和扩散通量研究", 《中国农业科学》, vol. 49, no. 20, pages 3968 - 3980 * |
| 胡保安 等: "不同水分条件下巴音布鲁克天鹅湖高寒湿地夏季N2O日排放特征", 《生态环境学报》, vol. 24, no. 5, pages 811 - 817 * |
| 谢昱: "中国自然源排放及对对流层光化学特性和臭氧污染控制的影响研究", 《万方博士论文数据库》, pages 24 - 38 * |
| 鲁易 等: "大气甲烷浓度变化的源汇因素模拟研究进展", 《地球科学进展》, vol. 30, no. 7, pages 763 - 772 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113704954A (zh) * | 2021-04-13 | 2021-11-26 | 北京师范大学 | 化学物质环境系统行为模拟中大气行为的简化模拟方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Czader et al. | Modeling nitrous acid and its impact on ozone and hydroxyl radical during the Texas Air Quality Study 2006 | |
| Fowler et al. | Modelling photochemical oxidant formation, transport, deposition and exposure of terrestrial ecosystems | |
| Karset et al. | Strong impacts on aerosol indirect effects from historical oxidant changes | |
| CN106975311A (zh) | 烟气除霾减碳微藻养殖系统及烟气处理方法 | |
| CN110020448A (zh) | 边界层气象参数改进型中尺度大气光化学污染模拟预测算法 | |
| Paladino et al. | Scale-up of photo-bioreactors for microalgae cultivation by π-theorem | |
| CN207385503U (zh) | 一种光化学烟雾箱 | |
| CN108121892A (zh) | 增加甲烷化学动态的中尺度大气光化学污染模拟预测算法 | |
| CN206996160U (zh) | 烟气除霾减碳微藻养殖系统 | |
| Abdi et al. | Prediction compost criteria of organic wastes with biochar additive in in-vessel composting machine using ANFIS and ANN methods | |
| Zhao et al. | Linking soil dissolved organic matter characteristics and the temperature sensitivity of soil organic carbon decomposition in the riparian zone of the Three Gorges Reservoir | |
| Liu et al. | Investigating the impact of pretreatment strategies on photocatalyst for accurate CO2RR productivity quantification: a machine learning approach | |
| Platt et al. | Fundamental issues in measurement of primary production | |
| CN108170882B (zh) | 基于查算表的中尺度大气光化学模拟分段线性回归算法 | |
| CN102921418B (zh) | 一种立方块状氧化亚铜可见光催化剂的合成方法 | |
| Pareek et al. | Mathematical modeling and simulation of methane gas production in simulated landfill column reactors under sulfidogenic and methanogenic environments | |
| Aggarwal et al. | Application of native mix algal strain for gray water treatment and biofuel production: Preliminary study | |
| Rahman et al. | Multi-objective genetic algorithm optimization of energy efficiency and biomass concentration of Synechococcus HS-9 cultivation for third-generation biodiesel feedstock | |
| CN116840002B (zh) | 红树林幼苗沉积物-大气界面二氧化碳排放来源的解析方法 | |
| Sugata et al. | Simulation of Sulfate Aerosol in Ease ASIA using Models-3/CMAQ with RAMS Meteorological Data | |
| CN110003248A (zh) | 一种多羧酸铀酰配合物光催化剂及其制备方法 | |
| Huebert et al. | Marine aerosol and gas exchange and global atmospheric effects | |
| CN208159511U (zh) | 一种便于携带的植物碳同位素标记仪 | |
| Tudoroiu et al. | Fuzzy UML and Petri nets modeling investigations on the pollution impact on the air quality in the vicinity of the Black Sea Constanta Romanian resort | |
| Seakins et al. | Analysis of Chamber Data |
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 |