CN112700089A - 一种草原灌溉地水资源消耗量的计算方法 - Google Patents
一种草原灌溉地水资源消耗量的计算方法 Download PDFInfo
- Publication number
- CN112700089A CN112700089A CN202011469032.8A CN202011469032A CN112700089A CN 112700089 A CN112700089 A CN 112700089A CN 202011469032 A CN202011469032 A CN 202011469032A CN 112700089 A CN112700089 A CN 112700089A
- Authority
- CN
- China
- Prior art keywords
- grassland
- representing
- area
- land
- irrigation
- 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.)
- Granted
Links
- 238000003973 irrigation Methods 0.000 title claims abstract description 103
- 230000002262 irrigation Effects 0.000 title claims abstract description 102
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 230000008859 change Effects 0.000 claims abstract description 15
- 230000004907 flux Effects 0.000 claims description 40
- 230000005855 radiation Effects 0.000 claims description 27
- 239000002689 soil Substances 0.000 claims description 21
- 238000001556 precipitation Methods 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 12
- 238000001704 evaporation Methods 0.000 claims description 11
- 238000007781 pre-processing Methods 0.000 claims description 9
- 230000008020 evaporation Effects 0.000 claims description 8
- 230000012010 growth Effects 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 238000009834 vaporization Methods 0.000 claims description 3
- 230000008016 vaporization Effects 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 238000011160 research Methods 0.000 abstract description 4
- 238000010276 construction Methods 0.000 abstract description 3
- 238000011161 development Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000012546 transfer Methods 0.000 description 6
- 238000012217 deletion Methods 0.000 description 4
- 230000037430 deletion Effects 0.000 description 4
- 238000002310 reflectometry Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 244000025254 Cannabis sativa Species 0.000 description 1
- 208000027066 STING-associated vasculopathy with onset in infancy Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000005273 aeration Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004459 forage Substances 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 241000411851 herbal medicine Species 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000008635 plant growth Effects 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 230000005068 transpiration Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Primary Health Care (AREA)
- Entrepreneurship & Innovation (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种草原灌溉地水资源消耗量的计算方法,利用草原区域参考作物系数在时间尺度上稳定渐变特性,等价分解区域蒸散发ET长序列扩展,利用改进参考作物系数法,耦合区域参考作物蒸散发,实现了遥感影像缺失特别是非晴日干扰下草原区域蒸散发长序列扩展,并解决了遥感区域蒸散发ET计算困难问题,提高草原区域蒸散发ET长序列计算的精度。本发明针对草原灌溉地水资源消耗量计算问题,通过灌溉地水资源消耗计算模型,量化确定了草原地区灌溉地作物生长季的水资源消耗量GW值,界定了草原区地下水资源实时监控和管理的重点区域,为灌溉地植被建设对草原区地下水资源开发利用管理和草原生态系统研究提供了典型案例分析参考。
Description
技术领域
本发明属于牧区水资源管理和用水监控领域,具体涉及一种草原灌溉地水资源消耗量的计算方法。
背景技术
对于生态结构相对脆弱的草原地区,灌溉在增加饲草料和粮食经济作物产量、提高牧区抗灾保畜能力和草地承载能力的同时,也会改变草原生态系统的水循环过程,使其呈现“自然-社会”二元特性。同时草原会消耗当地水资源直接影响生态环境变迁,合理计算灌溉地耗水量对草原生态系统和牧区经济社会稳定具有举足轻重的作用。
草原灌溉地耗水量计算过程中,灌溉产生的蒸散发(ET)和降水(Pr)是影响草原水资源系统的重要因素。如何合理计算这两个参数,并定量描述草原灌溉地耗水特征是牧区水资源可持续利用和草原生态系统稳定的重要基础。相比Pr,ET由于需要依据相对复杂的物理原理,并采用系统技术测量获得,该因素成为实时监测草原灌溉地耗水量的关键。传统点尺度监测方法并不能完全真实反映区域尺度上的ET特征,遥感技术利用其空间连续性和大跨度特点,可直接获得作为下垫面消耗的水量,解决了区域尺度特别是缺乏观测资料地区的ET定量表征问题,为区域ET计算提供了一种新途径。但是,受卫星过境时间、非晴日天气的干扰,遥感技术只能获取有限、离散的区域ET数据,而对草原地区水资源综合利用应用价值具有重要作用的长序列(月、季、年等)ET数据,很难获得。目前,现有区域ET长序列扩展是在有限、离散的晴日ET基础上实现的,当遥感影像过境时间不匹配、降水阴云非晴日出现时,天气因素会通过改变下垫面地表温度和土壤湿度等直接影响ET变化。现有诸如蒸发比不变法、正弦曲线法、冠层阻力法、恒定参考作物系数法等开展时间尺度长序列扩展时,利用两幅区域ET遥感计算结果,通过插补遥感影像缺失的区域ET,除主观弱化了地表温度、土壤湿度等改变对ET的干扰外,更容易弱化天气特别是非晴日所带来的气象参数变化影响、忽视下垫面植被生长带来的渐变,这种忽视天气和植被参数变化得到的长序列ET结果精度不能保证,很难有效控制误差的聚散,严重影响计算结果和精度。
发明内容
针对现有技术中的上述不足,本发明提供的一种草原灌溉地水资源消耗量的计算方法解决了现有技术中存在的问题。
为了达到上述发明目的,本发明采用的技术方案为:一种草原灌溉地水资源消耗量的计算方法,包括以下步骤:
S1、采集草原灌溉地的基础数据并对基础数据进行预处理,根据预处理后的数据,使用METRIC模型反演得到地表反照率、地表比辐射率、地表温度、归一化植被指数和地表植被粗糙度;
S2、根据地表反照率、地表比辐射率、地表温度、归一化植被指数和地表植被粗糙度,构建草原灌溉地的ET遥感计算模型,并通过ET遥感计算模型获取有效晴日的草原灌溉地区域蒸散发值ETd;
S3、通过草原灌溉地区域蒸散发值ETd计算草原灌溉地作物生长季的区域蒸散发值ET';
S4、采集草原灌溉地所在地域的有效降水数据,并利用空间插值法计算草原灌溉地的有效降水数据Pa;
S5、构建草原灌溉地水资源消耗量计算模型,并通过草原灌溉地作物生长季的区域蒸散发值ET'和有效降水数据Pa获取草原灌溉地水资源消耗量;
所述步骤S3具体为:
S3.1、利用FAO P-M公式空间插值获取草原灌溉地的逐日参考作物蒸散发值ET0;
S3.2、将草原灌溉地区域蒸散发值ETd与逐日参考作物蒸散发值ET0进行数据耦合,获取不同时段区域参考作物系数值ETrF的离散数据集;
S3.3、针对区域参考作物系数值ETrF在时间尺度具有稳定渐变特性,利用区域参考作物系数值ETrF的离散数据集,分段建立相邻晴日之间的区域参考作物系数值ETrF随儒略日DOY变化的函数关系式ETrF=f(DOY)具体为:
ETd,DOY=ETrF·ET0,DOY
ETrF=c·DOY+d
d=ETrFn-c·DOYn
S3.4、根据函数关系式ETrF=f(DOY),分析区域参考作物系数值ETrF在时间尺度长序列上的变化特征;
S3.5、根据区域参考作物系数值ETrF在时间尺度长序列上的变化特征,通过参考作物系数法获取草原灌溉地ET时间尺度长序列的扩展,得到草原灌溉地作物生长季的区域蒸散发值ET';
其中,ETd,DOY表示儒略日DOY对应当天的日蒸散发值,ET0,DOY表示儒略日DOY对应当天的日参考作物蒸散发值,c和d均为中间系数,ETrFn和ETrFm分别表示相邻两个晴日DOYn和DOYm的区域参考作物系数值。
进一步地,所述步骤S1中基础数据包括数字高程模型DEM、气象数据和遥感数据;所述步骤S1中对基础数据进行预处理的具体方法为:采用地理信息系统ARCGIS对数字高程模型DEM、气象数据和遥感数据进行预处理。
进一步地,所述步骤S2具体为:
S2.1、通过地表反照率和地表比辐射率获取地表净辐射量Rn,通过地表反照率、归一化植被指数、地表温度和地表净辐射量获取土壤热通量G,通过地表温度和地表植被粗糙度获取显热通量H;
S2.2、根据地表净辐射量Rn、土壤热通量G和显热通量H,构建草原灌溉地的ET遥感计算模型,获取潜热通量LE,所述潜热通量LE具体为:
LE=Rn-G-H
S2.3、以潜热通量LE为基础,并通过ET遥感计算模型获取有效晴日的草原灌溉地区域蒸散发值ETd具体为:
Rn24=(1-α)·Kin24-110·τsw
其中,Λ表示蒸发比,λ表示水的汽化热,Rn24表示日地表净辐射量,G24表示日土壤热通量,LE24表示日潜热通量,α表示地表反照率,Kin24表示日短波辐射,τsw表示大气透过率,表示相对日地距离订正系数,π表示圆周率,表示像元的地理纬度,δ表示太阳赤纬,Gsw表示太阳常数,ω2表示日落时的太阳时角。
进一步地,所述步骤S2.1中地表净辐射量Rn具体为:
Rn=(1-α)·Kin+(Lin-Lout)-(1-ε)·Lin
其中,α表示地表反照率,ε表示地表比辐射率,Kin表示短波辐射,Lin表示下行长波辐射,Lout表示上行长波辐射。
进一步地,所述步骤S2.1中土壤热通量G具体为:
其中,Ts表示地表温度,NDVI表示归一化植被指数。
进一步地,所述步骤S2.1中显热通量H具体为:
其中,ρair表示草原灌溉地空气密度,cp表示空气定压比热容,T1表示零平面位移以上高度z1处的温度,T2表示零平面位移以上高度z2处的温度,ra表示空气动力学阻抗,Pair表示空气压强。
进一步地,所述空气动力学阻抗ra具体为:
其中,z1表示草原灌溉地略高于植被冠层的平均高度,z2表示略低于边界层的参考高度,u*表示摩擦风速,κ表示常数;
所述摩擦风速u*具体为:
其中,ux表示高度x处的风速,zx表示高度为x的值,z0m表示地表植被粗糙度;
所述T1-T2具体为:
dT=a·Ts+b
其中,dT表示T1-T2的值,a和b均为中间参数,Rn干、G干、ra干、Ts干和ρ干分别表示在草原灌溉地遥感图像上选定的干点对应像元点的地表净辐射量、土壤热通量、空气动力学阻抗、地表温度和空气密度,Rn湿、G湿、ra湿、Ts湿和ρ湿分别表示在草原灌溉地遥感图像上选定的湿点对应像元点的地表净辐射量、土壤热通量、空气动力学阻抗、地表温度和空气密度,β湿=H/LE,β湿表示湿点波文比。
进一步地,所述步骤S5中草原灌溉地水资源消耗量GW具体为:
GW=ET'-Pa
其中,ET'表示草原灌溉地作物生长季的区域蒸散发值。
本发明的有益效果为:
(1)本发明针对遥感影像缺失条件下的遥感区域蒸散发ET扩展问题,利用区域参考作物系数在时间尺度上稳定渐变特性,通过对区域ET长序列扩展等价分解,提出改进参考作物系数法,耦合区域参考作物蒸散发,实现了遥感影像缺失特别是非晴日干扰下的草原区域蒸散发长序列扩展,提高了草原区域ET长序列计算的精度。
(2)本发明针对草原灌溉地水资源消耗量化问题,通过提出灌溉地水资源消耗计算模型,量化确定了草原地区灌溉地作物生长季的水资源消耗量GW值,界定了草原区地下水资源实时监控和管理的重点,为灌溉地植被建设对草原区地下水资源开发利用管理和草原生态系统影响研究提供了基础典型案例参考。
附图说明
图1为本发明提出的一种草原灌溉地水资源消耗量的计算方法流程图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
下面结合附图详细说明本发明的实施例。
如图1所示,一种草原灌溉地水资源消耗量的计算方法,包括以下步骤:
S1、采集草原灌溉地的基础数据并对基础数据进行预处理,根据预处理后的数据,使用METRIC模型反演得到地表反照率、地表比辐射率、地表温度、归一化植被指数和地表植被粗糙度;
S2、根据地表反照率、地表比辐射率、地表温度、归一化植被指数和地表植被粗糙度,构建草原灌溉地的ET遥感计算模型,并通过ET遥感计算模型获取有效晴日的草原灌溉地区域蒸散发值ETd;
S3、通过草原灌溉地区域蒸散发值ETd计算草原灌溉地作物生长季的区域蒸散发值ET';
S4、采集草原灌溉地所在地域的有效降水数据,并利用空间插值法计算草原灌溉地的有效降水数据Pa;
S5、构建草原灌溉地水资源消耗量计算模型,并通过草原灌溉地作物生长季的区域蒸散发值ET'和有效降水数据Pa获取草原灌溉地水资源消耗量;
所述步骤S3具体为:
S3.1、利用FAO P-M公式空间插值获取草原灌溉地的逐日参考作物蒸散发值ET0;
S3.2、将草原灌溉地区域蒸散发值ETd与逐日参考作物蒸散发值ET0进行数据耦合,获取不同时段区域参考作物系数值ETrF的离散数据集;
S3.3、针对区域参考作物系数值ETrF在时间尺度具有稳定渐变特性,利用区域参考作物系数值ETrF的离散数据集,分段建立相邻晴日之间的区域参考作物系数值ETrF随儒略日DOY变化的函数关系式ETrF=f(DOY)具体为:
ETd,DOY=ETrF·ET0,DOY
ETrF=c·DOY+d
d=ETrFn-c·DOYn
S3.4、根据函数关系式ETrF=f(DOY),分析区域参考作物系数值ETrF在时间尺度长序列上的变化特征;
S3.5、根据区域参考作物系数值ETrF在时间尺度长序列上的变化特征,通过参考作物系数法获取草原灌溉地ET时间尺度长序列的扩展,得到草原灌溉地作物生长季的区域蒸散发值ET';
其中,ETd,DOY表示儒略日DOY对应当天的日蒸散发值,ET0,DOY表示儒略日DOY对应当天的日参考作物蒸散发值,c和d均为中间系数,ETrFn和ETrFm分别表示相邻两个晴日DOYn和DOYm的区域参考作物系数值。
所述步骤S1中基础数据包括数字高程模型DEM、气象数据和遥感数据;所述步骤S1中对基础数据进行预处理的具体方法为:采用地理信息系统ARCGIS对数字高程模型DEM、气象数据和遥感数据进行预处理。
在本实施例中,所述步骤S1中地表反照率α具体为:
其中,b表示波段,ρs,b表示b波段的反射率,ωb表示b波段的权重系数,ρt,b表示b波段的大气外光谱反射率,ρa,b表示b波段的调整系数,τin,b表示b波段进入的大气透过率,τout,b表示b波段散射的大气透过率,π表示圆周率,Lb表示b波段的辐射亮度(W/m2·μm),表示相对日地距离订正系数,ESUNb表示大气外光谱辐照度,θref表示太阳天顶角,θhor表示太阳高度角,LMAXb表示b波段的最大辐射率,LMINb表示b波段的最小辐射率,Qcalmax表示最大像元灰度值,Qcalmin表示最小像元灰度值,Qcal表示影像的像元灰度值,C1表示第一中间系数,C2表示第二中间系数,C3表示第三中间系数,C4表示第四中间系数,C5表示第五中间系数,Pair表示空气压强,W表示水汽参数,K1表示0~1的常数,cosη=1。不同波段最大和最小辐射率如表1所示,Landsat卫星参数Ci和Wb如表2所示。
表1不同波段最大和最小辐射率单位:W/m2·μm·sr
表2 Landsat卫星参数Ci和Wb查询表
所述步骤S1中地表比辐射率ε具体为:
其中,LAI表示单位面积上植物叶面面积占土地面积的比值,SAVI表示土壤调整植被指数,ρi表示热红外波段i的反射率,ρi-1表示红外波段i-1的反射率,L表示常数;
所述步骤S1中地表温度Ts具体为:
C=ε·τsw
D=(1-τsw)·[1+(1-ε)·τsw]
其中,a6表示中间系数,C表示第一中间变量,D表示第二中间变量,b6表示中间系数,T6表示卫星高度上遥感器所观测到的亮度温度,Ta表示大气平均作用温度,τsw表示大气透过率,K1和K2均为常数,L6表示6波段的辐射亮度。
所述步骤S1中归一化植被指数NDVI具体为:
其中,ρi+1表示红外波段i+1的反射率;
所述步骤S1中地表植被粗糙度Z0m具体为:
Z0m=exp[(a1NDVI)+b1]
当坡度大于5度时,对地表植被粗糙度进行调整,调整如下:
其中,Z0m,ad表示调整后的地表植被粗糙度,a1表示第一经验系数,b1表示第二经验系数,s表示草原灌溉地坡度值。
所述步骤S2具体为:
S2.1、通过地表反照率和地表比辐射率获取地表净辐射量Rn,通过地表反照率、归一化植被指数、地表温度和地表净辐射量获取土壤热通量G,通过地表温度和地表植被粗糙度获取显热通量H;
S2.2、根据地表净辐射量Rn、土壤热通量G和显热通量H,构建草原灌溉地的ET遥感计算模型,获取潜热通量LE,所述潜热通量LE具体为:
LE=Rn-G-H
S2.3、以潜热通量LE为基础,并通过ET遥感计算模型获取有效晴日的草原灌溉地区域蒸散发值ETd具体为:
Rn24=(1-α)·Kin24-110·τsw
其中,Λ表示蒸发比,λ表示水的汽化热,Rn24表示日地表净辐射量,G24表示日土壤热通量,LE24表示日潜热通量,α表示地表反照率,Kin24表示日短波辐射,τsw表示大气透过率,表示相对日地距离订正系数,π表示圆周率,表示像元的地理纬度,δ表示太阳赤纬,Gsw表示太阳常数,ω2表示日落时的太阳时角。
所述步骤S2.1中地表净辐射量Rn具体为:
Rn=(1-α)·Kin+(Lin-Lout)-(1-ε)·Lin
其中,α表示地表反照率,ε表示地表比辐射率,Kin表示短波辐射,Lin表示下行长波辐射,Lout表示上行长波辐射。
在本实施例中,短波辐射Kin具体为:
W=0.14eaPair+2.1
其中,GSW表示太阳常数,取1367W/m2;dr2表示相对日地距离订正系数(无量纲);DOY表示儒略日,即遥感图像获取日期在太阳历中排列序号,例如1月1日的排列序号为1,2月1日的排列序号为32;δ表示太阳赤纬,表示像元的地理纬度,s表示研究区坡度,γ表示坡向(南为0,西为π/2,东为-π/2,北±π),ω表示太阳时角,Pair表示空气压强(kPa),ea表示水汽压(kPa),W表示水汽参数,τSW表示大气透过率。
所述太阳时角ω具体为:
ω=π/12[(t+0.0667·(Lz-Lm)+Sc-12)]
Sc=0.1645sin(2b)-0.1255cosb-0.025sinb
其中,t表示时段中点的时刻(hour),例如时段为14.00~15.00,t=14.5。Lz表示当地时区中心的经度(格林威治以西的度数),例如,Lz=75、90、105、120°分别为东、中心、落基山和太平洋的时区(美国),Lz=0°为格林威治,Lz=330°为开罗(埃及),Lz=225°为曼谷(泰国)。Lm表示测点经度(格林威治以西度数);Sc表示日照时间的季节修正(小时)当然,式(28)中的ω<-ωs或ω>ωs是表示太阳在地平线以下,定义Ra值为0。
所述下行长波辐射Lin和上行长波辐射Lout具体为:
Ta=16.011+0.92621·T0
其中,σ表示史蒂芬-波尔兹曼常数,取5.67*10-8W/m2·K4;Ta表示大气平均作用温度(K),T0表示是高度为2m的地面温度(K)。
所述步骤S2.1中土壤热通量G具体为:
其中,Ts表示地表温度,NDVI表示归一化植被指数。
所述步骤S2.1中显热通量H具体为:
其中,ρair表示草原灌溉地空气密度,cp表示空气定压比热容,T1表示零平面位移以上高度z1处的温度,T2表示零平面位移以上高度z2处的温度,ra表示空气动力学阻抗,Pair表示空气压强。
所述空气动力学阻抗ra具体为:
其中,z1表示草原灌溉地略高于植被冠层的平均高度(取0.01m),z2表示略低于边界层的参考高度(2m),u*表示摩擦风速,κ表示常数;
所述摩擦风速u*具体为:
其中,ux表示高度x处的风速,zx表示高度为x的值,z0m表示地表植被粗糙度;
所述T1-T2具体为:
dT=a·Ts+b
其中,dT表示T1-T2的值,a和b均为中间参数,Rn干、G干、ra干、Ts干和ρ干分别表示在草原灌溉地遥感图像上选定的干点对应像元点的地表净辐射量、土壤热通量、空气动力学阻抗、地表温度和空气密度,Rn湿、G湿、ra湿、Ts湿和ρ湿分别表示在草原灌溉地遥感图像上选定的湿点对应像元点的地表净辐射量、土壤热通量、空气动力学阻抗、地表温度和空气密度,β湿=H/LE,β湿表示湿点波文比。
在本实施例中,需要对摩擦风速u*和空气动力学阻抗ra进行更新,具体更新步骤为:
考虑表面受热导致近地层大气处于不稳定状态,模型应用Monin-Obukhov相似理论,引入大气热量传输与动量传输的稳定度订正因子Ψh、Ψm和L,对空气动力学阻抗进行校正后,迭代求解显热通量。
a.中性状态(L=0):
ψh=ψm=0
b.稳定状态(L>0):
c.不稳定状态(L<0):
求得L后,依据L的大小,根据稳定度订正因子Ψh和Ψm等参数,求得新的摩擦风速u*和空气动力学阻抗ra具体为:
根据新的空气动力学阻抗值,计算出的ai-1和bi-1值,获取Hi,进而得到新的参数ai和bi,比较ai和ai-1的相对误差,不满足精度,继续执行迭代计算,直到上一次求得的ai-1和本次求得的ai相对误差满足精度要求。
其中,g表示重力加速度,ψh表示大气热量传输与动量传输的第一稳定度订正因子,ψm表示大气热量传输与动量传输的第二稳定度订正因子,L表示第三大气热量传输与动量传输的稳定度订正因子,z表示高度,ψh(z)表示高度为z的第一稳定度订正因子,ψm(z)表示高度为z的第二稳定度订正因子,x(z)表示高度为z的中间变量。
在本实施例中,每一块灌溉用地作为独立的土地单元,将灌溉用地所处的空间垂向结构划分为大气层、灌溉用地层和地下水层,利用草原区域ET长序列扩展成果开展水资源消耗量化分析。结合灌溉用地水量均衡进行分析,灌溉用地在植物生长季的来水包括有效降水,以及地下水对灌溉用地的天然补给、灌溉抽取地下水的人工补给;灌溉用地的耗水包括作物蒸腾、土地蒸发、植物自身生长汲取消耗的水分(这部分水分相比前两者可忽略)、灌溉用地向地下水的渗漏量。因此,地下水资源补充了灌溉用地的包气带,但这部分水量未脱离灌溉用地进入大气中,此部分水量仍储存在灌溉用地所处的下垫面当中,属于当地水资源的内部消耗;除此之外,其余部分通过ET形式脱离了灌溉用地进入大气当中的水量,这部分水量是一种资源净消耗,定义此部分水量为地下水消耗量(以下简称GW)。换种角度解释,在没有外来水补给情况下,下垫面灌溉用地消耗的水资源量除了大气有效降水补给外,其余部分均来自地下水部分,区域ET扣除有效降水的部分定义为GW值。
所述步骤S5中草原灌溉地水资源消耗量GW具体为:
GW=ET'-Pa
其中,ET'表示草原灌溉地作物生长季的区域蒸散发值。
本发明的有益效果为:
(1)本发明针对遥感影像缺失条件下的遥感区域蒸散发ET扩展问题,利用区域参考作物系数在时间尺度上稳定渐变特性,通过对区域ET长序列扩展等价分解,提出改进参考作物系数法,耦合区域参考作物蒸散发,实现了遥感影像缺失特别是非晴日干扰下的草原区域蒸散发长序列扩展,提高了草原区域ET长序列计算的精度。
(2)本发明针对草原灌溉地水资源消耗量化问题,通过提出灌溉地水资源消耗计算模型,量化确定了草原地区灌溉地作物生长季的水资源消耗量GW值,界定了草原区地下水资源实时监控和管理的重点,为灌溉地植被建设对草原区地下水资源开发利用管理和草原生态系统影响研究提供了基础典型案例参考。
Claims (8)
1.一种草原灌溉地水资源消耗量的计算方法,其特征在于,包括以下步骤:
S1、采集草原灌溉地的基础数据并对基础数据进行预处理,根据预处理后的数据,使用METRIC模型反演得到地表反照率、地表比辐射率、地表温度、归一化植被指数和地表植被粗糙度;
S2、根据地表反照率、地表比辐射率、地表温度、归一化植被指数和地表植被粗糙度,构建草原灌溉地的ET遥感计算模型,并通过ET遥感计算模型获取有效晴日的草原灌溉地区域蒸散发值ETd;
S3、通过草原灌溉地区域蒸散发值ETd计算草原灌溉地作物生长季的区域蒸散发值ET';
S4、采集草原灌溉地所在地域的有效降水数据,并利用空间插值法计算草原灌溉地的有效降水数据Pa;
S5、构建草原灌溉地水资源消耗量计算模型,并通过草原灌溉地作物生长季的区域蒸散发值ET'和有效降水数据Pa获取草原灌溉地水资源消耗量;
所述步骤S3具体为:
S3.1、利用FAO P-M公式空间插值获取草原灌溉地的逐日参考作物蒸散发值ET0;
S3.2、将草原灌溉地区域蒸散发值ETd与逐日参考作物蒸散发值ET0进行数据耦合,获取不同时段区域参考作物系数值ETrF的离散数据集;
S3.3、针对区域参考作物系数值ETrF在时间尺度具有稳定渐变特性,利用区域参考作物系数值ETrF的离散数据集,分段建立相邻晴日之间的区域参考作物系数值ETrF随儒略日DOY变化的函数关系式ETrF=f(DOY)具体为:
ETd,DOY=ETrF·ET0,DOY
ETrF=c·DOY+d
d=ETrFn-c·DOYn
S3.4、根据函数关系式ETrF=f(DOY),分析区域参考作物系数值ETrF在时间尺度长序列上的变化特征;
S3.5、根据区域参考作物系数值ETrF在时间尺度长序列上的变化特征,通过参考作物系数法获取草原灌溉地ET时间尺度长序列的扩展,得到草原灌溉地作物生长季的区域蒸散发值ET';
其中,ETd,DOY表示儒略日DOY对应当天的日蒸散发值,ET0,DOY表示儒略日DOY对应当天的日参考作物蒸散发值,c和d均为中间系数,ETrFn和ETrFm分别表示相邻两个晴日DOYn和DOYm的区域参考作物系数值。
2.根据权利要求1所述的草原灌溉地水资源消耗量的计算方法,其特征在于,所述步骤S1中基础数据包括数字高程模型DEM、气象数据和遥感数据;所述步骤S1中对基础数据进行预处理的具体方法为:采用地理信息系统ARCGIS对数字高程模型DEM、气象数据和遥感数据进行预处理。
3.根据权利要求1所述的草原灌溉地水资源消耗量的计算方法,其特征在于,所述步骤S2具体为:
S2.1、通过地表反照率和地表比辐射率获取地表净辐射量Rn,通过地表反照率、归一化植被指数、地表温度和地表净辐射量获取土壤热通量G,通过地表温度和地表植被粗糙度获取显热通量H;
S2.2、根据地表净辐射量Rn、土壤热通量G和显热通量H,构建草原灌溉地的ET遥感计算模型,获取潜热通量LE,所述潜热通量LE具体为:
LE=Rn-G-H
S2.3、以潜热通量LE为基础,并通过ET遥感计算模型获取有效晴日的草原灌溉地区域蒸散发值ETd具体为:
Rn24=(1-α)·Kin24-110·τsw
4.根据权利要求3所述的草原灌溉地水资源消耗量的计算方法,其特征在于,所述步骤S2.1中地表净辐射量Rn具体为:
Rn=(1-α)·Kin+(Lin-Lout)-(1-ε)·Lin
其中,α表示地表反照率,ε表示地表比辐射率,Kin表示短波辐射,Lin表示下行长波辐射,Lout表示上行长波辐射。
7.根据权利要求6所述的草原灌溉地水资源消耗量的计算方法,其特征在于,所述空气动力学阻抗ra具体为:
其中,z1表示草原灌溉地略高于植被冠层的平均高度,z2表示略低于边界层的参考高度,u*表示摩擦风速,κ表示常数;
所述摩擦风速u*具体为:
其中,ux表示高度x处的风速,zx表示高度为x的值,z0m表示地表植被粗糙度;
所述T1-T2具体为:
dT=a·Ts+b
其中,dT表示T1-T2的值,a和b均为中间参数,Rn干、G干、ra干、Ts干和ρ干分别表示在草原灌溉地遥感图像上选定的干点对应像元点的地表净辐射量、土壤热通量、空气动力学阻抗、地表温度和空气密度,Rn湿、G湿、ra湿、Ts湿和ρ湿分别表示在草原灌溉地遥感图像上选定的湿点对应像元点的地表净辐射量、土壤热通量、空气动力学阻抗、地表温度和空气密度,β湿=H/LE,β湿表示湿点波文比。
8.根据权利要求7所述的草原灌溉地水资源消耗量的计算方法,其特征在于,所述步骤S5中草原灌溉地水资源消耗量GW具体为:
GW=ET'-Pa
其中,ET'表示草原灌溉地作物生长季的区域蒸散发值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011469032.8A CN112700089B (zh) | 2020-12-15 | 2020-12-15 | 一种草原灌溉地水资源消耗量的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011469032.8A CN112700089B (zh) | 2020-12-15 | 2020-12-15 | 一种草原灌溉地水资源消耗量的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112700089A true CN112700089A (zh) | 2021-04-23 |
CN112700089B CN112700089B (zh) | 2022-05-17 |
Family
ID=75509107
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011469032.8A Active CN112700089B (zh) | 2020-12-15 | 2020-12-15 | 一种草原灌溉地水资源消耗量的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112700089B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113836728A (zh) * | 2021-09-27 | 2021-12-24 | 广州极飞科技股份有限公司 | 灌水量计算方法及装置、电子设备及存储介质 |
CN115841628A (zh) * | 2022-12-12 | 2023-03-24 | 中国水利水电科学研究院 | 一种高精度灌区耗水信息自动解译系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110599360A (zh) * | 2019-09-20 | 2019-12-20 | 福州大学 | 一种干旱区农作物生长季蒸散发高分辨率遥感估算方法 |
US20200214231A1 (en) * | 2017-07-02 | 2020-07-09 | Manna Irrigation Ltd. | Methods and systems for irrigation guidance |
CN111553459A (zh) * | 2020-03-26 | 2020-08-18 | 北京大学 | 一种基于遥感信息的灌溉地区逐日实际蒸散发估算方法 |
-
2020
- 2020-12-15 CN CN202011469032.8A patent/CN112700089B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200214231A1 (en) * | 2017-07-02 | 2020-07-09 | Manna Irrigation Ltd. | Methods and systems for irrigation guidance |
CN110599360A (zh) * | 2019-09-20 | 2019-12-20 | 福州大学 | 一种干旱区农作物生长季蒸散发高分辨率遥感估算方法 |
CN111553459A (zh) * | 2020-03-26 | 2020-08-18 | 北京大学 | 一种基于遥感信息的灌溉地区逐日实际蒸散发估算方法 |
Non-Patent Citations (8)
Title |
---|
JUN WANG 等: "An estimation of the evapotranspiration of typical steppe areas using Landsat images and the METRIC model", 《JOURNAL OF WATER AND CLIMATE CHANGE》 * |
周剑 等: "改进SEBS模型评价黑河中游灌溉水资源利用效率", 《水利学报》 * |
夏浩铭 等: "遥感反演蒸散发时间尺度拓展方法研究进展", 《农业工程学报》 * |
张振宇 等: "干旱区膜下滴灌向日葵农田蒸散发特征", 《中国生态农业学报(中英文)》 * |
曾丽红 等: "2000年至2008年松嫩平原生长季蒸散量时空格局及影响因素分析", 《资源科学》 * |
曾丽红 等: "松嫩平原不同地表覆盖蒸散特征的遥感研究", 《农业工程学报》 * |
朱明承 等: "基于遥感双层模型的宝鸡峡灌区日蒸散发估算", 《北京师范大学学报(自然科学版)》 * |
王军 等: "基于遥感技术的区域蒸散发计算方法综述", 《节水灌溉》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113836728A (zh) * | 2021-09-27 | 2021-12-24 | 广州极飞科技股份有限公司 | 灌水量计算方法及装置、电子设备及存储介质 |
CN113836728B (zh) * | 2021-09-27 | 2024-04-26 | 广州极飞科技股份有限公司 | 灌水量计算方法及装置、电子设备及存储介质 |
CN115841628A (zh) * | 2022-12-12 | 2023-03-24 | 中国水利水电科学研究院 | 一种高精度灌区耗水信息自动解译系统 |
CN115841628B (zh) * | 2022-12-12 | 2023-09-19 | 中国水利水电科学研究院 | 一种高精度灌区耗水信息自动解译系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112700089B (zh) | 2022-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Los et al. | A global 9-yr biophysical land surface dataset from NOAA AVHRR data | |
Xu et al. | Assessing the impact of urbanization on regional net primary productivity in Jiangyin County, China | |
Ayenew | Evapotranspiration estimation using thematic mapper spectral satellite data in the Ethiopian rift and adjacent highlands | |
CN102650587B (zh) | 基于sebal-hj模型的农作物生物量反演方法 | |
Liu et al. | Estimation of evapotranspiration in the Mu Us Sandland of China | |
CN112700089B (zh) | 一种草原灌溉地水资源消耗量的计算方法 | |
Liu et al. | Mapping incident photosynthetically active radiation from MODIS data over China | |
Samuel et al. | Spatio-temporal variability of evapotranspiration and crop water requirement from space | |
Chatterjee et al. | Characterization of land surface energy fluxes in a tropical lowland rice paddy | |
CN110599360A (zh) | 一种干旱区农作物生长季蒸散发高分辨率遥感估算方法 | |
CN111553459B (zh) | 一种基于遥感信息的灌溉地区逐日实际蒸散发估算方法 | |
Wang et al. | Estimation of hourly actual evapotranspiration over the Tibetan Plateau from multi-source data | |
Chen et al. | Quantitative identification of temporal-spatial variations of urban heat island (UHI) effects in Changchun, China | |
CN110658325B (zh) | 一种确定sm监测中lst/fvc空间的方法 | |
CN113642191B (zh) | 一种基于短波红外的遥感蒸散模型构建方法 | |
CN115481943A (zh) | 一种生态系统恢复力的评价方法 | |
Parmar et al. | Estimation of Evapotranspiration using Remote Sensing and Surface Energy Balance Algorithm for Land (SEBAL) in Canal Command | |
Haugland et al. | The diurnal cycle of land–atmosphere interactions across Oklahoma’s winter wheat belt | |
Gao et al. | A coupled remote sensing and the Surface Energy Balance with Topography Algorithm (SEBTA) to estimate actual evapotranspiration under complex terrain | |
Li et al. | Surface temperature adjustment in METRIC model for monitoring crop water consumption in North China Plain | |
Zaeen | Using remote sensing techniques to monitoring and evaluate the water cover in AL_Razzaza lake: Iraq at deferent periods | |
Kumar et al. | Estimation of evapotranspiration using MODIS sensor data in Udupi District of Karnataka, India | |
Lamers | Soil slaking and the possibilities to record with infrared line scanning | |
Zhong | Land-atmosphere interactions and their relationships to the Asian monsoon in the Tibetan Plateau | |
Fu et al. | Land-atmosphere Interaction Observational Experiment of in China’s Summer Monsoon Transition Zone: a Case Study in Dingxi |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |