We investigate an energy-based formulation of the two-field poroelasticity model and related multiple-network as they appear in geosciences or medical applications. propose a port-Hamiltonian system equations, which is beneficial for preserving important properties after discretization model-order reduction. For this, we include commonly omitted second-order term consider corresponding first-or...