Почва — не сплошная губка, а лабиринт из пор разного калибра. После дождя вода заполняет их не все сразу: сначала самые узкие, потом те, что пошире. Физики построили теорию этого процесса, взяв за основу простой образ — множество сообщающихся сосудов. Перетоки между порами подчиняются законам течения вязкой жидкости, а движущей силой служит разность химических потенциалов. Грубо говоря, это «желание» воды занять ту или иную пору с учётом капиллярных сил и тяжести.
Ключевым оказалось безразмерное число Дамкёлера для каждой поры. Когда оно мало, вода успевает равномерно распределиться — это классическое медленное впитывание, которое описывает уравнение Ричардса. Когда же число велико, вода проскакивает по избранным путям, не успевая «осмотреться». Так возникают быстрые протоки, знакомые каждому, кто видел, как после ливня лужа исчезает за минуты, а не за часы. Теория показывает, что это не два разных явления, а плавный переход, управляемый одним параметром.
Такой взгляд делает ненужными десятки частных моделей. Вся необходимая информация о почве — геометрия поровой сети, которую можно измерить на микротомографе. Не нужно подгонять параметры под полевые данные. Уравнение само решает, где воде течь медленно, а где — быстро.
как её написали авторы, на английском arXiv:2607.09416
We develop a statistical-mechanical theory of water in unsaturated soil whose outcome is a continuum field equation for the pore-occupancy g(r,x,t), the fraction of pores of radius r that are water-filled at position x and time t. The theory is built across three scales: microscopic inter-pore transfers set by Hagen-Poiseuille rates and a driving potential (the difference of pore-class chemical potentials, taken in capillary-gravitational form but open to adsorptive, osmotic, or thermal refinement); a mesoscale master equation relaxing the occupancy toward the equilibrium step g_eq=H(r*-r); and, on contracting the averaging volume to a point, the continuum balance d_t g + div F = C[g] - E - T, of which everything else is a limit, a moment, or a boundary resolution. The kinetic equation is an Onsager gradient flow descending the Gibbs free energy, with an H-theorem for the isothermal unforced system and mass conservation as its zeroth moment. A single dimensionless group, the pore-resolved Damkohler number Da(r,x), organizes the behavior and unifies phenomenologies long modelled separately. A Chapman-Enskog reduction identifies Richards' equation as the quasi-static (Da->0) limit, with matric potential and hydraulic conductivity K emerging only there and K vanishing below the percolation threshold; capillary-bundle and critical-path models are its diagonal and spectral limits. Hysteresis is the holonomy of a forcing bundle, a geometric phase rather than per-pore bistability, with a falsifiable loop-area law H ~ I^2. Preferential flow is what the same equation does where Da>1, so the Richards/preferential-flow dichotomy becomes a continuous Da-controlled crossover. Out of the quasi-static limit g(r) is the irreducible state variable. All inputs are geometric properties of the pore network, measurable from micro-CT and calibrated against no macroscopic data.
Число Дамкёлера, которое управляет поведением воды в почве, изначально пришло из химической инженерии и описывало, насколько быстро химическая реакция успевает пройти по сравнению с потоком вещества. Теперь оно «работает» ещё и под ногами.