User

 Guide    
Version  1.6a,  November  2016    
FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Introduc3on    
FracPaQ  is  an  open  source  toolbox  wri;en  in  
MATLAB™  and  publicly  available  on  GitHub  and  
Mathworks  FileExchange.  FracPaQ  runs  on  any  
computer  with  MATLAB  installed  –  Windows,  OS  X  
or  Linux.  The  iniRal  release  version  is  v1.6a.  
FracPaQ  is  designed  to  quanRfy  fracture  pa;erns  
in  rock.  The  user  supplies  either  an  image  file  of  
fractured  rock  or  a  text  file  of  traced  fractures  and  
their  (x,y)  coordinates.  From  either  kind  of  input,  
the  code  calculates  the  fracture  lengths,  angles  
and  connecRvity.  These  are  displayed  as  maps  and  
graphs,  and  saved  as  *.Rf  files.  EsRmates  are  made  
of  fracture  intensity  (P21),  fracture  density  (P20),  
and  permeability  in  2D  using  a  simple  parallel-­‐
plate  model  and  assuming  constant  aperture.      
  Lengths  
We  hope  that  the  code  will  enable  researchers  to  
quanRfy  fracture  pa;erns  in  an  open,  objecRve,  
and  consistent  way.    We  also  hope  that  people  
will  contribute  new  func3ons  and  new  tools,  and   Angles  
report  any  bugs  J    

User  Guide  Version  1.6a,  November  2016    

FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Introduc3on    
FracPaQ  is  wri;en  in  MATLAB™  and  totals  over  
4,000  lines  of  code.    As  in  any  so]ware  project  of  
this  scale,  there  will  be  ‘bugs’  –  i.e.  coding  errors.    
If  you  encounter  a  bug,  please  let  us  know  –  
through  GitHub,  Mathworks  FileExchange  or  by  e-­‐
mail  ([email protected]);  please  provide  as  
many  details  as  you  can,  including  (where  
possible)  a  screen  shot  of  the  error,  the  input  file  
you  were  using  at  the  Rme,  and  the  MATLAB  
version.    
 
FracPaQ  uses  code  wri;en  by  others:  
lineSegmentIntersect.m  by  U.  Murat  
Erdem  and  readtext.m  by  Peder  Axensten  
(both  available  on  Mathworks  FileExchange).    
Permeability  
 
The  handling  of  image  file  input  in  FracPaQ  also  
uses  funcRons  from  the  Image  Processing  Toolbox  
(version  9.4),  a  MATLAB™  add-­‐on.        
ConnecRvity  

User  Guide  Version  1.6a,  November  2016    

FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Installing  FracPaQ  
You  can  get  all  the  source  code  and  run  it  directly  from  the  MATLAB  command  window.      
 
The  source  code  and  this  User  Guide  are  available  on  GitHub:  
h;p://davehealy-­‐aberdeen.github.io/FracPaQ/  
 
and  on  the  Mathworks  FileExchange:  
h;ps://uk.mathworks.com/matlabcentral/fileexchange/58860-­‐davehealy-­‐aberdeen-­‐fracpaq  
 
Steps:  
1.  Download  the  source  code  as  a  .zip  or  .tar.gz  file,  and  extract  all  of  the  files  
2.  Put  all  of  these  files  into  a  single  folder    
3.  Start  MATLAB  (NB.  we  developed  FracPaQ  using  MATLAB  R2016a)  
4.  Set  the  current  working  folder  in  MATLAB  to  the  folder  you  installed  the  code  in  
5.  At  the  MATLAB  command  prompt  type  ‘guiFracPaQ2D’    and  hit  Enter    
 
Output  files  (graphs,  maps  and  data)  will  appear  as  *.Rf  and  *.txt  files  in  the  same  folder.      

User  Guide  Version  1.6a,  November  2016    

FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Star3ng  FracPaQ  
There  is  only  one  window  in  
the  applicaRon  (see  right).  The  
window  can  be  minimised  or  
closed  using  the  standard  GUI  
controls  (top  le]  in  the  Mac  OS  
X  version  shown).      
 
Input  parameters  are  on  the  
le]  side,  and  output  opRons  
are  on  the  right.    The  central  
area  is  for  previewing  the  input  
data  as  a  fracture  trace  map.        

To  get  started,  click  Browse…  to  see  a  list  of  possible  input  files  in  the  current  working  folder.    Then  
click  Preview  to  open  the  data  file  and  view  the  traces  in  the  main  window;  basic  staRsRcs  on  the  
pa;ern  are  also  shown  (lower  le]  box).    Select  the  Outputs  you  need  (right  hand  side),  and  then  click  
Run  to  produce  the  selected  maps  and  graphs.  Each  output  is  shown  in  a  separate  figure  window,  
and  saved  as  a  separate  .jpeg  in  the  current  folder.        
 
Click  Exit  to  quit  FracPaQ.    
User  Guide  Version  1.6a,  November  2016    
FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Choosing  an  input  file  
*.txt  file  example  

Click  Browse…  to  select  an  input  file.    The  file  

types  shown  are  filtered  to  those  with  file   Shown  above  is  an  example  of  a  *.txt  file  for  
extensions  *.txt,  *.jpeg,  *.jpg,  *.Rff  and  *.Rf.       input.    The  file  contains  (x,y)  nodes  along  each  
  fracture  trace.    The  (x,y)  data  are  tab-­‐delimited,  
If  you  select  a  text  (*.txt)  file  of  (x,y)  fracture   and  each  fracture  trace  is  on  one  line.    There  is  
trace  nodes,  the  Input  file  type  changes   a  minimum  of  4  columns  for  each  line  (i.e.  each  
automaRcally  to  ‘Node  file’.       fracture  trace)  –  (x1,  y1)  for  node  1  and  (x2,  y2)  
  for  node  2.      
If  you  choose  a  graphic  file  format  (*.jpeg,    
*.jpg,  *.Rff,  *.Rf),  the  Input  file  type  changes   A  fracture  trace  can  be  made  of  many  
to  ‘Image  file’.       segments;  in  the  example  shown  above  the  
  longest  trace  (line  7)  has  5  segments,  bounded  
You  can  manually  override  these  defaults  if   by  6  nodes  –  i.e.  12  columns  =  6  x  2  (x,y)  pairs.      
you  wish.      
User  Guide  Version  1.6a,  November  2016    
FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Image  file  input  –  the  Hough  transform  
FracPaQ  can  also  read  binary  image  files  as  input.    
Two  examples  are  provided,  
MacduffBinary.tif  and  
OrkneyBinary.tif.  The  code  uses  a  Hough  
transform  to  find  straight  lines  within  the  pixels.  
The  user  can  choose  the  values  at  which  these  lines  
are  merged  (if  they  are  close  enough  together)  or  
discarded  (if  they  are  too  short).  More  details  on  
the  Hough  transform  method  can  be  found  in  the  
MATLAB  Help  documentaRon.      
 
The  key  parameters  that  affect  the  number  and  
length  of  the  fracture  traces  found  are  the  number  
of  peaks  and  the  threshold.  Using  the  Hough  
transform  is  a  ‘trial  and  error’  procedure:  load  an  
image  file  and  click  Preview.  Then  use  the  graphs   NB:    the  Hough  transform  method  can  
(shown  on  the  right,  created  as  Figure  1)  and  the   only  find  straight  line  fractures  each  
trace  map  displayed  in  the  main  window  to  tune   made  of  a  single  trace;  mulR-­‐segment  
the  Hough  transform  parameters.  Set  the  threshold   traces  cannot  be  found.  Ntraces  =  Nsegments  
to  a  value  between  0.0  and  1.0  that  captures  only   for  this  method.    
the  significant  peaks  (in  orange/yellow  on  the  
graphs).  0.33  is  a  good  start.  
User  Guide  Version  1.6a,  November  2016    
FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Text  file  input  –  conversion  from  .svg  format  
Also  included  with  FracPaQ  is  a  Unix  shell  script  
to  convert  scaleable  vector  graphics  (*.svg)  files  
into  the  correct  tab-­‐delimited  text  (*.txt)  file  
format.  This  file  is  called  svg2fracpaq.csh.  
 
Standard  graphics  packages,  such  as  Adobe  
Illustrator,  CorelDraw  and  Inkscape,  can  all  
produce  *.svg  files.  Typically,  the  user  imports  
an  image  or  map  of  fractured  rock,  and  then  
traces  the  fractures  onto  a  new  layer  in  the  
so]ware.  A]er  deleRng  the  original  image,  the  
traces  can  be  saved  as  *.svg  format.    We  
recommend  using  *.svg  format  version  1.1.        
 
The  shell  script  reads  in  the  *.svg  file  of  your  
choice  and  produces  a  new  *.txt  file  with  each  
fracture  on  a  separate  line,  and  all  the  (x,y)  
nodes  separated  by  tabs.  For  Windows  users,  
there  may  be  support  for  running  Unix  shell  
scripts  under  Windows  10.          

User  Guide  Version  1.6a,  November  2016    

FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Maps  –  Traces,  segments  

Click  on  the  checkbox  “Traces,  segments”  

in  the  Maps  output  panel  to  see  a  map  of  
all  fracture  traces  and  segments.    The  
numbers  of  traces,  segments  and  nodes  
(i.e.  the  end  points  of  traces  and  
segments)  are  shown  in  the  figure  Rtle.      
 
This  plot  is  saved  in  the  current  folder  as  
FracPaQ2D_tracemap.tif,  at  a  
default  resoluRon  of  300  dpi.      
 
The  figure  is  produced  by  the  script  
guiFracPaQ2Dtracemap.m.  Edit  this  
script  file  to  change  any  default  sesngs,  
such  as  the  plot  line  colour  or  the  print  
file  resoluRon.    

User  Guide  Version  1.6a,  November  2016    

 FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Maps  –  Traces,  segments;  Show  nodes  

Click  on  the  checkbox  “Traces,  segments”  

and  also  check  “Show  nodes”  to  see  a  map  of  
all  fracture  traces  and  segments  with  all  the  
nodes  displayed.  The  numbers  of  traces,  
segments  and  nodes  (i.e.  the  end  points  of  
traces  and  segments)  are  shown  in  the  figure  
Rtle.  Trace  and  segment  end  nodes  are  
shown  with  black  circles,  and  trace  centres  
are  shown  with  red  circles.        
 
This  plot  is  saved  in  the  current  folder  as  
FracPaQ2D_tracemap.tif,  at  a  default  
resoluRon  of  300  dpi.  The  figure  is  produced  
by  the  script  
guiFracPaQ2Dtracemap.m.  Edit  this  
script  file  to  change  any  default  sesngs,  
such  as  the  plot  line  colour  or  the  print  file  
resoluRon.    

User  Guide  Version  1.6a,  November  2016    

 FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Maps  –  Intensity  (P21),  Density  (P20)  

Click  on  the  checkbox  “EsRmated  Intensity,  P21”  and/or  

“EsRmated  Density,  P20”  to  see  maps  of  esRmated  fracture  
intensity  and/or  density.    These  contour  maps  are  produced  
from  the  fracture  segment  data  using  the  circular  scan  line  
method  of  Mauldon  et  al.  (2001).    
 
The  code  places  scan  circles  (using  the  number  in  the  text  box)  in  
the  x-­‐  and  y-­‐direcRons,  and  counts  the  intersecRons  of  fracture  
segments  with  the  circle  perimeter  (n),  and  the  number  of  
segments  that  terminate  inside  the  circle  (m).  These  measures  
are  used  to  esRmate  Intensity  (units  of  L-­‐1)  and  Density  (units  of  
L-­‐2)  from  the  equaRons  in  Mauldon  et  al.  (2001).  P21  and  P20  
refer  to  the  measures  used  by  Dershowitz  &  Herda  (1992).  
 
These  plots  are  saved  in  the  current  folder  as  
FracPaQ2D_intensityP21.tif  and  
FracPaQ2D_densityP20.tif,  at  a  default  resoluRon  of  
300  dpi.  The  figure  is  produced  by  the  script  
guiFracPaQ2Dpattern.m.  Edit  this  script  file  to  change  any  
default  sesngs,  such  as  the  plot  line  colour  or  the  print  file  
resoluRon.     User  Guide  Version  1.6a,  November  2016    
 FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Graphs  –  Lengths  

Click  on  the  checkbox  “Histogram  of  lengths”  and/or  “Log-­‐log  

lengths”  to  see  a  histogram  and/or  log-­‐log  plots  of  of  fracture  
lengths.  Separate  plots  are  produced  for  fracture  traces  (trace  
length  =  sum  of  segment  lengths)  and  fracture  segments.  Select  
the  number  of  bins  for  the  histogram  from  the  drop-­‐down  list.    
 
Check  “Censor  trace  lengths  at  edges”  to  remove  from  the  
graphs  any  fractures  with  nodes  on  the  edges  of  the  trace  map  
area.  The  red  lines  on  the  graphs  show  the  minimum,  maximum  
and  maximum  possible  (do;ed)  lengths.    
 
The  plots  are  saved  in  the  current  folder  as  
FracPaQ2D_*length.tif,  at  a  default  resoluRon  of  300  
dpi.  The  figures  are  produced  by  the  script  
guiFracPaQ2Dlength_new.m.  Edit  this  script  file  to  change  
any  default  sesngs,  such  as  the  plot  line  colour  or  the  print  file  
resoluRon.    
User  Guide  Version  1.6a,  November  2016    
FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Graphs  –  MLE  Sta3s3cal  Analysis  

Power  Law  

Click  on  the  checkbox  “MLE  analysis”  to  perform  

a  Maximum  Likelihood  EsRmaRon  (MLE)  for  
fracture  trace  and  segment  lengths.  Separate  
plots  are  produced  for  three  possible  underlying  
staRsRcal  distribuRons:  Power  Law,  Log-­‐Normal  
and  ExponenRal.  See  Clauset  et  al.  (2007)  and  
Rizzo  et  al.  (in  review)  for  further  details  on  
MLE.    
 
The  plots  for  fracture  segments  are  saved  in  the  
current  folder  as  
FracPaQ2Dlengths_Fitting*.tif,  at  a  
default  resoluRon  of  300  dpi.  The  figures  are   Log-­‐Normal   ExponenRal  
produced  by  the  script  
guiFracPaQ2Dlength_new.m.  Edit  this  
script  file  to  change  any  default  sesngs,  such  as  
the  plot  line  colour  or  the  print  file  resoluRon.    
User  Guide  Version  1.6a,  November  2016    
FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Graphs  –  MLE  Sta3s3cal  Analysis  

The  MLE  plots  are  supplemented  with  text  outputs  

in  the  MATLAB™  Command  Window.  For  each  of  
the  3  distribuRons,  the  corresponding  staRsRcal  
parameters  are  given  and  the  probability  –  
expressed  in  %  and  calculated  based  on  a  
Kolmogorov-­‐Smirnoff  (KS)  test  –  that  the  analysed  
fracture  lengths  are  distributed  according  to  one  
of  the  distribuRons.  
 
More  details  and  relaRve  references,  can  be  found  
in  the  funcRons  fitting*.m.  called  by  
guiFracPaQ2Dlength_new.m.  Edit  this  script  
to  change  the  variables  ‘uc’  and  ‘lc’    to  cut  off  a  
certain  percentage  from  the  beginnig  (‘uc’)  or  the  
end  ('lc')  of  the  data  set.  Default  values  are  0%  cut-­‐
off.  

User  Guide  Version  1.6a,  November  2016    

 FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Graphs  –  Angles  (Orienta3ons)  

FracPaQ  provides  3  plots  of  orientaRon  data,  derived  from  the  

angles  of  the  trace  segments.  Segment  angles  are  calculated  with  
respect  to  the  Y-­‐axis,  taking  clockwise  as  posiRve.  To  change  this  in  
the  event  that  the  Y-­‐axis  is  not  aligned  due  North,  enter  a  number  in  
the  text  box  labelled  ‘RotaRon  of  Y-­‐axis  from  N’.  This  number  is  then  
subtracted  from  each  of  the  segment  angles  to  produce  a  
‘corrected’  plot.  The  number  of  bins  in  the  histogram  and  the  bin  
size  in  the  rose  diagram  can  be  selected  with  the  dropdown  list  
boxes.  The  rose  diagram  is  plo;ed  using  area-­‐weighRng  as  this  
provides  a  more  robust  measure  of  orientaRon  distribuRons  in  2D  
(Nemec,  1988).                    
 
The  plots  are  saved  in  the  current  folder  as  
FracPaQ2D_*angle.tif,  at  a  default  resoluRon  of  300  dpi.  The  
figures  are  produced  by  the  script  guiFracPaQ2Dangle.m.  Edit  
this  script  file  to  change  any  default  sesngs,  such  as  the  plot  line  
colour  or  the  print  file  resoluRon.    
User  Guide  Version  1.6a,  November  2016    
 FracPaQ  –  Fluid  Flow  
 
Fluid  flow  –  Connec3vity    

Click  on  the  checkbox  “I-­‐Y-­‐X  connecRvity”  to  see  a  

ternary  plot  of  relaRve  proporRons  of  I,  Y  and  X  nodes  
in  the  fracture  network  (Manzocchi,  2002).  I  nodes  are  
the  isolated  ends  of  traces,  Y  nodes  are  splays  or  
abutments,  and  X  nodes  are  intersecRons.      
 
Also  shown  on  this  plot  are  two  contour  lines  for  CL  –  
the  number  of  connecRons  per  line  (or  trace  in  our  
terminology)  from  Sanderson  &  Nixon  (2015).  Higher  
values  of  CL  are  indicaRve  of  be;er  connecRvity,  
although  the  contours  depend  on  other  a;ributes  of   The  figure  is  saved  in  the  current  folder  as  
the  fracture  network,  such  as  the  length  distribuRon   FracPaQ2D_IYXtriangle.tif,  at  a  default  
and  the  spaRal  distribuRon.     resoluRon  of  300  dpi.  The  figures  are  produced  by  the  
  script  guiFracPaQ2Dpattern.m.  Edit  this  script  
Be;er  connected  networks  tend  to  plot  towards  the   file  to  change  any  default  sesngs,  such  as  the  colour  
bo;om  of  this  diagram  (higher  proporRons  of  Y  and  X   or  the  print  file  resoluRon.    
nodes).      
User  Guide  Version  1.6a,  November  2016    
 FracPaQ  –  Fluid  Flow  
 
Fluid  flow  –  Permeability    

Click  on  the  checkbox  “Permeability  ellipse”  to  see  the  

2D  permeability  ellipse  for  the  fracture  network.  This  
esRmate  is  produced  from  a  combinaRon  of  lengths,  
orientaRons  and  spaRal  densiRes  using  the  tensorial  
method  of  Suzuki  et  al.  (1998).  The  long  axis  of  the  
ellipse  is  proporRonal  to  the  maximum  permeability  
(√k1),  and  the  short  axis  is  proporRonal  to  the  
minimum  permeability  (√k2)  (this  is  for  permeability  in  
the  direc3on  of  flow;  see  Long  et  al.,  1982  for  details).    
 
The  user  can  make  changes  to  the  Aperture  (default  is   The  figure  is  saved  in  the  current  folder  as  
1,  pixels  or  metres)  and  the  Lambda  factor,  which  lies   FracPaQ2D_permtensor.tif,  at  a  default  
between  0  and  1.  If  all  fractures  are  fully  connected,   resoluRon  of  300  dpi.  The  figures  are  produced  by  the  
Lambda  can  be  set  to  1  (Suzuki  et  al.,  1998).  The  user   script  guiFracPaQ2Dtensor.m.  Edit  this  script  file  
can  adjust  this  variable  based  on  the  results  obtained   to  change  any  default  sesngs,  such  as  the  colour  or  
from  the  connecRvity  plot.   the  print  file  resoluRon.    

User  Guide  Version  1.6a,  November  2016    

 FracPaQ  –  Fluid  Flow  
 
Permeability  Tensor  Ellipse  

The  permeability  ellipse  is  supplemented  with  text  

outputs  in  the  MATLAB™  Command  Window.  
 
The  numerical  values  of  k1  (maximum),  k2  (minimum),  
and  the  azimuth  of  k1  (theta,  in  degrees  measured  from  
the  Y-­‐axis),  are  printed  in  the  Command  Window.      
 
The  units  of  k  default  to  pixels2,  unless  a  Scaling  Factor  
(pixels/metre)  has  been  entered  on  the  le]  hand  side  of  
the  main  window,  in  which  case  the  units  are  metres2.        

User  Guide  Version  1.6a,  November  2016    

 FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Bibliography  
Bonnet,  E.,  Bour,  O.,  Odling,  N.E.,  Davy,  P.,  Main,  I.,  Cowie,  P.  &   Mauldon,  M.,  Dunne,  W.M.  &  Rohrbaugh,  M.B.  (2001).  Circular  
Berkowitz,  B.  (2001).  Scaling  of  fracture  systems  in  geological   scanlines  and  circular  windows:  new  tools  for  characterizing  the  
media.  Reviews  of  geophysics,  39(3),  pp.347-­‐383.   geometry  of  fracture  traces.  Journal  of  Structural  Geology,  23(2),  
  pp.247-­‐258.  
Clauset,  A.,  Shalizi,  C.  R.  &  Newman,  M.  E.  J.  (2009).  Power-­‐law    
distribuRons  in  empirical  data.  SIAM  Review,  51(4),  661.     Nemec,  W.  (1988).  The  shape  of  the  rose.  Sedimentary  Geology,  
  59(1-­‐2),  pp.149-­‐152.  
Dershowitz,  W.S.  &  Herda,  H.H.  (1992).  InterpretaRon  of  fracture    
spacing  and  intensity.  In:  The  33th  US  Symposium  on  Rock   Newman,  M.  E.  J.  (2005).  Power  laws,  Pareto  distribuRons  and  
Mechanics  (USRMS).  American  Rock  Mechanics  AssociaRon.   Zipf’s  law.  Contemporary  Physics,  46(5),  323–351.  
   
Farrell,  N.J.C.  (2015).  Unpublished  PhD  thesis.    University  of   Rizzo,  R.E.,  Healy,  D.  &  De  Siena,  L.  (in  revision).  The  Benefits  of  a  
Aberdeen.     Maximum  Likelihood  EsRmator.  
   
Gillespie,  P.A.,  Howard,  C.B.,  Walsh,  J.J.  &  Wa;erson,  J.  (1993).   Sanderson,  D.J.  &  Nixon,  C.W.  (2015).  The  use  of  topology  in  
Measurement  and  characterisaRon  of  spaRal  distribuRons  of   fracture  network  characterizaRon.  Journal  of  Structural  Geology,  
fractures.  Tectonophysics,  226(1-­‐4),  pp.113-­‐141.   72,  pp.55-­‐66.  
   
Healy,  D.,  Rizzo,  R.E.,  Cornwell,  D.C.,  Farrell,  N.J.C.,  Watkins,  H.,   Smith,  M.  (2016).  Unpublished  MSc  thesis.  University  of  
Timms,  N.E.,  Gomez-­‐Rivas,  E.  &  Smith,  M.  (in  revision).    FracPaQ:   Aberdeen.    
a  MATLAB™  toolbox  for  the  quanRficaRon  of  fracture  pa;erns.        
  Suzuki,  K.,  Oda,  M.,  Yamazaki,  M.  &  Kuwahara,  T.  (1998).  
Long,  J.C.S.,  Remer,  J.S.,  Wilson,  C.R.  &  Witherspoon,  P.A.  (1982).   Permeability  changes  in  granite  with  crack  growth  during  
Porous  media  equivalents  for  networks  of  disconRnuous   immersion  in  hot  water.  InternaLonal  Journal  of  Rock  Mechanics  
fractures.  Water  Resources  Research,  18(3).     and  Mining  Sciences,  35(7),  pp.907-­‐921.  
   
Manzocchi,  T.  (2002).  The  connecRvity  of  two-­‐dimensional   Watkins,  H.  (2015).    Unpublished  PhD  thesis.    University  of  
networks  of  spaRally  correlated  fractures.  Water  Resources   Aberdeen.  
Research,  38(9).  
 
  User  Guide  Version  1.6a,  November  2016    
FracPaQ  –  Fracture  Pa1ern  Quan3fica3on  
 
Acknowledgements  
lineSegmentIntersect.m  wri;en  by  U.  Murat  Erdem;  available  on  Mathworks  FileExchange:  
h;p://uk.mathworks.com/matlabcentral/fileexchange/27205-­‐fast-­‐line-­‐segment-­‐intersecRon  
 
readtext.m  wri;en  by  Peder  Axensten;  available  on  Mathworks  FileExchange:  
h;p://uk.mathworks.com/matlabcentral/fileexchange/10946-­‐readtext  
 
We  thank  Paul  Gillespie  (Statoil),  Mark  Fischer  (Northern  Illinois  University)  and  Bill  Dunne  (University  of  
Tennessee;  and  Editor  of  the  Journal  of  Structural  Geology)  for  comments  on  our  submi;ed  manuscript,  which  
have  led  to  major  improvements  in  the  code,  the  paper  and  the  User  Guide.      

User  Guide  Version  1.6a,  November  2016