基于协同区域化变量的非饱和带平均粒径估算.docxVIP

基于协同区域化变量的非饱和带平均粒径估算 1 非饱和带岩性的空间变异特征及其定量测量 地质统计学是研究区域变量的空间变异性的有力工具。Aboufirassi M和Martin M A、Shakeel Ahmed和Ghislain de Marsily用协克立格法最优估计了含水层导水系数;Clifton Peter M和Neuman Shlomo P用克立格法和反演模型对区域含水层进行了条件模拟;Neuman Shlomo P将地质统计学用于地下水水文学的反演问题;朝伦巴根、刘廷玺等用泛克立格模型研究了区域含水层渗透系数空间变异性、用时间因子克立格模型、多时段泛克立格模型和协克立格模型研究了水文随机变量插值与预报。 非饱和带岩性直接影响土壤各种理化特性及水分保持和运移情况,因此在大气降雨、地表水、土壤水、地下水(四水)在非饱和带相互转化关系中的起到媒介和传播作用。因而,充分地认识非饱和带岩性及其空间变异性对揭示非饱和带中四水转换机理具有重要意义。 非饱和带岩性通常用土壤的粒径级配来表示。获得土壤级配数据,要实地采样,通过筛分试验做出粒径级配曲线。在较大尺度研究区通过采样-筛分-做粒径级配曲线的方法研究非饱和带岩性及其空间变异性能获得理想结果,但耗费的人力物力相当可观。为此,我们在研究区通过踏勘选一定数量的采样点,并按采样-筛分-做粒径级配曲线获得岩性平均粒径信息。与此同时有经验的野外工作人员又在这些采样点进行了干测、湿测法并用的岩性定性测定工作。对采样点上同时进行的颗分实验和干测、湿测法并用的岩性定性测定法所得平均粒径数据进行对比、分析和统计计算,证明颗分实验所得非饱和带岩性平均粒径和定性测定法得到的平均粒径可作为协同区域化变量。在此基础上有经验的野外工作人员,在研究区以按一定比例尺设定的定性采样点上用采样钻分层采样,通过干测、湿测并用的定性测定法判定出各层岩性,依据国际制土壤质地分类标准计算出定性测定法平均粒径。研究区非饱和带岩性平均粒径的空间变异性就利用定性测定法得到的平均粒径数据和与它有协同区域化关系的颗分试验法得到的平均粒径数据用协同克立格法进行了分析与评价,不仅取得满意结果,而且也提供了一种省力省钱的方法。 2 非饱和带岩性网络结构 研究区地处西辽河冲淤积平原中部(图1),面积为1 137.0 km2,该区地形开阔平缓,地面坡降为7/10000～8/10000。 本区非饱和带岩性主要受近代河流冲淤积及风蚀、剥蚀影响,地层结构在垂向上自下而上颗粒由粗变细;水平方向上自北向南东,土壤质地由砂质土逐渐向黏质土转变。 在1 137 km2的研究区,以3 km×3 km的正方形网格布设了147个非饱和带岩性采样点(图1),其中的44个测点上开展了两种试验,并计算了非饱和带岩性平均粒径:(1)颗分试验。分层取样、筛分、绘制粒径级配曲线、计算非饱和带岩性平均粒径,(2)干-湿法定性判定土壤质地,根据国际制土壤质地分类标准计算非饱和带岩性平均粒径。其余103个采样点分层取土样,用干-湿法定性测定法判定了土壤质地,计算了非饱和带岩性平均粒径。 2.1 采用试验法和定性测量法获得的平均粒径，构成区域化变量的横向分析 由于协同区域化变量必须在空间上相关,因此利用44个采样点颗分试验法、定性测定法所得平均粒径进行了统计分析。 2.1.1 davdsyj—颗分试验法平均粒径的计算 令DAVGSYjiSYji为通过颗分试验得到的j采样点处i层中各粒组的平均粒径(mm),DSΤjSTj为颗分试验的j采样点处非饱和带岩性平均粒径。平均粒径的计算式如式(1)、式(2)。 DAVGSYji=1LL∑l=1ˉdjilΡSΤjil;J=1,2,…,44;i=1,2,…,Ij(1) DSYj=Ιi∑i=1hjiDAVGSYjiΙj∑i=1hji;j=1,2,…,44 (2) 式中:DAVGSYji——通过颗分试验得到的j采样点处i层中各粒组的平均粒径(mm);DSYj——颗分试验的j采样点处非饱和带岩性平均粒径;ˉdjil——j点i层l种粒径级配的粒级;PSYjil——j采样点处i层l种粒径级配;hji——j点处i层岩性厚度;L——j采样点处i层的粒组数;Ij——j采样点处非饱和带岩性分层数。 2.1.2 岩性层数计算 令定性测定的k种岩性(k=1表示粗砂,k=2表示中粗砂,…)的统计平均粒径为DRΜGk,那么 DRΜGk=1Νk∑Yji=kDAVGSYji;?k(3) 式中:Nk——所有定性测定法采样点中判定为k种岩性的层数;Yji——定性测定法测定的j处i层岩性序号。事实上DRΜGk等于定性测定法在所有采样点上判定为k种岩性层颗分试验平均粒径DAVGSYji的平均值。 经定性测定的n点处i层土壤质地的平均粒径为DAVGΜGni,当层序号与岩性种