搜索
1,1,1,2(1,710048;2,100050):,Matlab,Viney1981~1990,,.:,,,,;30%,,.,,.:;Viney;;;:X131.2:A:02503301(2006)07129207:20051122;:20060310:(05JS35);(2001282):(1960~),,,,,Email:Lhuaien@mail.xaut.edu.cnContinuouslyModelingResearchontheTransferandTransformProcessesofPhosphorusinHeiheRiverWatershedLIHuaien1,XUYanling1,ZHANGQiang1,NIYongming2(1.InstituteofEnvironmentScienceofXianUniversityofTechnology,Xian710048,China;2.BeijingNaturalHistoryMuseum,Beijing100050,China)Abstract:TakingsupplywatersourceofXiancityasanexample,bythehelpofMatlabsoftware,usingfillandyieldmodel,theInverseGaussflowconcentrationmodel,runoffsedimentrelationmodelandVineymodel,thetransferandtransformprocessesofphosphorusfrom1981to1990weremodeledcontinuously.PhosphoruspollutionloadsofHeiheRiverwerealsoestimated.Theresultshowsthatthecontinuousmodelingresultaccordedwiththegenericlawofphosphorusloss.Therelativeerrorbetweenmodelingresultandmonitorinterpolationresultwasnotmorethan30%.SothemodelingmethodproposedinthispapercouldbeusedtomodelcontinuouslythetransferandtransformprocessesofphosphorusinHeiheriverwatershed.Thisstudywasonlyprimaryattempttononpointsourcepollutioncontinousmodel,therefore,themodelingmethodneededtobeimprovedandperfectedfurther.Keywords:phosphorus;Vineymodel;continuouslymodeling;HeiheRiver;watershed,,.,.[1~4],.,1481km2,45km2,,,,.,,Viney1981~1990,,,,.1,.,Matlab,27720067ENVIRONMENTALSCIENCEVol.27,No.7Jul.,2006,.,,1981~1990.,8,3h.22.1,,,,,()[5~8].2.2,.,.2.2.1[9~11],(InverseGaussianDistribution)Tweedie[12].h(t),Q(t)h(t)f(t):Q(t)=t0f(t-)h()d(1),,.,[9~11].2.2.2(),,,.,.(1)QG(t)=QG(t-1)KKG+RG(t)(1-KKG)U(2),QG;KKG;RG;U.(2)QSS(t)=QSS(t-1)KKSS+RSS(t)(1-KKSS)U(3),QSS;KKSS;RSS.,,5mm,..t,Rt,:M=24tN=24tintR5+1(4)(KKSSD)(KKSSD)(KKSS)(KKG):KKSSD=KKSS1M(5)KKGD=KKG1M(6)(KSSD)(KGD)(KSS)(KG):KSSD=1-[1-(KG+KSS)]1N1+KGKSS(7)KGD=KSSDKGKSS(8),,,,().19880303~031323m3/s10a,23m3/s.,(2)KKG1h,..,KKSS.,2.1,K,,1~4,,12937,,KSSKGKKSS:Min[FUN(KG,KSS,KKSS)]=Min[ER]=Min|Ri-ri|2,KSSKGKKSS,3.198919902a,K.:Min[FUN(K)]=Min[ER]=Min|Ri-ri|,Ri;ri.,,.3,.,,.[13],.:,1981~19906~9,,1(a).1,3,(168,1985)(242,145)(234,128).,,3,1(b).:S=0.000001Q2+0.0022Q-0.0222(9),S(kg/m3),Q(m3/s),R2=06114.(9),,..113,,0,4~100,,(9)2:1,1130,4~10;2,102m3/s(9),,,0.110aFig.1RelationbetweenmonitorrunoffandsedimentofHeiyukoustation4Viney[14],2.4.1(1)UP=kUPddt(10),UP;kUP,;ddt.ddt,ddt0,,.A()U(),PIAPIUPB.(2)1294272Fig.2StructuralmapofphosphoruscirculationHP=kHPPB(11),HP;kHP;PB.(3);,21:MP=KMP1r+1(PO-rPI)(12),MP;r;KMP,1.2,.4.2(1),,cP.,,:cPg=cPp/pg(13),cPg,;p;pg.,qS,cPgqS.(2)A,,,:SP=kSPqSPI(14),SPA;kSP;qS;PI.(3)UA(U)(A),:TP=kTPqssPIU(15),TPUA;kTP;qss;PIUU.(4)A,A,A,:AP=PIAqA/A(16),AP;PIAA;qA;AA.,:API=cPgqS+SP+AP(17),API.m3/s,1000t,g/s.4.3(1)Viney,1.XP=kXP/YXPE(18),XP,1XP10;kXP;YE;XP,0XP1.XP,:EP=kEPXPYE(19),EP;kEP.(2),,,12957:RP=POC(YR/SC))RP(20),RP;POC;YR;SC;RP1.(3),:WP={[Z/(YE+YI)]DP-1}(EP+IP)(21),Z;WP;YE;YI;EP;IP;DP.(4),.:BP=kBPPOAYB(22),BP;kBP;POAA;YB.(5):DP=(1-DP)(EP+IP+BP-WP)(23),DP;;DP.,:QPO=EP+IP+BP-WP-DP(24):OC=WP+DP-RP(25)4.42,,,,,.(12)()().Viney,KMP1,r.,3h01~05kg/s,r49.POPI,.,.1981,1995,1.9408%,224%()[15],6mg/kg.Viney5cm,(5cm),103t.10a,,,10a.(18)(19)Viney,.,.kXPXPkEP.Viney0XP1,1XP10,,1(3h)10~100t,(18)kXP=100XP=03,kEP.1Table1MonitorresultofsoiltotalphosphorusinHeiheRiverwatershed/mgkg-1/mgkg-1165778982113773671421417275079309261041,6,KMPkSPkTPkXPXPkEP.KMP=1kXP=100XP=03,kSPkTPkEP.(1981~1990),199519966,199519966(Pi=0.6198Q-17.396,R=09997;Q,104m3),1981~1990,,129627Pi,,pi.,1,,.1981~1983,1984~1990.:Min[FUN(kSPkTPkEP)]=Min[ER]=Min[|Pi-pi|](26),.,.2.2Table2PreferencesresultofnutrientmodelkSPkTPkEP060131645,1981~1990,199519966,1981~1990,3.1981~19903.3Table3Totalphosphorusloadsofmodelcalculationandmonitorinterpolation1981198219831984198519861987198819891990108/m312581122883521511701646631669/t744360756547323316434400391415/t89627599625269305415134256/%170-309236125166-38-592206926031981~1990Fig.3Calculationresultofnutrimentmodelfrom1981to19903,30%,1982198319881990;19811983896t99t,1983,;19821985198610a,275t269t305t;,6~9108m3,40~65t.,:,,,.31981,1986,4,1988,5.41986Fig.4Calculationresultofnutrimentmodelin198651988Fig.5Calculationresultofnutrimentmodelin1988129776,,,Viney.,,.,.,,,.:[1]OECD.Eutrophicationofwaters:monitoring,assessmentandcontrol[R].Paris:OECD,1982.154.[2]FreedmanB.Theimpactsofpollutionandotherstressesonecosystemstructureandfunction[M].San.Diego:AcademicPressInc,1989.[3]SchindlerDW.Eutrophicationandrecoveryinexperimentallakes:implicationoflakemanagement[J].Science,1974,184:897~899.[4]SchindlerDW.Theevaluationofphosphoruslimitationinlakes[J].Science,1977,195:260~262.[5],,,.[J].(),2003,36(5):38~42.[6].[M].:,1984.106~130.[7],.()[J].,2000,14(1):7~10.[8],,.[J].,2002,21(3):32~35.[9],,.[J].,1994,(2):147~151.[10],.[M].:,1996.43~52.[11],.[J].,1997,17(2):141~147.[12]TweedieMCK.StatisticalPropertiesofInverseGaussianDistributions:,[J].Ann.Math.Statist,1957.28.[13].[R].:,1998.[14]VineyNR,SivapalanM,DeeleyD.Aconceptualmodelofnutrientmobilizationandtransportapplicableatlargecatchmentscales[J].JournalofHydrology,2000,240: