保定平原区地下水水位变幅带自由孔隙率空间变异研究  

Spatial variability of free porosity in the groundwater level fluctuation zone in the Baoding Plain area

作　　者：徐步云 杨会峰[1,2] 白华 宋博 孟瑞芳[1,2] XU Buyun;YANG Huifeng;BAI Hua;SONG Bo;MENG Ruifang(Institute of Hydrogeology and Environmental Geology,Chinese Academy of Geological Sciences,Shijiazhuang,Hebei050061,China;Hebei Cangzhou Groundwater and Land Subsidence National Observation and Research Station,Shijiazhuang,Hebei050061,China)

机构地区：[1]中国地质科学院水文地质环境地质研究所,河北石家庄050061 [2]河北沧州平原区地下水与地面沉降国家野外科学观测研究站,河北石家庄050061

出　　处：《水文地质工程地质》2023年第3期23-33,共11页Hydrogeology & Engineering Geology

基　　金：国家自然科学基金地质联合基金项目(U2244214);河北省创新能力提升计划高水平人才团队建设专项(225A4204D);中国地质调查局地质调查项目(DD20221752);中国地质科学院科研业务费项目(SK202118)。

摘　　要：科学确定地下水水位变幅带的地学参数是水资源评价与管理中的重要环节。保定平原近40年来强烈开采地下水,水位持续下降形成规模巨大的厚包气带层,南水北调工程通水后,随着河湖生态补水与地下水压采工作的推进,保定平原局部地区地下水水位止跌回升。回补水量与水位变化的定量关系成为超采治理的一个重要科学问题,但水位回升条件下的计算中仍用表示释水过程的给水度参数将造成结果偏差,因此表述水位恢复过程的自由孔隙率参数研究是解决这一问题的重点所在。文章基于保定平原区67个工程地质钻孔采样数据,确定水位变幅带的综合自由孔隙率,通过趋势分析和结果交叉验证方法遴选半变异函数模型,结合普通克里金插值方法对未知点进行无偏最优估计。结果显示:(1)保定平原区水位变幅带综合自由孔隙率的最优半变异函数模型为一阶趋势效应指数模型。数据具有强空间自相关性,主要受变幅带空间位置、地层类型等结构性因素影响。(2)综合自由孔隙率分布表现为西南、西北为高值区,极值可达0.25,数值向中部及东部逐步降低,最小降至0.02。(3)与惯用给水度值进行对比,自由孔隙率值在南北部地区整体升高约0.03,约为惯用给水度的1.2倍。在中部地区降低了约0.06,变为惯用给水度值的一半左右。研究成果对南水北调受水区生态补水与水资源调控具有重要研究价值和借鉴意义。Scientific determination of geological parameters of the groundwater level fluctuation zone is an important step for water resources evaluation and management.In the past 40 years,over-exploitation of groundwater in the Baoding Plain led to a serious deficit of aquifers,thus forming a huge deep vadose zone.After the South-to-North Water Diversion Project was put into operation,with the progress of ecological water supplement of rivers and lakes and groundwater limited-over-exploitation,the groundwater levels in some areas of the Baoding Plain stopped falling and rose.The quantitative relationship between the amount of water supplement and the change of groundwater levels has become an important scientific issue in the management of over-exploitation of groundwater.However,under the condition of groundwater level recovery,the calculation results will be deviated if the specific yield parameter of water release process is still used.Therefore,the free porosity parameter study of groundwater level recovery process is the key to solve this problem.Based on 67 engineering geological boreholes in the Baoding Plain,the comprehensive free porosity during groundwater level recovery is determined according to the lithologic characteristic parameters of the groundwater level fluctuation zone.The semi-variogram model is selected by the trend analysis and cross validation,and the spatial unknown points are interpolated by the Ordinary Kriging interpolation.The results show that(1)the best semi-variogram model of the comprehensive free porosity in the groundwater level fluctuation zone of the Baoding Plain is the 1 order index model.The spatial autocorrelation of the data is obvious,which is mainly affected by structural factors such as spatial location of the groundwater level fluctuation zone and stratigraphic type.(2)The comprehensive free porosity distribution shows that the southwest and northwest are the high value areas,and the extreme value can reach 0.25.The parameter gradually decreases to the central and eastern reg

关 键 词：自由孔隙率 给水度 非饱和带 地质统计 空间变异 保定平原区 

分 类 号：P641.2[天文地球—地质矿产勘探]