专利名称：一种计算历史时期湖泊水位的技术的制作方法
技术领域：
本发明属于水生态系统监测领域，涉及一种以湖泊沉积物中孢粉组合带、典型理 化指标为依据，利用数学建模的方法计算历史时期湖泊水位的技术。
背景技术：
孢粉具有质轻量多，散布极远的特点。作为生物化石在各沉积层中均有保存，对恢 复古环境状况具有很高的参考价值，尤其在古气候研究中已被广泛应用。历史时期湖泊受 到的人为干扰因素较少，湖泊的水位水量与植物生长基本处于一种动态平衡状态，同时，生 物指标对水位水量的变化比较敏感，如生物量、多度、丰富度等。计算历史时期湖泊水位可 推测湖泊生态系统处于健康状况所需的水文条件，通过监测现代湖泊水位了解生物的生长 状况，及时调水补水，对维持湖泊生态系统的健康起着重要作用。近十几年来许多国内外学者利用孢粉组合带反演湖泊历史时期的环境状况。国 内，有学者利用封闭流域水量平衡原理，结合各降雨参数和孢粉组合带信息计算湖泊的古 降雨量，但该方法有一定的局限性，不适用于半封闭和开放流域降雨量时计算(邵兆刚等， 2004 ；吴敬禄等，1993)。有学者将孢粉组合带与沉积物的磁学指标相结合，定性分析湖泊面 积变化、湖水动力特征及气候状况，此种方法需要结合沉积物粒径、T0C、TN等常规指标综合 分析，同时需要有少量的湖泊历史监测资料做参考(曹希强等，2004)。有学者探索性地采 用有序聚类、回归、主成分分析等数学方法量化孢粉与沉积物理化指标的关系。这些方法需 要大量的样品原始统计数据，可以对湖泊流域环境变迁及人类活动记录等提供解释(张佳 华等，1997 ；王立新等，1995)。国外，在研究沉积物孢粉方面，有一套比较系统的方法，通过 常年监测建立了欧洲、北美、拉丁美洲、非洲、西伯利亚孢粉数据库，使用RCA、DCA、MAT等线 性或单峰模型建立相似集，运用WA、PLS、ANN等数学方法建立校正集，并进行误差估计，分 析物种的丰度、均度，种群结构、时间空间序列，从而建立与历史时期湖泊的气候、物种多样 性、降雨量、地表植被类型等的关系(H. S印pa，2007 ；S. Brewer et al.，2007 ；H. J. B. Birks, 2007 Joanna C. Ellison, 2008 ；N. John Anderson, 1995)，方法具有系统性和连贯性。检索发现预测河川水位的发明专利，中请号为03120092. 3，名称为河川水位预测 装置，公开号为CN1438386，它提供了一种能高精度预测河川水位的预测装置，包括三部分 计量并存储河川水位的水位计量部，计量并存储流域降雨量的降雨量计量部和根据这些计 量值预测河川水位的水位预测模型，水位预测模型包括自回归部分和FIR模型部分，两部 分的参数通过模型鉴别部计算得到。检索发现计算现代湖泊生态需水的发明专利，申请号为200910170085. 7，名称为 基于功能设定法的湖泊生态需水分析技术，公开号为CN101650762，它提供了一种综合考虑 湖泊生态系统不同生态服务功能对淡水输入要求的分析技术，将湖泊生态需水划分为蒸散 需水量、渗漏需水量、水生生物及其栖息地需水、湖泊出湖地表径流需水、湖泊出湖地下径 流需水、能源生产需水、自净需水、航运需水和景观娱乐需水等九种类型，根据最大值原则 和湖泊管理目标差异，确定水资源功能系数，最终确定多目标要求下的湖泊生态需水量。
3
由此可见，我们可以借助数学方法建立孢粉组合带和历史时期湖泊水位的关系， 了解湖泊在低干扰度下的植被类型和水文情况，对于湖泊生态系统健康的维持具有非常重 要的意义。
发明内容
历史时期湖泊水位的计算包括三部分筛选部分，分析部分和合成部分。筛选部分 是运用主成分分析法将相关环境因子转化为几个相互独立的综合指标，分析部分是运用灰 色关联法选取与湖泊水位相关性较大的水生植物作为指示生物，合成部分是将综合指标、 湖泊水位作为自变量，指示生物作为因变量用回归的方法拟合。以近代沉积物样品指标、湖 泊监测水位和孢粉组成数据为校正集，对拟合方程变量的系数进行调整，提高方程推算的 精度。
历史时期湖泊水位计算流程图
具体实施例方式(1)沉积物定年运用放射性元素21°Pb和137Cs，结合沉积模型对沉积物剖面进行定年，确定每层沉 积物对应的年代。①21°Pb定年模型At = A0e A关系式中At为在深度Z处21°Pbex的活度A0为在表层的21°Pbex的活度入是210Pb 的衰变常数(0. 693/22. 3a-1)②137Cs定年模型s = Z/(T0_Tz)关系式中Z为沉积柱深度I；为表层沉积时间Tz为深度Z处沉积时间(2)环境综合指标的构建运用实验方法测得各沉积剖面层的孢粉组成和物理化学指标，包括沉积物pH、盐 度、电导率、粒径、T0C、TN等。主成分分析法(PCA)，是将多个变量通过线性变换以选出较少个数重要信息的一 种多元统计分析方法，他们彼此之间相互独立，又能综合反映原来多个指标的大部分信息， 为使不同量纲数据可比，在计算过程中由软件对数据进行标准化，该计算过程由SPSS16. 0 实现。计算原理A、数据的标准化处理x*ij = Xij－Xj/Sj
其中，i = 1，2，...，11，11为样本点数；」=l，2，...，p，p为样本原变量数目。B、计算数据的协方差矩阵RC、求R的前m个特征值人工彡A2 ^ A3 ^ ,
彡Xm，以及对应的特征向量…，
叫，
Um D、求m个变量的因子载荷矩阵 A =
MP
121，
122，
Mm
pm
U12"/^2 5 …，u U2l、f^, U22 -> …，U
lmVAn 2m VAn
Upl V^", Up2 V^",…’ Upm JK
V p2' ”(3)指示生物的选取湖泊中植物的生长及生存状态与湖泊水位(水量)相关，由于湖泊生态系统的复 杂性和植物的多样性，分析所有的植物花粉特征是不现实的，只有选取部分指示生物作为 样本，指示生物的选取标准是生长状况(沉积层中的孢粉丰度)与湖泊水位的相关性较大。灰色系统理论是一种研究“少数据、贫信息”不确定性问题的新方法，以“小样本、 贫信息”的不确定系统为研究对象，主要通过对“部分”已知信息的生成、开发，提取有价值 的信息，实现对系统运行行为、演化规律的正确描述和有效监控。灰色关联分析是灰色系统 理论的重要内容，它通过计算系统特征序列和相关序列之间的灰色关联度来判断变量之间 关联的强弱，效果较好。本过程运用灰色关联法分析物种与湖泊水位之间的关系。计算原理①求各序列的初值像，即对数据进行标准化。现有数列Xi = [x“l)，x“2)， ，Xi(n)] (i = 1，2， ，1)初值像序列= [x； (1), x: (2), ...,x； ( )J，初值像‘ (k) = xt (k) / Xi ⑴1为系统相关序列的个数，n为各序列的数值。②求初值像的差序列== [6,(1), 6,(2),... , 6 , (n) ] (i = 1, 2, . . . , 1)幻为系统特征序列的初值像。
③计算两极最大差M和最小差mM = max max Sl (k)m = min min (k)
i k i k④求关联系数
,,、 m + iM广側+讽式中，k=1，2，⑤计算关联度
1 n7i
，n;i = 1，2，...，1 ；分辨系数 I G (0,1)
k=l
(4)孢粉_水位_综合指标方程的建立
运用多元回归的方法建立孢粉_水位_综合指标之间的关系Y = f (Z，&，X2，,
Xn)，设定水位为主变量，该部分用非线性拟合，环境综合指标为协变量，该部分用线性拟
1=1 o 多元变量回归
①指数拟合Y = b0+b1Zn+b2X1+b3X2+. . . +bpXp，Y为孢粉组合带数据，X为环境综合指标，Z为水 位。②对数拟合Y = bo+log^+b^^b^^. . . +bpXp，Y为孢粉组合带数据，X为环境综合指标，Z为水 位。(5)方程的验证将湖泊近代沉积物样品的物理化学指标，孢粉组合带数据代到方程中计算湖泊水 位，与湖泊监测水位对比计算误差，一定条件下对方程的变量系数进行调整，提高方程的精 度和准确性。实施例将本文的湖泊水位计算方法应用于中国北方某湖泊，利用21°Pb和137Cs测定沉积 物的年代范围是1830年-2006年。①应用主成分分析法对沉积物的pH、T0C、T0N、盐度、电导率、土壤湿度、燃烧值、黏 粒和沙粒构成作相关性分析，该过程借助SPSS16. 0实现。表1 KM0 和 Bartlett 检验 根据统计学家Kaiser给出的标准，KM0取值大于0. 6，适合因子分析，本例中KM0 为0. 63，符合标准，Bartlett球度检验给出的相伴概率为0. 000，小于显著性水平0. 05，认 为适合于因子分析。根据方差极大法对因子载荷矩阵旋转后，第一因子主要由TOC、TON、LOI、M0S解 释，主要反映土壤中有机质的含量，第二因子由盐度、电导率解释，主要反映了土壤中的无 机盐情况，第三部分主要由黏土、沙土比例和PH解释，反映了土壤所处的水力状况。得到的新因子函数为Fi = 0. 209x^0. 239x2_0. 195x3_0. 018x4_0. 078x5+0. 239x6+0. 2287+0. 076x8+0. 017
x9F2 = -0. OQSXi-O. 016x2_0. 025x3+0. 484x4+0. 516x5_0. 021x6_0. 028x7_0. 100x8+0. 05
7x9F3 = -0. 084x^0. 044x2_0. 028x3_0. 014x4_0. 034x5+0. 141x6_0. 027x7+0. 550x8+0. 53
5x9将F2，F3作为反映沉积物环境特征的综合指标。②孢粉是植物散落的花粉，植物开花期是孢粉形成的主要时期，3-9月是植物生长 和繁殖最旺盛的季节，选择该时间段的月均水位作为研究对象。同时孢粉有两个来源水域 和陆地，陆生植物主要靠风力把花粉吹到水面上，再通过扩散沉降作用进入到沉积物中，受 湖泊水位影响较小。因此，本文主要考虑水生植物孢粉与湖泊水位的关系，从所有检测的孢
6粉中选取人工禾本、自然禾本、香蒲、莎草科、狐尾藻属、蕨类孢子作为研究对象。运用灰色关联法选取与该湖泊水位相关性较大的孢粉作为指示生物。湖泊水位采 用1924年-1974年的记录水位。由于沉积层反映的是年代范围，需要对水位数据进行处理。处理方法第i年的水位
，代表3-9月的平均水位 第p层沉积剖面对应的年数为k年，该沉积层对应的水位为 取分辨系数I =0.8，得到灰色关联计算结果
表2灰色关联计算结果 根据灰色关联的结果，将关联系数大于0. 75的物种作为指示生物，依次为香蒲、 狐尾藻属、蕨类孢子、自然禾本。四种孢粉组成新的物种函数，关联系数作为孢粉变量的系 数，函数式为Y = 0. 975Xl+0. 845x2+0. 792x3+0. 783x4式中Y为不同孢粉组合量，x为各沉积层不同孢粉的含量(丰度表示)③用多元回归的方法建立孢粉_水位_综合指标方程根据F” F2，F3，Y，湖泊水位值X，通过多元回归得到最终方程，本过程通过 Origins. 0实现，二次、三次、指数方程均为自己定义。(1)指数拟合一次方程为Y =22: 72x-31,.95卩「15. 48F2-■21. 85F3-122.66
二次方程为Y =1.32x2-22.62F「3. 74F2-9..04F3-40.,26
三次方程为Y =0.11x3-13.81FA7. 54F2+3..54F3-21.,99
(2)对数拟合
自然对数拟合Y=194. 651nx-41. 82Ff27.59F2-34.89F3-330. 07
5为底对数拟合Y ==313. 281og5x-41. 82Ff-27. 59F2--34.89Fs-33010 为底对数拟合:Y = 448. 211gx_41. 82F「27. 59F2_34. 89Fs-330. 07当用对数拟合时，孢粉组合带结果主要取决于非线性部(水位)，线性部(环境综 合指标)对结果影响不大。④方程的验证将2个近代沉积物样品环境综合指标、孢粉组合带数据代入到方程中，计算湖泊 水位，与监测水位对比表3计算水位与监测水位的对比单位米
经比较，三次方程拟合值与监测值最接近，误差不超过0.3米。当然，方程还需要 更多的数据进行验证，一定条件下需要对变量的系数进行调整。
权利要求
一种利用沉积物孢粉组合带推算历史时期湖泊水位的技术，其特征在于运用主成分分析法将相关环境因子转化为独立的环境综合指标，通过灰色关联分析选取与湖泊水位相关性较大的指示生物，利用回归方程建立指示生物 水位 综合指标间的关系。经近代沉积物指标验证，三次方程拟合值与监测值最接近，水位误差小于0.3m。本技术还可以计算历史时期不同季节湖泊水位的变化，结合湖泊水力参数可以计算不同时期湖泊的面积 。
全文摘要
本发明针对湖泊历史水文数据短缺、不连续等问题，发明了一种以湖泊沉积物中孢粉组合带、典型理化指标为依据，运用主成分分析、灰色关联、回归的数学方法建立孢粉、环境指标与湖泊水位之间关系的技术，可以推算历史时期的湖泊水位，从而了解湖泊生态系统处于低人为干扰度下的物种状况和水文条件，对维持现代湖泊生态系统的健康起着重要作用。
文档编号G06F19/00GK101908104SQ20101027150
公开日2010年12月8日 申请日期2010年9月3日 优先权日2010年9月3日
发明者杨志峰, 郭通, 陈贺 申请人:北京师范大学