一种基于公里网格的浅层地下水蓄变量计算方法技术

本发明专利技术公开了一种基于公里网格的浅层地下水蓄变量计算方法，首先建立面向地下水需求应用的标准化公里网格，再基于标准化公里网格构建网格数据集，最后根据标准化公里网格和网格数据集计算研究区域的浅层地下水蓄变量。本发明专利技术在对区域浅层地下水蓄变量进行计算时既采用了等值线法，又考虑了水位站的实际控制面积，具有较高的计算精度；本发明专利技术通过建立标准化公里网格计算研究区域的浅层地下水蓄变量，由于公里网格的面积确定为1，因此大大减少了计算量；本发明专利技术计算得到研究区域的浅层地下水蓄变量后，研究区域内任意范围的浅层地下水蓄变量可通过查询对应的网格插值简单计算得到，便于共享，具有计算一次后可多区域查询分析，重复使用的优点。重复使用的优点。重复使用的优点。

【技术实现步骤摘要】
一种基于公里网格的浅层地下水蓄变量计算方法


[0001]本专利技术属于地下水资源管理
，具体涉及一种基于公里网格的浅层地下水蓄变量计算方法的设计。

技术介绍

[0002]地下水蓄变量是指均衡区计算时段初地下水储存量与计算时段末地下水储存量的差值。均衡区可利用一年或数年中地下水位的变化计算该时段含水层中地下水的盈亏情况，当水位上升时，蓄变量增加，反之蓄变量减少。因为深层地下水蓄变量变化一般比较平缓，因此针对浅层地下水蓄变量的研究和分析更具实用性，现有常规的区域浅层地下水蓄变量的计算方法主要包括：等值线法、泰森多边形法和算术平均法。
[0003]等值线法首先在地图上绘出等值线，用求积仪量出各相邻蓄变量等值线间的面积，以A1,A2,
…
,A
n
表示相邻两等值线间的面积，以B1,B2,
…
,B
n
作为该面积上的平均蓄变量(即“给水度
×
变幅”值)，计算区域浅层地下水蓄变量的公式为：R＝A1×
B1+A2×
B2+
…
+A
n
×
B
n
。该方法精度较高，但是每次计算都需要重绘等值线图，工作量大且繁琐，不方便数据成果共享。
[0004]泰森多边形法又称为垂直平分法，首先从地图上绘制三角网格，绘出各站点泰森多边形所占的区域面积A1,A2,
…
,A
n
，量算出各个地下水站点所对应的面积权重(即“给水度
×
变幅”值)K1,K2,
…
,K
n
，以此来计算出泰森多边形法中各个站点所对应的区域浅层地下水蓄变量R＝A1×
K1+A2×
K2+
…
+A
n
×
K
n
。该方法虽然考虑了权重，但权重是根据水位站的控制几何面积确定的，与实际控制面积不符合，因此计算精度较低。
[0005]算术平均法是求取蓄变量数据的算术平均值其中S为区域总面积，K1,K2,
…
,K
n
为个站点地下水的“给水度
×
变幅”值。该方法计算简单，但没有考虑各水位站的权重或看作等权重，只能粗略估计流域特定时段内的平均变幅量，计算精度同样较低。

技术实现思路

[0006]本专利技术的目的是为了解决现有的区域浅层地下水蓄变量的计算方法无法兼顾高计算精度和低计算量的问题，提出了一种基于公里网格的浅层地下水蓄变量计算方法。
[0007]本专利技术的技术方案为：一种基于公里网格的浅层地下水蓄变量计算方法，包括以下步骤：
[0008]S1、建立面向地下水需求应用的标准化公里网格。
[0009]S2、基于标准化公里网格构建网格数据集。
[0010]S3、根据标准化公里网格和网格数据集计算研究区域的浅层地下水蓄变量。
[0011]进一步地，步骤S1包括以下分步骤：
[0012]S11、确定公里网格的范围和投影。
[0013]S12、对公里网格进行编码。
[0014]S13、在公里网格尺度下，提取网格中心点。
[0015]S14、对研究区域边界处的网格进行处理，得到面向地下水需求应用的标准化公里网格。
[0016]进一步地，步骤S12中对公里网格进行9位编码，并将编码分为A部分和B部分，其中A部分代表省份，由两位数字表示；B部分代表网格序号，由7位数字组成，网格序号按照网格空间位置从左到右，从上到下的顺序进行编码，没有达到7位的数用补0处理。
[0017]进一步地，步骤S14中对研究区域边界处的网格进行处理的具体方法为：当研究区域边界处的网格有大于或等于75％的面积在研究区域内时，将该网格纳入研究区域进行计算，否则该网格不纳入研究区域计算。
[0018]进一步地，步骤S2中构建的网格数据集包括水位/埋深网格数据集、水位/埋深变幅网格数据集和给水度网格数据集。
[0019]进一步地，水位/埋深网格数据集的构建方法为：采用算术平均法计算研究区域内每个测站水位/埋深的日均值和月均值，并采用克里金插值方法对计算得到的日均值和月均值在空间上离散，形成水位/埋深网格数据集。
[0020]进一步地，水位/埋深变幅网格数据集的构建方法为：采用网格数据空间计算的方式计算相邻两个水位/埋深网格数据集中的数据差值，得到水位/埋深变幅网格数据集。
[0021]进一步地，给水度网格数据集的构建方法为：基于国家地下水监测一期工程采集的给水度参数矢量数据，将给水度与公里网格中心点相连接，并对研究区域边界部分网格采用相邻网格给水度进行补充，形成给水度网格数据集。
[0022]进一步地，步骤S3包括以下分步骤：
[0023]S31、获取研究区域的浅层统测水位，并利用克里金插值绘制研究区域地下水分层流场。
[0024]S32、对地下水分层流场进行网格离散化处理，并对每个标准化公里网格赋予一个水位变幅值。
[0025]S33、根据每个标准化公里网格的水位变幅值计算得到研究区域的浅层地下水蓄变量。
[0026]进一步地，步骤S33中研究区域的浅层地下水蓄变量的计算公式为：
[0027][0028]其中Q
s
表示研究区域的浅层地下水蓄变量，F
i
表示第i个公里网格的面积，数值上等于1，μ
i
表示第i个公里网格的给水度，通过给水度网格数据集获取，Δh
i
表示第i个公里网格的水位变幅值，通过水位/埋深变幅网格数据集获取，n表示研究区域的公里网格数量。
[0029]本专利技术的有益效果是：
[0030](1)本专利技术在对区域浅层地下水蓄变量进行计算时既采用了等值线法，又考虑了水位站的实际控制面积，具有较高的计算精度。
[0031](2)本专利技术通过建立标准化公里网格计算研究区域的浅层地下水蓄变量，由于公里网格的面积确定为1，因此大大减少了计算量。
[0032](3)本专利技术计算得到研究区域的浅层地下水蓄变量后，研究区域内任意范围的浅
层地下水蓄变量可通过查询对应的网格插值简单计算得到，便于共享，具有计算一次后可多区域查询分析，重复使用的优点。
附图说明
[0033]图1所示为本专利技术实施例提供的一种基于公里网格的浅层地下水蓄变量计算方法流程图。
具体实施方式
[0034]现在将参考附图来详细描述本专利技术的示例性实施方式。应当理解，附图中示出和描述的实施方式仅仅是示例性的，意在阐释本专利技术的原理和精神，而并非限制本专利技术的范围。
[0035]本专利技术实施例提供了一种基于公里网格的浅层地下水蓄变量计算方法，如图1所示，包括以下步骤S1～S3：
[0036]S1、建立面向地下水需求应用的标准化公里网格。
[0037]步骤S1包括以下分步骤S11～S14：本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于公里网格的浅层地下水蓄变量计算方法，其特征在于，包括以下步骤：S1、建立面向地下水需求应用的标准化公里网格；S2、基于标准化公里网格构建网格数据集；S3、根据标准化公里网格和网格数据集计算研究区域的浅层地下水蓄变量。2.根据权利要求1所述的浅层地下水蓄变量计算方法，其特征在于，所述步骤S1包括以下分步骤：S11、确定公里网格的范围和投影；S12、对公里网格进行编码；S13、在公里网格尺度下，提取网格中心点；S14、对研究区域边界处的网格进行处理，得到面向地下水需求应用的标准化公里网格。3.根据权利要求2所述的浅层地下水蓄变量计算方法，其特征在于，所述步骤S12中对公里网格进行9位编码，并将编码分为A部分和B部分，其中A部分代表省份，由两位数字表示；B部分代表网格序号，由7位数字组成，网格序号按照网格空间位置从左到右，从上到下的顺序进行编码，没有达到7位的数用补0处理。4.根据权利要求2所述的浅层地下水蓄变量计算方法，其特征在于，所述步骤S14中对研究区域边界处的网格进行处理的具体方法为：当研究区域边界处的网格有大于或等于75％的面积在研究区域内时，将该网格纳入研究区域进行计算，否则该网格不纳入研究区域计算。5.根据权利要求1所述的浅层地下水蓄变量计算方法，其特征在于，所述步骤S2中构建的网格数据集包括水位/埋深网格数据集、水位/埋深变幅网格数据集和给水度网格数据集。6.根据权利要求5所述的浅层地下水蓄变量计算方法，其特征在于，所述水位/埋深网格数据集的构建方法为：采用算术平均法计算研究区域内每个测站水位/埋深的日均值和月均值，并采用克里金...

【专利技术属性】
技术研发人员：刘杰，王卓然，孙龙，卢洪健，周扬，章雨乾，
申请(专利权)人：水利部信息中心，
类型：发明
国别省市：