TY - JOUR
T1 - Mathematical algorithm development and parametric studies with the GEOFRAC three-dimensional stochastic model of natural rock fracture systems
AU - Ivanova, Violeta M.
AU - Sousa, Rita
AU - Murrihy, Brian
AU - Einstein, Herbert H.
PY - 2014/6
Y1 - 2014/6
N2 - This paper presents results from research conducted at MIT during 2010-2012 on modeling of natural rock fracture systems with the GEOFRAC three-dimensional stochastic model. Following a background summary of discrete fracture network models and a brief introduction of GEOFRAC, the paper provides a thorough description of the newly developed mathematical and computer algorithms for fracture intensity, aperture, and intersection representation, which have been implemented in MATLAB. The new methods optimize, in particular, the representation of fracture intensity in terms of cumulative fracture area per unit volume, P32, via the Poisson-Voronoi Tessellation of planes into polygonal fracture shapes. In addition, fracture apertures now can be represented probabilistically or deterministically whereas the newly implemented intersection algorithms allow for computing discrete pathways of interconnected fractures. In conclusion, results from a statistical parametric study, which was conducted with the enhanced GEOFRAC model and the new MATLAB-based Monte Carlo simulation program FRACSIM, demonstrate how fracture intensity, size, and orientations influence fracture connectivity.
AB - This paper presents results from research conducted at MIT during 2010-2012 on modeling of natural rock fracture systems with the GEOFRAC three-dimensional stochastic model. Following a background summary of discrete fracture network models and a brief introduction of GEOFRAC, the paper provides a thorough description of the newly developed mathematical and computer algorithms for fracture intensity, aperture, and intersection representation, which have been implemented in MATLAB. The new methods optimize, in particular, the representation of fracture intensity in terms of cumulative fracture area per unit volume, P32, via the Poisson-Voronoi Tessellation of planes into polygonal fracture shapes. In addition, fracture apertures now can be represented probabilistically or deterministically whereas the newly implemented intersection algorithms allow for computing discrete pathways of interconnected fractures. In conclusion, results from a statistical parametric study, which was conducted with the enhanced GEOFRAC model and the new MATLAB-based Monte Carlo simulation program FRACSIM, demonstrate how fracture intensity, size, and orientations influence fracture connectivity.
KW - Connectivity
KW - Fractures
KW - Model
KW - Spatial
KW - Stochastic
UR - http://www.scopus.com/inward/record.url?scp=84897480144&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84897480144&partnerID=8YFLogxK
U2 - 10.1016/j.cageo.2013.12.004
DO - 10.1016/j.cageo.2013.12.004
M3 - Article
AN - SCOPUS:84897480144
SN - 0098-3004
VL - 67
SP - 100
EP - 109
JO - Computers and Geosciences
JF - Computers and Geosciences
ER -