Laboratory-based mathematical modelling of fines migration allows predicting well productivity reduction during the geothermal reservoir exploitation. The analytical model for one-dimensional flow with ionic strength alteration has been derived. Good adjustment of the permeability and breakthrough concentration data from coreflood test by the analytical model has been achieved, and the tuned model coefficients fall in the common ranges. The obtained maximum retention function of multi-sized fines allows calculating their size distribution. During the temperature rise, weakening of electrostatic attraction on fines attached to rock surface overwhelms the reduction of detaching drag force due to water viscosity decrease. It leads to increased fines detachment and more severe permeability decline at elevated temperatures, typical for geothermal fields. (C) 2015 Elsevier Ltd. All rights reserved.