本发明专利技术公开了一种确定稻田温室气体通量与风速风向倾斜关系的方法，包括以下步骤：（1）同时搜集稻田温室气体通量及风速风向的观测值；（2）将风向观测值转换为弧度计量单位；（3）将温室气体通量与风向的倾斜周期函数进行回归（周期回归）；（4）计算温室气体通量与风向周期回归方程结果间的残差（周期残差）；（5）计算温室气体通量周期残差与风速间的线性回归方程；（6）整理得到稻田温室气体通量与风速风向倾斜关系方程；（7）评估稻田温室气体通量的风向效应与风速效应强度。本发明专利技术能够同时精确估计稻田温室气体通量与风速风向间的倾斜周期关系，对评估风速风向对稻田温室气体通量的影响有重要作用。
A method to determine the relationship between greenhouse gas flux and wind speed and wind direction tilt in rice field
【技术实现步骤摘要】
一种确定稻田温室气体通量与风速风向倾斜关系的方法
本专利技术涉及一种确定稻田温室气体通量与风速风向倾斜关系的方法，属于气候变化管理领域。
技术介绍
温室效应能够引起气候变暖，温度上升，生态系统退化等多种自然灾害，因此受到广泛的关注。温室效应主要是由人类生产活动过程中排放温室气体引起的。稻田是氧化亚氮和甲烷等多种温室气体的主要来源之一。稻田温室气体的通量既受到风速的影响也受到风向的影响，正确估计风速风向对稻田温室气体排放的影响是计算稻田温室气体排放通量，减少温室气体排放，延缓全球气候变暖趋势的前提；也是国际气候谈判的重要内容；对保护人类生存的环境具有重要的作用。
技术实现思路
专利技术目的是为了克服现有技术中存在的不足，提供一种确定稻田温室气体通量与风速风向倾斜关系的方法。本专利技术充分利用了稻田温室气体通量及其风速风向的观测值能够同时精确估计不同风速和风向条件下稻田温室气体通量，为进行国际气候谈判，减少温室气体排放奠定基础。本专利技术的目的是这样实现的，本专利技术的一种确定稻田温室气体通量与风速风向倾斜关系的方法，包括以下步骤：(1)同时搜集或测量所需研究的稻田温室气体通量及风速风向的观测值；记稻田温室气体通量数据为y1,y2,",yM,；确定某一方向为0度，某一旋转方向为正方向；一般设置正北方向为0度，以逆时针方向旋转为正方向；在观测温室气体同时，观测风向数据为s1,s2,",sM,风速数据为w1,w2,",wM，其中，M为通量数据的观测个数其不能小于50个，M≥50；<br>(2)将风向观测值转换为弧度计量单位；(3)将温室气体通量与风向的倾斜周期函数进行回归(周期回归)；(4)计算温室气体通量与风向周期回归方程结果间的残差(周期残差)；(5)计算温室气体通量周期残差与风速间的线性回归方程；(6)整理得到稻田温室气体通量与风速风向倾斜关系方程；(7)评估稻田温室气体通量的风向效应、风速效应以及倾斜效应强度。作为优选，所述步骤(2)具体包括以下步骤：a确保所有的风向数据为0度(包含)至360度(不含)之间，不在此范围内的风向角度通过加减360度的整数倍，将其转换到此范围内；经过此转换的风向角仍记为s1,s2,",sM；b.通过公式将风向角转为弧度制，转换结果记为d1,d2,",dM。作为优选，所述步骤(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度之间，用公式将所有观测到的风向值转为弧度制，其与甲烷通量散点图如图2所示。(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.一种确定稻田温室气体通量与风速风向倾斜关系的方法，其特征在于，包括以下步骤：/n(1)同时搜集或测量所需研究的稻田温室气体通量及风速风向的观测值,记稻田温室气体通量数据为y
【技术特征摘要】
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.通过公式将风向角转为弧度制，转换结果记为d1,d2,…,dM。
3.根据权利要求2所述一种确定稻田温室气体通量与风速风向倾斜关系的方法，其特征在于，所述步骤(3)具体包括以下步骤：
a.将观测到的稻田温室气体通量值y1,y2,…,yM,与其同时观测到的风向角d1,d2,…,dM，代入下面的倾斜回归方程，估计其中的未知参数λ0，λ1，β，α
Y＝λ0+λ1cos[(ωd-β)+αcos(ωd-β)](1)
其中...
【专利技术属性】
技术研发人员：王玉琳，汪靓，程浩淼，何成达，程吉林，
申请(专利权)人：扬州大学，
类型：发明
国别省市：江苏;32
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。