CN107861914A - 一种统计分析大气甲烷长期变化趋势和季节变化规律的方法 - Google Patents
一种统计分析大气甲烷长期变化趋势和季节变化规律的方法 Download PDFInfo
- Publication number
- CN107861914A CN107861914A CN201711276691.8A CN201711276691A CN107861914A CN 107861914 A CN107861914 A CN 107861914A CN 201711276691 A CN201711276691 A CN 201711276691A CN 107861914 A CN107861914 A CN 107861914A
- Authority
- CN
- China
- Prior art keywords
- methane
- data
- variation tendency
- mrow
- seasonal
- 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
Classifications
-
- 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/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Discrete Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
本发明公开了一种统计分析大气甲烷长期变化趋势和季节变化规律的方法,包括:(一)建立长期变化趋势和季节变化规律模型V:V(t,b)=b0+b1cos2πt+b2sin2πt+b3cos4πt+b4sin4πt+...,然后建立地基、机载或星载观测的主被动大气甲烷总柱浓度数据用函数模型F:F(t,a,b)=at+V(t,b);(二)采用自举重采样技术对非高斯分布的数据进行严格拟合分析。本发明的方法通过建立低阶傅里叶级数模型以及更严格的统计分析,可获得给定数据集的平均趋势以及严格的季节变化估计值,准确量化甲烷的长期变化趋势和季节变化规律,在大气探测以及气候变化领域具有重要意义,可用于地基、机载、星载主被动大气甲烷观测数据分析,为大气甲烷源汇分析及其动态变化研究和模拟提供有价值的数据源。
Description
技术领域
本发明涉及一种统计分析气体变化趋势和规律的方法,具体涉及一种统计分析大气甲烷长期变化趋势和季节变化规律的方法,属于大气探测以及气候变化技术领域。
背景技术
全球气候变化是当今世界面临的重要环境问题之一,这一问题的关键要素是理解辐射活性气体以及参与大气化学反应的温室气体在大气中的行为特性。长期观测这些气体为研究这些气体的演变以及源汇变化提供了实验数据。
甲烷作为最主要的温室气体之一,其在大气中的含量远低于二氧化碳,由于长期甲烷观测数据中存在大量季节变化和其他因素影响,使得准确量化甲烷的长期变化趋势和季节变化规律成为一种挑战。
现有的趋势分析技术,在估计甲烷长期变化趋势和季节变化规律方面,主要基于不严格的数据统计分布假设,通常分析结果不稳定,并且容易给出错误的分析结果,实用性尚未被证明。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种统计分析大气甲烷长期变化趋势和季节变化规律的方法,该方法通过更严格的统计分析,可获得给定数据集的平均趋势(大气甲烷年增长速度及其不确定度)以及严格的季节变化估计值,准确量化甲烷的长期变化趋势和季节变化规律。
为了实现上述目标,本发明采用如下的技术方案:
一种统计分析大气甲烷长期变化趋势和季节变化规律的方法,其特征在于,包括以下步骤:
一、建立长期变化趋势和季节变化规律模型
长期变化趋势和季节变化规律模型基于低阶傅里叶级数,该低阶傅里叶级数用函数模型V表示:
V(t,b)=b0+b1cos2πt+b2sin2πt+b3cos4πt+b4sin4πt+...
式中,t为年时间序列,b为季节变化拟合系数,季节变化为自然年的函数;
地基、机载或星载观测的主被动大气甲烷总柱浓度数据用函数模型F表示:
F(t,a,b)=at+V(t,b)
式中,t为年时间序列,a为年变化趋势系数,b为季节变化拟合系数;
二、采用自举重采样技术对非高斯分布的数据进行严格拟合分析自举重采样技术基于与随机影响相关的数据分布,可由模型拟合残差来表示,拟合残差用Q表示:
式中,M为地基、机载或星载观测的主被动大气甲烷总柱浓度数据集,F为数据集M的拟合模型,i=1,2,.....,m,m为数据集M中数据个数;
通过最小二乘拟合算法获得该式最小值,此时可得拟合系数a0和b0,此时的拟合残差用R表示,即:
Ri,0=M-F(ti,a0,b0);
假定Ri,q是数据集{Ri,0}的随机采样,可得下式:
Mi,q=Ri,q+F(ti,a0,b0)
式中,q=1,2,…,N,N为重复迭代次数;
通过重复迭代取样后,得到新的数据集{(ti,Mi,q)};
再次用模型F拟合这一新的数据集{(ti,Mi,q)}获得参数aq和bq,这一过程重复q=1,…,N次生成向量A={aq}和矩阵B={bq}。
本发明的有益之处在于:
本发明的方法通过建立低阶傅里叶级数模型以及更严格的统计分析,可获得给定数据集的平均趋势(大气甲烷年增长速度及其不确定度)以及严格的季节变化估计值,准确量化甲烷的长期变化趋势和季节变化规律,在大气探测以及气候变化领域具有重要意义,可用于地基、机载、星载主被动大气甲烷观测数据分析,为大气甲烷源汇分析及其动态变化研究和模拟提供有价值的数据源。
附图说明
图1是地基傅里叶变换光谱仪长期甲烷时间序列趋势拟合结果图;
图2(a)是地基傅里叶变换光谱仪长期甲烷时间序列与过境卫星GOSAT甲烷数据对比结果图;
图2(b)是地基傅里叶变换光谱仪长期甲烷时间序列与过境卫星GOSAT甲烷数据相关性分析结果图。
具体实施方式
以下结合附图和具体实施例对本发明作具体的介绍。
第一部分:统计分析大气甲烷长期变化趋势和季节变化规律的方法
本发明提供的统计分析大气甲烷长期变化趋势和季节变化规律的方法主要包括以下两步:
一、建立长期变化趋势和季节变化规律模型
长期变化趋势和季节变化规律模型基于低阶傅里叶级数,在长期观测的甲烷序列和年份之间建立一定关系。
低阶傅里叶级数用函数模型V表示:
V(t,b)=b0+b1cos2πt+b2sin2πt+b3cos4πt+b4sin4πt+...
式中,t为年时间序列,b为季节变化拟合系数,季节变化为自然年的函数。
地基、机载或星载观测的主被动大气甲烷总柱浓度数据用函数模型F表示:
F(t,a,b)=at+V(t,b)
式中,t为年时间序列,a为年变化趋势系数,b为季节变化拟合系数。
函数模型F能够模拟出数据潜在的变化周期。
二、采用自举重采样技术对非高斯分布的数据进行严格拟合分析自举重采样技术基于与随机影响相关的数据分布,可由模型拟合残差来表示。
拟合残差用Q表示:
式中,M为地基、机载或星载观测的主被动大气甲烷总柱浓度数据集,F为数据集M的拟合模型,i=1,2,.....,m,m为数据集M中数据个数。
通过最小二乘拟合算法获得该式最小值,此时可得拟合系数a0和b0,此时的拟合残差用R表示,即:
Ri,0=M-F(ti,a0,b0)
假定Ri,q是数据集{Ri,0}的随机采样,可得下式:
Mi,q=Ri,q+F(ti,a0,b0)
式中,q=1,2,…,N,N为重复迭代次数。
通过重复迭代取样后,得到新的数据集{(ti,Mi,q)}。
再次用模型F拟合这一新的数据集{(ti,Mi,q)}获得参数aq和bq,这一过程重复q=1,…,N次生成向量A={aq}和矩阵B={bq}。
向量A由趋势参数分布的采样组成,可以选择2.5和97.5经验分布百分位数值来表示95%置信区间。
第二部分:验证本发明提供的方法及其置信区间的有效性
在此,我们采用偏差校正公式Z0来量化本发明方法的偏差:
Z0=C-1(r)
式中,C-1表示标准正态累积分布函数的反函数,r为自举重采样比重。
当r=1/2时,Z0=0,这意味着本发明的方法没有偏差。
将此偏差校正公式应用到自举重采样分析大气甲烷长期观测时间序列中,r值范围一般处于0.43至0.54之间,通过实际数据多次分析发现均接近0.5,可见,本发明的方法对甲烷长期变化趋势和季节变化规律的分析结果不会引入误差。
第三部分:应用案例
为了让大家能够更直观、更透彻的理解本发明的统计分析方法,下面我们结合具体的案例来对本发明的统计分析方法做进一步说明。
地点:安徽合肥地区。
时间:2014年7月至2017年7月。
应用本发明提供的统计分析方法,我们得到该地区的大气甲烷年增长趋势如图1所示。
经拟合,该地区的大气甲烷年增长趋势为14.08±1.33ppb。
从图1中我们还可以看出明显的季节变化规律:冬季甲烷柱浓度较低,夏季甲烷柱浓度较高。
我们将采用本发明的统计分析方法得到的地基傅里叶变换光谱仪长期甲烷时间序列与过境卫星GOSAT甲烷数据进行了对比和相关性分析。对比结果见图2(a),相关性分析结果见图2(b)。
由图2(a)和图2(b)可知,应用本发明的统计分析方法,我们获得了与过境卫星GOSAT甲烷数据一致的趋势拟合结果。
需要说明的是,上述实施例不以任何形式限制本发明,凡采用等同替换或等效变换的方式所获得的技术方案,均落在本发明的保护范围内。
Claims (1)
1.一种统计分析大气甲烷长期变化趋势和季节变化规律的方法,其特征在于,包括以下步骤:
一、建立长期变化趋势和季节变化规律模型
长期变化趋势和季节变化规律模型基于低阶傅里叶级数,该低阶傅里叶级数用函数模型V表示:
V(t,b)=b0+b1cos2πt+b2sin2πt+b3cos4πt+b4sin4πt+...
式中,t为年时间序列,b为季节变化拟合系数,季节变化为自然年的函数;
地基、机载或星载观测的主被动大气甲烷总柱浓度数据用函数模型F表示:
F(t,a,b)=at+V(t,b)
式中,t为年时间序列,a为年变化趋势系数,b为季节变化拟合系数;
二、采用自举重采样技术对非高斯分布的数据进行严格拟合分析
自举重采样技术基于与随机影响相关的数据分布,可由模型拟合残差来表示,拟合残差用Q表示:
<mrow>
<mi>Q</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>M</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>F</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
式中,M为地基、机载或星载观测的主被动大气甲烷总柱浓度数据集,F为数据集M的拟合模型,i=1,2,.....,m,m为数据集M中数据个数;
通过最小二乘拟合算法获得该式最小值,此时可得拟合系数a0和b0,此时的拟合残差用R表示,即:
Ri,0=M-F(ti,a0,b0);
假定Ri,q是数据集{Ri,0}的随机采样,可得下式:
Mi,q=Ri,q+F(ti,a0,b0)
式中,q=1,2,…,N,N为重复迭代次数;
通过重复迭代取样后,得到新的数据集{(ti,Mi,q)};
再次用模型F拟合这一新的数据集{(ti,Mi,q)}获得参数aq和bq,这一过程重复q=1,…,N次生成向量A={aq}和矩阵B={bq}。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711276691.8A CN107861914A (zh) | 2017-12-06 | 2017-12-06 | 一种统计分析大气甲烷长期变化趋势和季节变化规律的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711276691.8A CN107861914A (zh) | 2017-12-06 | 2017-12-06 | 一种统计分析大气甲烷长期变化趋势和季节变化规律的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107861914A true CN107861914A (zh) | 2018-03-30 |
Family
ID=61705293
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711276691.8A Pending CN107861914A (zh) | 2017-12-06 | 2017-12-06 | 一种统计分析大气甲烷长期变化趋势和季节变化规律的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107861914A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112668166A (zh) * | 2020-12-21 | 2021-04-16 | 西安理工大学 | 一种建立随时空变化大气温度模型的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155404A (zh) * | 2014-09-03 | 2014-11-19 | 中国科学院广州地球化学研究所 | 大气甲烷和一氧化碳分析系统和方法 |
-
2017
- 2017-12-06 CN CN201711276691.8A patent/CN107861914A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155404A (zh) * | 2014-09-03 | 2014-11-19 | 中国科学院广州地球化学研究所 | 大气甲烷和一氧化碳分析系统和方法 |
Non-Patent Citations (1)
Title |
---|
T.GARDINER ET AL.: "Trend analysis of greenhouse gases over Europe measured by a network of ground-based remote FTIR instruments", 《ATMOSPHERIC CHEMISTRY AND PHYSICS》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112668166A (zh) * | 2020-12-21 | 2021-04-16 | 西安理工大学 | 一种建立随时空变化大气温度模型的方法 |
CN112668166B (zh) * | 2020-12-21 | 2023-03-14 | 西安理工大学 | 一种建立随时空变化大气温度模型的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bao et al. | The flexible global ocean-atmosphere-land system model, spectral version 2: FGOALS-s2 | |
Baek et al. | Climate change in the 21st century simulated by HadGEM2-AO under representative concentration pathways | |
Bugmann | Sensitivity of forests in the European Alps to future climatic change | |
Haerter et al. | Climate model bias correction and the role of timescales | |
Masson-Delmotte et al. | Recent changes in north-west Greenland climate documented by NEEM shallow ice core data and simulations, and implications for past-temperature reconstructions | |
Zhang et al. | Runoff sensitivity to global mean temperature change in the CMIP5 Models | |
Schilt et al. | Isotopic constraints on marine and terrestrial N2O emissions during the last deglaciation | |
Kerkhoff et al. | Assessment of bias assumptions for climate models | |
M’Po et al. | Comparison of daily precipitation bias correction methods based on four regional climate model outputs in Ouémé Basin, Benin | |
Wang et al. | Preliminary evaluations of FGOALS-g2 for decadal predictions | |
Graf et al. | Intercomparison of methods for the simultaneous estimation of zero-plane displacement and aerodynamic roughness length from single-level eddy-covariance data | |
Mitchell et al. | Accounting for model error due to unresolved scales within ensemble Kalman filtering | |
Gray | Climate change 2007: the physical science basis summary for policymakers | |
Wei et al. | The effect of precipitation on hydropower generation capacity: a perspective of climate change | |
CN107861914A (zh) | 一种统计分析大气甲烷长期变化趋势和季节变化规律的方法 | |
Kuhn et al. | Use of field experiments in soil erosion research | |
Kim et al. | Spatial analysis of future East Asian seasonal temperature using two regional climate model simulations | |
Liang et al. | Asymmetrical precipitation sensitivity to temperature across global dry and wet regions | |
CN104156559A (zh) | 确定含有重现期修正系数的暴雨强度公式的一种方法 | |
Poulain et al. | Evaluation of the inter-annual variability of stratospheric chemical composition in chemistry-climate models using ground-based multi species time series | |
Le Lay et al. | How changing rainfall regimes may affect the water balance: a modelling approach in West Africa | |
Fanego et al. | A novel procedure for generating solar irradiance TSYs | |
Nyeko-Ogiramoi et al. | Nonparametric statistical downscaling of precipitation from global climate models | |
CN116681315B (zh) | 温室气体净排放的估算方法、装置、设备和存储介质 | |
Rai et al. | The dichotomy of wet and dry trends over India by aerosol indirect effects in CMIP5 models |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180330 |
|
WD01 | Invention patent application deemed withdrawn after publication |