An exemplary embodiment of the invention relates a method for determining hydraulic permeability of rocks in a subsurface region, the method comprising: in-situ measuring and/or calculating at least one of independent seismic velocities of rocks at different locations in said subsurface region; determining at least one lithological unit in said subsurface region based on the measured seismic velocities; for the at least one lithological unit, acquiring at least one rock sample at an at least one location of said subsurface region; measuring the rock-sample porosity and permeability as functions of stress; measuring and/or calculating of at least one of independent seismic velocities of said at least one rock sample as functions of stress; computing the stiff and compliant porosity of said at least one rock sample; computing numerical coefficients of a given analytical permeability model based on the stiff and compliant porosities; computing coefficients of a given analytical model of an elastic-compliance characteristic of the rock based on the stiff and compliant porosity; computing the stiff and compliant porosity for a plurality of other locations in said subsurface region; and computing the permeability of rocks belonging to said at least one lithological unit, for said at least one location in said subsurface region and for said plurality of other locations in said subsurface region.