一种农田净灌溉用水模型及灌溉用水量估算方法
阅读：862发布：2020-06-08
专利汇可以提供一种农田净灌溉用水模型及灌溉用水量估算方法专利检索，专利查询，专利分析的服务。并且本 发明 涉及生态 水 文学科学领域，尤其涉及一种农田净 灌溉 用水模型。模型包括中央 数据处理 器、植被截留模 块 、径流模块、蒸散发模块和 土壤 含水量 模块，其中，中央数据处理器分别与植被截留模块、径流模块、蒸散发模块和土壤含水量模块有数据传输。本发明充分考虑大尺度的地表不均一性,将空间离散化为像元，在像元尺度上反演农田净灌溉用水模型中各子模块所需的部分参数，耦合生态水文过程模型，进行降水植被截留、土壤含水量、径流和蒸散发全过程模拟，实现了对农田土壤含水量的实时监控，并能及时了解农田是否需要灌溉，或灌溉所需用水量，避免农田得不到及时灌溉或者灌溉量过大的情况，从而造成农业生产损失和水资源浪费。，下面是一种农田净灌溉用水模型及灌溉用水量估算方法专利的具体信息内容。
1.一种农田净灌溉用水模型，其特征在于，包括中央数据处理器，处理各部分数据；植被截留模块，计算植被截留量；径流模块，计算径流量；蒸散发模块，计算蒸散发量；土壤含水量模块，模拟土壤含水量，其中，中央数据处理器分别与植被截留模块、径流模块、蒸散发模块和土壤含水量模块有数据传输。
2.根据权利要求1所述的农田净灌溉用水模型，其特征在于，所述的植被截留模块采用半经验半理论模型进行计算。
3.根据权利要求1所述的农田净灌溉用水模型，其特征在于，所述的蒸散发模块包括实际蒸散发模块和潜在蒸散发模块。
4.根据权利要求1所述的农田净灌溉用水模型，其特征在于，所述的土壤含水量模块包括根系吸水模块、层间土壤水分增量模块和根系层补给量模块。
5.一种利用农田净灌溉用水模型实现估算农田灌溉用水量的方法，其特征在于，包括以下步骤：
S1：获取下垫面参数和气象参数；
S2：分别向植被截留模块、蒸散发模块和径流模块中输入下垫面参数和气象参数；
S3：将计算出的植被截留量、蒸散发量和径流量传输给中央数据处理器；
S4：中央数据处理器根据各模块传输来的数据控制土壤含水量模块进行土壤含水量模拟，
S5：根据模拟结果计算出农田灌溉用水量。
6.根据权利要求5所述的方法，其特征在于，所述的步骤S1中的下垫面参数包括植被类型、植被盖度、叶面积指数、田间持水量Wf、萎蔫系数Wm、各层土壤含水量和表层土壤厚度，气象参数包括降水量、气温、长波辐射、短波辐射、地表反照率和大气压。
7.根据权利要求5所述的方法，其特征在于，所述的步骤S4在计算时需向中央数据处理器输入降水量。
一种农田净灌溉用水模型及灌溉用水量估算方法
技术领域
背景技术
[0003] 净灌溉用水是指在一定时间内，为满足作物正常生长，除降水补给外需要的灌溉水量。当土壤含水量不能满足作物正常生长的需要，需要进行灌溉补给。按照水量平衡原理，观测期内植被降水截留量、有效作物蒸散发、径流量和观测始末土壤水分差值之和为观测期内降水总量和灌溉总量之和，进一步可得到观测期内需要净灌溉用水量。
[0004] 下垫面是指地球表面，包括海洋、陆地、陆地上的高原、山地、平原、森林、草原以及城市等。下垫面各部分温度、水分以及表面形状等参数均有较大差异，从而导致下垫面具有非均一性。
[0005] 现有的净灌溉用水量估算方法中，大多根据区域类型估算农田蒸散发量，并结合降水、径流、土壤含水量等实地观测与统计资料得到，或者直接根据观测期内一定面积农田统计的灌溉用水量得到，简单、实用。
[0006] 通常情况下，区域观测尺度较小时，认为气象参数、下垫面均一、农田管理方式等相差不明显。但是对于大的区域或流域，上述方法由于未充分考虑地表的非均一性而造成不能准确地估算特定区域空间分布状况，具有较大的局限性。
[0007] 综合上述，提出一种既考虑了下垫面的非均一性且精确度又高的大尺度像元尺度的农田净灌溉用水量估算方法是急需解决的问题。
发明内容
[0008] （一）要解决的技术问题
[0009] 本发明要解决的技术问题是：提供一种考虑了下垫面的非均一性且精确度又高的大尺度像元尺度的农田净灌溉用水模型，以及利用该模型实现的农田净灌溉用水量估算方法。
[0010] （二）技术方案
[0011] 为了解决上述技术问题，本发明提供了一种农田净灌溉用水模型，其包括中央数据处理器，处理各部分数据；植被截留模块，计算植被截留量；径流模块，计算径流量；蒸散发模块，计算蒸散发量；土壤含水量模块，模拟土壤含水量及灌溉量。
[0012] 其中，中央数据处理器分别与植被截留模块、径流模块、蒸散发模块和土壤含水量模块有数据传输。
[0013] 进一步的，植被截留模块采用半经验半理论模型进行计算。
[0014] 进一步的，蒸散发模块包括实际蒸散发模块和潜在蒸散发模块。
[0015] 进一步的，土壤含水量模块包括根系吸水模块、层间土壤水分增量模块和根系层补给量模块。
[0016] 本发明还提供了一种利用农田净灌溉用水模型实现估算农田灌溉用水量的方法，其包括以下步骤：
[0017] S1：获取下垫面参数和气象参数；
[0018] S2：分别向植被截留模块、蒸散发模块、径流模块和土壤含水量模块中输入下垫面参数和气象参数；
[0019] S3：将计算出的植被截留量、蒸散发量和径流量传输给中央数据处理器；
[0020] S4：中央数据处理器根据各模块传输来的数据控制土壤含水量模块进行土壤含水量模拟，
[0021] S5：根据模拟结果计算出农田灌溉用水量。
[0022] 其中，步骤S1中的下垫面参数包括植被类型、植被盖度、叶面积指数、田间持水量Wf、萎蔫系数Wm、各层土壤含水量和表层土壤厚度，气象参数包括降水量、气温、长波辐射、短波辐射、地表反照率和大气压。
[0023] 其中，步骤S4在计算时需向中央数据处理器输入降水量。
[0024] （三）有益效果
[0025] 本发明的上述技术方案具有如下优点：充分考虑大尺度的地表不均一性,将空间离散化为像元，在像元尺度上反演农田净灌溉用水模型中各子模块所需的部分参数，耦合生态水文过程模型，进行降水植被截留、土壤含水量、径流和蒸散发全过程模拟，实现了对农田土壤含水量的实时监控，并能及时了解农田是否需要灌溉，或灌溉所需用水量，避免农田得不到及时灌溉或者灌溉量过大的情况，从而造成农业生产损失和水资源浪费。附图说明
[0026] 图1为本发明模型结构示意图。
[0027] 图2为本发明净灌溉用水量推导过程图。
[0028] 图3为本发明植被截留量退出过程图。
[0029] 图4为本发明径流量推导过程图。
[0030] 图5为本发明蒸散发量推导过程图。
具体实施方式
[0031] 下面结合附图和实施例对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
[0032] 如图1所示，本发明提供了一种农田净灌溉用水模型，其包括中央数据处理器、植被截留模块、径流模块、蒸散发模块和土壤含水量模块。
[0033] 其中，中央数据处理器分别与植被截留模块、径流模块、蒸散发模块和土壤含水量模块有数据传输。
[0034] 净灌溉用水是指在一定时间内，为满足作物正常生长，除降水补给外需要的灌溉水量。水量交换侧重垂直方向，根据北方流域地下水位情况，在研究中忽略地下水对农田根系层的补给作用。农田净灌溉用水量I用如下公式表示：
[0035] I=Sv+R+ETa+ΔW-P
[0036] 式中，Sv为期内累计植被截留量，R为期内径流量，ETa为期内累计蒸散发量，ΔW为期内末期和起始期土壤含水量差值，P为期内累计降水量。
[0037] 由公式可知，如图2所示，要计算出净灌溉用水量就须知道植被截留量、径流量、蒸散发量、土壤含水量和降水量。除降水量可以通过遥感或者通过监测站获取以外，其他数据均需要计算才能得出。
[0038] 植被截留量可以通过植被截留模块计算得出。
[0040] 植被截留模块可以为经验模型、半经验半理论模型或理论模型，其中以Rutter模型和Gash解析模型较为完善且被广泛应用，但是模型的参数较难获取，因此，本发明采用半经验半理论模型进行计算。
[0041] 降雨过程中植被对降水的截留主要有枝叶作用，因此截留量主要与植被特征有关，如植被盖度、叶面积、叶面粗糙度，随植被类型和生长时段而变化。植被截留量Sv的计算方程为：
[0042]
[0043] 式中：cv为植被盖度，反映覆被空间分布情况；P为累积降水量，Smax为植被最大截留量；η为校正系数。
[0044] 植被最大截留量主要取决于叶面特征，即叶面积和叶面粗糙度。基于叶面积指数LAI可以计算出植被最大截留量，其公式为：
[0045] Smax=0.935+0.498×LAI-0.00575×LAI2
[0046] η=0.046×LAI
[0047] 综上所述，可得出植被截留量在植被截留模块中输入参数降雨量、叶面积指数和植被盖度后可以计算出，如图3所示。
[0048] 遥感技术在这些参数的获取上均已有较为成熟的应用，可以提供较高时间分辨率的LAI产品，能为模型提供驱动数据。
[0049] 径流量可以通过径流模块进行计算得出。
[0051] 降雨经林冠截留后到达表层土壤，一部分在表面发生地表径流，剩余部分进入表层土壤补充土壤含水量。地表径流采用VGTM模型计算，模型认为地表径流与有效降雨呈非线性关系，这种关系通过时变增益因子表示，时变增益因子通过表层土壤水分求算，同时增加下垫面覆被状况对产流状况的影响通过地表覆被因子C来表达，其主要与植被生长状况有关。地表径流单元计算公式如下：
[0052]
[0053] 式中： 为地表径流量，AWu为表层土壤含水量，WMu为表层土壤饱和含水量，(P-Sv)为有效降雨量，为降雨量扣除植被截留后到达地面的净雨量；g1与g2是时变增益因子的有关参数（0
[0054] 扣除地表径流外的降雨以一定的速率下渗进入土壤，一部分保存在土壤中，使土壤含水量增大；一部分沿着土壤空隙流动，从坡侧土壤空隙流出，转换为地表径流，注入河槽的部分径流称为表层壤中流，采用自由水蓄水库线性出流计算，表层壤中流单元计算公式如下：
[0055]
[0056] 式中： 为表层壤中流量，AWu为表层土壤含水量，Kr为土壤水出流系数；Thick为表层土壤厚度。
[0057] 表层土壤水分在重力和水势作用下向深层下渗，给定表层到深层的下渗率为fc，即可求算上层土壤渗入到下层的水量，地下径流亦相同。深层壤中流和地下径流计算方式同于表层壤中流，即
[0058]
[0059]
[0060] 式中： 和 分别为深层壤中流量和地下径流量，AWd和AWg分别为深层土壤含水量和地下层含水量，Kd和Kg分别深层土壤流出流系数和地下径流出流系数，地下径流一般出流较为稳定，在冰期，径流的主要供给为地下流，因此可根据冬季径流量确定出流系数。
[0061] 由图4和上述公式可知，在径流模块中输入表层土壤含水量AWu、表层土壤饱和含水量MWu、表层土壤厚度Thick、深层土壤含水量AWd、地下层含水量AWg、降水量P以及在植被截留模块得出的植被截留量Sv就可以得出径流量R，这些数据均可以通过遥感方式获取，如降水量P等还可以通过其他方式获取。
[0062] 蒸散发量ETa的计算是通过蒸散发模块完成的，如图5所示为蒸散发量ETa推导过程图。
[0064] ETa=Ks×ETp
[0065] 土壤水分胁迫系数Ks可以利用下列公式确定：
[0066] Ks=ln(Av+1)/ln(101)
[0067] Av=[(W-Wm)/(Wf-Wm)]×100
[0068] 式中，W为即时含水量，Wf为田间持水量，Wm为萎蔫系数。
[0069] 其中，田间持水量Wf与萎蔫系数Wm均与土壤的类别有关，其均可以用常规方法进行估算。
[0070] 在蒸散发模块中，潜在蒸散发量ETP由Priestley-Taylor公式计算得出，其公式为：
[0071]
[0073] 由上述公式可以看出，在计算潜在蒸散发量ETP时，需知道地表净辐射Rn、土壤热通量G、汽化潜热λ、饱和水气压-温度曲线斜率Δ和干湿表常数γ。
[0074] 地表净辐射Rn又称辐射平衡或辐射差额，是指地表面净得的短波辐射RS与长波辐射RL之和，即地表辐射能量收支的差额。它是地表面能量、动量、水分输送与交换过程中的主要能源，是地气能量交换中最重要的组分，也是驱动地表蒸散和显热通量的主要来源。
[0075] 在蒸散发模块中，地表净辐射采用MODIS数据估算卫星过境时刻的瞬时地表净辐射Rn：
[0076] Rn=（1-α）RS+RL
[0077] 式中，α为地表反照率，RS为短波辐射，RL为长波辐射。
[0078] 地表反照率α、短波辐射RS和长波辐射RL均可以通过遥感方式获取数据，故蒸散发模块可以通过上述公司计算出瞬时地表净辐射。
[0079] 瞬时地表净辐射Rn计算出以后，蒸散发模块再对瞬时地表净辐射Rn进行积分运算，就可以得出所需要的地表净辐射Rn。
[0080] 地表净辐射能量一部分被用作蒸散发，一部分被用来加热地表的大气，剩下的则被储藏在土壤或水体之中，这部分土壤或水体的热交换能量即土壤热通量。目前土壤热通量G还不能直接通过遥感技术获取，一般通过它与地表净辐射Rn及植被盖度cv的关系来确定，在蒸散发模块中的计算公式如下：
[0081] G=Rn[Γc+(1-cv)(Γs+Γc)]
[0082] 式中：Γc为土壤热通量与地表净辐射的比值，在全植被覆盖下Γc=0.05，Γs为土壤热通量与地表净辐射比值，在裸地情况下Γs=0.315，cv为植被盖度。
[0083] 对于水体和冰雪，土壤热通量G的计算多采用与地表净辐射Rn取比值的方法，本发明中，水体、积雪与冰川的土壤热通量G与地表净辐射Rn的比值取0.5，即认为进入水体和冰雪的地表净辐射Rn有一半被吸收成为冰雪的G：
[0084] G=0.5Rn
[0085] 由上述可知，土壤热通量G的计算是通过地表净辐射Rn和植被盖度cv完成的。
[0086] 汽化潜热λ在蒸散发模块中是通过气温Ta进行估算：
[0087] λ=2.5-0.0022Ta
[0088] 而饱和水气压-温度曲线斜率Δ在蒸散发模块中是通过气温Ta进行计算：
[0089]
[0090] 干湿表常数γ在蒸散发模块中是通过下列公式进行计算：
[0091]
[0092] 式中，Cp为空气定压比热，取值为1.013×10-3，Pr为大气压，ε是水汽分子量和干空气分子量的比值，取值为0.622；λ为汽化潜热。
[0093] 以上可以得出瞬时的蒸散发量ETa，但是蒸散发量ETa在每一时间是不一样的，在得出瞬时蒸散发量ETa以后，再进行积分运算，就可以得出所需要的蒸散发量。
[0094] 由上述可知，在蒸散发模块中输入田间持水量Wf、萎蔫系数Wm、长波辐射RL、短波辐射RS、地表反照率α、植被盖度cv、大气压Pr和气温Ta就可以得出蒸散发量ETa，而上述数据均可以通过遥感进行获取，部分数据还可以通过其他方式获取。
[0095] 土壤含水量模块可以进行土壤含水量的模拟，以确定土壤内的水分是否对足够生物生长需要。
[0096] 利用遥感获取初始土壤含水量和检测时的土壤含水量之后，就可以对农田净灌溉用水进行计算了。
[0097] 综上所述，可以得出利用农田净灌溉用水模型实现估算农田灌溉用水量的方法，其包括以下步骤：
[0098] S1：获取下垫面参数和气象参数；
[0099] S2：分别向植被截留模块、蒸散发模块、径流模块和土壤含水量模块中输入下垫面参数和气象参数；
[0100] S3：将计算出的植被截留量SV、蒸散发量ETa和径流量R传输给中央数据处理器；
[0101] S4：中央数据处理器根据各模块传输来的数据控制土壤含水量模块进行土壤含水量模拟，
[0102] S5：根据模拟结果计算出农田灌溉用水量I。
[0103] 其中，步骤S1中的下垫面参数包括植被类型、植被盖度cv、叶面积指数LAI、田间持水量Wf、萎蔫系数Wm、各层土壤含水量W和表层土壤厚度Thick，气象参数包括降水量P、气温Ta、长波辐射RL、短波辐射RS、地表反照率α和大气压Pr。
[0104] 其中，步骤S4在计算时需向中央数据处理器输入降水量P。
[0105] 需要用遥感获取的数据，如果有其他方式可以获取，也可以用其他方式，如气温Ta可以用温度计，降水量P可以在气象监测站获取等。
[0106] 本发明充分考虑大尺度的地表不均一性,将空间离散化为像元，通过逐个像元反演水平衡方程各子模型所需的部分参数，耦合生态水文过程模型，并进行降水植被截留、土壤含水量、径流和蒸散发全过程模拟，实现了对农田土壤的含水量进行实时监控，能及时知道农田是否需要灌溉，或需要灌溉时的用水量，从而避免农田没有及时灌溉而造成损失或灌溉用水量太大而受到损失和水资源浪费。
[0107] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变型，这些改进和变型也应视为本发明的保护范围。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 基于人水耦合互馈系统动力学演化模型的预测方法和系统 | 2020-05-23 | 4 |
| 一种小区降雨实验径流与水质综合指标的测量方法 | 2020-05-23 | 21 |
| 一种城市与山地水文学实验模型装置 | 2020-05-11 | 142 |
| 一种旱涝急转标准的计算方法 | 2020-05-26 | 471 |
| 一种基于下垫面特征的重配置产流模拟方法 | 2020-05-27 | 118 |
| 一种冰封期河流生态基流计算方法 | 2020-05-20 | 88 |
| 一种时空动态优化的嵌套式自稳定汇流演算方法 | 2020-05-22 | 522 |
| 一种干旱半干旱区高山融雪径流观测方法 | 2020-05-15 | 795 |
| 一种自动控制闸门泄流量的方法和系统 | 2020-05-12 | 489 |
| 一种城市与山地水文学实验模型装置 | 2020-05-11 | 582 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
分析报告
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
QQ群二维码
意见反馈
意见反馈