基于地质统计学理论的海拉尔河流域降水时空变异性研究..docVIP

基于地质统计学理论的海拉尔河流域降水时空变异性研究..doc

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
基于地质统计学理论的海拉尔河流域降水时空变异性研究..doc

基于地质统计学理论的海拉尔河流域降水 时空变异性研究 流域平均降水量的计算方法主要有泰森多边形法、算术 平均法、距离反比法或距离平方反比法、等降水量线法,以 及上世纪70年代发展起来的地质统计学方法,,也简称克立 格法克立格法对于平稳和非平稳地质信息的处理,是较为有 效而精确的方法[1 3],它包括简单克立格法、泛克立格 法、泛克立格法、协同克立格法4种方法。对于水文变量而 言,大多数学者多采用泛克立格法进行计算,但泛克立格法 方法存在缺陷:其一,由于我国水文测站在大江大河上分布 较多,而在中小河流上则稀疏分布,为了更好地研究局部地 区和流域的降水量或其他水文变量的空间变异性,必须通过 周边大量水文信息插补研究区域的数据,工作量较大;其二, 对于小流域特别是水文测站较少或无水文测站的小流域,该 方法在利用邻域水文资料时,由于各变量成因类型不同,破 坏了变量一致性假设。 为了克服克立格法在水文领域应用的局限性,本文在计 算海拉尔河流域上各流域平均降水量时,将多时段泛克立格 空间分析理论[4 6]进一步延伸,采用几何异向性的变 差函数来模拟将多时段水文变量的结构函数,从而考虑了不 同方向上水文数据的分布特点与变异程度,以解决小流域或 局部区域内少站点情况下降水量的最优估计问题。 1多时段泛克立格空间估计理论 首先假定一时空随机函数Z,式中X代表空间坐标,t 为时间维。对流域内的任一测站xi,Z为一随机过程,而Z 在各个时刻的实现则组成一个时间序列Z;多个测站、在 各时刻,Z的观测值则组成一个随机时间序列簇。在有限邻 域内空间上任意两站的时间序列都是相关的,测站距离越大, 其相关程度逐渐减弱。 基本假设 X, y随机过程Z是一阶平稳的,即: E [Z] =mx 本 x、y,随机过程Z、Z的互方差函数C是平稳的, 即: E [ZZ] =mxmy+C t,根据随机过程Z的性质,确定其增量[Z-Z]具 有非平稳的方差函数和非平稳的数学期望[rn-m],且它们都 具有一阶平稳的特性。对 x、y、t: E [Z-Z] =m-m=mx-my 1 2E [Z-Z] 2=r=r 式中:r为变差函数。 设Z可以分解为飘移mx和剩余R两部分,即Z=mx+R, 其中,mx代表长历时现象;R代表短周期内变化。 变差函数的计算 在确定上述假设的前提下,结合各水文测站的时间序列 {Zi; i=l,2,…,N; j=l,2,…,T},推算出实验变差函 数,其公式如下: mi=lTZTj=lZi (tj) Ri(tj)=Zi(tj)-mi 空间点i (xi)处剩余Ri (t)的一个实现用Ri(tj)代替, 则不同两点i、k在时段tj内的剩余Ri (tj)、Rk(tj)可形 成一组数据对,如时段为tj(j=l,2,…,T),则可获得T个 数据对Ri(tj)、Rk(tj),即可估计i、k两点剩余的变差函 数: r(i,k)=12TZTM :Ri(tj)-Rk(tj)] 2 如果结构函数r(i,k)在各向都是同性的,则r(i,k)与i、 k两点的方向没有关系,与两点的距离存在一定关系,则上 式变形为: r(i,k)=r(hik)=T2TEHT [Ri -Rk] 2 式中:r代表空间上每两站的年降水量时间序列。如果 有N个测站,即可获得[N/2]个r。以各水文测站之间的距 离作h为横坐标,以其对应的r为纵坐标,绘制实验变差函 根据对不同水文测站所获取的资料分析可知,受地形、 地貌、地质成因条件等影响,大多数水文变量都呈现出较为 明显的异向性。根据流域或地区内水文测站的分布情况,可 按照前文所列的方法,做出各向的变差函数图,再将各向变 差函数图套合,最终形成套合变差函数。 多时段泛克立格方程组 前文方法求得的变差函数在各个时段内是相同的,但可 根据站年降水量序列{Zi(tj);i=l,2,…,N ; J=l,2,…, T},利用泛克立格空间估计理论,估计任意时段tO、空间任 一点xO处的Z值。其估计式为: Z *=ZNd=l A tOdZd 式中:X tOd为to时段第个测站观测值Zd的估计权重。 假设区域化变量Z的飘移式为: m=kl=0alf1 式中:fl为第1个飘移基函数;al为其系数。 通常,飘移式多采用线性飘移式或二次飘移式。在二维 空间下,也就是: 线性式:m=a0+alxl+a2x2 二次式:m=a0+alxl+a2x2+a3x21+a4x22+a5xlx2 式中:为空间点x的坐标;xl、x2、x21、x22、xlx2均 为基函数;aO为飘移常数因子。 可依据地区实际情况,初步确定飘移式形式。然后根据 数据的空间变化特点或趋势面分析法进一步确定。确定的形 式要通过模型检验,即给定多个飘移式形式,分别代入泛克 立格方程组: ZN3 =1 A tO r+ZKl =

文档评论(0)

ggkkppp + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档