Determination of the throat size distribution of a porous medium as an inverse optimization problem combining pore network modeling and genetic and hill climbing algorithms

The pore size distribution of a porous medium is often estimated from the retention curve or the invading fluid flow rate curve using simple relationships more or less explicitly based on the consideration that the porous medium is made of a bundle of cylindrical parallel tubes. This type of determination is tested using pore network simulations. Starting from two- or three-dimensional networks, the characteristics of which are known a priori, the estimation of the throat size distribution (TSD) is performed using the standard methods in the case of drainage. Results show a significant discrepancy with the input data. The disagreement is more pronounced when the fluid flow rate curve is employed together with the parallel tubes assumption. The physical origins of these shortcomings are identified. A method, based on pore network simulations combined with a genetic algorithm and the hill climbing algorithm, is then designed, which makes simultaneous use of the nonwetting fluid flow rate curve and the retention curve of the medium. Very significant improvement is achieved in the estimation of the TSD using this procedure.