气象干旱潜在蒸散量的计算方法、使用参数法计算短序列气象观测站气象干旱综合指数的方法.docx

附录A（规经性附录）降水量距平百分率的计算方法降水量距平百分率反映某时段降水It与同期平均状态的儡离程度,按式（A,1）计算:P-PPA=-J=-X100JV（.1）式中：PA某时段降水股距平百分率,用先表示：P某时段降水量，值位为墨米（丽）：P计匏时段同期气候平均降水量，舱位为:至米（ma），按式（A.2）计算.W（a-2）式中：n般取3。，指3011（月或年）：一一计算时段第i日（月或年降水电，单位为亳米（mm）.附录B（规范性附录）相对湿泄度指数的计算方法相对湿润度指数为某段时间的降水M与同时段内潜在然敌峡之差再除以同时段内潜在蒸敢呆得到的指数，按式（BJ）计算；式中：MI某时段相对湿润度：P某时段的降水盘，单位奏米<11m）:PET某时段的潜在热散HIFAOPenman-Mon1.eiIh或ThornthWaite方法计算,单位/米（rr<n）.附录C(规经性附录)气象干旱潜在蒸散量的计算方法C.1滞在IKtUt的计算本标准推荐两种方法计算潜在然放依，即Thornthwaite方法和FAoPenman-NJonteith方法.FAOPenman-Monteith方法计算误差小，但需要的气象要素多.Thornthuaite方法计算相对简单,需要的气象要常少,但有一定的局限性.使用者请根据资料条件选择合适的计算方法.C.2TbornthwaiteThornthwaite方法求算潜在蒸散玳是以月平均温慢为主要依据，并考虑纬度因子<11照长度建立的经验公式，需要输入的因子少，计算方法简单：PET=16.0X(i-/<C.1)式中：PET潜在蒸敢量，此处是指月的潜在蒸散地，单位为亳米每月(mm/Bon)；Ti月的平均气温，单位为摄氏度CC);H年热冰指数：A常数.各月热At指数,由下式汁豫：”,=百严(C.2)年热量指数H为：/=9,=*（"产<C.3）ZI-I5常数A由下式计算:A=6.75XIO77'-7.71XIO52+1.792X10:H+0.49（C.1>当月平均气温TWOC时，月热此指数，=0,潜在然敌依用r=6>mmon).C.3FAOPenaan-Monteith方法C.3.1FAOPenBanTonteith方法介绍FOPenmanUonteith方法是世界粮农组织(FAO)推荐计算潜在凝欣出的方法.这里.定义潜在蒸散量为一种假想参照作物冠层的蒸散速率,假设作勒植株高度为0.12，固定的作物表面阻力为70ms,反射率为0.23,非'荒类似F衣面开阔、诲度一致、生长旺盛、完全遮盖地面而水分充分适宜的绿色草地的蒸放;丸FAOPeruKin-Mo1Heithf修正公式表达如下：408Ar尼-6J+y亭X。.-QPET=(C.5>+f1.+0.34w2)式中：PET潜在然散Sih单位为富米昨天(mud')；&一一地表冷福射，单位为兆焦每米每天(MJm'd')：G-土壤热通盘.单位为兆焦每平方米每天(町m?d')；Tman-H平均气讯,第位为摄氏度(C):U2一2米高处风速，单位为米每杪(Ins):饱和水汽压，单位为帕(kPa：J实际水汽压，单位为帕(kPa)；一一饱和水汽压曲规斜率，通位为千帕每排氏度(kPa-'O;/干湿表常数,单位为千帕每摄氏度(kPaX?).C.3.2FAOPganTonteith方法潜在蒸«AiF算步C.3.2.1计算日平均气温(7；M)I1.IjFP0Pen三nMonteith公式中温度资料的非线性分布,某时段的平均气温以该时段的日最高气温、日能低气部计算得来。月、季、年的日最育气混、日最低气温为月、季、年日最海气温、日最低气温的总和除以月、季、年的总11数得到，FAOPenman-Monteith公式中用到的11平均气温(，3)，建议由日最高气温(»“)和日最低气温(<n,n)的平均值计算知到，而不毡当日24小时逐时(或一日4次、8次观测气温的平均(ft.T+TJ=3yEC6)式中：Twm一一口平均气温，单位为摄氏度(C)：7；八11城海气温，单位为摄氏度(C)；Tma日最低气温,单位为摄氏度(C).C.3.2.2计算2米高处风速(U)在计。潜在蒸散时,衢要2米高处测址的风速.其他高度观测到的风速可以根据下式进行订正，Hi=it,<C.7>211.n(67.8z-5.42)式中：“22米诲处的风速，单位为米砥抄（ms'>n:一-:米高处测量的风速，单位为米每秒（ms'）：Z风速计仪器安放的离地面离度，单位为米（m）.C.3.2.3计算平均博和水汽压（e,）饱和水汽压与气温相关，计算公式如下：1727Te0（T）=0.6108×cxp-<C.8）7'+237.3式中：ei,（T）一一气温为7时的饱和水汽压，单位为千帕（kPa）：T空气温度，单位为摄氏度（£）.由于饱和水汽压方程（C8）的非战性，11,旬、月等时间段的平均饱和水汽压应当以该时段的日最庙气淑、日最低气温计算出来的饱和水汽压的平均1来计算：e_%（二N）+%（7iB）（Cg）,一2式中：e,平均饱和水汽压.电位为千帕（kPa）：en（TmJ一一为口鼓高气温7；八时的饱和水汽压，单位为帕（kPa）：。（，/一为日最低气温时的泡和水汽JA胞位为千帕（kPa）.如果用平均气温代替日最高气温和日最低气温会造成偏低估计饱和水汽压e1.的值，相应的饭和水汽压与实际水汽压的差M少.最终的潜在蒸散仪的计算结果也会减少.C.3.2.4计算实际水汽压（Q）实际水汽压e“就是露点温度工*”下的饱和水汽压，单位为千帕（kPa）。实际水汽压计算公式；%=备（J）=060SXeXP二：晶（C.10）式中：e“实际水汽压，单位为千柏（kPa）：T1.itv露点温度,单位为摄氏度（C）：e1（乙“）一一为露点温度Ae下的饱和水汽的位为千帕（kPa）.C.3.2.5计算愎和水汽压曲线斜率（）饱和水汽压与温度的斜率的计算公式:172774O98O.6IO8cxp()(C.11)A7+2,37.3(7+237.3)2式中：一一在气温为T时的饱和水汽压斜率，单位为帕姆摄氏度(kPaC');T一一空气温度，雉位为接氏度(C),C.3.2.6计算土集编JM<G)运用复杂模式可以计算土壤热通量.相对于冷幅射RC来说,土地热通JAG是很小的fit特别是当地表被植被滑盖、计算时间尺度是24小时或更长时，当计算较长的时间尺度时，下面的简化公式可以用来计算土塘热通法:G二式中：G一-土塘热通吊：，单位为兆焦每平方米旬天(MJ.jday,>sCf土堆热容依，单位为兆焦每立方米每度(MJ'C')：i晡刻i时的空气温度，维位为摄氏度<r)：ii时刻i-时的空气阻度，单位为摄氏度<r)：AI-时间步长.单位为天(day)：z为有效士耀深度.单位为米(11>).土壤热容量与土堰如成成分和水分含量有关.一天至十天的时间尺慢，斟考隼地的土壤热容业相当小，可以忽略不计;GNo(C.13)月时间尺度,转改在适当的土塔深度、土塘热容城为常数2.1MJm'C时，由方程(C1.4)可以用来估算月土地热通JftG:G=GVz=吟(iw11,一，M1.IbX)=OJ-(皿_。CM)JOJ式中：ZEW第"1时的平均气温，单位为摄氏度(C):Aewt一一上月平均气温单位为摄氏度(X?)：C.3.2.7计算干衰常数(/)干湿表常数y由下式计算得到：(C.15)=0.665X10'P式中:Y一一干湿表常数，单位为千柏埒摄氏度kPaC'；一一蒸发潜热，2.，15兆焦每克(JUk8,);J空气定压比热，1.013X】。'兆侔每千克每摄氏度(MJkg"O:£一一水与空气的分子度之比，取0.622：z当地的海拔拓度,单.位为米(m)：P一大气压.胞位为千帕(kPa).无观测值时.可由公式(C.16)计驾。Eo""(C.16)I293)C.3.2.8计算JHi射(/?.,)净辎射凡是收入的短波净辐射aU和支出的长波净辐射R2之差：R,a<ci7>式中：Ru一一净辎射，单位为兆您年平方米每天<MJm-day,);R1.U一太阳净幅射或短波净辐射,单位为兆烬每平方米每天(MJ.,):R1.t1.一一长波净辎射，单位为兆焦好平方米每天<MJmay1.>.计算地衰净辐射的步如下：第一步计*i"外射(R1).不同纬度一年中每日的地球外辐射R1.I可以由太阳常数、太阳磁偏角和这一天在一年中位置来fAih04×6()Ra-Gdr,sin(f7)sin()+cos>)cos()sin(<u1.)(C.18)11式中:R11一一地球外相射，单位为兆焦每平方米每天(MJmiday'：GK一一太阳常数,取值0.0820兆供每平方米每分钟(MJtnmin')：dr反转日地平均用离，出方程(C.19)计算；,日落时加，单位为孤度(rad，由方程<C.21>和(C.22)计尊：玮慢，单位为孤度(rad)；S太阳磁偏角,单位为弧度(rad),由方程(C.20)计算.日地平均距窗4和太阳磁偏角6由下式汁分：(1.r=I+0.033CoS巨7)(C.19)3652=0.408sin(一J-1.39)(C.20)365式中：J日序，取值范附为1到365或366,1月1日取日序为1日落时角叫由下式计算:t=arccos-tan()tan()<C.21)如果在所使用的计算机沿言中没有反余弦函数，日落时角叱也可以用反正切函数计算:K.r-tan(p)tiin(<)1“介八W,=-arctan1<c-22>式中：X=1.-tan()2tan(J)f(C.23)如果X0,X=O1(XXX)1.第二步计算可日惠时数(N),由下式计制：74N=一<,(C.24)r式中：N可11照时角：牝一一为方程(C,2I>或(C.22)计算的日落时痢.第三步计算太阳射(Ri)t如果没有太阳辐射R,的观测侪.可以由太阳相射与地球外辐射和相对日照的关系公式求得:Rt(at+-(C.25)N式中：R,太阳幅射或短波辐射,单位为兆焦每平方米每天(Nj三May')：w-实际H照时数,单位为小时(hour)：N一一最大可能日照时数，单位为小时(hour)：n/N一一相对11照：Ru一一地球外辎射，单位为兆焦好平方米i天<MJm1<tay')；a,一一表示阴天(=0)时到达地域我面的地球外辐射的透过系数：凡十一暗天(=N)时到达地球表面的地球外相射透过率.%和”的大气状况湿度、尘埃)和太阳磁偏角纬度和月份)而变化。当没有实际的太阳辐射资料和经验参数UJ以利用时，推荐使用a,=025,¼=0.50.第四步计算太阳净射或短波净制射(RIU)地表短波净辎射由接收和反射的太阳辐射的平衡来计算：R,u=(I-Ct)/?,<C.26)式中：凡,一一太阳净幅射或短波净辐射,单位为兆焦每平方米柘天(MJ三'):反照率，此处取绿色草地参考作物的反照率0.23：凡一一接收的太阳辎射，单位为兆焦PF方米饵大(MJ.2day,>.第五步计算空太阳辑射(/?“,)在接近海平面或者凡和"有经验参数可以利用时，踏空太阳辐射由下式计修；&,=(。,+")凡<C.27)式中：R11,一一暗空太阳播射.单位为亮焦每平方米每天<M.Jm2day,):4+2一晴天(A=N)时到达地球表面的地球外辐射透过率：R11一一地球外辐射，单位为兆便短平方米郎大(MJmsday,>.在没有。，和"的势验信可以利用时,以下式计算晴空太阳辐射；RS=(O75+2XK)TZ)&(C.28>式中：2为站点海拔高度，单位为米()，第六步计算长波削射(R1.1.1.)长波辐射与地我绝对沿度的4次事成比例关系,这种关系可以由斯蒂芬-波尔兹曼定律(Stefan-Bo1.tzmann1.aw)定圻发示.然而，由于大气的吸收和向下相射,地表的冷能量通Ift要少于用斯蒂芬-波尔兹亚定律计算出来的值.水汽、云、二氧化碳和尘埃都吸收和锋放长波辐射，在估算净支出轴射通量时应当知道它们的浓度.由于湿度和云俄的影响大，所以在使用斯蒂芬-波尔拉亚定律时估算长波轴射净支出通收时，用这两个因子进行订正，并假设其他的吸收体的浓度为常数：r。+7'RRn1.=53,i1(0.34-0.1AM)(1.35J-0.35)(C.29)24式中：RZ长波冷辐射,单位为亮焦每平方米每天(MJmda屋)：斯蒂芬-波尔兹晚常数，数值为4.903X10"(MJK,.jday'>：TZ天(24小时)中最i绝对温度，取位为开尔文<K>(K=C+273.16：Tmmx天(24小时)中最低绝对温度，单位为开尔文<K>(K=r+273,16):Q一实际水汽压，单位为千帕(kPa)；氏-太阳幅射.单位为兆焦每平方米每天(MJm；day)：K,.,晴空相射单位为亮保每平方米每大<MJrn-day')：RjR1.o一一相对短波辐射(W1.0：(0.34-0.14)一一空气湿度的订正攻，如果空气湿度增加,该项的值将变小：(135-J-0.35)云的订正项，如果云量埴加，兄将及少，该项的值也相应减少。R,、附录D(规范性附录)标准化降水指数的计算方法由于降水量的分布一般不是正态分布，而是一种偏态分布所以在进行降水分析和干旱监测、评估中，采用分布概率来描述降水fit的变化.标准化降水指标(简称SPD就是在计算出某时段内降水班的分布概率后，再进行正态标准化处理,最终用标准化降水累枳颓率分布来划分干旱等找.标准化降水指数(简称SP1.)的计算步骤为：a)假设某时段降水用为陵机变域x则其分布的概率密度函数为：八幻=A>°-°其中:/7>O./>0分别为尺度和形状参数./?和/可用极大蚁然估计方法求得：-1+I+4A3、4Ax(D.3)其中：InA=Inx-q2Inxi1-1(D.4)式中：-V1一一为降水量资料样本：X一-为降水的气候平均值.确定概率需度函数中的答数后，对于某一年的降水玳MV可求出随机变量X小于公事件的概率为：尸(X<X0)=,'"(x)dx(D.5)利用数值积分可以计算用(DJ)式代入(D5)式后的事件概率近似估计值.b)降水量:为0时的犷件概率由下式估计；F(x=0)=m/n(D.6)式中：w降水技为O的样本数:总样本数，c）对分布概率进行正态标准化处理，即将<D.5>,（D.6）式求得的概率值代入标准化正态分布函数,即：(D.7)对（D,7）式进行近似求解可得:(D.8)其中t=,E为（D5）式或（D.6）式求汨的概率；并当F>0.5时，F值取1.0尸.S=I：当尸"0.5时.S=-I.CQ=2.515517,C1=0.802853,<?2=0.010328.J1=1.432788,4=01X92691d、=().(X)13O8III（D.8）式求得的Z值就是此标准化降水指数SP1.附录E(规范性附录)标准化降水蒸散指数的计算方法E.1标准化降水蒸微指数JK理干旱不仅受到降水的影响,而且与蒸敢密切相关.2010年ViCenIC-Serran。采用降水与蒸散的差值构建了SPI-I指数，并采用了3个参数的1.og-1.ogistic概率分布函数来描述其变化，通过正态标准化处理，最终用标掂化降水与蒸敝差俏的黑枳频率分布来划分干旱等级.E.2标准化降水IMt指数计算步第一捧计算潜在然散(PET)0ViCeme-SeMw。推荐的是Thom1.hwai【e方法，该方法的优点是考虑了温度变化，能较好反映地表潜在蒸散。第二步用公式(E.I)计算逐月降水I1.t与潜在蒸散量的基tt1.：D.三P.-PET.(E.1)式中：D1一一降水吊马潜在泰放吊的差值：P1月降水量：PETi月潜在蒸故也.计算方法见CJ.第三步如同SPI方法，对D1数据序列进行正态化处理，计算每个数做对应的SPE1.指数“I1.1.于原始数据序列。中可能存在负值.所以SPEI指数采用了3个参数的1.og-1.ogistic概率分布”1.og-1.ogistic概率分布的累积函数为：上式中的参数a、。、分别采用线性电的方法拟合获得:(E.3)aW.-2犷"，rd+)/z?)r(i-)(E.4)<E.5)=w.-<z(1.+1)(1-)式中，为阶乘函数,W1,用门W,为原始数据序列。的概率加权矩,计算方法如卜：MYir'J"(E-6)-O35<E.7)匕-V式中,N为参与计算的月份数.然后对累积概率密度进行标准化：F=I-F()<E.8)当累枳概率<0.5时：w=J-2In(P)(E.9)2SPEJ=C“-Ctr+CW(E.IO)1+diw+dzw+dtw式中，<o=2.515517.c,=0.802853,C2=0.010328.rf,=1.432788.<6=0.189269,rff=O-OOI3O8-当尸>0.5时,产值取I-P：SPE1.=-(r-C即+?r)<E11)1+dw+dr+dsw附录F(规范性附录)帕默尔干旱指数的计算方法F.1帕默尔干旱指数原理帕默尔干旱指数PDSI(DJ1.aer匕OUghtSeVerity1.ndeX)是表征在一段时间内,该地区实际水分供应持续地少于当地气候适宜水分供应的水分亏缺情况。域本原咫是土壤水分平衡原理。该指数是她于月做资料来设计的，指数经标准化处理，指数值-股在-6(干)和-6(湿之间变化，可以对不同地区、不同时间的土壤水分状况诳行比较。PKI在计算水分收支平衡时，考虑了前期降水麻和水分供儒，物埋意义明晰.F.2帕默尔干度指数计算方法R2.1第一步计*水分异常指敷Z水分供需达到气候适应的水平衡方程友示如下：AAAAAP=ET+RtRO-1.<F.1)式中：P气候适宜降水%ET一一气候适宜蒸散量：R气候适宜补水室：Z"气候适宜失水量：RO-气候适宜径流Jft.上述气候适宜值分别由下列方程计匏：ETaPE(F.2)R=PR(F.3)RO=yPRO(F.4)1.=P1.<R5>式中：PET潜在蒸侬依由FAoPennIan-MOnteith或ThOrnt.huaite方法计究.计算方法见附录C：PR一一土堞可能水分供给量，计算方程:(F.6)PR=AWC-(Si+Su)PRO一一可能径流，计算方程:PRO=AWC-PR=S1+Sn(F.7)PI.土城可能水分损失量.计算方程:P1.=P1.s+P1.v(F.8)P1.s=11in(PE,SJ；即：PE和S,两者选小的(F.9)P1.c=(PE-P1.s)SnAWC(F.10)式F.6-F.10中：AIC整层土壤H1.间有效持水埴:S,初始上层土壤有效含水量:SII一初始下层土壤疗效含水最。a、fi、7、b分别为蒸散系数、土填水供给系数、径流系数和土堞损失系数，等站年月分别有四个鸟的常系数值，计蚱如Ra=叵<FU>(3ueo)_(F,13)(PfiO)6=色<F.H)国)各段上面的横浅代去其多年平均值。ET一一平均实际蒸散信：7½7平均潜在然Ifiift：R平均土壤实际水分供给最；Tr一一平均土墙可能水分供给;匕RO一一平均实际径流埴；PRO平均可能径流量；1.一一平均实际上块水分损失量；P1.一一平均可能上堞水分损失I限PaImCr假定土填为上下两层模式，当上层土堞中的水分全部丧失，下层土壤才开始失去水分，而且下层土塘的水分不可能全部失去,在计算然放量、径流量、土堆水分交换城的可能伯与实际值时，需要遵循系列的规则和假定。另外，土墙有效持水汆AWCSvai1.ab1.eWaterHo1.dingCapacity)也作为初始蛤入f1.t.在计算PDSI过程中,实际伯与正常值相比的水分距平d表示为实际降水量(P)与气候适宜降水肽(户)的基：(I=P-P(F.15)式中：d一一水分距平：P实际降水累：P气候适宜降水玳，为了使PDSI成为一个标准化的指数,水分距平d求出后，又拘其与指定地点给定月份的气候权重系数K相乘,得出水分异常指数Z,也称Pa1.三erZ指数.表示给定地点给定月份.实际气候干湿状况与其多年平均水分状态的偏禹程度.Z=(IK(F.16)式中：Z水分异常指数：d水分矩平：K一气候权重系数.气候权重系数K的取值也月份和地理位置决定,由公式(V.17)计算.F.2.2第二步计算修正的气候特征系数K公式(F.16)的气候特征系数K，根据中国气候特点进行修正得：y12八，j,三西+瓦+而-÷X.o其中：<=1.6-Iog10+0.4(F.18)式中：K一一气快特征系数说权曳因子：f-表示第i个月.i=1.2,12：DfK.一一多年平均年绝对水分异出，j去示1至12月：/-1Di一一第i月的水分距平d的绝时值的多年平均值；PETi第i月的平均潜在蒸侬fit：R一一第i月的平均土塔实际水分供绐ht：ROt第j月的平均实际径流中；Pi一一第i月的平均实际降水早:：1.1一一第i月的平均实际土壤水分损失鼠：F.2.3第三步建立修正帕默尔干旱指数根据帕默尔早度模式的思路.利川我国气象站资料对帕默尔早度模式,进行修正得:PDS1.i=0.755PoS/,+Z1(F.19)''1.63'式中：PDSI1第i月PI)S1.干下指数值：Z1.，月水分异常指数：PDS1.1.第i-1月的PDSI干旱指数(ft.公式(F.中的0.755和一为持续因子，起始月份的PDW,的计算公式为1.63PDS1.i=Z.(F.20)1.63附录G(规范性附录)标准化权重降水指蛾的讨算方法标准化权受降水指数(SP1.W)首先对某一时段内逐日降水量进行加权累积，然后对权更累积的降水量(M'AP.WeightedAverageOrPreCiPiIa1.ion)进行标准化处理而得到的指数，标准化处理方法见附录D.标准化降水权里指数的计算步骤为：a)第一步计算权施累积降水盘,按式<G.1)计算:=£a11Pn(G.1)»-0式中；K伊权重累枳降水量，单位毫米(mm)：,V-某一时段的长度，单位天：«贡献多数，当N为60天时，则a取0.85：A距离当天前第n天的降水fit单位亳米(11m>.b)第二步计算标准化权正降水指数，按式(G.2)计算：Sm=SPIW)(G.2)式中：SPIW-标准化权重降水指数；SPI一一标准化处理,计算方法见附录C：VAP-一权求累枳降水指数，单位耻米(伽)，计律见式(GJ，附录H（资料性附录）气象干旱综合指散ZC1.中单项指侬的权重系数取值通过计算和比较,确定了各省（区、市）ZCI各参数的推荐取值（表H.D.各省（区、市）n取值为150天或120大，a、b和C的取值为0.20.8之间。表H.1各省（区、市）气象干旱综合指数ZC1.的参数推荐取值省（区、市）ZCI的套致取值nabC北京1500.20.20.8天津1200.20.20.8河北1500.20.20.8山西1500.20.80.2内蒙古1500.20.20.8辽宁1500.20.20.8Ji林1500.30.30.6黑龙江1500.20.20.8上海1200.20.20.8江苏1500.20.80.2浙江1500.80.20.2安徽1500.20.20.8福建1500.20.20.8江西ISO0.30.40.5山东1500.20.20.8河南1500.20.20.8湖北1500.30.40.5湖南1500.30.40.5广东1200.30.40.5广西1500.30.10.5海南重庆四川贵州云雨西藏陕西甘京出海宁更新源0.20.20.80.30.40.50.20.80.20.20.20.80.30.40.50.20.80.20.20.80.20.20.20.80.30.-10.50.20.20.81501501501501201501501501501501500.20.20.8附录I（规范性附录）使用参效法计算短序列气象观测站气象干旱综合指数的方法使用参数法计算短序列气观测站气象干公综合指数的方法和步骤如下.1.1 短序列站及其周边多源数据的收集及预处理。包括短序列孙（拟算站）及周边具有50年以上观测资料长序列站的逐日降水资料的数据收集及预处理，具体包括以下步骤。（I）首先.收集拟獴站及周边的多派数据,所述多源数据至少包括：拟算站的经度、端度、海拔i度及站间IH地形、地貌和下垫面状况资料，所需统计气象干旱缘合指数时段及前期近36S天以上逐日降水/、日平均气温资料.（2）选取拟能站周边100km半径内具干j30年以上长序列降水资料站点的经度、纬度、海拔高度及站周用地形、地貌和卜垫面状况资料，以及近30年以上逐日降水做资料，且用列越长越好.（3）针对所收集的拟驾站和周边长序列气象观测数据进行数据的质型控制，包括气候阈值及时空一致性检查的方式.引除原始气飘观测数据中的假值和虚假值.对于连续多个时次的观测值固定不变的做他,通过“a-.=。进行筛选和剔除，其中，M,.表示时刻的观测值，表示I时刻的观测值,1.2 选取拟豫站周边长序列站，计翼其与短序列站的距席和海拔高度差.选取与拟算站降水概率分布相似的周边站点.基于I中值处理后的气&观测数据.对拟算站周边100km半径内选取出的长序列站进行进一步暗选，选取与拟算站近一年降水和气温进行相关性发好的3S个站点，同时保证其相关显著性检脸均超过9S%信度，然后再注曾选取的长序列站与拟算站的距离和海拔高度差.长序列站与拟算站之间距离的具体计算过程为：基于1.1.中模处理后的自站经度、纬度.及1.2中挑选得到的周边长序列站A,对应台站经度、纬度假设为1.onA,1.aiA;拟算站设为M,对应台站羟度、纬度假设为1.onM,1.xitM;则M与A之间的用禹DMA计算公式如下：MfA=R*arccos(sin(Zarf)*sinZrJ)+COS(IafH)*cos(Z<)CoS(1.omI-1.onff1*/180(1.(1)式中:R为地球平均半径，取6371km：11=3.14s1.onA,1.a1.A和1.onM,1.atM分别为A站、M站的经度和绊度.1.3 统计长序列站多时间尺度降水分布概率密度函数的特征参数收.荔于I中预处理后的气象观涮数据，统计计算由1.2选取得到的周边长序站的逐11演动60天权里降水麻、90天降水盘、150天降水盘尺度的分布概率密侬函数的尺度参数.形状参数丫和零降水欢概率FO特征卷数状，并构建其特征多数域的标准气候伯.不同时间尺度降水量的分布概率密度函数的特征参数量的计算步骤见附录D.权重降水量的分布概率密度函数的特征参数值的计算方法见附录G.1.4 利用插依方法计算拟算站降水厂分布概率密度函数的特征参数Jih¼1.2中得到的拟算站与周边长序列降水站点的距施，将1.3构建的多尺度的分布概率密度函数的特征冬数量的标准气候值,通过反距离加权插值法得到拟算站的多尺度r分布柢率密度函数的特征参数!ftB、y和H),从而建立拟算站的6()大权重降水出、90天降水依、150天降水玳的多时间尺度厚水求分布概率密度函数。反距起加权插伯法的计算公式为：z(.j/)=yj1./y-1.d.2>式中：Z(W)一一拟算站的特征多数出B、丫和F0：Zj-周边第/站的特征参数量B、Y和F0：明一周边第1站与拟笄站之间的距两：选取得到的站数。1.5 计算拟算站多时间尺度标准化降水指数.利用1.1中处理后好的拟算站的降水量资料，分别计算出近60天权重降水量、90天降水业、150天降水量,代入1.4得到的拟算站相应尺度降水法的分布概率密度函数，进而分别求得近60天标准化权IR降水指数Sw%,近90天标准化降水指数SP90,近150天标准化降水指数.%7闾.多时间尺度的SW的计律方法见附录D.标准化权重降水指数9小的计算方法见附录G1.6 计算拟算咕近30年热欣依和相对湿润度指数，利用1.1预处理后好的“拟算站”的逐H降水盘和平均气温资料,计算出近30天得到相对湿润度指数J：,。.相对湿润度指数Mfia的计算方法见附录B.1.7 计算拟算站气象F早综合指数.利用I.5得到的“拟制站”近60天标准化权电降水指数SV%,近90天标准化降水Iifft577,0,近150天标准化降水指数S/m以及1.6得到的拟算站的相对湿制度指数if,n,就可以计算出拟我站气象干？综合指数MC1.和ZCI.