CN112529723A - 一种基于像元尺度的地下水补给量估算方法及系统 - Google Patents
一种基于像元尺度的地下水补给量估算方法及系统 Download PDFInfo
- Publication number
- CN112529723A CN112529723A CN202010359817.3A CN202010359817A CN112529723A CN 112529723 A CN112529723 A CN 112529723A CN 202010359817 A CN202010359817 A CN 202010359817A CN 112529723 A CN112529723 A CN 112529723A
- Authority
- CN
- China
- Prior art keywords
- unit
- monthly
- area
- data set
- pixel area
- 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.)
- Pending
Links
- 239000003673 groundwater Substances 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000011160 research Methods 0.000 claims abstract description 98
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 85
- 239000002689 soil Substances 0.000 claims abstract description 31
- 230000008018 melting Effects 0.000 claims abstract description 24
- 238000002844 melting Methods 0.000 claims abstract description 24
- 238000004519 manufacturing process Methods 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims description 64
- 238000006243 chemical reaction Methods 0.000 claims description 14
- 238000007781 pre-processing Methods 0.000 claims description 14
- 238000005520 cutting process Methods 0.000 claims description 7
- 238000003702 image correction Methods 0.000 claims description 7
- 238000007689 inspection Methods 0.000 claims description 7
- 230000007547 defect Effects 0.000 abstract description 3
- 238000011156 evaluation Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 9
- 238000012544 monitoring process Methods 0.000 description 8
- 239000005413 snowmelt Substances 0.000 description 4
- 238000001556 precipitation Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
- Y02A40/22—Improving land use; Improving water use or availability; Controlling erosion
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于像元尺度的地下水补给量估算方法及系统,包括获取研究区域的气候要素原始遥感数据集以及研究区域像元面积;所述气候要素原始遥感数据集包括单位像元面积月降水量数据、单位像元面积月实际蒸散发量数据、单位像元面积月融雪量数据、单位像元面积月土壤水分变化值数据以及月径流量数据;根据所述气候要素原始遥感数据集和所述研究区域像元面积,利用水量平衡方程,计算研究区域水资源总产量;根据所述研究区域水资源总产量和所述月径流量数据,估算研究区域的地下水补给量。采用本发明提供的方法或系统,克服了现有技术和现有模型难以实现基于空间像元尺度的地下水赋存量评估技术缺陷,弥补了现有模型的空白。
Description
技术领域
本发明涉及地下水补给量估算领域,特别是涉及一种基于像元尺度的地下水补给量估算方法及系统。
背景技术
地下水是与地表水同等重要的水资源,赋存量大,具有极高的开采价值。然而,由于气候变化和人类活动的原因,地下水补给量变化动态较大,且有逐年衰竭的趋势。当前地下水赋存量的估算方法主要是基于地下水出口或泉点进行实时监测数据进行估算,在时空尺度上具有很大的不确定性,这主要是由于地下水埋藏较深,监测难度较大,时间长,监测点位上的数据不充分,难以表征真实的地下水变化特征。由于监测数据的缺乏,监测区域小,而地下水资源在空间大尺度上具有流动性和补给性能,因此,地下水资源评估往往基于跨省跨国尺度。此外,地下水补给量赋存具有极其显著的时空异质性,亟需快速、高效、精准地实现监测,并获取数据进行分析,而现阶段技术和方法难以实现空间上的监测任务,给地下水资源的实时评估带来了很大的困难。因此,亟需一套能够快速、高效、精准估算地下水补给量的新方法,以实现地下水补给量评估研究,为社会经济发展和生态恢复建设提供新的技术支撑。
发明内容
本发明的目的是提供一种基于像元尺度的地下水补给量估算方法及系统,克服现有技术和现有模型难以实现基于空间像元尺度的地下水赋存量评估技术缺陷,弥补现有模型的空白。
为实现上述目的,本发明提供了如下方案:
一种基于像元尺度的地下水补给量估算方法,包括:
获取研究区域的气候要素原始遥感数据集以及研究区域像元面积;所述气候要素原始遥感数据集包括单位像元面积月降水量数据、单位像元面积月实际蒸散发量数据、单位像元面积月融雪量数据、单位像元面积月土壤水分变化值数据以及月径流量数据;
根据所述气候要素原始遥感数据集和所述研究区域像元面积,利用水量平衡方程,计算研究区域水资源总产量;
根据所述研究区域水资源总产量和所述月径流量数据,估算研究区域的地下水补给量。
可选的,在执行根据所述气候要素原始遥感数据集和所述研究区域像元面积,利用水量平衡方程,计算研究区域水资源总产量之前,还包括:
对所述气候要素原始遥感数据集中的数据进行预处理,得到研究区域的气候要素处理数据集;预处理的操作步骤依次包括格式转换、影像校正、裁剪、配准、质量检查和投影转换。
可选的,所述水量平衡方程为S(QSN+P)=S(ET+ΔS)+R+G;
其中,S为像元面积,单位为m2;QSN为单位像元面积月融雪量,单位为mm;P为单位像元面积月降水量,单位为mm;ET为单位像元面积月实际蒸散发量,单位为mm;ΔS为单位像元面积月土壤水分变化值,单位为mm;R为月径流量,单位为m3;G为地下水补给量,单位为m3。
可选的,所述根据所述气候要素原始遥感数据集和所述研究区域像元面积,利用水量平衡方程,计算研究区域水资源总产量,具体包括:
根据公式W=R+G=S(QSN+P-ET-ΔS)计算研究区域水资源总产量;
其中,W为研究区域水资源总产量,单位为m3;R为气候要素处理数据集中的月径流量,单位为m3;G为地下水补给量,单位为m3;S为像元面积,单位为m2;QSN为气候要素处理数据集中的单位像元面积月融雪量,单位为mm;P为气候要素处理数据集中的单位像元面积月降水量,单位为mm;ET为气候要素处理数据集中的单位像元面积月实际蒸散发量,单位为mm;ΔS为气候要素处理数据集中的单位像元面积月土壤水分变化值,单位为mm。
可选的,所述根据所述研究区域水资源总产量和所述月径流量数据,估算研究区域的地下水补给量,具体包括:
根据公式G=W-R=S(QS+P-ET-ΔS)-R估算研究区域的地下水补给量;
其中,G为研究区域的地下水补给量,单位为m3;W为研究区域水资源总产量,单位为m3;R为气候要素处理数据集中的月径流量,单位为m3;S为像元面积,单位为m2;QSN为气候要素处理数据集中的单位像元面积月融雪量,单位为mm;P为气候要素处理数据集中的单位像元面积月降水量,单位为mm;ET为气候要素处理数据集中的单位像元面积月实际蒸散发量,单位为mm;ΔS为气候要素处理数据集中的单位像元面积月土壤水分变化值,单位为mm。
一种基于像元尺度的地下水补给量估算系统,包括:
信息获取模块,用于获取研究区域的气候要素原始遥感数据集以及研究区域像元面积;所述气候要素原始遥感数据集包括单位像元面积月降水量数据、单位像元面积月实际蒸散发量数据、单位像元面积月融雪量数据、单位像元面积月土壤水分变化值数据以及月径流量数据;
研究区域水资源总产量计算模块,用于根据所述气候要素原始遥感数据集和所述研究区域像元面积,利用水量平衡方程,计算研究区域水资源总产量;
地下水补给量估算模块,用于根据所述研究区域水资源总产量和所述月径流量数据,估算研究区域的地下水补给量。
可选的,还包括:
预处理模块,用于对所述气候要素原始遥感数据集中的数据进行预处理,得到研究区域的气候要素处理数据集;预处理的操作步骤依次包括格式转换、影像校正、裁剪、配准、质量检查和投影转换。
可选的,所述研究区域水资源总产量计算模块,具体包括:
研究区域水资源总产量计算单元,用于根据公式W=R+G=S(QSN+P-ET-ΔS)计算研究区域水资源总产量;
其中,W为研究区域水资源总产量,单位为m3;R为气候要素处理数据集中的月径流量,单位为m3;G为地下水补给量,单位为m3;S为像元面积,单位为m2;QSN为气候要素处理数据集中的单位像元面积月融雪量,单位为mm;P为气候要素处理数据集中的单位像元面积月降水量,单位为mm;ET为气候要素处理数据集中的单位像元面积月实际蒸散发量,单位为mm;ΔS为气候要素处理数据集中的单位像元面积月土壤水分变化值,单位为mm。
可选的,所述地下水补给量估算模块,具体包括:
地下水补给量估算单元,用于根据公式G=W-R=S(QS+P-ET-ΔS)-R估算研究区域的地下水补给量;
其中,G为研究区域的地下水补给量,单位为m3;W为研究区域水资源总产量,单位为m3;R为气候要素处理数据集中的月径流量,单位为m3;S为像元面积,单位为m2;QSN为气候要素处理数据集中的单位像元面积月融雪量,单位为mm;P为气候要素处理数据集中的单位像元面积月降水量,单位为mm;ET为气候要素处理数据集中的单位像元面积月实际蒸散发量,单位为mm;ΔS为气候要素处理数据集中的单位像元面积月土壤水分变化值,单位为mm。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种基于像元尺度的地下水补给量估算方法及系统,利用月径流量数据、月降水量数据、月实际蒸散发量数据以及月融雪量数据,根据水量平衡方程推导出地下水补给量估算模型,为地下水储量监测评估提供了一种新方法。该地下水补给量估算模型具有快速、高效、适用于全球大尺度的优势,解决了地下水补给量监测难度大、时间久、精度不高的问题,为生态恢复和社会经济发展相关研究提供了新的技术支撑和理论基础。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明基于像元尺度的地下水补给量估算方法的流程图;
图2为本发明基于像元尺度的地下水补给量估算系统的结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种基于像元尺度的地下水补给量估算方法及系统,克服现有技术和现有模型难以实现基于空间像元尺度的地下水赋存量评估技术缺陷,弥补现有模型的空白。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示,本发明提供了一种基于像元尺度的地下水补给量估算方法,具体包括以下步骤。
步骤101:获取研究区域的气候要素原始遥感数据集以及研究区域像元面积;所述气候要素原始遥感数据集包括全球陆地总水储量变化数据(GRACE)、单位像元面积月降水量数据、单位像元面积月实际蒸散发量数据、单位像元面积月融雪量数据、单位像元面积月土壤水分变化值数据以及月径流量数据等。
气候要素原始遥感数据集包括各月尺度的单位像元面积的降水量、实际蒸散发量、厚度为0-10cm、10-40cm、40-100cm和100-200cm的土壤水含量数据、融雪量数据,并合并成年值数据。这些数据均来源于美国国家航空航天局(NASA)(网址:https://www.nasa.gov/)的FLDAS数据集(FLDAS Noah Land Surface Model L4 Global Monthly0.1x 0.1degree(MERRA-2and CHIRPS)V001(FLDAS_NOAH01_C_GL_M)at GES DISC(网址:https://ldas.gsfc.nasa.gov/FLDAS/)),FLDAS数据集的空间分辨率为0.1°×0.1°,或GLDAS数据集。气候要素原始遥感数据集的时间分辨率为每月,空间范围覆盖全球(60S,180W,90N,180E)。
此外,用全球土壤深度数据来计算月尺度的土层厚度数据,全球土壤深度数据来源于网址:https://daac.ornl.gov/,空间分辨率为0.1°×0.1°,和网址:https://www.isric.org/explore/soilgrids,也可以根据研究尺度选择不同空间分辨率的的土壤深度数据,例如250m×250m、1km×1km、5km×5km和10km×10km。2015年最新行政区划矢量数据均来源于中国科学院资源环境科学数据网(网址:http://www.resdc.cn/)和国家测绘地理信息局(网址:http://www.sbsm.gov.cn/article/zxbs/dtfw/)。
全球陆地融雪量及地表径流量数据(GLDAS)来源于全球陆地数据同化系统(Goddard Earth Sciences Data and Information Services Center,简称GES DISC,(GLDAS Noah Land Surface Model L4 Monthly 0.25x 0.25degree),网址:https://mirador.gsfc.nasa.gov/)。
步骤102:对所述气候要素原始遥感数据集中的数据进行预处理,得到研究区域的气候要素处理数据集;预处理的操作步骤依次包括格式转换、影像校正、裁剪、配准、质量检查和投影转换等。
本发明利用数据同化方法,将气候要素原始遥感数据集中所有栅格数据的栅格单元大小均转为相同的空间分辨率,投影方式全部采用ALBERS等积圆锥投影(Krasovsky-1940-Albers,为一种投影坐标系)方式。
上述全球尺度的栅格数据通过格式转换、影像校正、裁剪、配准、质量检查、投影转换后,最终得到研究区域的气候要素处理数据集。
步骤103:根据所述气候要素原始遥感数据集和所述研究区域像元面积,利用水量平衡方程,计算研究区域水资源总产量。具体包括:
首先,基于水量平衡方程推导出研究区域水资源总产量的计算公式。
其次,将研究区域像元面积和气候要素处理数据集中的单位像元面积月融雪量、单位像元面积月降水量、单位像元面积月实际蒸散发量、单位像元面积月土壤水分变化值,输入到研究区域水资源总产量的计算公式中,确定研究区域水资源总产量。
其中,水量平衡方程为
S(QSN+P)=S(ET+ΔS)+R+G (1);
公式(1)中,S为像元面积,单位为m2;QSN为单位像元面积月融雪量,单位为mm;P为单位像元面积月降水量,单位为mm;ET为单位像元面积月实际蒸散发量,单位为mm;ΔS为单位像元面积月土壤水分变化值,单位为mm;R为月径流量,单位为m3;G为地下水补给量,单位为m3。
研究区域水资源总产量的计算公式为
W=R+G=S(QSN+P-ET-ΔS) (2);
公式(2)中,W为研究区域水资源总产量,单位为m3;R为气候要素处理数据集中的月径流量,单位为m3;G为地下水补给量,单位为m3;S为像元面积,单位为m2;QSN为气候要素处理数据集中的单位像元面积月融雪量,单位为mm;P为气候要素处理数据集中的单位像元面积月降水量,单位为mm;ET为气候要素处理数据集中的单位像元面积月实际蒸散发量,单位为mm;ΔS为气候要素处理数据集中的单位像元面积月土壤水分变化值,单位为mm。
融雪量数据来源于空间分辨率为0.125°×0.125°,覆盖于全球的(60S,180W,90N,180E)单位像元面积月尺度融雪量数据的重采样结果。上述融雪量数据来源于NLDAS NoahLand Surface Model L4 Monthly 0.125°×0.125°数据集(网址:https://mirador.gsfc.nasa.gov/)。
土壤水分变化值来源于利用美国国家航空航天局(NASA)(网址:https://www.nasa.gov/)的FLDAS数据集(FLDAS Noah Land Surface Model L4 Global Monthly0.1x 0.1degree(MERRA-2and CHIRPS)V001(FLDAS_NOAH01_C_GL_M)at GES DISC(网址:https://ldas.gsfc.nasa.gov/FLDAS/)),空间分辨率为0.1°×0.1°的土壤含水量栅格数据经重采样后作为土壤水分变化值,该土壤含水量栅格数据的土壤深度为2m。
步骤104:根据所述研究区域水资源总产量和所述月径流量数据,估算研究区域的地下水补给量。具体为:
首先,基于研究区域水资源总产量的计算公式推导出研究区域的地下水补给量的计算公式。
其次,将研究区域水资源总产量和气候要素处理数据集中的月径流量,输入到研究区域的地下水补给量的计算公式中,确定研究区域的地下水补给量。
其中,研究区域的地下水补给量的计算公式为
G=W-R=S(QS+P-ET-ΔS)-R (3);
公式(3)中,G为研究区域的地下水补给量,单位为m3;W为研究区域水资源总产量,单位为m3;R为气候要素处理数据集中的月径流量,单位为m3。
为实现上述目的,本发明还提供了一种基于像元尺度的地下水补给量估算系统,如图2所示,包括:
信息获取模块201,用于获取研究区域的气候要素原始遥感数据集以及研究区域像元面积;所述气候要素原始遥感数据集包括单位像元面积月降水量数据、单位像元面积月实际蒸散发量数据、单位像元面积月融雪量数据、单位像元面积月土壤水分变化值数据以及月径流量数据。
预处理模块202,用于对所述气候要素原始遥感数据集中的数据进行预处理,得到研究区域的气候要素处理数据集;预处理的操作步骤依次包括格式转换、影像校正、裁剪、配准、质量检查和投影转换。
研究区域水资源总产量计算模块203,用于根据所述气候要素原始遥感数据集和所述研究区域像元面积,利用水量平衡方程,计算研究区域水资源总产量。
地下水补给量估算模块204,用于根据所述研究区域水资源总产量和所述月径流量数据,估算研究区域的地下水补给量。
所述研究区域水资源总产量计算模块203,具体包括:
研究区域水资源总产量计算单元,用于根据公式W=R+G=S(QSN+P-ET-ΔS)计算研究区域水资源总产量。
其中,W为研究区域水资源总产量,单位为m3;R为气候要素处理数据集中的月径流量,单位为m3;G为地下水补给量,单位为m3;S为像元面积,单位为m2;QSN为气候要素处理数据集中的单位像元面积月融雪量,单位为mm;P为气候要素处理数据集中的单位像元面积月降水量,单位为mm;ET为气候要素处理数据集中的单位像元面积月实际蒸散发量,单位为mm;ΔS为气候要素处理数据集中的单位像元面积月土壤水分变化值,单位为mm。
所述地下水补给量估算模块204,具体包括:
地下水补给量估算单元,用于根据公式G=W-R=S(QS+P-ET-ΔS)-R估算研究区域的地下水补给量。
其中,G为研究区域的地下水补给量,单位为m3;W为研究区域水资源总产量,单位为m3;R为气候要素处理数据集中的月径流量,单位为m3;S为像元面积,单位为m2;QSN为气候要素处理数据集中的单位像元面积月融雪量,单位为mm;P为气候要素处理数据集中的单位像元面积月降水量,单位为mm;ET为气候要素处理数据集中的单位像元面积月实际蒸散发量,单位为mm;ΔS为气候要素处理数据集中的单位像元面积月土壤水分变化值,单位为mm。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (9)
1.一种基于像元尺度的地下水补给量估算方法,其特征在于,包括:
获取研究区域的气候要素原始遥感数据集以及研究区域像元面积;所述气候要素原始遥感数据集包括单位像元面积月降水量数据、单位像元面积月实际蒸散发量数据、单位像元面积月融雪量数据、单位像元面积月土壤水分变化值数据以及月径流量数据;
根据所述气候要素原始遥感数据集和所述研究区域像元面积,利用水量平衡方程,计算研究区域水资源总产量;
根据所述研究区域水资源总产量和所述月径流量数据,估算研究区域的地下水补给量。
2.根据权利要求1所述的一种基于像元尺度的地下水补给量估算方法,其特征在于,在执行根据所述气候要素原始遥感数据集和所述研究区域像元面积,利用水量平衡方程,计算研究区域水资源总产量之前,还包括:
对所述气候要素原始遥感数据集中的数据进行预处理,得到研究区域的气候要素处理数据集;预处理的操作步骤依次包括格式转换、影像校正、裁剪、配准、质量检查和投影转换。
3.根据权利要求1所述的一种基于像元尺度的地下水补给量估算方法,其特征在于,所述水量平衡方程为S(QSN+P)=S(ET+ΔS)+R+G;
其中,S为像元面积,单位为m2;QSN为单位像元面积月融雪量,单位为mm;P为单位像元面积月降水量,单位为mm;ET为单位像元面积月实际蒸散发量,单位为mm;ΔS为单位像元面积月土壤水分变化值,单位为mm;R为月径流量,单位为m3;G为地下水补给量,单位为m3。
4.根据权利要求2所述的一种基于像元尺度的地下水补给量估算方法,其特征在于,所述根据所述气候要素原始遥感数据集和所述研究区域像元面积,利用水量平衡方程,计算研究区域水资源总产量,具体包括:
根据公式W=R+G=S(QSN+P-ET-ΔS)计算研究区域水资源总产量;
其中,W为研究区域水资源总产量,单位为m3;R为气候要素处理数据集中的月径流量,单位为m3;G为地下水补给量,单位为m3;S为像元面积,单位为m2;QSN为气候要素处理数据集中的单位像元面积月融雪量,单位为mm;P为气候要素处理数据集中的单位像元面积月降水量,单位为mm;ET为气候要素处理数据集中的单位像元面积月实际蒸散发量,单位为mm;ΔS为气候要素处理数据集中的单位像元面积月土壤水分变化值,单位为mm。
5.根据权利要求2所述的一种基于像元尺度的地下水补给量估算方法,其特征在于,所述根据所述研究区域水资源总产量和所述月径流量数据,估算研究区域的地下水补给量,具体包括:
根据公式G=W-R=S(QS+P-ET-ΔS)-R估算研究区域的地下水补给量;
其中,G为研究区域的地下水补给量,单位为m3;W为研究区域水资源总产量,单位为m3;R为气候要素处理数据集中的月径流量,单位为m3;S为像元面积,单位为m2;QSN为气候要素处理数据集中的单位像元面积月融雪量,单位为mm;P为气候要素处理数据集中的单位像元面积月降水量,单位为mm;ET为气候要素处理数据集中的单位像元面积月实际蒸散发量,单位为mm;ΔS为气候要素处理数据集中的单位像元面积月土壤水分变化值,单位为mm。
6.一种基于像元尺度的地下水补给量估算系统,其特征在于,包括:
信息获取模块,用于获取研究区域的气候要素原始遥感数据集以及研究区域像元面积;所述气候要素原始遥感数据集包括单位像元面积月降水量数据、单位像元面积月实际蒸散发量数据、单位像元面积月融雪量数据、单位像元面积月土壤水分变化值数据以及月径流量数据;
研究区域水资源总产量计算模块,用于根据所述气候要素原始遥感数据集和所述研究区域像元面积,利用水量平衡方程,计算研究区域水资源总产量;
地下水补给量估算模块,用于根据所述研究区域水资源总产量和所述月径流量数据,估算研究区域的地下水补给量。
7.根据权利要求6所述的一种基于像元尺度的地下水补给量估算系统,其特征在于,还包括:
预处理模块,用于对所述气候要素原始遥感数据集中的数据进行预处理,得到研究区域的气候要素处理数据集;预处理的操作步骤依次包括格式转换、影像校正、裁剪、配准、质量检查和投影转换。
8.根据权利要求7所述的一种基于像元尺度的地下水补给量估算系统,其特征在于,所述研究区域水资源总产量计算模块,具体包括:
研究区域水资源总产量计算单元,用于根据公式W=R+G=S(QSN+P-ET-ΔS)计算研究区域水资源总产量;
其中,W为研究区域水资源总产量,单位为m3;R为气候要素处理数据集中的月径流量,单位为m3;G为地下水补给量,单位为m3;S为像元面积,单位为m2;QSN为气候要素处理数据集中的单位像元面积月融雪量,单位为mm;P为气候要素处理数据集中的单位像元面积月降水量,单位为mm;ET为气候要素处理数据集中的单位像元面积月实际蒸散发量,单位为mm;ΔS为气候要素处理数据集中的单位像元面积月土壤水分变化值,单位为mm。
9.根据权利要求7所述的一种基于像元尺度的地下水补给量估算系统,其特征在于,所述地下水补给量估算模块,具体包括:
地下水补给量估算单元,用于根据公式G=W-R=S(QS+P-ET-ΔS)-R估算研究区域的地下水补给量;
其中,G为研究区域的地下水补给量,单位为m3;W为研究区域水资源总产量,单位为m3;R为气候要素处理数据集中的月径流量,单位为m3;S为像元面积,单位为m2;QSN为气候要素处理数据集中的单位像元面积月融雪量,单位为mm;P为气候要素处理数据集中的单位像元面积月降水量,单位为mm;ET为气候要素处理数据集中的单位像元面积月实际蒸散发量,单位为mm;ΔS为气候要素处理数据集中的单位像元面积月土壤水分变化值,单位为mm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010359817.3A CN112529723A (zh) | 2020-04-30 | 2020-04-30 | 一种基于像元尺度的地下水补给量估算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010359817.3A CN112529723A (zh) | 2020-04-30 | 2020-04-30 | 一种基于像元尺度的地下水补给量估算方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112529723A true CN112529723A (zh) | 2021-03-19 |
Family
ID=74978654
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010359817.3A Pending CN112529723A (zh) | 2020-04-30 | 2020-04-30 | 一种基于像元尺度的地下水补给量估算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112529723A (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101248848B1 (ko) * | 2012-09-13 | 2013-04-02 | (주)한서엔지니어링 | 전력소모를 최소화하는 토양 및 지하수 환경의 실시간 모니터링 시스템 |
CN103344344A (zh) * | 2013-07-04 | 2013-10-09 | 武汉珈和科技有限公司 | 一种基于多源遥感数据的全球热异常点信息遥感快速监测及发布方法 |
JP2014062871A (ja) * | 2012-09-24 | 2014-04-10 | Shinshu Univ | 地下水資源の解析方法 |
CN104021283A (zh) * | 2014-05-23 | 2014-09-03 | 清华大学 | 融雪期的日径流量的预测方法和装置 |
CN105912836A (zh) * | 2016-03-25 | 2016-08-31 | 环境保护部卫星环境应用中心 | 一种纯遥感数据驱动的流域水循环模拟方法 |
CN106295194A (zh) * | 2016-08-12 | 2017-01-04 | 中国水利水电科学研究院 | 一种基于流域水量平衡原理的冰川蓄变量评估方法 |
CN108763797A (zh) * | 2018-06-04 | 2018-11-06 | 中国水利水电科学研究院 | 一种基于地下水模型的湖泊与地下水稳定流作用模拟方法 |
CN109035105A (zh) * | 2018-06-15 | 2018-12-18 | 河海大学 | 一种月尺度蒸散发量的定量估算方法 |
CN109872041A (zh) * | 2019-01-17 | 2019-06-11 | 北京师范大学 | 一种基于地质资源环境的评价指标筛选方法 |
-
2020
- 2020-04-30 CN CN202010359817.3A patent/CN112529723A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101248848B1 (ko) * | 2012-09-13 | 2013-04-02 | (주)한서엔지니어링 | 전력소모를 최소화하는 토양 및 지하수 환경의 실시간 모니터링 시스템 |
JP2014062871A (ja) * | 2012-09-24 | 2014-04-10 | Shinshu Univ | 地下水資源の解析方法 |
CN103344344A (zh) * | 2013-07-04 | 2013-10-09 | 武汉珈和科技有限公司 | 一种基于多源遥感数据的全球热异常点信息遥感快速监测及发布方法 |
CN104021283A (zh) * | 2014-05-23 | 2014-09-03 | 清华大学 | 融雪期的日径流量的预测方法和装置 |
CN105912836A (zh) * | 2016-03-25 | 2016-08-31 | 环境保护部卫星环境应用中心 | 一种纯遥感数据驱动的流域水循环模拟方法 |
CN106295194A (zh) * | 2016-08-12 | 2017-01-04 | 中国水利水电科学研究院 | 一种基于流域水量平衡原理的冰川蓄变量评估方法 |
CN108763797A (zh) * | 2018-06-04 | 2018-11-06 | 中国水利水电科学研究院 | 一种基于地下水模型的湖泊与地下水稳定流作用模拟方法 |
CN109035105A (zh) * | 2018-06-15 | 2018-12-18 | 河海大学 | 一种月尺度蒸散发量的定量估算方法 |
CN109872041A (zh) * | 2019-01-17 | 2019-06-11 | 北京师范大学 | 一种基于地质资源环境的评价指标筛选方法 |
Non-Patent Citations (1)
Title |
---|
李鹏;许海丽;潘云;孙颖;王新娟;: "北京市平原区地下水补给量计算方法对比研究", 水文, no. 02, pages 108 - 109 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kjeldsen et al. | Spatial and temporal distribution of mass loss from the Greenland Ice Sheet since AD 1900 | |
US11300709B2 (en) | Method for determining surface runoff yield in vegetation-covered area | |
Jiang et al. | Comparative analysis of different environmental loading methods and their impacts on the GPS height time series | |
Dussaillant et al. | Two decades of glacier mass loss along the Andes | |
US11493496B2 (en) | Method and system for estimating surface runoff based on pixel scale | |
Woodworth et al. | Trends in UK mean sea level revisited | |
Feng et al. | Groundwater storage variations in the North China Plain from GRACE with spatial constraints | |
Berthier et al. | Recent rapid thinning of the “Mer de Glace” glacier derived from satellite optical images | |
US20210341445A1 (en) | Method for monitoring change of vegetation water conservation | |
Spada et al. | Greenland uplift and regional sea level changes from ICESat observations and GIA modelling | |
CN106054283B (zh) | 一种反演上对流层与下平流层风场的方法及装置 | |
Navratil et al. | Requirements of 3D cadastres for height systems | |
CN105069295A (zh) | 基于卡尔曼滤波的卫星以及地面降水测量值同化方法 | |
Sharma et al. | Vertical accuracy of digital elevation model from Shuttle Radar Topographic Mission–a case study | |
Belart et al. | Mass balance of 14 Icelandic glaciers, 1945–2017: spatial variations and links with climate | |
US20210341453A1 (en) | Method and System for Estimating Groundwater Recharge Based on Pixel Scale | |
Ramsankaran et al. | Decadal estimates of surface mass balance for glaciers in Chandra Basin, Western Himalayas, India—A geodetic approach | |
Li et al. | Changes in the glacier extent and surface elevation in Xiongcaigangri region, Southern Karakoram Mountains, China | |
Gärtner-Roer | Sediment transfer rates of two active rockglaciers in the Swiss Alps | |
Kierulf et al. | Time-varying uplift in Svalbard—an effect of glacial changes | |
Suwandana et al. | Thematic information content assessment of the ASTER GDEM: A case study of watershed delineation in West Java, Indonesia | |
CN112529723A (zh) | 一种基于像元尺度的地下水补给量估算方法及系统 | |
Üstün et al. | Height biases of SRTM DEM related to EGM96: from a global perspective to regional practice | |
Kumari et al. | Modelling ice thickness distribution and volume of Patsio Glacier in Western Himalayas | |
CN112529722A (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 |