逐日降水的随机模拟模型及其在黑龙江省的应用.pdfVIP

2O11年9月 农 机 化 研 究 第 9期 逐 日降水 的随机模 拟模 型及 其在 黑龙江 省 的应 用 王 斌h ，付 强h ，张金萍 ，徐文仲 (1．东北农业大学 a．水利与建筑学院；b．理学院，哈尔滨 150030；2．黑龙江省高校节水农业重点实验室， 哈尔滨 150030；3．水利部 松辽水利委员会 ，长春 130021) 摘 要：在模拟逐 日降水数据方面 ，近些年来 的一种通用方法是联合应用一阶马尔科夫链和伽玛分布函数建立 随机模拟模型，这种方法在 国外很大的环境范 围内被证 明是有效的。为此 ，应用黑龙江省 14个站点 56年 的实 测降水数据，针对马尔科夫两状态转移概率和伽玛分布参数的不同来源，即统计分析实测降水数据直接获得和 利用实测降水数据建立 回归模型间接推求 ，分别模拟 了各站点的 1000年逐 日降水过程 。结果表明：应用一阶马 尔科夫链和伽玛分布函数可以很好地模拟黑龙江省 14个站点的逐 日降水过程 ，模型的模拟精度很高，证明了这 种方法是有效的；但 国外推求转移概率和伽玛分布参数的经验模型和长系列法在黑龙江省不完全适用，应用重 新建立的回归模型推求的转移概率和伽玛分布参数模拟的逐 日降水数据的精度 同样很高。 关键词：降水事件；随机模型；转移概率；伽玛分布 中图分类号：P426．6 文献标识码 ：A 文章编号：1003-188X【2011)09-0065—05 日降水数据验证这种方法的有效性，以期为评估降水 0 引言 对农业影响、雨水高效利用、降水量频率计算、制定合 降水是农田土壤水分的主要来源之一 ，其数量和 理的灌溉制度和抗旱措施等方面提供科学依据。 分配过程不仅直接影响土壤水分平衡，也间接影响气 1 逐 日降水的随机模拟方法 温、湿度等天气变量，同时也密切影响作物、杂草、害 虫等的生长发育过程。可见 ，降水量及其时空分布对 1．1 马尔科夫链的两状态转移概率 于农业生产具有重要意义。在雨量站的空间布设既 大量研究表明，采用一阶马尔科夫链模型即可很 定的前提下，研究单站降水量的时间分布规律显得尤 好地描述逐 日降水过程-4一 。该模型的两个状态为 为重要。由于降水现象具有随机性特点，应用数理统 降水 日(wetday)和非降水 日(dryday)，状态的转换可 计方法分析降水现象时需要长系列的观测资料。然 以用两个转移概率描述，即由降水日到降水 日的转移 而，我国大部分地区的降水观测数据只有几十年，且 概率P(W／W)和由非降水 日到降水 日的转移概率P 多以逐El数据为主。基于这种短系列的降水观测资 (W／D)。ShuGeng等对7个 (美国5个 ，菲律宾2个) 料，企图解决诸如一年中某时段降水量超过 (或小于) 站点降水数据的统计分析结果表明，在每个站点，各 某数值的概率、一年中某干旱期 (如连续 30天无雨) 月份的P(W／D)，P(w／w)均与同月降水 日出现的频 发生的概率等问题是很难的。 率间存在良好的线性关系，尤其 P(W／D)与，之间相 在模拟产生逐 日降水数据方面，近些年来较通用 关性更强，可用公式表示如下 J 的一种方法是联合应用马尔科夫链 (Markovchain)和 P(W／D)=n+ (1) 伽玛分布函数 (gammadistributionfunction)建立随机 P(W／W)=(1一b)+P(W／D) (2) 模拟模型，即用马尔科夫链描述降水 日的发生，再用 其中，厂为某月降水 日出现频