CN105930420A - 适用于青藏高原东南缘地区的山区降水测算方法及应用 - Google Patents

适用于青藏高原东南缘地区的山区降水测算方法及应用 Download PDF

Info

Publication number
CN105930420A
CN105930420A CN201610242788.6A CN201610242788A CN105930420A CN 105930420 A CN105930420 A CN 105930420A CN 201610242788 A CN201610242788 A CN 201610242788A CN 105930420 A CN105930420 A CN 105930420A
Authority
CN
China
Prior art keywords
precipitation
river
area
calculating
measuring
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
Application number
CN201610242788.6A
Other languages
English (en)
Inventor
苏鹏程
韦方强
董丹丹
刘晶晶
孙铭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Mountain Hazards and Environment IMHE of CAS
Original Assignee
Institute of Mountain Hazards and Environment IMHE of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Mountain Hazards and Environment IMHE of CAS filed Critical Institute of Mountain Hazards and Environment IMHE of CAS
Priority to CN201610242788.6A priority Critical patent/CN105930420A/zh
Publication of CN105930420A publication Critical patent/CN105930420A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29—Geographical information databases
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • 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
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Algebra (AREA)
  • Remote Sensing (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种适用于青藏高原东南缘地区的山区降水测算方法及应用。所述山区降水测算方法首先获取拟进行降水测算山区的海拔、坡向和坡度，然后根据拟进行降水测算山区的地理位置确定其所在分区，最后应用所属分区的降水垂直分布关系式计算得到拟进行降水测算山区的降水量。本发明方法原理可靠，计算过程科学简便，能够普遍适用于获取位于青藏高原东南缘地区内缺少气象观测站的高海拔地区的降水数据，对泥石流等山地灾害形成区的降水量估算和汇流计算，具有重要的参考价值。

Description

适用于青藏高原东南缘地区的山区降水测算方法及应用
技术领域
本发明涉及一种缺气象观测站点的山区降水计算方法，特别是涉及一种适用于青藏高原东南缘地区、缺气象观测站点的高海拔山区降水测算方法及其应用。
背景技术
降水是大气中的水的相变，水汽凝结成雨雪等过程。从降水的形成条件可以看出，除大气环流外，降水的垂直分布受地理环境影响很大，影响山区降水的地形因子主要是：高大山脉、地形、海拔高度、坡地方位等。多年来，山区降水一直是研究的难点。
针对山区降水，目前众多学者主要采用：①空间插值法。常用的有反距离加权插值法、全局多项式插值法、局部多项式插值法、径向基函数插值法、克里金插值等方法。插值法不能充分考虑气象要素空间分布与诸多地理环境要素之间的复杂函数关系，其插值能力差，误差较大。②统计模型法。即根据实测站点信息，建立降水量与地理位置、地形及气象等因子间的关系，即研究降水与经度、纬度、海拔高度、坡度、坡向等因子的关系，建立相应的空间分布方程，能定量反映地形影响下的降水空间分布规律。但目前文献中的统计模型方法主要是针对3000m以下的降水进行分析，而泥石流等灾害的形成区海拔高度常远超这个范围，因此现有的统计模型方法在使用范围方面具有一定的局限，并且由于需要的因子较多，基础数据的收集也相对困难。
青藏高原东南缘地区分布着近于平行的、近南北走向的山脉和水系，在这样的地形条件下，影响山区降水的垂直分布特征的地形因子当中，海拔高度、季风活动所形成的迎风坡和背风坡、坡度是最主要因素。此区域是泥石流等山地灾害的高发区，要对泥石流等山地灾害形成区的汇流进行计算，就需要充足的降水数据提供参考。但目前我国山区的气象站和雨量站多设在河谷低处，在海拔较高的地方雨量测点很少，急需适用于青藏高原东南缘地区、缺气象观测站点的高海拔山区降水测算方法。
发明内容
本发明的目的就是针对青藏高原东南缘高海拔地区缺少气象观测站、无法获得降水数据的不足，提供一种适用于青藏高原东南缘地区的山区降水测算方法及应用，该方法根据缺气象观测站点的山区所在地理位置、海拔高度、坡向和坡度特征用降水垂直分布关系式计算降水，能够实现获取缺气象观测站点的山区降水数据。
为实现上述目的，本发明的技术方案是：
本发明提出一种适用于青藏高原东南缘地区的山区降水测算方法，其技术思路是：在同一个大山脉内和相似地形形态下，山区降水主要受海拔高度、坡向和坡度的影响，因此，将青藏高原东南缘地区分为6个分区，同一个分区内(即在同一个大山脉内和相似地形形态下)，存在一个比较适用的降水垂直分布关系式P＝ao+a1×H+a2×H2+a3×A+a4×A2+a5×B+a6×B2，该关系式根据已有气象观测站点多年年平均降水量与海拔、坡向、坡度的统计关系确定；式中，P为降水量(mm)、由气象观测站点降雨资料得到，H为海拔(m)、由实际测量得到，A为坡向(rad)、由DEM提取得到，B为坡度(rad)、由DEM提取得到，a0、a1、a2、a3、a4、a5、a6为系数(不同分区内，其关系式的系数不同)。通过已有气象观测站点的降水量与海拔高度、坡向、坡度建立的统计关系式(6个分区所对应的降水垂直分布关系式有所不同)，应用于青藏高原东南缘地区、缺气象观测站点的山区计算降水；对于缺气象观测站点的山区，分析其基本特征，确定其在青藏高原东南缘的分区，根据建立的相应分区的降水垂直分布关系式计算得到其降水数据。
上述技术思路的技术依据主要基于：第一，前人建立的降水随海拔高度变化的数学模型中，具有代表性的是傅抱璞根据半经验理论建立的高斯模式，其本质依然为类抛物线模式；在类抛物线模式的基础上，海拔高度、坡向、坡度为主因子，建立降水垂直分布的非线性回归模型。第二，降水垂直分布关系式只适用于青藏高原东南缘的缺测地区，此区的主体为横断山脉中北段，分布着彼此平行的、近南北走向的山脉和水系，根据山脉和流域进行分区，同一区域内，能够建立降水垂直分布关系式。第三，根据缺气象观测站点的山区的基本特征，确定其所在的分区，根据降水垂直分布关系式，能够得到其降水数据。
所述降水垂直分布关系式的具体推导过程如下：
第一步，采用滑动平均法滤掉数据中频繁起伏的随机误差，探究青藏高原东南缘降水随海拔高度的变化趋势。即对研究区范围内所有站点的降水数据，按海拔高度从高到低的变化排序，进行11步滑动平均。用平滑后的降水与海拔高度数据序列进行作图。结果显示，青藏高原东南缘降水垂直变化并不是简单的随高度线性减少，而是为抛物线型的变化趋势，即考虑降水P与海拔高度H的关系式为P＝ao+a1×H+a2×H2，a0、a1、a2为系数。
第二步，在确定青藏高原东南缘降水垂直变化为抛物线型的变化趋势的基础上，考虑季风活动所形成的迎风坡、背风坡和坡度对降水的影响，即引入坡向A和坡度B这两个影响因子。
第三步，在前两步的基础上，将青藏高原东南缘地区分为6个区域，同一个分区内(即在同一个大山脉内和相似地形形态下)，经相关性检验，存在一个比较适用的降水垂直分布关系式，即P＝ao+a1×H+a2×H2+a3×A+a4×A2+a5×B+a6×B2。式中，P为降水量(mm)、由气象观测站点降雨资料得到，H为海拔(m)、由实际测量得到，A为坡向(rad)、由DEM提取得到，B为坡度(rad)、由DEM提取得到，a0、a1、a2、a3、a4、a5、a6为系数。
第四步，利用spss软件，将各个分区的降水数据P与相应的海拔高度H、坡向A和坡度B数据，建立分区降水垂直分布关系式。
所述青藏高原东南缘地区的边界为：西与雅鲁藏布江流域的林芝、波密一带相连，东部边界是北起龙门山，往南经大渡河、小相岭、锦屏山，到高黎贡山一线，北部边界为曲麻莱—石渠—甘德一线。在青藏高原东南缘地区，根据空间分布连续性原则，在保证区域气候特征相对相似性基础上，参照地形(分水岭)划定界线，从西到东，进行分区，具体分为6个区域，即怒江流域、澜沧江流域、金沙江流域、雅砻江流域、大渡河流域和岷江流域。具体而言，所述适用于青藏高原东南缘地区的山区降水测算方法步骤如下：
(一)通过现场GPS(Global Positioning System，全球定位系统)实地测量或利用GIS(Geographic Information System，地理信息系统)技术从数字高程模型(DEM)上提取(即先建立DEM，然后利用GIS技术从DEM上提取)，确定拟进行降水测算的山区的海拔高度H，单位m；通过建立数字高程模型(DEM)，确定拟进行降水测算的山区海拔高度H处的坡向A和坡度B，单位均为rad。
(二)根据拟进行降水测算的山区的所在地理位置、通过河流水系基本单元的地貌特征，结合利用GIS技术从数字高程模型(DEM)提取的地形、流域形态与分水岭的空间分布，确定该山区是属于怒江流域，或澜沧江流域，或金沙江流域，或雅砻江流域，或大渡河流域，或岷江流域。所述怒江流域为米林-洛隆-福贡-泸水一线，所述澜沧江流域为杂多-昌都-兰坪一线，所述金沙江流域为治多-德格-中甸一线，所述雅砻江流域为甘孜-雅江-盐源一线，所述大渡河流域为甘德-壤塘-泸定一线，所述岷江流域为若尔盖-黑水-理县一线。
(三)如果拟进行降水测算的山区属于怒江流域，则将步骤(一)中得到的海拔高度H(单位m)、坡向A(单位rad)和坡度B(单位rad)代入公式P＝1443.187-0.263×H+9.79×10-6×H2+361.33×A-47.297×A2-3745.491×B+4039.651×B2计算得到拟进行降水测算的山区的降水量P，单位mm；
如果拟进行降水测算的山区属于澜沧江流域，则将步骤(一)中得到的海拔高度H(单位m)、坡向A(单位rad)和坡度B(单位rad)代入公式
P＝4242.051-2.215×H+2.22×10-4×H2+
1949.504×A-352.076×A2-9216.931×B+12378.636×B2
计算得到拟进行降水测算的山区的降水量P，单位mm；
如果拟进行降水测算的山区属于金沙江流域，则将步骤(一)中得到的海拔高度H(单位m)、坡向A(单位rad)和坡度B(单位rad)代入公式P＝-1250.575+0.982×H-1.41×10-4×H2+7.465×A+0.479×A2+484.695×B-488.259×B2计算得到拟进行降水测算的山区的降水量P，单位mm；
如果拟进行降水测算的山区属于雅砻江流域，则将步骤(一)中得到的海拔高度H(单位m)、坡向A(单位rad)和坡度B(单位rad)代入公式P＝3063.093-1.391×H+2.06×10-4×H2-24.288×A+1.996×A2-236.562×B+713.528×B2计算得到拟进行降水测算的山区的降水量P，单位mm；
如果拟进行降水测算的山区属于大渡河流域，则将步骤(一)中得到的海拔高度H(单位m)、坡向A(单位rad)和坡度B(单位rad)代入公式P＝55.214+0.579×H-1.09×10-4×H2-87.343×A+13.935×A2+437.648×B-637.494×B2计算得到拟进行降水测算的山区的降水量P，单位mm；
如果拟进行降水测算的山区属于岷江流域，则将步骤(一)中得到的海拔高度H(单位m)、坡向A(单位rad)和坡度B(单位rad)代入公式P＝440.453+0.027×H-8.437×10-7×H2+187.412×A-32.139×A2+558.147×B-535.386×B2计算得到拟进行降水测算的山区的降水量P，单位mm。
本发明方法主要针对青藏高原东南缘内高海拔地区缺少气象观测站点导致基础气候特征缺乏数据支撑的现状，以年降雨量为突破口，利用海拔高度、坡向和坡度等指标，寻找青藏高原东南缘年降水量垂直分布的计算方法，为提高该地区降水量时空分布的精度，更好的服务于该区域防灾减灾工作提供数据支撑，特别是对泥石流等山地灾害形成区降水梯度分布的计算尤为有效。
与现有技术相比，本发明的有益效果是：在青藏高原东南缘，以海拔、坡向、坡度为降水的主因子，通过历史降水资料与海拔、坡向、坡度的统计关系，建立非线性回归模型，进而得到比较适用的降水垂直分布关系式(不同分区内降水垂直分布关系式有所不同)，通过计算得到缺气象观测站点的山区降水数据，对泥石流等山地灾害形成区的降水量估算和汇流计算，具有重要的参考价值。本发明方法原理可靠，计算过程科学简便，能够普遍适用于获取位于青藏高原东南缘地区内缺少气象观测站的高海拔地区的降水数据，且计算结果代表的是整个区域的年降水量平均值水平。
附图说明
图1是实施例一中得到的江咀沟降水梯度分布测算结果图。其中，加正方形标记点的实线表示主沟降水梯度分布，加圆形标记点的虚线表示支沟降水梯度分布。
图2是实施例二中得到的拉哇沟降水梯度分布测算结果图。
具体实施方式
下面结合附图，对本发明的优选实施例作进一步的描述。
实施例一
如图1所示。江咀沟位于大渡河左岸,沟口距离大渡河黄金坪水电站5km,四川甘孜州康定县管辖。江咀沟源头海拔约3800m,沟口海拔约1420m,主沟长约19.37km，具有流域面积大、纵坡降缓、径流区长、支沟发育等特点,具备典型的泥石流沟地形地貌特征。江咀沟主沟为瓜达沟，在其流域范围内沟道较为曲折，主沟流向变动很大，在下火地附近由南东向转为近东西向，流经一段距离后，转为南西向，汇入大渡河。清水汇流区为3800m-2800m，物源区为2800-1800m，1800m以下为流通区和堆积区，麦笨乡上游的两河口至沟源植被发育。支沟磨子沟沟道顺直,大体为北西流向,发育3条支沟,该沟长约11.64km,沟源海拔约4264m,于高程1905m汇入主沟。根据卫星遥感图像显示，主沟瓜达沟植被覆盖度较高的地带位于海拔2700-3000m左右，支沟磨子沟植被覆盖度较高的地带位于海拔2600-3000m左右，大体上较好反应当地降水最大高度带的分布规律。
为了对江咀沟流域泥石流形成区的降水梯度分布进行测算，采用本发明的适用于青藏高原东南缘地区的山区降水测算方法对江咀沟流域进行模型验证，计算其降水量的梯度分布，并验证其最大降水高度带分布的区间。所述山区降水测算方法步骤如下：
第一步，利用GIS技术从数字高程模型上提取，确定拟进行降水测算的江咀沟的主沟和支沟海拔高度H为2000-4000m，每隔100m取一个测算点；通过建立数字高程模型，确定江咀沟相应测算点处的坡向A和坡度B，分别如下表1(主沟)和表2(支沟)所示。
第二步，根据江咀沟所在地理位置、通过河流水系基本单元的地貌特征，结合GIS技术从数字高程模型提取的地形、流域形态与分水岭的空间分布，确定江咀沟属于大渡河流域。
第三步，将第一步中得到的海拔高度H、坡向A和坡度B代入公式
P＝55.214+0.579×H-1.09×10-4×H2-87.343×A+13.935×A2+437.648×B-637.494×B2
计算得到江咀沟主沟和支沟海拔2000-4000m之间相应测算点处的年降水量P，具体计算结果分别如下表1(主沟)和表2(支沟)所示。
表1江咀沟主沟测算点
如上表1所示，即绘制得到江咀沟主沟的降水梯度分布测算结果图(图1)。可以看出，最大降水高度带位于2700m左右，与实际植被分布总体上比较吻合。且由于所取测算点的具体位置不同，降水量变化不同，如2800m处所取测算点处于沟内拐角较深位置，被高山遮挡较为严重，所以降水量急剧减少，如3700m、3800m处所取测算点植被状况比3600m好，降水量比3600m略微增加是符合事实的。
海拔高度H(m) 经度(°) 纬度(°) 坡向A(rad) 坡度B(rad) 降水量P(mm)
2000 102.226 30.221 6.223 0.296 847.029
2100 102.236 30.221 4.910 0.412 769.636
2200 102.240 30.219 6.262 0.175 857.976
2300 102.247 30.219 4.257 0.194 751.974
2400 102.250 30.221 3.611 0.176 740.651
2500 102.254 30.221 4.985 0.255 802.516
2600 102.258 30.220 6.264 0.116 865.669
2700 102.262 30.219 3.746 0.314 766.805
2800 102.266 30.217 6.048 0.347 878.439
2900 102.269 30.214 4.948 0.364 801.434
3000 102.272 30.213 6.251 0.590 846.076
3100 102.275 30.211 4.991 0.469 778.910
3200 102.277 30.210 4.300 0.297 747.661
3300 102.282 30.208 5.928 0.391 824.431
3400 102.285 30.206 0.245 0.249 812.642
3500 102.286 30.202 1.284 0.332 732.349
3600 102.290 30.202 4.864 0.313 706.340
3700 102.293 30.201 5.521 0.263 718.790
3800 102.296 30.200 5.572 0.269 699.048
3900 102.297 30.197 0.322 0.438 698.102
4000 102.303 30.196 5.594 0.230 641.535
表2江咀沟支沟测算点
如上表2所示，即绘制得到江咀沟支沟的降水梯度分布测算结果图(图1)。可以看出，最大降水高度带位于2800m左右，与实际植被分布密集区域大体上吻合。且由于所取测算点的具体位置不同，降水量变化不同，如支沟2000m处测算点处于沟内平坦地区，四周无高山遮挡，降水量较大，且对比主沟2000m测算点，其位于高山的阴影内，两者的海拔高度一致，坡度相差较小，但坡向相差较多，降水量相差较大。
实施例二
如图2所示。拉哇沟是金沙江左岸的Ι级支流。该沟口坐标为30°5′2.24″N，99°2′40.09″E，流域总体呈EW流向，并以73°夹角汇入金沙江。主沟长度8.13km，流域面积27.77km2。整个流域上游的地形坡度较为平缓，中、下游地形较为陡峭，尤其是沟道下游左岸地形十分陡峭，为崩塌等的发育提供了有利条件。拉哇沟源头高程4374m，属当地Ⅱ级夷平面(海拔4200-4400m)，沟口海拔2532m，流域内大部分海拔高度在3000m以上，相对高差1842m。谷地2800m以下为干热河谷小叶灌丛，2800-3100m为山地稀树小叶灌丛，再向上逐步进入寒温性常绿针叶林、高山杜鹃灌丛与高山草甸，3000-4000m的阳坡或半阳坡上以灌丛为主，阴坡或半阴坡上为亚高山针叶林为主，3800-4400m处为高山草甸植被，4700m左右为流石滩稀疏植被，4900-5000m为冰雪带。
为了对拉哇沟泥石流形成区的降水梯度分布进行测算，采用本发明的适用于青藏高原东南缘地区的山区降水测算方法对拉哇沟流域进行模型验证，计算其降水量的梯度分布，并根据植被的分布验证其最大降水高度带分布的区间。所述山区降水测算方法步骤如下：
第一步，通过现场GPS实地测量，确定拟进行降水测算的拉哇沟海拔高度H为2700-4300m，每个100m取一个测算点；通过建立数字高程模型，确定拉哇沟相应测算点处的坡向A和坡度B，分别如下表3所示。
第二步，根据拉哇沟所在地理位置、通过河流水系基本单元的地貌特征，结合GIS技术从数字高程模型提取的地形、流域形态与分水岭的空间分布，确定拉哇沟属于金沙江流域。
第三步，将第一步中得到的海拔高度H、坡向A和坡度B代入公式
P＝-1250.575+0.982×H-1.41×10-4×H2+7.465×A+0.479×A2+484.695×B-488.259×B2
计算得到拉哇沟海拔2700-4300m之间相应测算点处的年降水量P，具体计算结果如下表3所示。
海拔高度H(m) 经度(°) 纬度(°) 坡向A(rad) 坡度B(rad) 降水量P(mm)
2700 99.048 30.086 3.990 0.324 515.720
2800 99.051 30.086 5.268 0.310 549.065
2900 99.055 30.088 4.231 0.277 547.933
3000 99.055 30.088 4.213 0.277 562.722
3100 99.061 30.088 5.808 0.485 617.922
3200 99.064 30.088 0.043 0.333 555.163
3300 99.067 30.088 4.633 0.188 572.893
3400 99.070 30.088 5.464 0.454 632.294
3500 99.074 30.089 3.644 0.170 560.524
3600 99.077 30.087 0.111 0.201 535.425
3700 99.089 30.083 3.592 0.253 576.472
3800 99.097 30.084 3.387 0.301 576.805
3900 99.103 30.088 3.954 0.379 584.696
4000 99.107 30.090 3.182 0.237 536.801
4100 99.110 30.092 3.907 0.314 545.362
4200 99.112 30.094 3.661 0.334 527.235
4300 99.115 30.096 2.745 0.306 491.076
表3拉哇沟测算点
如上表3所示，即绘制得到拉哇沟的降水梯度分布测算结果图(图2)。可以看出，总体降水的分布随着海拔高度的增加在3400m左右达到最大值，然后又呈下降趋势，这与植被分布总体上比较吻合。且由于测算点的选取具有一定的随机性，在相同海拔高度上其坡度、坡向对降水量也有一定影响。

Claims (4)

1.一种适用于青藏高原东南缘地区的山区降水测算方法，其特征在于：所述山区降水测算方法步骤如下：
(一)通过现场GPS实地测量或利用GIS技术从数字高程模型上提取，确定拟进行降水测算的山区的海拔高度H，单位m；通过建立数字高程模型，确定拟进行降水测算的山区海拔高度H处的坡向A、单位rad，和坡度B、单位rad；
(二)根据拟进行降水测算的山区的所在地理位置、通过河流水系基本单元的地貌特征，结合GIS技术从数字高程模型提取的地形、流域形态与分水岭的空间分布，确定该山区是属于怒江流域，或澜沧江流域，或金沙江流域，或雅砻江流域，或大渡河流域，或岷江流域；
(三)如果拟进行降水测算的山区属于怒江流域，则将步骤(一)中得到的海拔高度H、坡向A和坡度B代入公式
P＝1443.187-0.263×H+9.79×10-6×H2+361.33×A-47.297×A2-3745.491×B+4039.651×B2计算得到拟进行降水测算的山区的降水量P，单位mm；
如果拟进行降水测算的山区属于澜沧江流域，则将步骤(一)中得到的海拔高度H、坡向A和坡度B代入公式
P＝4242.051-2.215×H+2.22×10-4×H2+
1949.504×A-352.076×A2-9216.931×B+12378.636×B2
计算得到拟进行降水测算的山区的降水量P，单位mm；
如果拟进行降水测算的山区属于金沙江流域，则将步骤(一)中得到的海拔高度H、坡向A和坡度B代入公式
P＝-1250.575+0.982×H-1.41×10-4×H2+7.465×A+0.479×A2+484.695×B-488.259×B2
计算得到拟进行降水测算的山区的降水量P，单位mm；
如果拟进行降水测算的山区属于雅砻江流域，则将步骤(一)中得到的海拔高度H、坡向A和坡度B代入公式
P＝3063.093-1.391×H+2.06×10-4×H2-24.288×A+1.996×A2-236.562×B+713.528×B2
计算得到拟进行降水测算的山区的降水量P，单位mm；
如果拟进行降水测算的山区属于大渡河流域，则将步骤(一)中得到的海拔高度H、坡向A和坡度B代入公式
P＝55.214+0.579×H-1.09×10-4×H2-87.343×A+13.935×A2+437.648×B-637.494×B2
计算得到拟进行降水测算的山区的降水量P，单位mm；
如果拟进行降水测算的山区属于岷江流域，则将步骤(一)中得到的海拔高度H、坡向A和坡度B代入公式
P＝440.453+0.027×H-8.437×10-7×H2+187.412×A-32.139×A2+558.147×B-535.386×B2
计算得到拟进行降水测算的山区的降水量P，单位mm。
2.根据权利要求1所述适用于青藏高原东南缘地区的山区降水测算方法，其特征在于：所述青藏高原东南缘地区的边界为：西与雅鲁藏布江流域的林芝、波密一带相连，东部边界是北起龙门山，往南经大渡河、小相岭、锦屏山，到高黎贡山一线，北部边界为曲麻莱—石渠—甘德一线。
3.根据权利要求1所述适用于青藏高原东南缘地区的山区降水测算方法，其特征在于：所述怒江流域为米林—洛隆—福贡—泸水一线，所述澜沧江流域为杂多—昌都—兰坪一线，所述金沙江流域为治多—德格—中甸一线，所述雅砻江流域为甘孜—雅江—盐源一线，所述大渡河流域为甘德—壤塘—泸定一线，所述岷江流域为若尔盖—黑水—理县一线。
4.如权利要求1所述适用于青藏高原东南缘地区的山区降水测算方法的应用，其特征在于：适用于山地灾害形成区降水梯度分布的测算。
CN201610242788.6A 2016-04-19 2016-04-19 适用于青藏高原东南缘地区的山区降水测算方法及应用 Pending CN105930420A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610242788.6A CN105930420A (zh) 2016-04-19 2016-04-19 适用于青藏高原东南缘地区的山区降水测算方法及应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610242788.6A CN105930420A (zh) 2016-04-19 2016-04-19 适用于青藏高原东南缘地区的山区降水测算方法及应用

Publications (1)

Publication Number Publication Date
CN105930420A true CN105930420A (zh) 2016-09-07

Family

ID=56838412

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610242788.6A Pending CN105930420A (zh) 2016-04-19 2016-04-19 适用于青藏高原东南缘地区的山区降水测算方法及应用

Country Status (1)

Country Link
CN (1) CN105930420A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109829027A (zh) * 2019-01-29 2019-05-31 中国科学院、水利部成都山地灾害与环境研究所 区域泥石流前期有效降水量测算方法
CN110619109A (zh) * 2019-08-01 2019-12-27 中国林业科学研究院 一种高寒地区降水相态识别方法
CN111123408A (zh) * 2019-11-26 2020-05-08 深圳震有科技股份有限公司 一种基于gis预测降水分布的方法、系统及存储介质
CN112729768A (zh) * 2021-01-20 2021-04-30 中国科学院、水利部成都山地灾害与环境研究所 一种泥石流入汇主河运动过程实验测量系统及测算方法
CN113536521A (zh) * 2020-09-27 2021-10-22 石河子大学 一种天山北坡中低海拔地区的年平均降雨量推算方法
CN116049341A (zh) * 2023-03-08 2023-05-02 北京七兆科技有限公司 一种水文数据标准化方法、装置、设备及存储介质
CN116150570A (zh) * 2022-11-07 2023-05-23 浙江省水文管理中心 一种高精度网格降雨的存储和叠加分析计算方法及装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101256590A (zh) * 2008-04-03 2008-09-03 北京艺龙天地文化传播有限公司 结合WebGis的三维在线虚拟现实生态环境仿真系统及其方法
CN102955863A (zh) * 2011-08-17 2013-03-06 长江水利委员会长江科学院 一种基于分布式水文模拟的旱情评估和预报模型方法
CN103218485A (zh) * 2013-04-03 2013-07-24 郑州大学 一种gis技术支持下的小流域环境容量计算方法
CN103364532A (zh) * 2012-03-28 2013-10-23 中国环境科学研究院 一种定量评估区域水源涵养量的实用方法
CN104537083A (zh) * 2014-12-31 2015-04-22 北京奇虎科技有限公司 一种基于电子地图的降水预测方法及装置
CN104851048A (zh) * 2015-05-28 2015-08-19 四川农业大学 一种复杂山区烤烟适宜移栽期空间分布的确定方法
CN105160465A (zh) * 2015-08-25 2015-12-16 四川省宇行规划设计有限公司 一种基于3s技术的退耕还林方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101256590A (zh) * 2008-04-03 2008-09-03 北京艺龙天地文化传播有限公司 结合WebGis的三维在线虚拟现实生态环境仿真系统及其方法
CN102955863A (zh) * 2011-08-17 2013-03-06 长江水利委员会长江科学院 一种基于分布式水文模拟的旱情评估和预报模型方法
CN103364532A (zh) * 2012-03-28 2013-10-23 中国环境科学研究院 一种定量评估区域水源涵养量的实用方法
CN103218485A (zh) * 2013-04-03 2013-07-24 郑州大学 一种gis技术支持下的小流域环境容量计算方法
CN104537083A (zh) * 2014-12-31 2015-04-22 北京奇虎科技有限公司 一种基于电子地图的降水预测方法及装置
CN104851048A (zh) * 2015-05-28 2015-08-19 四川农业大学 一种复杂山区烤烟适宜移栽期空间分布的确定方法
CN105160465A (zh) * 2015-08-25 2015-12-16 四川省宇行规划设计有限公司 一种基于3s技术的退耕还林方法

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109829027A (zh) * 2019-01-29 2019-05-31 中国科学院、水利部成都山地灾害与环境研究所 区域泥石流前期有效降水量测算方法
CN110619109A (zh) * 2019-08-01 2019-12-27 中国林业科学研究院 一种高寒地区降水相态识别方法
CN111123408A (zh) * 2019-11-26 2020-05-08 深圳震有科技股份有限公司 一种基于gis预测降水分布的方法、系统及存储介质
CN111123408B (zh) * 2019-11-26 2022-02-18 深圳震有科技股份有限公司 一种基于gis预测降水分布的方法、系统及存储介质
CN113536521A (zh) * 2020-09-27 2021-10-22 石河子大学 一种天山北坡中低海拔地区的年平均降雨量推算方法
CN112729768A (zh) * 2021-01-20 2021-04-30 中国科学院、水利部成都山地灾害与环境研究所 一种泥石流入汇主河运动过程实验测量系统及测算方法
CN112729768B (zh) * 2021-01-20 2023-08-15 中国科学院、水利部成都山地灾害与环境研究所 一种泥石流入汇主河运动过程实验测量系统及测算方法
CN116150570A (zh) * 2022-11-07 2023-05-23 浙江省水文管理中心 一种高精度网格降雨的存储和叠加分析计算方法及装置
CN116150570B (zh) * 2022-11-07 2024-04-23 浙江省水文管理中心 一种高精度网格降雨的存储和叠加分析计算方法及装置
CN116049341A (zh) * 2023-03-08 2023-05-02 北京七兆科技有限公司 一种水文数据标准化方法、装置、设备及存储介质
CN116049341B (zh) * 2023-03-08 2023-08-15 北京七兆科技有限公司 一种水文数据标准化方法、装置、设备及存储介质

Similar Documents

Publication Publication Date Title
CN105930420A (zh) 适用于青藏高原东南缘地区的山区降水测算方法及应用
Singh et al. Topographical influence on precipitation distribution in different ranges of western Himalayas
Bendix et al. Cloud occurrence and cloud properties in Ecuador
Yao et al. The mass elevation effect of the Tibetan Plateau and its implications for Alpine treelines.
Bañares et al. Seasonal and diurnal variations of observed convective rain events in Metro Manila, Philippines
Li et al. Glacier extent changes and possible causes in the Hala Lake Basin of Qinghai-Tibet Plateau
Mätlik et al. Synoptic weather types that have caused heavy precipitation in Estonia in the period 1961–2005
Thomas The climate of the Gongga Shan range, Sichuan province, PR China
CN105912860A (zh) 基于gis的林冠截留模型
Loboda et al. Maximum flow of rivers of the Ukrainian Carpathians (in the upper Dniester) in the climate change conditions
CN116542009A (zh) 一种大尺度陆面水文过程模拟系统及方法
Tarasova et al. Impact of air transport on seasonal variations and trends of surface ozone at Kislovodsk high mountain station
Yembuu Climate and climate change of Mongolia
Kanno Classification of the Yamase (cold northeasterly wind around northeastern Japan) based upon its air-mass vertical structures
Li et al. Mountain effect and differences in storm floods between northern and southern sources of the Songhua River Basin
Zhang et al. Why are glacial lakes in the eastern Tianshan Mountains expanding at an accelerated rate?
CN110543542B (zh) 用于变电设备的暴雨淹没灾害分布图的绘制方法及系统
Cogliati et al. Air flow analysis in the upper Rio Negro Valley (Argentina)
Kejna Topoclimatic conditions in the vicinity of the Arctowski Station (King George Island, Antarctica) during the summer season of 2006/2007
Nojarov et al. Influence of atmospheric circulation on the spatial distribution of precipitation in the area of Sofia city
Kim et al. Mapping monthly temperature normals across North Korea at a landscape scale
Enea et al. Multi-criterial gis analysis for identifying optimum location for vineyard placement. Case Study: Moldova Region, Romania
Korjenić Climate as spatial planning factor of the Una Sana Canton, Bosnia and Herzegovina
Yang et al. The Short-Duration Heavy Rainfall in Different Quadrants of Northeast China Cold Vortices
Planchon A study of the coastal climates in France using temperature and precipitation data (1961–1990)

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20160907

RJ01 Rejection of invention patent application after publication