一种降雨入渗补给地下水临界埋深计算方法技术

本发明专利技术涉及一种降雨入渗补给地下水临界埋深计算方法，包括：饱和入渗深度计算；饱和入渗深度之下非饱和入渗深度计算；水力传导计算分层内的湿润锋推移速度；判断推移速度是否大于零；计算土壤含水量；计算非饱和入渗深度；计算入渗深度。本发明专利技术通过分析降雨入渗补给地下水的作用机理，阐述降雨入渗补给地下水的物理过程，并定义了降雨入渗补给地下水的临界埋深；采用饱和与非饱和下渗理论，以湿润锋推移距离为计算目标，构建了降雨入渗深度计算模型确定地下水临界埋深。所述方法物理过程与作用机制明确，对地下水超采引发的可再生机制的破坏具有通用性，为指导地下水资源合理开发利用与地下水含水层保护提供可靠的理论与技术支撑。

【技术实现步骤摘要】

本专利技术涉及，是一种生态水文计算方 法，一种以降雨入渗补给深度为控制性指标的地下水临界埋深计算方法。
技术介绍
对半干旱区与半湿润区而言，确定入渗补给的地下水位，是维持地下水含水层稳 定的关键，对于保护地下水含水层再生机制意义重大。许多农区以井灌为主，而地下水补给 主要来源于降雨和灌溉回归水。由于地下水长期开发利用，地下水位下降明显，埋深增大， 地下含水层与农作物直接的水力联系已被割断。为了保障地下水含水层稳定，需要维持一 个科学合理的入渗补给地下水埋深，也即采补平衡水位，这是地下水可持续利用的核心。目 前多利用水井动态观测以及进行打井实验，通过经验分析方法确定采补平衡水位。这种方 法简单直观，但因为缺乏作用机理和物理过程，且带有主观偏好，难于推广应用。
技术实现思路
为了进行理性客观的技术分析，本专利技术提出了一种降雨入渗补给地下水临界埋深 计算方法。所述的方法以湿润锋的推进速度为核心，以湿润锋运动速率趋近于零的位置定 义为入渗临界埋深，通过构建降雨入渗深度计算模型，确定地下水临界埋深。 本专利技术的目的是这样实现的：，所述 方法的步骤如下： 饱和入渗深度计算的步骤：用于对降雨入渗阶段土壤中的水分饱和入渗深度进行计 算，计算公式为：其中： HP 为饱和入渗深度，单位：米， F--为降雨入渗量，单位：晕米， 一一土壤饱和含水量，cm3/cm3或用百分比表示， 0i--土壤初始含水量，cm3/cm3或用百分比表示； 饱和入渗深度之下非饱和入渗深度计算的步骤：对土壤分层，根据湿润锋推进速度，对 土壤水再分配阶段非饱和入渗深度进行计算； 水力传导计算分层内的湿润锋推移速度的步骤：用于通过公式计算第i层的湿润锋推 移速度：其中： K{^ :)一一为第i减1层的土壤含水量为奋时的非饱和水力传导度， K( 9,)一一为第i层的土壤含水量为#时的非饱和水力传导度， 一一为第i层的基质势对应的当量水头，hm>1t一一为第i减1层的基质势％对应的当量水头， Az--为各层土壤的厚度； 判断推移速度是否大于零的步骤：用于判断湿润锋推移速度是否大于零，如果大于零 则进入"计算土壤含水量的步骤"，如果小于等于零则进入"计算非饱和入渗深度的步骤"； 计算土壤含水量的步骤：用于根据公式计算当前时刻湿润锋上部土壤含水量：其中： 0t--甜寸刻湿润锋上部的土壤含水量， 〇tl--卜1时刻湿润锋上部的土壤含水量， Htj 为卜1时刻湿润层厚度， 计算结束后，将i加1，并回到"饱和入渗深度之下非饱和入渗深度计算的步骤"； 计算非饱和入渗深度的步骤：用于将各分层土壤厚度相加SAz，形成非饱和入渗深 度； 计算入渗深度的步骤：用于将饱和入渗深度尽与非饱和入渗深度SAz相加，即为降 雨入渗补给地下水临界埋深。 进一步的，所述的"计算分层内的湿润锋推移速度的步骤"中的"湿润锋面的土壤 水势"计算经验公式为：其中： h-为土壤吸力，cmH20 ; 9s--为土壤饱和含水率； K一一为土壤凋萎含水率； 9--土壤实际含水率； a、in、n--为参数，其中，'猶:政:1:-:#。 .将. 进一步的，所述的基质势与非饱和土壤水力传导度的经验关系为： K(e) =Ks(e/f)b 其中： K( 0)一一为非饱和水力传导度， 9一一为土壤含水量， Ks一一为饱和水力传导度， f-一为土壤孔隙度，b-一为经验系数。 本专利技术产生的有益效果是：本专利技术通过分析降雨入渗补给地下水的作用机理，阐 述入渗补给地下水的物理过程，并以机理描述定义降雨入渗补给地下水的临界埋深；采用 饱和与非饱和下渗理论，以湿润锋推移距离为计算目标，构建了降雨入渗深度计算模型确 定地下水临界埋深。所述方法物理过程与作用机制明确，为地下水含水层保护进而指导地 下水资源合理开发利用提供可靠的理论与技术支撑。【附图说明】 下面结合附图和实施例对本专利技术作进一步说明。 图1是地表水入渗过程与土壤水分分布和本专利技术实施例一的土壤分层示意图； 图2是本专利技术的实施例一所述计算方法的流程图。【具体实施方式】 实施例一： 本实施例是。 根据水量下渗原理和补给地下水临界定义，当下渗过程完成，湿润锋消失所处的 位置为下渗补给地下水的临界埋深。实际发生的下渗补给地下水位应该在临界埋深之上。 如图1所示，入渗可能达到的入渗最大深度为临界埋深，当潜水埋深大于最大入渗深度时， 入渗达不到含水层即终止，不能实现对地下水的补给。而在最大深度之内任何地下水埋深 都可以完成入渗补给。确定临界埋深的最重要的两个物理量是供水量(降雨量+回归水量） 和土壤初始含水量。 根据土壤水下渗理论，确定降雨入渗补给地下水的临界埋深，关键是分析降雨入 渗和入渗后土壤水的再分配过程。本实施例结合饱和与非饱和入渗理论，将入渗划分为降 雨入渗和土壤水再分配两个阶段，构建降雨入渗深度的计算模型。在降雨入渗阶段，采用饱 和土壤水运动方法；在土壤水再分配阶段，采用非饱和土壤水运动方法。再分配阶段，通过 土壤分层，计算土壤水再分配阶段的入渗深度。计算流程如图2所示，具体计算步骤如下： 饱和入渗深度计算的步骤：用于对降雨入渗阶段土壤中的水分饱和入渗深度进行计 算，计算公式为：其中： HP一一为饱和入渗过程中湿润锋位置，单位：米， F--为降雨入渗量，单位：晕米， 〇s-一土壤饱和含水量，cm3/cm3或用百分比表示， 6,--土壤初始含水量，cm3/cm3或用百分比表示。 本步骤计算的是土壤在降雨入渗阶段，采用饱和入渗理论，计算降雨期间土壤含 水量的改变，进而推断入渗土壤层的厚度，考虑西辽河平原土壤质地，饱和入渗的边界条件 取土壤的田间持水量，即降雨自表层向下逐层补充土壤含水量，当上层土壤含水量达到田 间持水量时，则继续向下层补充，直至实现上层均达田间持水量的深度。 饱和入渗深度之下非饱和入渗深度计算的步骤：对土壤分层，用于计算湿润锋推 进速度，据此推求非饱和入渗深度。分层只是人为对土壤的概化，通过假设土壤分层，便于 湿润锋的计算，原理类似于数学中分段积分的概念。分层的厚度影响计算结构的精度，其误 差也即为单层的厚度。一般分层越薄，计算精度越高，但计算量越大。因此，需要通过调试计 算，确定分层厚度。一般分层厚度0.01m基本可以满足计算需求。图1是土壤分层示意图， 图1中的数字表示分层的层数，最终分层的层数由湿润锋推移速度趋近于〇的位置决定。 水力传导计算分层内的湿润锋推移速度的步骤：用于通过公式计算第i层的湿润 锋推移速度：其中： i一一为分层的某个层数，i是自然数， K{ 〇x ,)一一为第i减1层的土壤含水量为故时的非饱和水力传导度， K(当前第1页1 2 本文档来自技高网...

【技术保护点】
一种降雨入渗补给地下水临界埋深计算方法，其特征在于，所述方法的步骤如下：饱和入渗深度计算的步骤：用于对降雨入渗阶段土壤中的水分饱和入渗深度进行计算，计算公式为：，其中：——为饱和入渗过程中湿润锋位置，单位：米，——为降雨入渗量，单位：毫米，——土壤饱和含水量，cm??/cm??或用百分比表示，θi——土壤初始含水量，cm??/cm??或用百分比表示；饱和入渗深度之下非饱和入渗深度计算的步骤：用于对土壤分层，根据湿润锋推进速度，对土壤水再分配阶段非饱和入渗深度进行计算；水力传导计算分层内的湿润锋推移速度的步骤：用于通过公式计算第i层的湿润锋推移速度：，其中：K(θi‑1)——为第i减1层的土壤含水量为时的非饱和水力传导度，K(θi)——为第i层的土壤含水量为时的非饱和水力传导度，hm,i——为第i层的基质势Ψm,i对应的当量水头,hm,i‑1——为第i减1层的基质势Ψm,i‑1对应的当量水头,Δz——为各层土壤的厚度；判断推移速度是否大于零的步骤：用于判断湿润锋推移速度是否大于零，如果大于零则进入“计算土壤含水量的步骤”，如果小于等于零则进入“计算非饱和入渗深度的步骤”；计算土壤含水量的步骤：用于根据公式计算当前时刻湿润锋上部土壤含水量：，其中：θt——t时刻湿润锋上部的土壤含水量，θt‑1——t‑1时刻湿润锋上部的土壤含水量，Ht‑1——为t‑1时刻湿润层厚度，计算结束后，将i加1，并回到“饱和入渗深度之下非饱和入渗深度计算的步骤”；计算非饱和入渗深度的步骤：用于将各分层土壤厚度相加ΣΔz，形成非饱和入渗深度；计算入渗深度的步骤：用于将饱和入渗深度HP与非饱和入渗深度ΣΔz相加，即为降雨入渗补给地下水临界埋深。...

【技术特征摘要】

【专利技术属性】
技术研发人员：陈敏建，张秋霞，马静，汪林，曾向辉，陈根发，王鹏，刘盈斐，闫龙，
申请(专利权)人：中国水利水电科学研究院，
类型：发明
国别省市：北京;11