CN115687883A - 一种基于公里网格的浅层地下水蓄变量计算方法 - Google Patents
一种基于公里网格的浅层地下水蓄变量计算方法 Download PDFInfo
- Publication number
- CN115687883A CN115687883A CN202211446533.3A CN202211446533A CN115687883A CN 115687883 A CN115687883 A CN 115687883A CN 202211446533 A CN202211446533 A CN 202211446533A CN 115687883 A CN115687883 A CN 115687883A
- Authority
- CN
- China
- Prior art keywords
- grid
- kilometer
- data set
- water level
- shallow groundwater
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 53
- 239000003673 groundwater Substances 0.000 title claims abstract description 43
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 113
- 238000000034 method Methods 0.000 claims abstract description 38
- 238000011160 research Methods 0.000 claims abstract description 34
- 238000009933 burial Methods 0.000 claims description 27
- 238000009825 accumulation Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 6
- 238000011835 investigation Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 230000001502 supplementing effect Effects 0.000 claims description 3
- 239000002689 soil Substances 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 229920006395 saturated elastomer Polymers 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 240000004282 Grewia occidentalis Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000036571 hydration Effects 0.000 description 1
- 238000006703 hydration reaction Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于公里网格的浅层地下水蓄变量计算方法,首先建立面向地下水需求应用的标准化公里网格,再基于标准化公里网格构建网格数据集,最后根据标准化公里网格和网格数据集计算研究区域的浅层地下水蓄变量。本发明在对区域浅层地下水蓄变量进行计算时既采用了等值线法,又考虑了水位站的实际控制面积,具有较高的计算精度;本发明通过建立标准化公里网格计算研究区域的浅层地下水蓄变量,由于公里网格的面积确定为1,因此大大减少了计算量;本发明计算得到研究区域的浅层地下水蓄变量后,研究区域内任意范围的浅层地下水蓄变量可通过查询对应的网格插值简单计算得到,便于共享,具有计算一次后可多区域查询分析,重复使用的优点。
Description
技术领域
本发明属于地下水资源管理技术领域,具体涉及一种基于公里网格的浅层地下水蓄变量计算方法的设计。
背景技术
地下水蓄变量是指均衡区计算时段初地下水储存量与计算时段末地下水储存量的差值。均衡区可利用一年或数年中地下水位的变化计算该时段含水层中地下水的盈亏情况,当水位上升时,蓄变量增加,反之蓄变量减少。因为深层地下水蓄变量变化一般比较平缓,因此针对浅层地下水蓄变量的研究和分析更具实用性,现有常规的区域浅层地下水蓄变量的计算方法主要包括:等值线法、泰森多边形法和算术平均法。
等值线法首先在地图上绘出等值线,用求积仪量出各相邻蓄变量等值线间的面积,以A1,A2,…,An表示相邻两等值线间的面积,以B1,B2,…,Bn作为该面积上的平均蓄变量(即“给水度×变幅”值),计算区域浅层地下水蓄变量的公式为:R=A1×B1+A2×B2+…+An×Bn。该方法精度较高,但是每次计算都需要重绘等值线图,工作量大且繁琐,不方便数据成果共享。
泰森多边形法又称为垂直平分法,首先从地图上绘制三角网格,绘出各站点泰森多边形所占的区域面积A1,A2,…,An,量算出各个地下水站点所对应的面积权重(即“给水度×变幅”值)K1,K2,…,Kn,以此来计算出泰森多边形法中各个站点所对应的区域浅层地下水蓄变量R=A1×K1+A2×K2+…+An×Kn。该方法虽然考虑了权重,但权重是根据水位站的控制几何面积确定的,与实际控制面积不符合,因此计算精度较低。
发明内容
本发明的目的是为了解决现有的区域浅层地下水蓄变量的计算方法无法兼顾高计算精度和低计算量的问题,提出了一种基于公里网格的浅层地下水蓄变量计算方法。
本发明的技术方案为:一种基于公里网格的浅层地下水蓄变量计算方法,包括以下步骤:
S1、建立面向地下水需求应用的标准化公里网格。
S2、基于标准化公里网格构建网格数据集。
S3、根据标准化公里网格和网格数据集计算研究区域的浅层地下水蓄变量。
进一步地,步骤S1包括以下分步骤:
S11、确定公里网格的范围和投影。
S12、对公里网格进行编码。
S13、在公里网格尺度下,提取网格中心点。
S14、对研究区域边界处的网格进行处理,得到面向地下水需求应用的标准化公里网格。
进一步地,步骤S12中对公里网格进行9位编码,并将编码分为A部分和B部分,其中A部分代表省份,由两位数字表示;B部分代表网格序号,由7位数字组成,网格序号按照网格空间位置从左到右,从上到下的顺序进行编码,没有达到7位的数用补0处理。
进一步地,步骤S14中对研究区域边界处的网格进行处理的具体方法为:当研究区域边界处的网格有大于或等于75%的面积在研究区域内时,将该网格纳入研究区域进行计算,否则该网格不纳入研究区域计算。
进一步地,步骤S2中构建的网格数据集包括水位/埋深网格数据集、水位/埋深变幅网格数据集和给水度网格数据集。
进一步地,水位/埋深网格数据集的构建方法为:采用算术平均法计算研究区域内每个测站水位/埋深的日均值和月均值,并采用克里金插值方法对计算得到的日均值和月均值在空间上离散,形成水位/埋深网格数据集。
进一步地,水位/埋深变幅网格数据集的构建方法为:采用网格数据空间计算的方式计算相邻两个水位/埋深网格数据集中的数据差值,得到水位/埋深变幅网格数据集。
进一步地,给水度网格数据集的构建方法为:基于国家地下水监测一期工程采集的给水度参数矢量数据,将给水度与公里网格中心点相连接,并对研究区域边界部分网格采用相邻网格给水度进行补充,形成给水度网格数据集。
进一步地,步骤S3包括以下分步骤:
S31、获取研究区域的浅层统测水位,并利用克里金插值绘制研究区域地下水分层流场。
S32、对地下水分层流场进行网格离散化处理,并对每个标准化公里网格赋予一个水位变幅值。
S33、根据每个标准化公里网格的水位变幅值计算得到研究区域的浅层地下水蓄变量。
进一步地,步骤S33中研究区域的浅层地下水蓄变量的计算公式为:
其中Qs表示研究区域的浅层地下水蓄变量,Fi表示第i个公里网格的面积,数值上等于1,μi表示第i个公里网格的给水度,通过给水度网格数据集获取,Δhi表示第i个公里网格的水位变幅值,通过水位/埋深变幅网格数据集获取,n表示研究区域的公里网格数量。
本发明的有益效果是:
(1)本发明在对区域浅层地下水蓄变量进行计算时既采用了等值线法,又考虑了水位站的实际控制面积,具有较高的计算精度。
(2)本发明通过建立标准化公里网格计算研究区域的浅层地下水蓄变量,由于公里网格的面积确定为1,因此大大减少了计算量。
(3)本发明计算得到研究区域的浅层地下水蓄变量后,研究区域内任意范围的浅层地下水蓄变量可通过查询对应的网格插值简单计算得到,便于共享,具有计算一次后可多区域查询分析,重复使用的优点。
附图说明
图1所示为本发明实施例提供的一种基于公里网格的浅层地下水蓄变量计算方法流程图。
具体实施方式
现在将参考附图来详细描述本发明的示例性实施方式。应当理解,附图中示出和描述的实施方式仅仅是示例性的,意在阐释本发明的原理和精神,而并非限制本发明的范围。
本发明实施例提供了一种基于公里网格的浅层地下水蓄变量计算方法,如图1所示,包括以下步骤S1~S3:
S1、建立面向地下水需求应用的标准化公里网格。
步骤S1包括以下分步骤S11~S14:
S11、确定公里网格的范围和投影。
本发明实施例中,参考国家测绘地理信息局提供1:100万的标准地图,其四角坐标分别为:73.25°E,53.67°N、73.25°E,18°N、135°E,18°N、135°E,53.67°N。网格范围与其一致,同时参考《地下水动态月报》的投影参数,具体如下:
地理坐标系:China Geodetic Coordinate System 2000。
投影坐标系:阿伯斯投影(双标准纬线等积圆锥投影),其中中央经线105°,双纬线分别为25°N、47°N。
在此投影坐标系下,其标准地形图的四角坐标为分别为:-2664038.237862,5904091.998909、-2664038.237862,1850574.589658、2240815.194745,1850574.589658、2240815.194745,5904091.998909。
S12、对公里网格进行编码。
本发明实施例中,对公里网格进行9位编码,并将编码分为A部分和B部分,其中A部分代表省份,由两位数字表示;B部分代表网格序号,由7位数字组成,网格序号按照网格空间位置从左到右,从上到下的顺序进行编码,没有达到7位的数用补0处理,其取值范围为0000001-9999999。具体网格编码如表1和表2所示。
表1网格编码为规定
表2省级行政区划代码
代码 | 名称 | 代码 | 名称 |
11 | 北京市 | 42 | 湖北省 |
12 | 天津市 | 43 | 湖南省 |
13 | 河北省 | 44 | 广东省 |
14 | 山西省 | 45 | 广西壮族自治区 |
15 | 内蒙古自治区 | 46 | 海南省 |
21 | 辽宁省 | 50 | 重庆市 |
22 | 吉林省 | 51 | 四川省 |
23 | 黑龙江省 | 52 | 贵州省 |
31 | 上海市 | 53 | 云南省 |
32 | 江苏省 | 54 | 西藏自治区 |
33 | 浙江省 | 61 | 陕西省 |
34 | 安徽省 | 62 | 甘肃省 |
35 | 福建省 | 63 | 青海省 |
36 | 江西省 | 64 | 宁夏回族自治区 |
37 | 山东省 | 65 | 新疆维吾尔自治区 |
41 | 河南省 |
S13、在公里网格尺度下,提取网格中心点,作为后续网格数据集的赋值依据。
S14、对研究区域边界处的网格进行处理,得到面向地下水需求应用的标准化公里网格。
计算各平原区的浅层地下水蓄变量时,边界处容易出现网格边界与实际边界不贴合的情况,对蓄变量计算结果产生的影响不可忽略,因此需要合理的处理方案。本发明实施例中,以网格面积占比75%为界限,对研究区域边界处的网格进行处理。当研究区域边界处的网格有大于或等于75%的面积在研究区域内时,将该网格纳入研究区域进行计算,否则该网格不纳入研究区域计算。
S2、基于标准化公里网格构建网格数据集。
本发明实施例中,构建的网格数据集包括水位/埋深网格数据集、水位/埋深变幅网格数据集和给水度网格数据集。
其中,水位/埋深网格数据集的构建方法为:采用算术平均法计算研究区域内每个测站水位/埋深的日均值和月均值,并采用克里金插值方法对计算得到的日均值和月均值在空间上离散,形成水位/埋深网格数据集。本发明实施例中选取了京津冀地区2018年1月至2020年3月的水位/埋深月均值数据进行计算,形成了27个月的水位/埋深网格数据集。
水位/埋深变幅网格数据集的构建方法为:采用网格数据空间计算的方式计算相邻两个水位/埋深网格数据集中的数据差值,得到水位/埋深变幅网格数据集。本发明实施例中选取了2018年1月至2020年3月逐月与上月相比的水位/埋深变幅进行网格数据空间计算,形成26个月的水位/埋深变幅网格数据集。
给水度网格数据集的构建方法为:基于国家地下水监测一期工程采集的给水度参数矢量数据,将给水度与公里网格中心点相连接,并对研究区域边界部分网格采用相邻网格给水度进行补充,形成给水度网格数据集。
给水度是含水层的释水能力,其定义为饱和的土壤或岩层在重力作用下排出的水量与土壤或岩层体积的比值。它表示单位面积的含水层,当潜水面下降一个单位长度时在重力作用下所能释放出的水量。数值上给水度等于释出的水的体积与释水的饱和岩土总体积之比,等于孔隙度减去持水度。
给水度是地下水资源评价中最主要的参数,其计算公式为:
其中μ表示给水度,VW表示所研究土体经饱和后能由重力排出的水量体积,V表示所研究土体的总体积。
S3、根据标准化公里网格和网格数据集计算研究区域的浅层地下水蓄变量。
步骤S3包括以下分步骤S31~S33:
S31、获取研究区域的浅层统测水位,并利用克里金插值绘制研究区域地下水分层流场。
S32、对地下水分层流场进行网格离散化处理,并对每个标准化公里网格赋予一个水位变幅值。
S33、根据每个标准化公里网格的水位变幅值计算得到研究区域的浅层地下水蓄变量:
其中Qs表示研究区域的浅层地下水蓄变量,Fi表示第i个公里网格的面积,数值上等于1,μi表示第i个公里网格的给水度,通过给水度网格数据集获取,Δhi表示第i个公里网格的水位变幅值,通过水位/埋深变幅网格数据集获取,n表示研究区域的公里网格数量。
本发明实施例中,计算得到研究区域的浅层地下水蓄变量Qs后,研究区域内任意范围的浅层地下水蓄变量可通过查询对应的网格插值简单计算得到,便于共享,具有计算一次后可多区域查询分析,重复使用的优点。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (10)
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、根据每个标准化公里网格的水位变幅值计算得到研究区域的浅层地下水蓄变量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211446533.3A CN115687883B (zh) | 2022-11-18 | 2022-11-18 | 一种基于公里网格的浅层地下水蓄变量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211446533.3A CN115687883B (zh) | 2022-11-18 | 2022-11-18 | 一种基于公里网格的浅层地下水蓄变量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115687883A true CN115687883A (zh) | 2023-02-03 |
CN115687883B CN115687883B (zh) | 2024-04-30 |
Family
ID=85053029
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211446533.3A Active CN115687883B (zh) | 2022-11-18 | 2022-11-18 | 一种基于公里网格的浅层地下水蓄变量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115687883B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008084243A (ja) * | 2006-09-29 | 2008-04-10 | Hitachi Engineering & Services Co Ltd | 氾濫シミュレーション装置およびプログラム |
CN102156779A (zh) * | 2011-04-13 | 2011-08-17 | 北京石油化工学院 | 地下水流仿真与预测分析方法 |
CN109492299A (zh) * | 2018-11-07 | 2019-03-19 | 南开大学 | 基于swmm与modflow耦合的水资源模拟方法 |
CN112268595A (zh) * | 2020-10-29 | 2021-01-26 | 水利部水利水电规划设计总院 | 一种区域平均地下水水位监测方法 |
CN113139308A (zh) * | 2020-12-19 | 2021-07-20 | 武汉大学 | 空间地理数据与地下水数值模型的嵌套交互方法及装置 |
CN114357678A (zh) * | 2022-01-06 | 2022-04-15 | 余楚 | 一种新型区域地下水位监测网优化设计方法 |
CN114861502A (zh) * | 2022-05-27 | 2022-08-05 | 北京林业大学 | 基于Modflow模型的三维动态地下水污染模拟的安全饮用水区域确定方法 |
CN114970302A (zh) * | 2022-05-10 | 2022-08-30 | 沈阳工业大学 | 一种基于地下水监测系统的区域地下水情预测方法 |
-
2022
- 2022-11-18 CN CN202211446533.3A patent/CN115687883B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008084243A (ja) * | 2006-09-29 | 2008-04-10 | Hitachi Engineering & Services Co Ltd | 氾濫シミュレーション装置およびプログラム |
CN102156779A (zh) * | 2011-04-13 | 2011-08-17 | 北京石油化工学院 | 地下水流仿真与预测分析方法 |
CN109492299A (zh) * | 2018-11-07 | 2019-03-19 | 南开大学 | 基于swmm与modflow耦合的水资源模拟方法 |
CN112268595A (zh) * | 2020-10-29 | 2021-01-26 | 水利部水利水电规划设计总院 | 一种区域平均地下水水位监测方法 |
CN113139308A (zh) * | 2020-12-19 | 2021-07-20 | 武汉大学 | 空间地理数据与地下水数值模型的嵌套交互方法及装置 |
CN114357678A (zh) * | 2022-01-06 | 2022-04-15 | 余楚 | 一种新型区域地下水位监测网优化设计方法 |
CN114970302A (zh) * | 2022-05-10 | 2022-08-30 | 沈阳工业大学 | 一种基于地下水监测系统的区域地下水情预测方法 |
CN114861502A (zh) * | 2022-05-27 | 2022-08-05 | 北京林业大学 | 基于Modflow模型的三维动态地下水污染模拟的安全饮用水区域确定方法 |
Non-Patent Citations (2)
Title |
---|
BANGLE HE ET.AL: "Fault Identification of High-voltage Cable Sheath Grounding System Based on Ground Current Analysis", 《PREPRINTS OF THE 4TH IEEE CONFERENCE ON ENERGY INTERNET AND ENERGY SYSTEM INTEGRATION》, pages 3248 - 3251 * |
李冬田等: "安集海地区地下水空间信息系统", 《河海大学学报》, pages 8 - 13 * |
Also Published As
Publication number | Publication date |
---|---|
CN115687883B (zh) | 2024-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107180450B (zh) | 一种基于dem的河谷横断面形态的算法 | |
CN109815305B (zh) | 一种无资料地区场次洪水径流过程反演的方法 | |
CN103093114B (zh) | 一种基于地形和土壤特性的分布式流域缺水量测算方法 | |
CN104102845A (zh) | 尺度自适应的插值方法及尺度自适应的插值系统 | |
CN110852526B (zh) | 一种基于雨洪过程相似性判别的实时洪水预报方法 | |
CN104298841A (zh) | 一种基于历史数据的洪水预报方法和系统 | |
CN110501733B (zh) | 一种自适应格网vrs生成与服务方法 | |
CN115048354A (zh) | 一种水文模型的创建及径流预测方法、装置及计算机设备 | |
CN104392113B (zh) | 一种近海海面冷空气大风风速的估算方法 | |
Li et al. | A spatially explicit reconstruction of cropland based on expansion of polders in the Dongting Plain in China during 1750–1985 | |
CN110765677B (zh) | 一种高精度、快速三维地质模型有限元模型的建模方法 | |
Qie et al. | Method and its application of natural quality evaluation of arable land based on self-organizing feature map neural network | |
CN103422094A (zh) | 一种阴极保护浅埋式阳极地床接地电阻数值模拟方法 | |
CN115687883A (zh) | 一种基于公里网格的浅层地下水蓄变量计算方法 | |
Hui et al. | The Fractal Characteristics of Drainage Networks and Erosion Evolution Stages of Ten Kongduis in the Upper Reaches of the Yellow River, China | |
CN101957193A (zh) | 一种海岛礁高程传递的优化方法 | |
CN105588569A (zh) | 一种利用影子定位的方法 | |
CN110765298B (zh) | 矢量数据几何属性解耦的瓦片编码方法 | |
CN115858955B (zh) | 水电站入库流量及发电量时序分析方法、装置及存储介质 | |
CN111523077A (zh) | 一种海岸带陆地重力测量中区—远一区地形改正计算方法 | |
CN103903289A (zh) | 一种基于分形发生器的随机型沼泽地花纹图案生成方法 | |
CN113610436A (zh) | 一种承灾体动态脆弱性评估方法及系统 | |
CN101546262B (zh) | 测算地貌现象信息盒维数的方法 | |
CN113986906A (zh) | 一种基于态势目标的轨迹重构方法 | |
CN106709992A (zh) | 基于信息熵理论的海底地貌复杂度表示方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |