In the upper 5--10 km of the crust fractures may remain partially open, exerting a strong influence over the response of their host material to seismic waves due to their superior deformability. This raises the possibility of using seismic data to determine not only the location but also the properties of fracture networks. We present a particle-based numerical scheme for modeling rock deformation and wave propagation and show that it accurately captures wave propagation across fractures. We investigate the effect of fracture tensile properties (i.e. fracture cohesion) on the wavefield and demonstrate that seismic data can, in theory, be used to distinguish fracture tensile and compressional properties.