На счет Гиббса, Вы не правы, так как там стоит не скорость одной частицы (динамическая переменная), а аргумент (переменная) функции распределения.
Привожу цитату из:
Ландау Л.Д., Лифшиц Е.М. - Теоретическая физика. Том 05 из 10.
Статистическая физика. Часть 1, 2002г, 5-е изд.
Страница 23.
§ 3. Теорема Лиувилля
Вернемся к дальнейшему изучению свойств функции статистического распределения. Предположим, что мы наблюдаем в течение весьма длительного промежутка времени некоторую подсистему. Разделим этот промежуток времени на очень большое (в пределе — бесконечное) количество одинаковых малых интервалов, разделенных моментами времени t1, t2,... В каждый из этих моментов рассматриваемая подсистема изобразится в ее фазовом пространстве точкой (назовем эти точки А1, А2, А3...). Совокупность полученных точек распределится в фазовом пространстве с плотностью, в пределе пропорциональной в каждом данном месте значению функции распределения p(p,q), по самому смыслу последней, как определяющей вероятности различных состояний подсистемы.
Вместо того чтобы рассматривать точки, изображающие состояния одной подсистемы в различные моменты времени t1,t2,..., можно формальным образом ввести в рассмотрение одновременно очень большое (в пределе — бесконечное) число совершенно одинаковым образом устроенных
подсистем 1) , находящихся в некоторый момент времени (скажем, t = 0) в состояниях, изображающихся точками А1, А2,...
1) Такую воображаемую совокупность одинаковых систем обычно называют статистическим ансамблем.
Так что необязательно статистические законы рассматривать как функции распределения — можно и как динамические.