一种基于区域尺度的方法：利用流域尺度有效值识别浅层含水层岩性特异性的水力传导系数分布

《Journal of Hydrology》：A regionally based method to identify lithology-specific hydraulic conductivity distributions in shallow aquifers using catchment-scale effective values

　　

本研究提出了一种基于流域尺度有效值估算岩性特异性水力传导系数的新方法，将每个有效值解释为底层各岩性水力传导系数的线性组合。该方法特别适用于由径流量推导的有效水力传导系数估算，其中河流网络整合了来自相邻丘陵坡以浅地下水流的综合影响。研究人员利用岩性面积分数作为权重构建超定线性系统，并在假设水力传导系数服从对数正态分布以考虑空间变异性的前提下，采用改进的矩量法求解。该方法在合成案例中得到了验证，并应用于阿莫里坎地块（Armorican Massif，法国）的113个流域。研究表明该方法具有良好的收敛特性和高预测精度：在85%的情况下，模拟水力传导系数落在观测值的90%置信区间内。收敛分析表明，每个岩性使用5个流域即可达到75%的预测成功率，均值和标准差参数的估算误差分别控制在6%和30%以内。当流域数量增至80个时，预测精度提高至80%，参数估算误差降至3%和15%。该方法计算效率高、可推广性强，为区域尺度岩性特异性水力传导系数数据集的探索提供了有力工具，也可 adapted 于孔隙度等其他适用线性平均的变量，从而在水文地质建模中获得更广泛的应用。

该论文发表于《Journal of Hydrology》，旨在解决浅层含水层中岩性特异性水力传导系数难以获取的关键科学问题。水力传导系数（hydraulic conductivity, K）是描述地下水流动与能量梯度关系的核心参数，但其获取面临显著的尺度效应挑战。局部测量值（厘米至百米尺度）常受采样偏差影响，尤其在裂隙介质中可能因过度采样高渗透区或遗漏低渗透区而产生系统性偏差。全球数据库（如GLHYMPS）虽提供了大尺度参考数据，但难以反映区域变异性和特定岩性的局部特征。河流径流分析为千米尺度水力参数表征提供了补充途径，但既往研究仅能获得代表多种岩性综合贡献的有效水力传导系数（effective hydraulic conductivity, K_eff），无法分解至单一岩性。针对这一瓶颈，研究人员开发了水力传导系数分解模型（Hydraulic Conductivity Decomposition Model, HCDM），利用HS1D（one-equivalent hillslope model）丘陵坡水文模型从径流时间序列中提取流域尺度有效水力传导系数，进而反演各岩性类别的水力传导系数概率分布。

研究采用的核心技术方法包括：（1）基于超定线性系统的矩量法优化求解，将流域有效水力传导系数表示为各岩性贡献的加权线性组合，权重由岩性面积分数确定；（2）假设岩性水力传导系数服从对数正态分布，通过最小二乘估计求解均值参数μ_l和方差参数σ_l²；（3）合成数据集的蒙特卡洛验证与留一法（leave-one-out）交叉验证；（4）在法国阿莫里坎地块113个实测流域的应用验证，HS1D模型校准采用SURFEX水文气象模型的入渗补给时间序列，强迫数据为SAFRAN再分析资料。

**方法开发与验证**

**超定线性系统构建与求解**：HCDM将第i个流域的有效水力传导系数K_i建模为m种岩性水力传导系数的线性组合，即K_i = Σ_l=1^mα_ilK_l，其中α_il为岩性l在流域i中的面积分数，满足Σα_il=1。假设各岩性的水力传导系数K_l为独立对数正态分布随机变量，即log₁₀(K_l)~N(μ_l, σ_l<2)。通过n个流域（n>>m）构建超定线性系统，首先求解期望值的优化问题获得μ_l的估计，进而通过偏差分析估计模型诱导变异，最终求解方差σ_l²。优化问题在MATLAB中使用quadprog函数求解，单次运算耗时仅0.1-0.4秒。

**合成案例验证**：通过大量蒙特卡洛模拟验证HCDM的收敛性和预测能力。结果表明，参数相对误差随样本量增加呈幂律衰减（均值误差∝n^-0.87，标准差误差∝n^-0.86）。对于未测流域的预测，当每个岩性配置约5个流域（即n≈5m）时即可达到75%的预测成功率；当n增至50-100时，预测精度快速提升至75%-83%，此后继续增加样本量收益递减。岩性数量m对预测精度的影响相对较弱，当m从2增至10时，相同条件下的预测精度仅从83%微降至77%。

**阿莫里坎地块应用**

**地质与水文地质背景**：研究区为法国西北部阿莫里坎地块的前寒武纪-古生代结晶岩基底，主要由变质岩和深成岩组成。基于1:1,000,000地质图，将109种岩性归并为5个主导岩性类别：Brioverian页岩（占流域出现率100%，下同）、原生页岩（74%）、古生代深成岩（92%）、元古代深成岩（71%）和变质岩（100%）。从DREAL管理的HydroPortail数据库选取113个水文测站，采用25米分辨率BD ALTI数字高程模型进行流域划分，确保流域互不重叠。

**岩性水力传导系数估算**：HS1D模型校准获得113个流域的有效水力传导系数，进而应用HCDM求解。结果显示，不同岩性间水力传导系数均值存在显著差异（跨度达1.68个数量级）：Brioverian页岩最低（log₁₀(K)=-7.68 m/s），原生页岩次低（-6.92），古生代深成岩（-6.29）、元古代深成岩（-6.19）和变质岩（-6.01）较高。标准差σ_l普遍小于0.5个数量级，表明岩性分类对有效水力传导系数具有良好的约束作用。值得注意的是，基于岩石年龄的细分（元古代vs古生代深成岩）并未提供额外信息量，两类深成岩的水力传导系数分布高度相似，说明广义的岩性分类而非年龄或风化程度是更主要的控制因素。Brioverian页岩的水力传导系数显著低于原生页岩，这一反直觉结果可能源于其强烈风化导致黏土充填初始裂隙，而非增强渗透性。

**模型预测能力评估**：将HCDM估算的岩性水力传导系数用于重新计算113个流域的有效值，并与HS1D实测值比较。在96个流域（85%）中，实测值落在HCDM估计的90%置信区间内；剔除17个异常值后，预测值与实测值的相对误差从0.71降至0.51。随样本量增加，参数估计误差逐步降低：n=20时均值和标准差误差分别为7%和30%；n≥60后显著改善，n=80时预测精度达80%，此后边际增益有限。

**讨论与展望**

**方法适用性与推广**：HCDM适用于浅层、无压、地形驱动局部流的含水层系统，如结晶岩地区。类似地质背景还包括伊比利亚地块、波希米亚地块、斯堪的纳维亚地盾等。该方法特别适用于需要流域尺度有效水力传导系数的水文模拟、水资源评估和低流量分析，而非精细的钻孔尺度表征。与全球数据库（如GLHYMPS）相比，HCDM提供的是区域约束更严、变异性更窄（约2个数量级）的估算结果，更能反映特定地区的水文地质结构。

**岩性分类策略**：当前五类分类方案基于专家经验，兼顾地质特征与水文行为。虽然数据集信息量可能支持更细致的分类，但增加岩性类别需满足两个条件：（1）子类别须基于可识别的流域特征（如年龄、风化程度、矿物组成或裂隙密度）；（2）预测性能须改善或至少不下降。Brioverian页岩的高变异性提示未来可考虑将其进一步细分。

**线性组合假设的拓展**：HCDM假设有效水力传导系数为各岩性贡献的算术加权平均（对应平行于流向的岩性排列）。若岩性排列垂直于主流向，则需采用调和平均。通过对数转换X_l=1/K_l可实现线性化，但不同分布假设下均值与方差的转换关系各异，需针对具体分布形式推导。

**不确定性集成与参数拓展**：当前假设各流域仅具单一有效值，但野外估算通常伴随不确定性。将测量误差纳入模型可扩展至质量参差不齐的数据集，如结合Brutsaert法等方法从大量测站提取参数。此外，HCDM框架可拓展至孔隙度（porosity）、水力扩散系数（hydraulic diffusivity）等同样适用线性平均的其他水文地质参数，前提是有效参数可表示为各岩性参数的线性组合。

**研究结论**：水力传导系数分解模型（HCDM）是一种计算高效、可推广的方法，能够从流域尺度有效水力传导系数测量中反演岩性特异性的水力传导系数分布。该方法基于以下核心假设：有效水力传导呼和浩特市模型的变异主要来自不同流域间岩性组成的差异。合成验证表明参数估算误差随样本量增加呈幂律收敛，预测精度在50-100个流域区间快速提升至平台期。在阿莫里坎地块的应用揭示了岩性对水力传导系数的主导显控作用，五类岩性的水力传导系数差异达1.68个数量级，其中变质岩和深成岩显著高于页岩。模型对113个流域的预测中，85%的实测值落入90%置信区间，证实了方法的稳健性和区域尺度水文地质表征的可靠性。