




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
/长江水质的数学评价和预测朱金杨浩乐摘要:本文应用插值原理对长江水质现状进行综合评价,分析了各地区水质的污染状况,并应用状态空间模型和改进的kalmam递归式对长江未来水质状况进行预测分析。得出:如果不及时加以治理,长江污染状况将会随着工业的发展和时间的推移迅速恶化;并根据计算结果提出解决长江水质的建议和意见。关键词:水质,插值评判,状态空间模型,Kalmam递归式。问题复述:水是生命之源,是人类赖以生存的资源,保护水资源就是保护我们自己,是我们每个人应尽的责任。对于我国大江大河水资源的保护和治理应是重中之重.长江是我国第一、世界第三大河流,是华夏5000年文明的发祥地.如今,长江水质污染程度日趋严重已经引起相关政府部门和专家们的高度重视。根据国家环保局的政府网站〔〕2003年6月至2005年9月发布的长江沿线8个省、17个水文监测站的水质报告,长江水利委员会编辑出版的《长江年鉴》1995年至2004年长江流域水质报告表以及2004年4月至2005年4月发布的长江干流主要观测站点的基本数据表明,原本“一江清水向东流”的母亲河正在悄悄被污染,而且污染指数在急剧攀升.2004年10月10日,全国政协与中国发展研究院联合掀起了“保护长江万里行”活动。11月13日,南方新闻中心〔〕刊登了《专家称若不及时拯救,长江生态10年濒临崩溃》的文章.12月14日根据2002年6月国家重新修订颁布的《地表水环境质量标准》,长江水质中Ⅳ类,Ⅴ类,劣Ⅴ类的比例正在逐年增加。新标准采用单因子极值评价方式似乎过于严格。所谓“极值法”是指以其中最差的指标来评价水质等级,也就是说水质标准当中只要1项达不到Ⅴ类,那么整体水质就要被判为劣Ⅴ类。这种“一票否决”式的考评方式有一定的片面性,对水质的考评应该综合各项指标后得出结论,现在只考虑其中1项指标不是太科学。为此,应用插值原理,并采用国家水利部门提供的最新监测数据〔附件3〕,对长江水质进行多因子线性综合评价。对于长江未来水质的预测,根据资料〔附件4〕中数据,通过分析实验,长江水质污染状态具有很强的随机性,因此这种状态在时间状态的变化下具有“无后效性”的特点.为此,本文采用了状态空间转移模型和改进了的Kalmam递归式进行了较为合理、准确的预测分析.问题〔1〕:对长江近2年多的水质情况做出定量的综合评价,并分析各地区水质的污染状况。1插值评判原理1。1定义:已知n+1节点,j=0,1…,n,其中互不相同,不妨设,求任一插值点x*〔≠〕处的插值y*。可以看成是由某个函数y=g〔x〕产生的,g的解析表达式可能十分复杂,或不存在封闭形式,也可以未知.㈠1.2求解的基本思路是:构造一个相对简单的函数y=f﹙x﹚,使f通过全部节点,即f﹙﹚=﹙j=0,1,…,n﹚,再用f﹙x﹚计算插值,即y*=f﹙x*﹚。见参考文献1。2水质综合评价的SP插值模型设根据水质评价标准产生的某次水样的标准水质等级及其水质指标分别为y﹙i﹚及﹛x*﹙i,j﹚|j=1~m﹜,i=1~n。其中,n、m分别为样本容量和水质指标数目。污染越严重,水质等级就越高,最低水质等级设为1、最高水质等级设为N.为消除各水质指标的量纲效应,使建模具有一般性,对水质指标进行标准化处理=(1)式中:Ex(j)、Sx(j)分别为原第j个水质标准{x*(i.j)|I=1~n}的均值和标准差。当得到研究水体的水质指标值{x(n+1,j)|j=1~m}后,利用上述n个样本内插研究水体的水质等级yc(n+1),使下式:E=达到最小.式(2)中:(2)式中:di为第i个样本的水质指标与研究水体水质指标之间的距离:wi为权重,表示第i个样本对内插研究水体的水质等级yc(n+1)的贡献大小;b为待定参数,一般为大于1的常数,b取得越大,则在点{x(n+1,j)|j=1~m}附近的拟合曲面将变得越平坦,而使远离点{x(n+1,j)|j=1~m}处的拟合曲面将变得越陡峻.对式(2)求导数并令其为0,可解得式(2)的最小值为:(3)这就是所求的对应点{x(n+1,j)|j=1~m}的水质等级值。建立SP模型的步骤可归纳为如下3步:(1)根据水质评价标准随机生成水质等级样本系列x(i。j)及y(i),i=1~n,j=1~m。(2)根据样本系列对参数b进行优化估计.在样本系列中任取某样本I,由其它n-1个样本进行Shepard插值,得到相应于水质等级y(i)的插值记yc(i).可通过求解如下优化问题来优化估计参数b(4)s.t1≤b≤5式(6)的区间是根据笔者的经验确定的。这是一个一维非线性优化问题,模拟生物进化过程中优胜劣汰规则与群体内部染色体信息交换机制的加速遗传算法(AcceleratingGeneticAlgorithm,简称AGA)[7]是一种通用的全局性优化方法,用它来求解该问题显得十分简便而有效。AGA的详细算法可参考文献[7]。(3)进行水质综合评价.当得到研究水体各水质指标值{x(n+1,j)|j=1~m}后,与样本系列一起代入式(3)、式(4),即可由n个样本插出研究水体的水质等级yc(n+1),作为该水质综合评价的结果。3模型求解现以长江水系四川攀枝花地区水质综合评价为例,应用SP插值模型求解。水质评价标准见附表:《地表水环境质量标准》附表:《地表水环境质量标准》(GB3838-2002)中4个主要项目标准限值单位:mg/L序号
分类
标准值
项目Ⅰ类Ⅱ类Ⅲ类Ⅳ类Ⅴ类劣Ⅴ类1溶解氧(DO)
≥7.5(或饱和率90%)653202高锰酸盐指数(CODMn)≤2461015∞3氨氮(NH3-N)
≤0。150.51.01。52.0∞4PH值(无量纲)6-——9根据附表和附件(3)利用4个主要项目标准限值对四川攀枝花地区28个月水质生成6个等级的样本系列,对应1、2、3、4、5、6,其中1、2、3类为可饮用水,4类为轻度污染,5类为中度污染,6类和6类以上为重度污染.利用SP插值模型求解可得表1-1(见附录),同理可得表1—2至表1-17(见附录),此为长江水系17个地区28个月内水质具体状况。根据表1中的17张表中每个地区在28个月内各种等级水样所占比例得表2。ⅠⅡⅢⅣⅤ劣Ⅴ四川攀枝花龙洞0.2857140。5357140.1071430。07142900四川乐山岷江大桥00.1428570。3214290.3214290。2142860四川宜宾凉姜沟0.0357140.6785710.1428570。0357140.1071430四川泸州沱江二桥0。1071430.4642860。1428570.1428570.0357140.107143重庆朱沱00。7142860。285714000湖北宜昌南津关00。8928570。107143000湖南岳阳城陵矶00。6428570.357143000湖南长沙新港00.1071430。50.39285700湖南岳阳岳阳楼00。5357140.3571430.10714300湖北丹江口胡家岭0。7142860.2857140000湖北武汉宗关00。6785710。2857140.03571400江西九江河西水厂0.0714290.9285710000江西南昌滁槎000。1071430.3571430.0357140.5江西九江蛤蟆石00。6785710。1785710.14285700安徽安庆皖河口0。0357140.9285710。035714000江苏南京林山0。0714290.8928570.035714000江苏扬州三江营00。7142860。250.03571400表2表2为17个地区近两年多水质污染状况。(例如:四川攀枝花地区,1级水比例为0.285714;2级水比例为0.535714;3级水比例为0。107143;4级水比例为0.071429;5级水、劣5级水比例均为0)由表1-1至表1—17(见附录)中SP计算值的平均值与GB平均值对长江近两年多水质情况定量的综合评价得表3:17地区两年水质综合评价sp均值GB均值1四川攀枝花龙洞1.9833181.9642862四川乐山岷江大桥3.9867463.6785713四川宜宾凉姜沟2。5833712.54四川泸州沱江二桥3。7658682。9642865重庆朱沱2.1460212.2857146湖北宜昌南津关2。2617362.1071437湖南岳阳城陵矶2。2779792.3571438湖南长沙新港3。2972753。2857149湖南岳阳岳阳楼2.3079892.57142910湖北丹江口胡家岭1。3779821。28571411湖北武汉宗关2。3643292.35714312江西九江河西水厂1.9814461。92857113江西南昌滁槎4.6827714.92857114江西九江蛤蟆石2.4261462.46428615安徽安庆皖河口2.068407216江苏南京林山1.9926641。96428617江苏扬州三江营2.5079822.321429表34误差分析利用误差绝对值落在绝对误差区间的百分比(%)得表4:长江水质等级的标准值与SP模型计算值之间误差分析结果。误差绝对值小于下列标准的百分比0.10.20.30。40.50。60.70。80.911441711872312803163563793893970.3010.3580.3920.4840.5870。6620.7460.7940.8150。8321.11.21.31。41.51。61.71。81。924204264354414474514554584624640。8800.8930.9110。9240.9370.9450。9530。9600.9680.9726由于SP模型对17个地区的综合评价结果绝大部分与国家标准上限值相近,且精确度较高,因此SP模型评价的长江全流域水质综合等级与各地区综合等级均是合理准确的。问题(2)研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?解答:在问题1中表1—1至表1-17里,依据高锰酸盐指数和氨氮含量取出列向量生成表5-1攀枝花时间DOCODMnNH3-N水质实际标准实际标准实际标准2003年6月6.8≥50。2≤60.1≤1。0II2003年7月8.1≥55.6≤60.1≤1.0III2003年8月8。1≥50。5≤60.15≤1。0I2003年9月8.8≥51。1≤60.11≤1。0I2003年10月8。6≥53.1≤60。17≤1.0Ⅱ2003年11月9.1≥51.5≤60.06≤1。0Ⅰ2003年12月9.39≥52.5≤60.07≤1.0Ⅱ2004年1月9.47≥52≤60。42≤1.0Ⅱ2004年2月9.49≥50.4≤61.22≤1.0Ⅳ2004年3月8.97≥52.3≤60.13≤1。0Ⅱ2004年4月9。04≥52.3≤60.15≤1。0Ⅱ2004年5月9.09≥54.3≤60.07≤1。0Ⅲ2004年6月8.79≥52。5≤60.04≤1.0Ⅱ2004年7月8.46≥52.4≤60。04≤1.0Ⅱ2004年8月13.9≥55.8≤61≤1。0Ⅲ2004年9月9.42≥56.1≤60.09≤1。0Ⅳ2004年10月9。72≥50.8≤60.08≤1.0Ⅰ2004年11月9。34≥52.8≤60.06≤1.0Ⅱ2004年12月10。3≥51。6≤60.08≤1。0Ⅰ2005年1月10.2≥51.2≤60。07≤1。0Ⅰ2005年2月9.38≥50。9≤60.15≤1。0Ⅰ2005年3月9。6≥51。1≤60.26≤1.0Ⅱ2005年4月9.03≥51。1≤60。1≤1。0Ⅰ2005年5月8。38≥52.5≤60.07≤1。0Ⅱ2005年6月8.8≥52。9≤60。08≤1.0Ⅱ2005年7月8.62≥54≤60.08≤1.0Ⅱ2005年8月8.93≥52。6≤60。09≤1.0Ⅱ2005年9月8。5≥54≤60.08≤1.0Ⅱ表5-1此表为四川攀枝花地区高锰酸盐和氨氮指数统计情况表。同理可得长江干流6个地区的高锰酸盐和氨氮指数统计情况表5-2至5—6(见附录)。根据表5的数据可知长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要来自支流地区,具体较重的是南昌、泸洲、乐山、长沙地区。问题(3)假如不采取更有效的治理措施,依照过去10年的主要统计数据,对长江未来水质污染的发展趋势作出预测分析,比如未来10年的情况。3.1状态空间转移3.1。1定义参考《时间序列的理论与方法》〔田铮译〕第二版第12章第一节3.1.2求解方法参考《时间序列的理论与方法》〔田铮译〕第二版第12章3.2。1由附件4中1995-—2004年长江流域水质分析可知无论评价河长为多少,五类水等级的百分比包含各种评价指标、长江水体的自我清洁能力等各方面因素,且这种百分比无论是在时间还是空间上均无后效性,呈现随机状态,因此以每种等级的百分比作为状态变量。通过状态转移基本模型多次实验分析知,此种状态在向下一状态转移的过程中符合Xt=Ft*Xt-1+Vt,依据表3—1中数据做状态预测,通过误差分析可知:Vt=0时预测结果更加接近监测值。因此以改进的Kalmam递归式:Xt+1=Ft*Xt—1为预测模型。3.2.2模型求解利用改进的Kalmam递归式和题目附件4里面对长江10年来主要统计数据进行状态预测,得表6(见附录)。对表6做归一化处理得表7(见附录),即是对长江未来10年的发展趋势的预测结果。同样,利用这个模型和10年来长江总流量和废水量来预测出未来10年中每年的总流量和废水量表8(见附录)。由表7分析可知,随着时间的推移,长江水质污染程度日趋加重,即1,2,3类水比例逐年减小,4,5,6类比例逐年增加。长江的总流量每年都成递减趋势,而废水量日趋增加.、3。2.3模型评价(1)此模型为随机状态下的状态转移模型,不具有连续性,因此比较符合长江未来水质状况的实际情况。(2)在对系数矩阵Ft可靠性检验的过程中,次误差范围控制在0.05%之内,因此此模型比较可靠,具有良好的适应范围和准确性。(3)此模型可用于符合假设条件的任意随机状态,能够对下一状态作出比较准确的预测分析,具有广泛推广的意义。问题4根据你的预测分析,如果未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内,且没有劣Ⅴ类水,那么每年需要处理多少污水?解答:根据问题3中对未来10年长江水质的预测数据,结合题设要求设:每年要处理的废水量Q;每年长江总流量AQ;4、5、6级水的百分比分别为、、,其中、、为水文年的数据; 则根据题设建立等式
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 甘泉县2024-2025学年三下数学期末监测试题含解析
- 提升保安证考试成绩的正确姿势试题及答案
- 湖北省襄阳五中2025年高三考前实战演练历史试题含解析
- 河南水利与环境职业学院《病原生物学与医学免疫学》2023-2024学年第二学期期末试卷
- 2025年保安证考试小技巧试题及答案
- 闽南师范大学《钢琴舞台表演》2023-2024学年第二学期期末试卷
- 广西生态工程职业技术学院《兽医药理学实验》2023-2024学年第二学期期末试卷
- - 2025年化学试题及答案必要症结分析
- 长春理工大学《影视文学赏析》2023-2024学年第一学期期末试卷
- 保安证考试全面复习试题及答案
- 2024年廊坊市第二人民医院招聘工作人员考试真题
- 【珍藏版】鲁教版初中英语全部单词表(带音标、按单元顺序)
- 第三单元分数乘法测试卷(单元测试)北师大版五年级下册数学
- 2025年安徽医学高等专科学校单招职业适应性测试题库完整版
- 医院知识产权培训课件
- 2025年度KTV娱乐门店转让协议
- 人教部编版道德与法治八年级下册:3.2 《依法行使权利》听课评课记录
- 机电一体化专业课程改革调研报告及改革建议
- 2025年甘肃甘南州国控资产投资管理集团有限公司招聘笔试参考题库附带答案详解
- 新生儿各种导管的护理
- 《天津天狮奖金制度》课件
评论
0/150
提交评论