[發明專利]一種基于公里網格的淺層地下水蓄變量計算方法

在審

• 申請號: 202211446533.3
• 申請日: 2022-11-18
• 公開（公告）號: CN115687883A
• 公開（公告）日: 2023-02-03
• 發明（設計）人: 劉杰;王卓然;孫龍;盧洪健;周揚;章雨乾
• 申請（專利權）人: 水利部信息中心
• 主分類號: G06F17/18
• 分類號: G06F17/18;G06Q50/06
• 代理公司: 北京正華智誠專利代理事務所(普通合伙) 11870
• 代理人: 劉方正

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表示研究區域的公里網格數量。