[发明专利]一种基于公里网格的浅层地下水蓄变量计算方法

权利要求书

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所述的浅层地下水蓄变量计算方法，其特征在于，所述水位/埋深网格数据集的构建方法为：采用算术平均法计算研究区域内每个测站水位/埋深的日均值和月均值，并采用克里金插值方法对计算得到的日均值和月均值在空间上离散，形成水位/埋深网格数据集。

7.根据权利要求5所述的浅层地下水蓄变量计算方法，其特征在于，所述水位/埋深变幅网格数据集的构建方法为：采用网格数据空间计算的方式计算相邻两个水位/埋深网格数据集中的数据差值，得到水位/埋深变幅网格数据集。

8.根据权利要求5所述的浅层地下水蓄变量计算方法，其特征在于，所述给水度网格数据集的构建方法为：基于国家地下水监测一期工程采集的给水度参数矢量数据，将给水度与公里网格中心点相连接，并对研究区域边界部分网格采用相邻网格给水度进行补充，形成给水度网格数据集。

9.根据权利要求5所述的浅层地下水蓄变量计算方法，其特征在于，所述步骤S3包括以下分步骤：

S31、获取研究区域的浅层统测水位，并利用克里金插值绘制研究区域地下水分层流场；

S32、对地下水分层流场进行网格离散化处理，并对每个标准化公里网格赋予一个水位变幅值；

S33、根据每个标准化公里网格的水位变幅值计算得到研究区域的浅层地下水蓄变量。

10.根据权利要求9所述的浅层地下水蓄变量计算方法，其特征在于，所述步骤S33中研究区域的浅层地下水蓄变量的计算公式为：

其中Q_s表示研究区域的浅层地下水蓄变量，F_i表示第i个公里网格的面积，数值上等于1，μ_i表示第i个公里网格的给水度，通过给水度网格数据集获取，Δh_i表示第i个公里网格的水位变幅值，通过水位/埋深变幅网格数据集获取，n表示研究区域的公里网格数量。