（计算机应用技术专业论文）基于数据挖掘技术的黄土高原降雨侵蚀预测模型的研究.pdf

摘要 降雨是引起土壤侵蚀的主要原因之一，为了减少降雨侵蚀造成的危害，构建降雨侵 蚀预测模型已经成为我国水土保持工作的一个研究重点。降雨侵蚀预测模型是土壤侵蚀 研究的重要部分，可以作为预测降雨引起的土壤侵蚀量的主要工具，为定量研究土壤侵 蚀提供依据。为了能够使模型具有实际应用价值，采用黄土高原( 黄河流域) 实际测量 得到的数据作为分析基础，使用数学建模方法得到降雨侵蚀预测模型的方程表达式。 降雨侵蚀预测模型属于土壤侵蚀模型的组成部分，在土壤侵蚀的建模理论中，因子 分析模型提供的方法为本文的研究提供了分析依据。因此，本文利用降雨侵蚀因子和土 壤侵蚀量的测量数据来构建降雨侵蚀预测模型。 在理解降雨侵蚀机理和土壤侵蚀相关理论的基础上，本文主要是以数据分析为手 段，以实际应用价值为目的构建降雨侵蚀预测模型。具体的，本文使用回归分析作为数 学建模的主要方法，并结合数据挖掘中的数据预处理技术对数据进行处理，从而达到优 化模型的效果。本文主要使用s p s s l 3 0 作为数据分析软件，应用环境是w i n d o w sx p 。 主要研究成果为： ( 1 ) 通过对降雨侵蚀数据的分析，结合数据预处理技术和数据的实际含义，提出 降雨侵蚀预测模型的数据预处理算法，结果表明，使用该预处理算法处理后的数据所建 立的模型要比原始数据建立的模型更好。 ( 2 ) 通过对降雨侵蚀数据的曲线拟合分析，结合对通用水土流失方程( u s l e ) 的 理解，提出降雨侵蚀预测模型的乘法模型。 ( 3 ) 通过乘法模型的线性化变换方法，将乘法模型转化为可以使用线性回归分析 进行求解的表达形式。将结果求解、变换，得到降雨侵蚀预测的乘法模型。 ( 4 ) 所使用的模型评价指标都是数值型，在模型的比较和评价中，这些都可以定 量地对结果进行评价。另外，还使用图形对模型的结果作了分析和比较，在模型求解中， 这种方法可以更直观地反映问题。 本文在最后总结了在研究过程中遇到的一些问题和解决方法，并对降雨侵蚀预测模 型的建立和应用作了展望。 关键词：土壤侵蚀，回归分析，数据挖掘，曲线拟合，数学模型 a b s t r a c t r a i n f a l li st h em a i nc a u s a t i o nf o rs o i le r o s i o n t om i t i g a t et h ej e o p a r d yf r o mr a i n f a l l e r o s i o n ，o u rg o v e r n m e n th a sl a i dt h ec r u c i a lp i v o to ft h en a t i o n sw a t e ra n ds o i lp r e s e r v a t i o n r e s e a r c ho nc o n s t r u c t i n gr a i n f a l le r o s i o np r e d i c t i o nm o d e l ，w h i c hi sai m p o r t a n tp a r to f r e s e a r c ho ns o i le r o s i o n ，a n di tc a l le s t a b l i s ht h ef o u n d a t i o nf o rr e s e a r c h i n gs o i le r o s i o ni n f i x e dq u a n t i t y , a st h em a i ni n s t r u m e n tt of o r e c a s tt h es o i le r o s i o nq u a n t u mc a u s e db yr a i n f a l l t h ef a c t u a lm e a s u r i n gd a t aa b o u tt h el o e s sp l a t e a u ( t h ey e l l o wr i v e rd r a i n a g ea r e a ) i s i n t r o d u c e da sa n a l y z i n ge l e m e n tt oe n d o wt h em o d e lw i t hp r a c t i c a la p p l i c a t i o nv a l u e ，a n d m a t h e m a t i c a lm o d e l i n gm e t h o di su s e dt of i g u r eo u tt h ee q u a t i o n e x p r e s s i o no fr a i n f a l l e r o s i o np r e d i c t i o nm o d e l t h er a i n f a l le r o s i o np r e d i c t i o nm o d e li sac o n t r i b u t i n gp a r to fs o i le r o s i o nm o d e l a c c o r d i n gt ot h et h e o r yo fm o d e l i n gt h es o i le r o s i o n ，t h er e s e a r c ho ft h i sd i s s e r t a t i o na c h i e v e s a n a l y z i n gb a s i su s i n gt h em e t h o dp r o v i d e db yg e n ea n a l y z i n gm o d e l t h u si nt h i sd i s s e r t a t i o n ， t h er a i n f a l le r o s i o np r e d i c t i o nm o d e li sc o n s t r u c t e di nt h el i g h to fr a i n f a l le r o s i o ng e n ea n dt h e m e a s u r i n gd a t ag a i n e df r o ms o i le r o s i o nq u a n t u m o nt h eb a s i so fu n d e r s t a n d i n gt h em e c h a n i s mo fr a i n f a l le r o s i o na n dr e l a t i v et h e o r y a b o u ts o i le r o s i o n ，t h i sd i s s e r t a t i o nc o n s t r u c t st h er a i n f a l le r o s i o np r e d i c t i o nm o d e lw i t hd a t a a n a l y z i n ga st h ei n s t r u m e n tt oa c h i e v et h e a i mo fp r a c t i c a la p p l i c a t i o n c o n c r e t e l y , t h e r e g r e s s i o na n a l y s i si si n t r o d u c e da st h em a i nm a t h e m a t i c a lm o d e l i n gm e t h o d ，c o m b i n e dw i t h t h ed a t ap r e t r e a t m e n tt e c h n i q u eo fd a t am i n i n g t h u st h ed a t ai sp r o c e e d e dt oa c c o m p l i s h e f f e c to fo p t i m i z i n gt h em o d e l s p s s1 3 0i ss e r v e da st h em a i nd a t aa n a l y z i n gs o f t w a r ei nt h e w i n d o w sx pa p p l i c a t i o ne n v i r o n m e n t t h em a i na c h i e v e m e n t so ft h i sd i s s e r t a t i o ni sa s f o l l o w s ： ( 1 ) b ya n a l y z i n gt h er a i n f a l le r o s i o nd a t a ，c o m b i n e dw i t hd a t ap r e t r e a t m e n tt e c h n i q u e a n da c t u a lm e a n i n g so ft h ed a t a ，t h ed a t ap r e t r e a t m e n ta l g o r i t h mi sp r o p o s e df o rr a i n f a l l e r o s i o np r e d i c t i o nm o d e l t h ea n a l y z i n gr e s u l t si n d i c a t et h a tt h em o d e lc o n s t r u c t e dw i t hd a t a p r o c e e d e db yt h i sa l g o r i t h mb e h a v e sb e t t e rt h a nt h eo n ec o n s t r u c t e dw i t ho r i g i n a ld a t a ( 2 ) t h em u l t i p l i c a t i v em o d e lf o rr a i n f a l le r o s i o np r e d i c t i o nm o d e li sp u tf o r w a r db y l u a n a l y z i n gt h er a i n f a l le r o s i o nd a t aw i t hc u r v ef i t t i n gt e c h n i q u ea n dc o m b i n i n gc o m p r e h e n s i o n o fu n i v e r s a ls o i ll o s se q u a t i o n ( u s l e ) ( 3 ) t h em u l t i p l i c a t i v em o d e l ，w i t ht h em e t h o do fl i n e a rt r a n s f o r m a t i o nf o ri t ， i s t r a n s l a t e di n t oe x p r e s s i o nf o r mw h i c hc a nb ef i g u r e do u tb yl i n e a rr e g r e s s i o na n a l y s i s t h e r e s o l v e dr e s u l ti st r a n s f o r m e dt oa c h i e v er a i n f a l le r o s i o np r e d i c t i o nm u l t i p l i c a t i v em o d e l ( 4 ) t h em o d e le v a l u a t i o nc r i t e r i o n su s e di nt h i sd i s s e r t a t i o na r ea l li nt h ef o r mo f n u m e r i c a lv a l u e s ot h e s ec r i t e r i o n sc a nb ea p p l i e df o rq u a n t i f i c a t i o n a le v a l u a t i o no ft h e r e s u l t si nc o m p a r i n ga n de v a l u a t i n gt h em o d e l i na d d i t i o n ，g r a p h i c sa r ee m p l o y e dt oa n a l y z e a n dc o m p a r er e s u l t so ft h em o d e l t h i sm e t h o dh e l p st of i n do u tp r o b l e m si nm o d e ls o l u t i o n a tl a s tt h i sd i s s e r t a t i o ns u m m a r i z e ds o m eq u e s t i o n sa n dm e t h o d si nt h ep r o c e s so f r e s e a r c ha n di n d i c a t e st h ef o r e g r o u n do fc o n s t r u c t i n ga n da p p l y i n gt h er a i n f a l le r o s i o n p r e d i c t i o nm o d e l k e y w o r d s ：s o i le r o s i o n ，r e g r e s s i o na n a l y s i s ，d a t am i n i n g , c u r v ef i t t i n g ，m a t h e m a t i c a l m o d e l 论文独创性声明 本人声明：本人所呈交的学位论文是在导师的指导下，独立进行研究工 作所取得的成果。除论文中已经注明引用的内容外，对论文的研究做出重 要贡献的个人和集体，均已在文中以明确方式标明。本论文中不包含任何 未加明确注明的其他个人或集体已经公开发表的成果。 本声明的法律责任由本人承担。 论文作者签名：乌交木自 功口- 7 年莎月占日 论文知识产权权属声明 本人在导师指导下所完成的论文及相关的职务作品，知识产权归属学 校。学校享有以任何方式发表、复制、公开阅览、借阅以及申请专利等权 利。本人离校后发表或使用学位论文或与该论文直接相关的学术论文或成 果时，署名单位仍然为长安大学。 ( 保密的论文在解密后应遵守此规定) 论文作者签名： 导师签名： 乌天桶 散嘲 2 口一年6 只ge l 匕7 年月寥e l 长安大学硕士学位论文 第一章绪论 水土流失治理关乎国计民生，造成水土流失主要原因是土壤侵蚀，土壤侵蚀的定 义是土壤及其母质在水力、风力、冻融、重力等外营力作用下，被破坏、剥蚀、搬运和 沉积的过程。在自然界，当携带和推移松散物质的力超过抵抗松散物质移动的力时，便 引起了土壤侵蚀。一般地，水力作用造成的土壤侵蚀称为水蚀，也就是水土流失。 由降雨产生的土壤侵蚀称为降雨侵蚀，在土壤侵蚀研究中，降雨侵蚀相对于其他 侵蚀而言，可以由测量得到较为可靠的数据，因此可以通过数学建模方法构建数学预测 模型。 降雨侵蚀预测模型在水土资源管理中有着重要作用，它是预报水土流失、指导水 土保持措施配置、优化水土资源利用的有效工具。同时，降雨侵蚀预测模型是土壤侵蚀 数学模型的重要组成部分，随着科学技术的不断进步和发展，建立模型的方法也层出不 穷，不断更新。 1 1 课题研究意义 1 1 1 背景介绍 水土流失是影响社会进步、经济发展和资源开发的重点问题之一，属于环境问题研 究范畴。防止水土流失从古至今都是人类社会发展一大研究领域，由于人类社会发展速 度不断提高，造成的土地资源破坏也日趋严重，水土流失问题对农业、交通、人类生存 环境等都构成了威胁，特别是在发展中国家，由于在早期对水土流失的危害缺乏认识， 已经构成了水土流失的恶性循环，在我国一些地区，如黄土高原黄河流域，水土流失尤 其严重。要建立可持续发展的社会，水土流失成为首要解决的问题之一。 在科学研究中，为了解决水土流失造成的问题，必须从土壤侵蚀的研究入手。研究 土壤侵蚀的成因、分析土壤侵蚀规律是防止水土流失的根本手段。在土壤侵蚀规律的研 究中，建立土壤侵蚀预报模型是现今研究的热点。构建土壤侵蚀预报模型主要是研究土 壤侵蚀过程中各种侵蚀因子在侵蚀过程中的作用。土壤侵蚀预报模型是解决水土流失危 害的有效工具，建立实用性强、应用范围广、结构简单的土壤侵蚀预报模型一直是研究 者们致力追求的目标。 降雨侵蚀是土壤侵蚀的主要侵蚀过程之一，在土壤侵蚀的众多因素中，降雨是引起 水土流失的直接原因之一。在多种土壤侵蚀因素综合作用下，土壤结构遭到破坏，而造 1 第一章绪论 成土壤侵蚀过程的一个环节，搬运过程是必不可少的，大多数情况下，降雨引起的土壤 搬运比其他的因素要大得多，所以对土壤侵蚀的研究大多是围绕降雨侵蚀进行展开的。 另外，在测量方面，相关的降雨数据的收集比较方便，也更可靠，使得构造出的降雨侵 蚀预测模型可信度和实用价值更高。 1 1 2 土壤侵蚀的危害 土壤侵蚀及其导致的土地退化和泥沙等是全球性的环境灾害问题，并且已经对人类 的生存和发展构成了威胁。土壤侵蚀的危害直接反映在对土壤的破坏上。土壤侵蚀引起 的水土流失带走大量土壤、养分和水分，破坏土地资源，危害生态环境，后果严重，影 响巨大。其危害表现在对当地粮食生产的威胁和对下游生态环境的破坏具体为： ( 1 ) 破坏土地资源、蚕食农田； ( 2 ) 削弱土地肥力、威胁粮食生产； ( 3 ) 水土流失与地区贫困恶性循环； ( 4 ) 泥沙淤积，影响水库湖泊及航运综合利用； ( 5 ) 影响现代城镇的发展； ( 6 ) 水土流失会加快温室气体的排放。 中国是世界上的人口大国、农业大国，也是世界上土壤侵蚀最严重的国家之一，各 种类型的侵蚀遍及全国，而且强度高，成因复杂，危害极大，特别是黄河流域的土壤侵 蚀及其造成的水土流失尤为严重。据水利部遥感中心1 9 9 0 年全国范围的土壤侵蚀调查， 全国轻度以上侵蚀的水蚀面积达1 7 9 万k m 2 、风蚀面积达到1 8 8 万k m 2 、融冻侵蚀面积 达到1 2 5 万k m 2 ，总计轻度以上侵蚀面积达到了4 9 2 万k m 2 、合5 1 的国土面积需要治 理。大部分水蚀区和农牧过渡带的风蚀区急需治理。 由于水土流失的危害和治理的需要，水土流失预报日益成为许多国家政府环境建设 和治理规划与决策的重要依据，也是水土保持学科领域国际发展的热点问题，因此开展 中国水土流失宏观动态预报研究具有重要的生产实践意义和科学研究意义【。 在黄土高原地区，由于植被破坏严重，加上气候、地质、地形和人类活动等多种因 素的影响，水土流失问题一直是国家重点治理的对象。黄土高原地表支离破碎，千沟万 壑，有7 0 是坡地，植被覆盖极差，每遇暴雨，泥沙与雨水齐下，水土流失量之大为世 界所罕见。现黄土高原的水土流失面积占高原总面积的9 0 ，其中水土流失严重的地区 占一半以上。据山西省大宁县县志记载，太德塬在清光绪年间，塬面面积约8 7 0 公顷， 2 长安大学硕士学位论文 现在只剩下了6 0 0 公顷，其余的都变成了沟壑。黄土高原的水土保持是一项十分艰巨的 工作，需要较长的时间才能完成。当前，水土流失仍然是我国国民经济和社会发展的严 重障碍。据遥感最新普查结果，现黄土高原每年流失掉的土壤就相当3 6 万公顷土地的 3 0 厘米耕作层，伴随失去的氮、磷、钾养分4 2 0 0 万吨，相当于1 9 8 9 年我国化肥使用量 的1 7 7 1 2 1 。 1 1 3 建立降雨侵蚀预测模型的意义 土壤侵蚀定量研究中主要是对土壤侵蚀量进行分析，土壤侵蚀量是流域内被侵蚀土 壤的数量表示。我们可以通过野外调查获得土壤侵蚀量，也可通过数学模型计算求得。 但是通过野外调查来测量土壤侵蚀过程和评价土壤侵蚀的影响是一项困难且需耗费大 量人力物力的工作，这是因为： ( 1 ) 在整个流域内土壤侵蚀的来源是巨大的，且流域范围小至小块农田，大至大 型流域集水区，因而要定量评价土壤侵蚀程度和数量是困难的； ( 2 ) 土壤侵蚀是多因素综合作用的结果，包括各种自然因素和社会经济人为因素， 难以全面考察土壤侵蚀的影响因素； ( 3 ) 从一个地方得出的土壤侵蚀量关系难以外推到其他无测站的流域，开发和评 价所有潜在侵蚀控制的方法是不可行的。 自从2 0 世纪6 0 7 0 年代，数学模型已在水文学和流域管理中广泛应用。同样在土 壤侵蚀研究中运用数学模型模拟流域内的土壤侵蚀过程是一项有意义的工作。同时计算 机技术也不断提高，更使数学模型模拟变得可行。它有助于克服野外调查方法的困难。 数学模型模拟的优点主要可概括为如下： ( 1 ) 模型是自然系统的抽象。许多不容易测量的过程可用具有一定假设的数学算 法进行描述。因此，它们有助于提高对真实系统的理解； ( 2 ) 模型可用于把个别过程放大到大比例，因此，它们是进行数据综合和研究整 个系统( 如流域等) 在多重时间尺度上对多重压力( 如土地利用变化、气候变化等) 响 应的有效工具； ( 3 ) 模型有益于绝对预测，有助于对自然和人文影响分析；模型也可重复运行以 模拟不同情况下的相对响应： ( 4 ) 模型可能是高成本，甚至是对培训模型构造者是耗费财力的，但数学模型被 认为是规划任何水土保持措施的影响、规划水资源、设计水质和水文系统的最有效工 3 第一章绪论 具。 综上所述，构造数学模型模拟土壤侵蚀过程，计算土壤侵蚀量具有比野外调查更多 的优点，而且它恰好也便于与计算机技术相结合，发挥计算机运算速度快的特长，以解 决数学模型的计算问题。另外，它便于和其它计算机软件，如地理信息系统相连，完成 更复杂的工作，因此其意义是重大的【3 1 。 构建黄土高原的降雨侵蚀预测模型可以改善黄土高原水土流失现状。由于降雨侵蚀 是引起黄土高原水土流失的一大因素，也是造成黄土高原土壤侵蚀的主要侵蚀过程，建 立降雨侵蚀预测模型对黄土高原土壤侵蚀定量分析有促进作用，并可以将建模过程的经 验总结推广，为降雨侵蚀严重地区提供有效的水土流失治理方案。 1 2 国内外土壤侵蚀研究进展 国外在土壤侵蚀研究方面起步较早，特别是在通用水土流失方程( u s l e ：u n i v e r s a l s o i ll o s se q u a t i o n ) 提出后，对土壤侵蚀的研究从概念定性研究提升到可以利用模型定 量研究，使通过测试数据构建数学模型成为可能。 我国土壤侵蚀定量评价的研究始于上世纪4 0 年代，上世纪8 0 年代初开始统计模型 的研究和开发，8 0 年代末开展了土壤侵蚀预报物理模型的探索，同时还在较大区域乃至 全国范围内进行了土壤侵蚀和水土保持的评价研究，特别在地貌复杂的黄土高原地区研 究较多。中国土壤侵蚀与水土保持定量评价研究可概括为三个层次和两种类型，即地块 小流域区域( 或国家) 三个层次，统计模型、物理模型两种类型。 土壤侵蚀预报模型研究是国际土壤学研究的前沿领域。由于我国土壤侵蚀环境的复 杂、侵蚀类型和过程的特殊及人类活动影响较大，使国际上比较成熟的侵蚀预报模型在 我国的应用受到极大的限制。研究开发中国土壤侵蚀预报模型，是水土保持实践的需求， 也是水土保持学科发展的需要。 中国土壤侵蚀预报模型的开发以已有的长期实验观测数据为基础，分析坡面、小流 域和区域等多种尺度土壤侵蚀发生发展的过程与机制；认识影响土壤侵蚀的主要因子， 并提取量化指标；从时间和空间上揭示径流形成、汇集及各种侵蚀方式中的土壤分散及 泥沙搬运与沉积的动态过程；并进行模拟与集成；吸收和借鉴国际土壤侵蚀预报的最新 研究成果，建立能预报复杂地形条件下的具有明确物理学意义的土壤侵蚀预报模型，为 我国水土保持与生态环境建设等决策工作提供科学依据，并带动我国土壤侵蚀及其相关 现代地表过程研究进入世界先进行列【4 1 。 4 长安大学硕士学位论文 对于水土流失严重的区域如长江中游，黄土高原地区等地，我国一直投入大量人力 物力，希望解决水土流失造成的危害。 作为我国土壤侵蚀预报模型开发的预研究，从2 0 0 1 年1 2 月立项以来，由水利部水 保监测中心和中科院水土保持与生态环境研究中心共同主持的“水土流失预报模型应用 研究及推广项目已完成了预定的研究任务，为国内开展土壤侵蚀预报模型研究和开发 提供了有益的借鉴和指导。该项目对国内外土壤侵蚀预报模型的研究成果进行了系统地 调研和总结，在土壤侵蚀预报模型评述、中国侵蚀模型开发战略、中国土壤侵蚀模型的 基本框架设计等方面取得了显著的进展。专家们认为，该项目针对国家水土保持和生态 环境建设的需求，就我国土壤侵蚀预报模型开发有关问题进行先期研究，对系统地开展 我国土壤侵蚀预报模型的开发研究工作具有重要意义；项目对国内外各种模型的开发背 景、基本结构、计算参数、使用范围、经验教训和推广应用的可能性等进行了系统分析 和评价，为国内开展土壤侵蚀预报模型研究提供了有益的借鉴；根据我国侵蚀环境的特 殊性，该研究提出了中国土壤侵蚀预报模型研究的初步框架，符合我国土壤侵蚀研究的 实际情况，为我国土壤侵蚀模型开发提供指导【5 1 。 在黄土高原土壤侵蚀数学模型的研究方面，我国政府给予了很高的重视，2 0 0 6 年9 月1 6 1 7 日，黄委组织在郑州召开会议，对黄土高原水土流失数学模型研发任务书 进行了审查，并通过了审查。 会议围绕着土壤侵蚀数学模型研发的目的、目标、主要类型与框架以及关键技术路 线等主要内容进行了讨论，并提出了审查意见，主要内容：为满足黄土高原水土保持 规划设计、效益评价和黄河水沙调控体系科学运行等工作的迫切需要，开展黄土高原土 壤侵蚀数学模型研究是非常必要的。模型研发要充分吸纳国内外现有模型研究成果的 精华。总体目标是以实测资料较为丰富的小流域为模拟对象，构建不同类型区的土壤 侵蚀模型，逐步建立中大尺度模型。近期目标是用4 年时间( 2 0 0 7 2 0 1 0 年) ，建立岔 巴沟小流域的机理模型、经验模型和中尺度典型支流土壤侵蚀经验模型。研究内容主 要包括：模型调研和实地查勘；论证确定模型类型与构架，编制概要设计报告；构建产 流产沙模型；构建信息支撑体系；可视化构件设计；模型率定与验证掣副。 1 3 本文研究内容 本文是利用观测到的黄土高原降雨相关数据和土壤侵蚀数据，研究黄土高原地区降 雨侵蚀预测模型，通过参考通用水土流失方程的框架，使用回归分析和数据挖掘相关的 5 第一章绪论 数据预处理方法，得出降雨侵蚀预测模型的数学方程式，并给出评价结论。 本文难点是降雨侵蚀数学方程式的框架无法预知，需要通过曲线拟合得出几个变量 之间的关系，通过试验比较得出最佳的模型。另外，由于土壤侵蚀并不是由降雨一个因 素决定，而是由很多其他侵蚀因素与之共同作用，所以在用降雨数据研究的时候，其他 因素产生的误差必然存在，如何分析数据的结构特征，并通过数据清洗技术除去噪声从 而减少误差，是构造精确数学方程的关键。 1 4 技术路线 本文以水土流失和降雨侵蚀主要理论为指导，结合实际数据，通过使用s p s s 统计 分析软件进行试验，主要的研究方法是使用回归分析建立数学模型，再通过使用数据挖 掘分析技术对模型进行优化。 本文构造降雨侵蚀的主要方法是线性回归分析，对数据的处理主要使用的是数据挖 掘的数据分析技术，辅助的分析方法还有曲线拟合、数据统计分析、方程线性变换方法 等。 主要分析过程和使用技术如下： ( 1 ) 利用线性回归技术和未通过处理的原始数据建立线性模型，求出模型的评价 指标，这个过程求出的结果是作为以后分析的基础，用于比较其后使用方法的好坏。 ( 2 ) 通过数据挖掘的数据预处理技术分析原始数据，并根据实际需要建立数据预 处理算法，将数据优化后再用线性回归进行分析。 ( 3 ) 通过曲线拟合方法建立乘法模型，通过线性变换和线性回归分析求出模型的 数学方程的表达式和评价指标，并与前边的分析结果比较。 ( 4 ) 在各个分析过程中，根据需要利用不同的图形分析，可以直观看出分析的好 坏。 1 5 论文结构 本论文的内容围绕降雨侵蚀预测模型的构建过程，共分为六章： 第一章绪论：阐述土壤侵蚀的相关概念，以及课题研究的意义，并对土壤侵蚀数 学模型，特别是降雨侵蚀预测模型的研究作了介绍。本章分析了研究的重点，给出了全 文研究的技术路线。 第二章土壤侵蚀理论与降雨侵蚀研究：本章从土壤侵蚀的分类和成因描述开始， 6 长安大学硕上学位论文 进而介绍了土壤侵蚀研究的主要方法，在此基础上对通用水土流失方程和相关概念作了 阐述。 第三章降雨侵蚀预测模型的建立方法和原理：首先从大的方面介绍了土壤侵蚀数 学模型的建立方法，进而介绍了降雨侵蚀中的几个因子，重点描述了本论文进行分析的 主要方法回归分析。 第四章数据预处理技术和模型的评价方法：简单介绍了数据挖掘中数据预处理的 相关技术，并对使用的模型分析方法做了阐述，最后利用s p s s 和降雨侵蚀原始数据进 行线性回归建模，用作分析比较的基础。 第五章结合数据预处理技术的线性模型研究：详细介绍了本论文使用的数据预处 理技术，建立降雨侵蚀数据预处理算法，并使用经过数据预处理的数据进行回归分析， 建立数学模型，并与第四章得出的模型比较得出结论，提出仍需解决的问题。 第六章降雨侵蚀预测模型的乘法模型研究：首先介绍了曲线拟合的概念，并介绍 使用曲线拟合的目的。通过分析给出降雨侵蚀预测模型的乘法模型，通过线性化变换和 回归分析得出乘法模型的方程表达式，并将其分析结果和以前的模型分析结果作比较， 得出结论。 总结与展望对降雨侵蚀预测模型的建模过程中遇到的问题作了总结，归纳了在研 究中尚未解决的问题，分析了模型的不足之处，最后展望了前景和相关的后续研究内容。 7 第二章降雨侵蚀理论与研究 第二章降雨侵蚀理论与研究 土壤侵蚀分类的内容包括侵蚀的类型、侵蚀强度、抗侵蚀程度等。根据外营力的不 同，土壤侵蚀可以分成水力侵蚀、重力侵蚀、风力侵蚀、冻融侵蚀、冰川侵蚀、混合侵 蚀、生物侵蚀等几类。按侵蚀成因上有无人类扰动，土壤侵蚀可分为地质侵蚀( g e o l o g i c e r o s i o n ) 和加速侵蚀( a c c e l e r a t e de r o s i o n ) 等。按土壤侵蚀发生的时代可分为古代侵蚀 ( a n c i e n te r o s i o n ) 和现代侵蚀( h u m a ne r o s i o n ) 。一般地，侵蚀强度和抗侵蚀程度以分 级评价的方法来研究。 2 1 土壤侵蚀的成因 土壤侵蚀是自然现象，其产生的原因既有自然因素，也有人为因素，黄河中游黄土 高原地区植被稀少土壤疏松、暴雨较多、地形破碎，易产生了强烈的土壤侵蚀。在研究 中，能够造成水土流失的土壤侵蚀研究因素是需要注意的。由于土壤侵蚀过程是持续的 自然过程，随时随地都会发生，土壤侵蚀到一定程度，加上外在一些条件如降雨等，从 而形成水土流失。在土壤侵蚀研究的实际意义中，形成水土流失的土壤侵蚀因素需要重 点研究，因此我们介绍的土壤侵蚀因素都是形成水土流失的主要成因，包括自然和人为 两方面。 2 1 1 自然因素 从自然方面考虑，主要有地形、降雨、土壤( 地面物质组成) 、植被四个方面。 ( 1 ) 地形。地面坡度越陡，地表径流的流速越快，对土壤的冲刷侵蚀力就越强。 坡面越长，汇集地表径流量越多，冲刷力也越强。黄土丘陵区、地面坡度大部在1 5 0 以 上，有的达3 0 0 ；坡长一般1 0 0 2 0 0 m ，甚至更长。每年每亩流失5 1 0 吨，甚至1 5 吨以 上。 - ； ( 2 ) 降雨。产生水土流失的降雨，一般是强度较大的暴雨，降雨强度超过土壤入 渗强度才会产生地表( 超渗) 径流，造成对地表的冲刷侵蚀。 ( 3 ) 地面物质组成。质地松软，遇水易蚀，抗蚀力很低的土壤，如黄土、粉沙壤 土等是产生水土流失的对象。 ( 4 ) 植被。达到一定郁闭度的林草植被有保护土壤不被侵蚀的作用。郁闭度越高， 保持水土越强。 长安大学硕士学位论文 2 1 2 人为因素 人类对土地不合理的利用、破坏了地面植被和稳定的地形，以致造成严重的土壤侵 蚀，引起大面积水土流失，最主要的有两个方面。 ( 1 ) 毁林毁草、陡坡开荒，破坏了地面植被。 ( 2 ) 开矿、修路等基本建设不注意水土保持，破坏了地面植被和稳定的地形，同 时，将废土弃石随意向河沟倾倒，造成大量新的水土流失。 2 2 土壤侵蚀研究方法 土壤侵蚀研究方法在本节特指降雨侵蚀研究方法，主要研究土壤侵蚀量和降雨、地 形等因素的关系。土壤侵蚀量是土壤侵蚀定量分析的指标之一，它是流域内被侵蚀土壤 的数量表示。 土壤侵蚀量是指土壤侵蚀作用的数量结果，通常把土壤、母质及地表松散物质在外 营力的破坏、剥蚀作用下产生分离和位移的物质量，称为土壤侵蚀量。单位时间：、单位 面积内产生的土壤侵蚀量( 单位通常为t k i n 2 a ) ，称为土壤侵蚀速率( 或速度) 。在特 定时间段内通过小流域出口断面的泥沙总量( 或径流总量) ，称为流域输沙量 或径流 输出量) ，一般以年平均输沙量( 单位通常为t a ) 计量。 从2 0 世纪2 0 3 0 年代起，一些学者注重研究土壤侵蚀的发生和发展规律，并相应 地产生了一些研究方法，在这些研究方法中有用于定性研究和半定量研究的田间观察法 和侵蚀针法，有用于定量研究的径流小区观测法、单源流域或汇水区的大面积试验区法、 核素示踪法和定量遥感法。 2 2 1 田问观测法 这种方法是指有经验的观测人员通过观测细沟、冲沟、底土裸落、沉积物在径流流 经平缓地段流速减慢时的沉积情况，将侵蚀土壤与同类林地或其他保护地土壤相比，用 表土减少的情况或其它情况来鉴别土壤侵蚀程度，属于定性研究，反映了土壤侵蚀研究 的初步阶段。 2 2 2 侵蚀针法 该法通过埋设特制有刻度的铁钉或水泥钉，观测土层在铁钉或水泥钉刻度上的位置 变化，根据量算侵蚀区或沉积区的面积和土壤容重来估算土壤侵蚀量或沉积量，该法只 能用于粗测土壤侵蚀量( b 兆宏，1 9 9 5 ) 。 9 第二章降雨侵蚀理论与研究 2 2 3 径流小区法 径流小区法是指用木板、石棉、混凝土、铁皮或其他隔湿材料围成矩形小区，在较 低的一段安装收集槽和测量设备，以确定每次降雨的径流量和土坡流失量。径流小区法 是定量研究土壤侵蚀最广泛采用的一种方法。它可以适应不同坡长、坡度、土坡类型、 土壤物理性状、作物轮作制度和植物类型的研究条件，并提供上述诸条件对径流和侵蚀 影响十分有用的资料。利用这些资料已建立了许多土壤侵蚀模型，其中最著名的是2 0 世纪5 0 年代，w i s c h m e i e r 和s m i t h 运用径流小区法建立的通用水土流失方程_ i i s l e ， 并于1 9 7 8 年由美国政府农业局和水土保持局对外公布，得到了广泛的应用。该法的特 点是准确、资料获取及时，但该法代表的面积较小，往往不能反映大范围空间的实际情 况，特别是对地形变化较大的坡面部位侵蚀强度垂直分布的测定无能为力，而侵蚀强度 的垂直分布规律是侵蚀区泥沙来源及侵蚀区土地整治，水土保持措施合理实施的重要依 据。对风蚀和水蚀都存在的地区不能获得风蚀结果，另外该法如大范围的布设小区，则 费用较高，且实地观测费时费力。 2 2 4 单源流域或汇水区的大面积试验区法 这种方法实际上是径流小区法的扩展，径流小区法因研究的范围较小，不能代表整 个坡面或整个单源流域或汇水区，因此水土保持研究工作者以单源流域或汇水区作为研 究对象，建立大面积的试验区。这种大面积的试验区可以是一个小的集水区或者至少大 到足以包括一条自然排水系统的小流域，当然这种大面积试验区的界限是以其分水岭为 界限的，在单源流域或汇水区的出水处设立沉沙收集池和测量设备，以测定每次降雨的 径流量和泥沙流失量。这种大面积试验区法的特点是可以同时观测所有的侵蚀过程和水 土保持措施的综合结果。其结果具有代表性和可信性，但试验结果重复性较差。 2 2 5 核元素示踪法 早在2 0 世纪6 0 年代初，m e n z e l ( 1 9 6 0 ) 运用沉降放射性核元素来研究土壤侵蚀规 律，为土壤侵蚀研究的新方法和新技术开辟了一条新路。2 0 世纪7 0 年代以来，美国 ( r i t c h i e e t a l 1 9 7 4 ：m c h e n r y r i t c h i e ，1 9 7 7 ) ，加拿大( d ej o n ge ta l ，1 9 8 6 ；k a c h a n o s k i ， 1 9 8 7 ) 、澳大利亚( l o n g m o r e e ta l ，1 9 8 3 ；c o m p b e l le ta l ，1 9 8 6 ) ，英国( w a l l i n g b r a d l e y ， 1 9 8 8 ；w a l l i n g & q u i n e ，1 9 9 0 ) 和中国( 张信宝等，1 9 8 9 ；汪阳春等，1 9 9 1 ) 等科学家 相继运用大气沉积的放射性元素来研究农林系统的土壤侵蚀和沉积的分布规律。2 0 世纪 8 0 年代美国学者k n a u s ( 1 9 8 9 ) 利用稳定性核素示踪法研究沼泽地土壤侵蚀和沉积分布 1 0 长安大学硕士学位论文 规律。9 0 年代澳大利亚科学家开展了7 b e 、1 3 7 c s 、2 1 0 p b 、2 2 6 r a 、2 3 2 t 1 l 等多种核素示踪 技术来研究小流域的泥沙来源( o l l c yj m ，1 9 9 3 ) 。我国学者从8 0 年代起利用核素示 踪法进行土壤侵蚀的研究，如张信宝( 1 9 8 9 ) 等利用”7 c s 来研究农耕地的土壤侵蚀； 李小龙( 1 9 9 5 ) 等利用2 2 6 r a 分析法研究t d , 流域的泥沙来源；马建国等( 1 9 9 0 ) 利用 土壤特征元素( 核素) 研究土壤中泥沙运移规律。核示踪技术应用促进了土壤侵蚀量化 研究进入了新的阶段。目前土壤侵蚀研究中的核素示踪技术的应用可分为放射性核素示 踪法和可活化稳定性核素示踪法。 2 2 6 定量遥感法 遥感作为一种不直接接触目标物体，通过远距离接受地物发射的电磁波信息而识别 物体的高新技术。我国2 0 世纪8 0 年代以来航天、遥感技术已广泛应用于土壤侵蚀的监 测，卜兆宏( 1 9 9 5 ) 等应用遥感技术结合g i s 方法以及结合w i s c h m e i e r 的通用水土流 失方程估算了我国南方红壤地区水土流失，并和实测资料进行了比较，表明应用遥感技 术结合g i s 方法能较好地预测水土流失量i r l 。 2 3 降雨侵蚀因子 2 3 1 含义 降雨侵蚀是降雨和土壤的相互作用。降雨的作用称为侵蚀力，土壤的作用称为可蚀 性。具体讲，降雨侵蚀力是指降雨引起土壤侵蚀的潜在能力，它是降雨的物理性质的函 数，在通用水土流失方程式中用r 表示。由降雨引起的水土流失大致可以分为两个过程： 一是土粒从土体中被分离；二是分离出的土粒被冲走，雨滴降落到地表，引起土壤击溅 分离，一部分雨水渗入地下，另一部分形成地表径流。美国水土保持学家维斯奇迈尔 ( w i s c h m e i e r ) 通过对各种降雨不同组合形式的土壤流失量进行相关分析，提出了降雨 的动能与其3 0 分钟最大强度的乘积是判断土壤流失的最好指标，其表达式为： 尺= 1 3 0 ( 2 1 ) 一 w 式( 2 1 ) 中：r _ 一次降雨的侵蚀力( j m 2 ) ； b 一一次降雨过程中某段降雨量所产生的动能( j m 2 ) ； 厶r 一次降雨过程中连续3 0 分钟的最大降雨强度( c m h ) 。 由于该方程是根据美国的资料求得的，而世界各国各地降雨、地形等自然地理因素 存在显著差异，所以许多国家的水土保持工作者根据本国本地区实际情况，经过多年研 第二章降雨侵蚀理论与研究 究，在通用方程基础上得出了许多计算尺的经验公式。诸如魏斯曼经验公式： 尺= 舡7 3 2 m 5 i g 知1 8 8 】 ( 2 2 ) 式( 2 2 ) 中：融各月平均降雨量( r a m ) ： p - 4 平均降雨量( m m ) 。 2 3 2 影响因素 由计算尺的通用方程可看出，影响降雨侵蚀因子尺的因素主要是暴雨动能和最大 降雨强度，而暴雨动能和最大降雨强度又受降雨量、降雨时段等因素的影响。为此，则 应从降雨量、降雨时段、降雨强度、雨滴能量大小几方面来阐述降雨侵蚀因子r 对土壤 流失强度的影响并找出引起土壤侵蚀的临界值或临界范围。 1 降雨量 降雨量是降雨特性指标之一，区域性降雨量级的划分见表2 1 。 一般l o m m 以下的小雨不产生径流，也不存在侵蚀；中雨能产生径流，但冲刷量 微弱，引起水土流失的强度不大。随着降雨量的增加，水土流失强度逐渐增大。当降雨 量达到5 0 r a m 时，将产生明显的水土流失，且强度较大、径流快，所以一般将降雨量 = 5 0 r a m 作为产生明显水土流失的临界值。 表2 1 区域性降水量级界定 量级 降雨量下限( m m ) 无雨 0 小雨0 2 ( 或0 5 ) 小中雨 5 中。大雨 1 7 大暴雨3 8 暴雨 5 0 大暴雨 1 0 0 特大暴雨 2 0 0 2 降雨时段 降雨因素中的暴雨强度，降雨时空分布是形成土壤侵蚀的动力条件。在一定降雨条 件下，降雨历时越短，强度越高，产生的土壤流失强度越大。据中国水土保持有关专家 1 2 长安大学硕士学位论文 的研究表明：一般高强度暴雨的降雨时间大多在半小时左右，即使是时间较长的一些暴 雨，其雨量也往往集中在1 小时或半小时左右( 据统计，在此时段的降雨可占总雨量的 7 0 ) 。其余时间的降雨为低强度的间歇性降雨，一般不产生水土流失。 由降雨侵蚀因子尺的计算方程式可看出，不同时期的降雨