本发明涉及一种确定稻田温室气体通量与风速风向倾斜关系的方法，属于气候变化管理领域。
背景技术：
温室效应能够引起气候变暖，温度上升，生态系统退化等多种自然灾害，因此受到广泛的关注。温室效应主要是由人类生产活动过程中排放温室气体引起的。稻田是氧化亚氮和甲烷等多种温室气体的主要来源之一。稻田温室气体的通量既受到风速的影响也受到风向的影响，正确估计风速风向对稻田温室气体排放的影响是计算稻田温室气体排放通量，减少温室气体排放，延缓全球气候变暖趋势的前提；也是国际气候谈判的重要内容；对保护人类生存的环境具有重要的作用。
技术实现要素：
发明目的是为了克服现有技术中存在的不足，提供一种确定稻田温室气体通量与风速风向倾斜关系的方法。本发明充分利用了稻田温室气体通量及其风速风向的观测值能够同时精确估计不同风速和风向条件下稻田温室气体通量，为进行国际气候谈判，减少温室气体排放奠定基础。
本发明的目的是这样实现的，本发明的一种确定稻田温室气体通量与风速风向倾斜关系的方法，包括以下步骤：
(1)同时搜集或测量所需研究的稻田温室气体通量及风速风向的观测值；记稻田温室气体通量数据为y1,y2,",ym,；
确定某一方向为0度，某一旋转方向为正方向；一般设置正北方向为0度，以逆时针方向旋转为正方向；在观测温室气体同时，观测风向数据为s1,s2,",sm,风速数据为w1,w2,",wm，其中，m为通量数据的观测个数其不能小于50个，m≥50；
(2)将风向观测值转换为弧度计量单位；
(3)将温室气体通量与风向的倾斜周期函数进行回归(周期回归)；
(4)计算温室气体通量与风向周期回归方程结果间的残差(周期残差)；
(5)计算温室气体通量周期残差与风速间的线性回归方程；
(6)整理得到稻田温室气体通量与风速风向倾斜关系方程；
(7)评估稻田温室气体通量的风向效应、风速效应以及倾斜效应强度。
作为优选，所述步骤(2)具体包括以下步骤：
a确保所有的风向数据为0度(包含)至360度(不含)之间，不在此范围内的风向角度通过加减360度的整数倍，将其转换到此范围内；经过此转换的风向角仍记为s1,s2,",sm；
b.通过公式
作为优选，所述步骤(3)具体包括以下步骤：
a.将观测到的稻田温室气体通量值y1,y2,",ym,与其同时观测到的风向角d1,d2,",dm，代入下面的倾斜回归方程，估计其中的未知参数λ0，λ1，β，α
y＝λ0 λ1cos[(ωd-β) αcos(ωd-β)](1)
其中，ω与风向角的周期有关，一般定为1；此外，α有一定的取值范围，
b.注意方程(1)是非线性的回归方程，且无法转化为线性的回归方程；求解该方程参数推荐用拟牛顿法。
作为优选，所述步骤(4)具体包括以下步骤：
a.将di,i＝1,2,3,",m分别代入方程(1)，得到温室气体通量倾斜周期回归方程的结果ri,i＝1,2,",m；
b.将yi-ri,i＝1,2,",m得到m个周期残差εi,i＝1,2,",m。
作为优选，所述步骤(5)具体包括以下步骤：
以εi,i＝1,2,",m为因变量，以εi,i＝1,2,",m对应风速w1,w2,",wm为自变量用极大似然法估计方程(2)中的系数c和λw：
ε＝cw λw(2)。
作为优选，所述步骤(6)具体包括以下步骤：
根据方程(1)和方程(2)的结果，整理得到对应于温室气体通量周期和趋势的统一方程为
yd,w＝λ1cos[(ωd-β) αcos(ωd-β)] cw λ0 λw(3)
yd,w代表风向为d风速为w的温室气体通量预测值。
作为优选，所述步骤(7)具体包括以下步骤：
a.c代表了风速对稻田温室气体通量的效应强度，即风速效应强度；
b.λ1是风向效应强度；
c.α是风向的倾斜效应强度，若α大于0则温室气体通量向左倾斜，即偏北风时通量较大；若α小于0则相反；若α等于0则温室气体通量无倾斜效应，即倾斜效应强度为0。
理论和实验研究都证明稻田排放是大气中温室气体的重要来源，估计稻田温室气体通量对控制其中温室气体的排放有着重要意义，其对减缓全球气候变暖有着重要意义。本发明利用稻田温室气体排放通量的观测数据，得到一定风速风向下计算公式，为进一步减少稻田温室气体的排放奠定基础。
本发明提供种确定稻田温室气体通量与风速风向倾斜关系的方法。本发明充分利用了稻田温室气体通量及其风速风向的观测值能够同时精确估计不同风速和风向条件下稻田温室气体通量，为进行国际气候谈判，减少温室气体排放奠定基础。
有益效果：本发明基于统计学理论通过两步估计法得到了精确估计一定风速和风向情况下稻田温室气体通量的计算方程；进而能够为国际气候变化谈判提供了依据，为评估及控制稻田温室气体排放指明方向。
附图说明
图1为本发明流程图。
图2为稻田甲烷通量及风向散点图。
具体实施方式
下面结合附图和我国某地区稻田甲烷的实际观测数据对本发明作更进一步的说明：
(1)按照图1所示流程，测量该农田甲烷通量观测值100个，并同时测量该通量值对应的风速和风向；其中，记录风向时以正北方为0度，以逆时针为正方向。
(2)搜集到所有数据的风向均在0到360度之间，用公式
(3)根据流程图和图2，将风向和甲烷通量代入式(1)的用拟牛顿法估计得到风向和甲烷通量的倾斜回归方程，
y＝1.80cos[(d-4.04)-0.38cos(d-4.04)] 1.93
即λ0＝1.93，λ1＝1.80，β＝4.04，α＝0.38
(4)根据上述方程依次计算出100个甲烷观测值对应该方程的100周期残差；
(5)以计算得到的甲烷通量残差为因变量，以对应的风速为自变量进行线性回归，得到风速与甲烷通量间的回归方程：
ε＝0.2w 0.47
(6)根据以上计算的结果，可得该稻田风向为d，风速为w的甲烷通量回归方程及预测值yd,w为
yd,w＝1.80cos[(d-4.04)-0.38cos(d-4.04)] 0.2w 2.4
(7)根据以上的结果得到该稻田甲烷通量的风速效应强度为0.2毫克/(天*10平方米*米)；风向效应强度为1.8毫克/(天*10平方米)；倾斜效应强度为0.38，该稻田在偏南风时甲烷通量较大。
1.一种确定稻田温室气体通量与风速风向倾斜关系的方法，其特征在于，包括以下步骤：
(1)同时搜集或测量所需研究的稻田温室气体通量及风速风向的观测值,记稻田温室气体通量数据为y1,y2,…,ym,；
确定某一方向为0度，某一旋转方向为正方向；一般设置正北方向为0度，以逆时针方向旋转为正方向；在观测温室气体同时，观测风向数据为s1,s2,…,sm,风速数据为w1,w2,…,wm，其中，m为通量数据的观测个数，其不能小于50个，m≥50；
(2)将风向观测值转换为弧度计量单位；
(3)将温室气体通量与风向的倾斜周期函数进行回归(周期回归)；
(4)计算温室气体通量与风向周期回归方程结果间的残差(周期残差)；
(5)计算温室气体通量周期残差与风速间的线性回归方程；
(6)整理得到稻田温室气体通量与风速风向倾斜关系方程；
(7)评估稻田温室气体通量的风向效应、风速效应以及倾斜效应强度。
2.根据权利要求1所述一种确定稻田温室气体通量与风速风向倾斜关系的方法，其特征在于，所述步骤(2)具体包括以下步骤：
a确保所有的风向数据为0度(包含)至360度(不含)之间，不在此范围内的风向角度通过加减360度的整数倍，将其转换到此范围内；经过此转换的风向角仍记为s1,s2,…,sm；
b.通过公式
3.根据权利要求2所述一种确定稻田温室气体通量与风速风向倾斜关系的方法，其特征在于，所述步骤(3)具体包括以下步骤：
a.将观测到的稻田温室气体通量值y1,y2,…,ym,与其同时观测到的风向角d1,d2,…,dm，代入下面的倾斜回归方程，估计其中的未知参数λ0，λ1，β，α
y＝λ0 λ1cos[(ωd-β) αcos(ωd-β)](1)
其中，ω与风向角的周期有关，一般定为1；此外，α有一定的取值范围，
b.注意方程(1)是非线性的回归方程，且无法转化为线性的回归方程；求解该方程参数用拟牛顿法。
4.根据权利要求3所述一种确定稻田温室气体通量与风速风向倾斜关系的方法，其特征在于，所述步骤(4)具体包括以下步骤：
a.将di,i＝1,2,3,…,m分别代入方程(1)，得到温室气体通量倾斜周期回归方程的结果ri,i＝1,2,…,m；
b.将yi-ri,i＝1,2,…,m得到m个周期残差εi,i＝1,2,…,m。
5.根据权利要求4所述一种确定稻田温室气体通量与风速风向倾斜关系的方法，其特征在于，所述步骤(5)具体包括以下步骤：
以εi,i＝1,2,…,m为因变量，以εi,i＝1,2,…,m对应风速w1,w2,…,wm为自变量用极大似然法估计方程(2)中的系数c和λw：
ε＝cw λw(2)。
6.根据权利要求5所述一种确定稻田温室气体通量与风速风向倾斜关系的方法，其特征在于，所述步骤(6)具体包括以下步骤：
根据方程(1)和方程(2)的结果，整理得到对应于温室气体通量周期和趋势的统一方程为
yd,w＝λ1cos[(ωd-β) αcos(ωd-β)] cw λ0 λw(3)
yd,w代表风向为d风速为w的温室气体通量预测值。
7.根据权利要求1所述一种确定稻田温室气体通量与风速风向倾斜关系的方法，其特征在于，所述步骤(7)具体包括以下步骤：
a.c代表了风速对稻田温室气体通量的效应强度，即风速效应强度；
b.λ1是风向效应强度；
c.α是风向的倾斜效应强度，若α大于0则温室气体通量向左倾斜，即偏北风时通量较大；若α小于0则相反；若α等于0则温室气体通量无倾斜效应，即倾斜效应强度为0。
技术总结