发明名称 一种测定流域蒸发量的方法
摘要 本发明涉及一种测定流域蒸发量的方法,其包括以下步骤:1)由降水和潜在蒸发数据统计得到流域多年平均月降水量和多年平均月潜在蒸发量,及日平均降水强度;并由土壤类型数据得到流域土壤饱和导水率。2)由多年平均月降水量与多年平均月潜在蒸发量计算得到气候季节性指数,并进一步结合日平均降水强度与土壤饱和导水率计算得到流域特性参数n。3)将参数n与流域多年平均年降水量和多年平均年潜在蒸发量代入流域水热耦合平衡方程,即可得到流域实际蒸发量。本发明所提出的根据流域气候条件及下垫面条件计算流域蒸发量的方法,相对于传统方法具有突出的优势,可广泛用于流域蒸发量的测定过程中。
申请公布号 CN102819050B 申请公布日期 2014.06.04
申请号 CN201210270001.9 申请日期 2012.07.31
申请人 清华大学 发明人 杨汉波;杨大文;吕华芳;胡庆芳
分类号 G01W1/00(2006.01)I 主分类号 G01W1/00(2006.01)I
代理机构 北京纪凯知识产权代理有限公司 11245 代理人 徐宁;宋少娜
主权项 1.一种测定流域蒸发量的方法,其包括以下步骤:1)确定需要进行测定的流域,获取该流域的以下数据:①利用该流域的地形数据生成该流域的边界;②根据该流域内及该流域周边气象站点观测到的逐月降水量及逐月潜在蒸发量,进行空间插值得到逐月降水量及逐月潜在蒸发量的空间分布,并进一步根据该流域范围统计得到该流域的面平均月降水量及面平均月潜在蒸发量,根据面平均月降水量及面平均月潜在蒸发量得到年降水量及年潜在蒸发量,根据多年的年降水量及年潜在蒸发量,进一步得到多年平均年降水量<img file="FDA00004777351200000111.GIF" wi="38" he="52" />和多年平均年潜在蒸发量<img file="FDA00004777351200000112.GIF" wi="48" he="63" />;根据多年平均年降水量<img file="FDA00004777351200000113.GIF" wi="44" he="51" />和多年平均年潜在蒸发量<img file="FDA00004777351200000114.GIF" wi="47" he="63" />进一步得到多年平均月降水量<img file="FDA00004777351200000115.GIF" wi="41" he="58" />和多年平均月潜在蒸发量<img file="FDA00004777351200000116.GIF" wi="84" he="60" />③由该流域气象站的降水数据统计得到日平均降水强度<img file="FDA00004777351200000117.GIF" wi="56" he="70" />④由土壤类型数据;得到该流域的土壤饱和导水率K<sub>s</sub>;2)计算该流域的气候季节性指数:①该流域的月降水量P和月潜在蒸发量E<sub>0</sub>的季节性变化与多年平均月降水量<img file="FDA00004777351200000118.GIF" wi="43" he="56" />和多年平均月潜在蒸发量<img file="FDA0000477735120000011.GIF" wi="70" he="71" />存在如下正弦关系,如公式(4)和公式(5)所示:<maths num="0001"><![CDATA[<math><mrow><mi>P</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mi>m</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>&delta;</mi><mi>P</mi></msub><mi>sin</mi><mi>&omega;t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>E</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mover><mi>E</mi><mo>&OverBar;</mo></mover><mrow><mn>0</mn><mo>,</mo><mi>m</mi></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>&delta;</mi><mi>E</mi></msub><mi>sin</mi><mi>&omega;t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>采用最小二乘法拟合得到δ<sub>P</sub>和δ<sub>E</sub>;其中,δ<sub>P</sub>和δ<sub>E</sub>分别是相对于多年平均月降水量<img file="FDA00004777351200000119.GIF" wi="47" he="58" />和多年平均月潜在蒸发量<img file="FDA0000477735120000014.GIF" wi="70" he="71" />的谐波放大率,t为时间,ω为变化周期,2π/ω=1年;②依据如下公式(6)计算该流域的气候季节性指数S:S=|δ<sub>P</sub>-δ<sub>E</sub>φ|,                    (6)其中<maths num="0003"><![CDATA[<math><mrow><mi>&phi;</mi><mo>=</mo><msub><mover><mi>E</mi><mo>&OverBar;</mo></mover><mrow><mn>0</mn><mo>,</mo><mi>m</mi></mrow></msub><mo>/</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mi>m</mi></msub><mo>;</mo></mrow></math>]]></maths>3)利用流域水热耦合平衡方程计算该流域的蒸发量:①根据如下公式(7)计算流域水热耦合平衡方程中的流域特性参数n:<maths num="0004"><![CDATA[<math><mrow><mi>n</mi><mo>=</mo><msup><mrow><mn>3.083</mn><mi>S</mi></mrow><mrow><mo>-</mo><mn>0.276</mn></mrow></msup><msup><mrow><mo>(</mo><msub><mi>K</mi><mi>s</mi></msub><mo>/</mo><msub><mover><mi>i</mi><mo>&OverBar;</mo></mover><mi>r</mi></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>0.240</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA0000477735120000017.GIF" wi="39" he="78" />为该流域的日平均降水强度;K<sub>s</sub>为该流域的土壤饱和导水率;S为该流域的气候季节性指数;②根据公式(8)计算该流域的蒸发量<img file="FDA00004777351200000120.GIF" wi="58" he="50" /><maths num="0005"><![CDATA[<math><mrow><mover><mi>E</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><mover><mi>P</mi><mo>&OverBar;</mo></mover><msub><mover><mi>E</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub></mrow><msup><mrow><mo>(</mo><msup><mover><mi>P</mi><mo>&OverBar;</mo></mover><mi>n</mi></msup><mo>+</mo><msup><msub><mover><mi>E</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mi>n</mi></msup><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mi>n</mi></mrow></msup></mfrac><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA0000477735120000019.GIF" wi="50" he="62" />为该流域多年平均年降水量;<img file="FDA00004777351200000110.GIF" wi="57" he="74" />为该流域多年平均年潜在蒸发量。
地址 100084 北京市海淀区清华大学水利系