本文主要记录利用不同方法估算水体体积的两篇文献笔记。
(一)
2019年发表在remote sensing for land & resource期刊的一篇介绍湖泊水储量算法的文章,传统利用观测水位与实测水下地形数据相结合的方法估算湖泊水储量,对于资料匮乏的区域很难获取湖泊水储量信息,所以提出一种利用湖泊水面以上的DEM构建高程与面积的关系以及面积与体积增量的关系,实验表明该方法与实测储水量相对误差为9.4%。
构建的纳木错水储量估算模型,可以用于纳木错水储量监测,只要获取湖泊水面面积即可估算水储量。
未来针对不同面积、不同类型的湖泊构建不同的水储量模型,形成一个体系。
原理:湖泊水面以下的湖盆和水面以上的湖盆环境是由相同的地质演化过程形成的,湖泊周围地形的高度的变化特征可能在其表面以下具有类似的变化特征。
(1)建立面积Si与高程Hi的函数关系;
(2)计算体积Vi和体积增量△Vi;构建Si与△Vi的关系;
(3)估算湖底高程;
(4)构建湖泊底下面积和体积的函数关系。
即通过数据构建湖泊上面水面面积与体积的函数关系,利用关系相似性,构建湖泊下面水面面积与体积的关系,进而可以得到水储量体积与面积的函数关系(模型)。
【注】图片源自论文基于地形参数的湖泊水储量估算方法。
数据:STRM DEM 分辨率 90m
MODIS MOD09Q1 时间分辨率8d
验证数据DEM
(二)
2013年发表在RSE期刊的一篇文章,Estimating water volume variations in lakes and reservoirs from four operational satellite altimetry databases and satellite imagery data。
目前可以从四个不同的卫星测高数据库(databases)获得湖泊和水库的水位:
(i)全球水库和湖泊监测(GRLM); http://www.pecad.fas.usda.gov/cropexplorer/global_reservoir/
(ii)河流湖泊水文(RLH);
http://tethys.eaprs.cse.dmu.ac.uk/RiverLake/shared/main
(iii) Hydroweb:提供了世界各地大型河流、大约150个湖泊/水库和湿地的时间序列的水位;
http://www.legos.obs-mip.fr/soa/hydrologie/hydroweb/
(iv) ICESat-GLAS level2 全球地表测高数据产品(ICESat-GLAS):提供土地的表面海拔,包括河流、湖泊和水库;激光足迹地理定位;距离测量;大地测量;仪器和大气校正参数。http://nsidc.org/data/icesat/
如ICESAT/GLAS参考系统为EGM2008 geoid.
【注】图片摘自论文。
本文主要利用databases卫星测高数据库数据结合测高影像数据,在不使用实测数据(in-situ)和深度测量地图(bathymetry maps)的前提下,估算湖泊和水库的储水量。
湖泊或水库中储存的水量取决于inflow流入(即降水precipitation、河流流入、社区和工业的排放和渗漏seepage)与流出outflow之间的平衡(即蒸发evaporation、地下水渗漏percolation、抽取withdraw和河流流出)。
储水量估算发展过程:
(1)传统方法:湖泊或水库的水量是根据现场水位和水深图来估算的。
(2)卫星雷达/激光测高:属于测绘工具,而不是成像设备,只能沿着地面轨迹记录测量结果,而不能真正覆盖全球。
卫星测高通过发射脉冲至地面,并且接收地面反射的回声。不同种类的脉冲产生两种高度计:雷达高度计和激光高度计。
- Radar雷达高度计采用微波脉冲,如ku波段、c波段和s波段;
- Laser激光高度计采用可见光、近红外的激光脉冲。
卫星雷达高度计(T/P,Jason-1,Jason-2,GFO,ENVISAT)地面轨道覆盖详情可见:https://www.aviso.altimetry.fr/en/home.html
界面如下:
估算水体体积的方法:
水体总体体积V = 某固定的最小水体体积Vcon + 随水位变化的水体体积Vvar,即:
- Vcon 可以从地形图中获得(在水库建成前);
- 卫星测高获取的最低水位值,设置为水位参考值以区分Vcon和Vvar;
- 一般研究相对水量变化,而不是绝对值,固不考虑Vcon。
- Water Level Above the Lowest Level (WLALL);
Water Volume Above the Lowest water Level (WVALL);
本文中估算水体体积的步骤:
(1)首先确定各卫星测高产品的最低水位;
(2)然后将各卫星测高产品所得的所有水位减去最低水位,得到最低水位以上的水位(WLALL);
(3)用回归分析方法建立了WLALL与相应表面积(area-WLALL)的关系;
(4)水量是水面面积与水位之间的函数关系的综合,假设面积与水位存在二次多项式的关系,如:
则体积与水位的函数关系为: