Simulation of pressure-induced fracture in two-dimensional (2D) and three-dimensional (3D) fully saturated porous media is presented together with some peculiar features. A cohesive fracture model is adopted together with a discrete crack and without predetermined fracture path. The fracture is filled with interface elements which in the 2D case are quadrangular and triangular elements and in the 3D case are either tetrahedral or wedge elements. The Rankine criterion is used for fracture nucleation and advancement. In a 2D setting the fracture follows directly the direction normal to the maximum principal stress while in the 3D case the fracture follows the face of the element around the fracture tip closest to the normal direction of the maximum principal stress at the tip. The procedure requires continuous updating of the mesh around the crack tip to take into account the evolving geometry. The updated mesh is obtained by means of an efficient mesh generator based on Delaunay tessellation. The governing equations are written in the framework of porous media mechanics and are solved numerically in a fully coupled manner. Numerical examples dealing with well injection (constant inflow) in a geological setting and hydraulic fracture in 2D and 3D concrete dams (increasing pressure) conclude the paper. A counter-example involving thermomecanically driven fracture, also a coupled problem, is included as well. The examples highlight some peculiar features of hydraulic fracture propagation. In particular the adopted method is able to capture the hints of Self-Organized Criticality featured by hydraulic fracturing.