摘要：亨利·达西（Henry Darcy，1803‑1858）是19世纪法国工程水力学的代表人物，也是定量水文地质学与现代多孔介质科学的奠基人。长期以来学界多聚焦于1856年提出的达西定律公式本身，却容易忽视该理论背后深厚的工程实践背景与知识演化脉络。达西的水力学思想根植于第戎城市供水工程建设，在水源勘察、管道输水、砂滤净水、井泉调查等工程问题驱动下逐步成型。1855‑1856年砂柱渗流实验将多孔介质内部水流转化为可控、可测量的物理过程，得到流量、水头差、介质特征之间的定量关系，即达西定律。该定律并非普适性物理定理，而是建立在饱和、单相、低速层流、连续介质假设之上的宏观本构关系。后世学者将达西定律与质量守恒原理结合，先后发展出Dupuit井流理论、Boussinesq自由面理论、Thiem稳定井流以及Theis非稳定井流理论，构建起地下水动力学核心框架。面对非达西流、裂隙‑岩溶介质、非饱和渗流、多相流以及热‑水‑力‑化（THMC）多场耦合等复杂条件，研究者不断界定其适用边界，发展各类修正与拓展模型。时至今日，达西的科学范式依旧支撑地下水数值模拟、地下资源开发、污染场地治理等领域，并开始与人工智能、物理约束数据模型相融合。达西最重要的科学遗产并非单一公式，而是建立起“工程问题‑受控实验‑参数化‑数学描述‑工程预测”的研究范式，为整个多孔介质科学开辟发展道路。

关键词：亨利·达西；第戎城市供水；达西定律；地下水动力学；多孔介质；定量水文地质学

1 引言

19世纪欧洲城市化与工业化快速推进，城市供水、矿山排水、隧道与铁路建设带来大量水工程难题，地下水从单纯的自然地质现象转变为需要勘察、计量、调控的工程对象。在这一时代背景下，工程师需要回答地下水补给来源、井的出水能力、地层透水阻力、砂滤层过水能力、管道水头损失等一系列现实问题。

达西因达西定律广为人知，1856年他在著作《第戎城市公共喷泉》（Les fontaines publiques de la ville de Dijon）附录公布砂柱实验结果，记录水在砂质多孔介质中的运动规律。但这部长达647页的著作是综合性工程文献，涵盖水源、自流井、含水层、管道、过滤工艺等诸多内容，渗流实验只是工程研究的衍生成果，并非专门的纯理论研究报告。现代教科书给出的达西定律矢量形式是后世整理推广得到，并非达西原始手稿中的数学表达。

以往部分研究将达西简单理解为“公式发现者”，割裂理论与其工程起源。本文结合科学史文献，沿着城市供水工程实践→地下水与含水层认知→管道水力学与砂滤技术→砂柱实验→达西定律→定量水文地质学→现代多孔介质科学这一历史链条，梳理达西水力学思想形成过程，剖析达西定律的科学突破、适用局限、后世理论谱系及其当代价值，重新认识达西在水文地质学史中的地位。

2 达西定律诞生的知识背景

达西定律的出现不是孤立的科学突破，而是建立在前人多领域知识积累之上，主要包含地下水工程经验、矿山排水实践、泉水调查传统以及泊肃叶毛细管黏性流动研究四大知识源流。

2.1 井泉与矿山排水的经验积累

人类很早就通过井、泉认识地下水，但早期认知以自然哲学与野外经验为主，缺少定量描述手段。随着欧洲采矿业发展，矿井持续涌水迫使工程师关注地下水来源、岩层透水性、水压力、坑道涌水量等实际问题，积累大量关于含水层、裂隙岩体的感性认识，但是始终没有形成描述地下水运动的定量规律。

2.2 法国泉水调查的经验传统

19世纪法国形成以Paramelle为代表的地下水调查传统，代表著作《发现泉水的艺术》尝试将地形、地质条件和泉水出露、水井取水能力建立联系。虽然这套体系属于经验归纳，缺少严密数学推导，但已经建立“地质‑地形‑地下水”的分析思路。达西在自己著作中专门讨论自流井成因，引用Paramelle的研究成果，说明这套经验知识已经被纳入当时工程界的知识体系之中。

2.3 泊肃叶毛细管黏性流动研究

泊肃叶开展毛细管内低速黏性流体实验，证明低速条件下流量与驱动压差呈线性关系，区别于高速流动下流速平方相关的惯性规律。达西注意到砂层渗流与毛细管黏性流动的相似性，砂粒之间复杂孔隙网络可以等效为大量曲折毛细管集合。泊肃叶的研究为达西理解砂柱实验结果提供重要的流体力学参照，达西定律可以看作管流、毛细管黏性理论向颗粒多孔介质的拓展与综合。

3 第戎供水工程：达西水力学思想的工程起源

3.1 工程背景

达西毕业于法国桥路工程学校，职业身份是桥路工程师，全部科学探索都服务于工程需求。19世纪第戎面临供水水量不足、水质恶化，霍乱等水媒传染病频繁暴发。达西主持第戎城市重力供水工程，将远处Rosoir泉水通过长距离暗渠引入城内，完成蓄水、配水系统建设。开展工程期间，他必须系统解决水源评价、泉水径流、长距离输水水头损失、管道阻力、砂滤净水等一系列问题，完整接触从地下水水源到城市末端供水的全链条水问题。

3.2 研究视角的转变：从“哪里有水”到“水如何流动”

此前地下水相关研究大多聚焦找水，即解决“水在哪里”。达西在工程实践中将问题向前推进：水为何赋存于地层之中？水如何穿过砂质地层？水流经受多大阻力？如何定量计算流量、预测供水能力？研究重心由静态资源勘察转向水流运动机理，这是地下水由经验描述走向定量科学的重要思想转折。

同时，砂滤净水工艺成为连接工程实践与地下水科学的关键桥梁：人工滤池内水穿过砂层，与天然含水层中地下水在颗粒介质内运动具备物理类比关系，人工过滤与天然地下渗流被联系在一起，为后续砂柱实验提出现实命题。

4 砂柱实验与达西定律及其科学突破

4.1 砂柱实验过程

达西的渗流思想并非一蹴而就，早在1830年Saint‑Michel井勘察时，他已经意识到含水层内部存在不可忽略的水流阻力；19世纪40年代积累大量自流井观测资料，观察到井的出流量与水位降深之间存在近似线性关系。经过长期问题积累，1855‑1856年开展系统砂柱实验，共计35组试验，使用索恩河石英砂，改变砂层厚度、压实程度，通过测压计控制水头差，计量渗流流量，完成控制变量实验研究。

实验得到核心结论：渗流流量与过水断面面积成正比、与上下游水头差成正比，与渗流路径长度成反比，得到原始实验关系：Q=KA\frac{\Delta h}{L}，后世改写为通量形式q=-K\nabla h，即现代通用的达西定律，负号代表水流沿水头降低方向运动。需要强调：达西本人只定义水力传导系数（渗透系数K）；表征介质本身属性的渗透率、“达西”单位，是后世石油行业研究者进一步提出的概念。

4.2 达西定律的科学意义

⑴实现地下水从定性观察向定量科学跨越。在达西之前，可以观测泉、井、自流井，但无法定量预测水头变化如何改变流量。砂柱实验实现现象可重复、条件可控制、结果可计算，地下水运动成为可参数化的物理过程。

⑵将地质介质纳入水力学体系。传统水力学主要研究管道内水流，达西把砂质地质介质作为水流系统的组成部分，揭示地层本身会产生渗流阻力，奠定含水层水力学的思想基础。

⑶打通水处理工程与地下水科学，证明人工砂滤池与天然含水层可以使用同一套宏观渗流规律描述，实现工程技术与地质科学的知识互通。

5 达西定律推动下定量水文地质学的理论演进

达西一维砂柱实验得到的是简单一维渗流关系，现代地下水动力学不是达西单独完成，而是几代学者将达西本构关系与质量守恒定律结合，从一维拓展到二维、三维，由稳定流拓展至非稳定流，逐步构建完整学科体系。

Dupuit与达西处于同一工程体系，率先将达西定律应用于抽水井径向汇流问题，建立Dupuit井流模型，实现实验室一维砂柱向野外真实含水层井流问题的跨越。之后Boussinesq建立潜水自由面地下水运动理论；Thiem完善稳定井流公式，奠定抽水试验求参的理论基础；1935年Theis提出非稳定井流理论，引入含水层储释水效应，描述水位随抽水时间动态变化规律，地下水动力学的理论框架正式成型。

整个理论发展链条为：达西一维渗流本构→Dupuit井流→Boussinesq潜水运动→Thiem稳定井流→Theis非稳定井流理论。在此基础之上，研究者继续发展溶质运移理论，解决地下水带动污染物迁移问题，形成环境水文地质学；后续数值计算方法出现，有限差分、有限元算法将上述理论转化为MODFLOW等数值模拟工具，广泛用于地下水资源评价、防渗计算、海水入侵、地面沉降模拟。

6 达西定律的适用边界与现代拓展

达西定律是宏观统计性本构关系，建立在代表性单元体REV尺度之上，存在严格适用条件：低速层流、单相饱和流体、介质结构稳定连续；当现实条件偏离假设，经典达西定律失效，需要开展理论修正与拓展。

6.1 高速非达西流

砾石、岩溶管道、宽大裂隙中渗流流速高，惯性力不可忽略，水头损失和流量不再保持线性，采用Forchheimer方程J=aq+bq^2，描述非线性渗流行为。

6.2 介质非均质与各向异性

天然含水层存在岩性变化、透镜体、断层、层理，渗透系数具备空间分布特征(K=K(x,y,z))；层理沉积岩体表现各向异性，渗透系数由标量升级为张量形式。同时渗透系数存在显著尺度效应，室内小试样测试参数不能直接等价于野外区域含水层等效参数，由此发展出随机场、均质化、有效介质等理论。

6.3 裂隙‑岩溶复杂介质

花岗岩裂隙、石灰岩岩溶管道系统，水流沿着优势通道快速运移，简单连续孔隙介质假设不再适用。学界发展双孔隙模型、双重渗透率模型、离散裂隙网络DFN、岩溶管道‑基质耦合模型。达西思想并未被抛弃，而是作为基质部分的控制方程嵌入更加复杂介质框架中。

6.4 非饱和渗流

包气带孔隙同时存在水与空气，水力传导系数随含水率变化，以达西通量框架结合质量守恒得到Richards方程，描述土壤‑地下水连续体系水分运移，广泛用于入渗、蒸发、农田排水计算。

6.5 多相渗流与多物理场耦合

油气开发、CO₂地质封存场景存在水‑油‑气多相流体，引入绝对渗透率与相对渗透率概念，建立多相达西方程。现代多孔介质科学进一步考虑温度场、应力场、化学反应，开展热‑水‑力‑化学THMC耦合模拟，用来解析地面沉降、地热开发、地质碳封存等复杂工程问题。

6.6 达西通量与真实孔隙流速区分

达西通量是按整个断面（固体颗粒+孔隙）计算的表观流速，不等于地下水在孔隙通道内实际平均流速，真实孔隙流速v=q/n_e，n_e为有效孔隙度。该区分对于地下水示踪、污染物运移、地下水年龄计算至关重要。

7 达西思想的当代价值与科学遗产

进入21世纪，传感器网络、遥感、地球物理勘探带来海量地下水监测数据，人工智能被引入地下水领域，用于参数反演、水位预测。但纯数据驱动模型容易出现物理不自洽问题，越来越多研究采用物理约束‑数据融合建模思路，把达西定律、质量守恒作为物理约束嵌入机器学习模型，实现传统物理理论与现代数据科学结合。

达西留给后世最重要的遗产，并不局限于达西定律这一公式，而是一套完整的科学范式：从真实工程中提炼科学问题，通过受控实验获取定量规律，借助数学完成参数化表达，再回归自然系统实现预测与工程管控。

达西定律不是多孔介质科学的终点，而是持续生长的理论起点。它本身存在明确适用边界，但具备极强可拓展性：研究者不去推翻原始定律，而是不断界定边界条件，针对复杂地质、流体、环境场景发展修正模型，把一维饱和砂柱的实验规律，拓展到多尺度、多相、多场耦合的复杂现实世界。

8 结论

亨利·达西作为一名法国桥路工程师，其水力学思想诞生于第戎城市供水工程。水源勘察、长距离输水、砂滤净水等工程需求推动他开展砂柱渗流实验，1856年得到描述饱和砂层水流的定量关系，即后世所称达西定律。该成果只是工程著作的附录，并非专门纯理论研究产物。

达西定律是建立在REV尺度之上的宏观本构关系，只适用于低速层流、单相饱和连续介质。Dupuit、Boussinesq、Thiem、Theis等学者将达西本构关系和质量守恒结合，完成从一维砂柱向野外含水层、从稳定流向非稳定流的拓展，奠定现代地下水动力学与定量水文地质学的理论根基。面对高速非达西流、裂隙岩溶介质、非饱和带、多相流体、热‑水‑力‑化耦合等复杂条件，学术界不断界定其适用边界，发展各类拓展模型。

历经近170年，达西的研究范式仍然支撑地下水数值模拟、地下资源开发、污染修复等诸多领域，并且与人工智能、物理信息模型相结合继续向前发展。达西的核心科学遗产，是建立“工程实际问题‑受控实验‑定量参数‑数学描述‑工程预测”的完整研究链条，为整个现代多孔介质科学开辟道路。

参考文献

（略）

支持单位：