Impact of mineralisation on digital coal properties - ACS Publications

fundamental requirement for enhanced CSG recovery because the mineralisation of coal is. 70 ... micro-CT data. We find that the characterisation of mi...
0 downloads 11 Views 2MB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

Article

Impact of mineralisation on digital coal properties Yu Jing, Ryan Troy Armstrong, and Peyman Mostaghimi Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.7b01401 • Publication Date (Web): 27 Aug 2017 Downloaded from http://pubs.acs.org on August 29, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

Energy & Fuels

Impact of mineralisation on digital coal properties

2

Yu Jing, Ryan T. Armstrong and Peyman Mostaghimi

3

School of Petroleum Engineering, The University of New South Wales, NSW 2052, Sydney,

4

Australia

5

Abstract

6

Coal seam gas (CSG) is an unconventional energy resource, whose production is mainly

7

controlled by the underlying fracture networks, known as “cleats”. The natural cleats are

8

generally mineralised during diagenesis, which significantly reduces fracture conductivity

9

because of the more tortuous flow pathways and smaller cross-sectional area perpendicular to the

10

flow. This paper aims at characterising the mineral fillings of the coal cleat network by utilising

11

X-ray micro-computed tomography (micro-CT) imaging to investigate the effect of minerals on

12

the petrophysical properties of coal. We apply a high-resolution micro-CT to obtain a 3D digital

13

representation of a mineralised coal sample collected from Moura mine of Bowen Basin. The

14

main components of the coal sample, including coal matrix, cleats and minerals, are individually

15

analysed to give the statistics of cleat orientation, length, and thickness using quantitative image

16

analysis. According to the measured statistical data, digital coal models with mineralisation are

17

stochastically constructed and then used for simulation to obtain petrophysical properties. Our

18

results show that minerals in cleats can significantly reduce the conductivity of the cleat network 1

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 31

19

by reducing the permeability by up to 75% and increasing the tortuosity by 21%. In addition, by

20

studying the deformed sample under external stresses, we find that porosity reduction due to

21

compression of the cleat network with minerals is significantly less than models without

22

minerals. In addition, mineralisation has a detrimental effect on the matrix-fracture contact area,

23

which is reduced by 17% with the addition of minerals. Lastly, by comparing digital coal models

24

to the original micro-CT images, the mineralised digital coal models are found to be more

25

representative in terms of permeability estimation, with an error of only 4.4%.

26

Keywords: Coal seam gas, mineralisation, digital coal model, coal permeability, pore-scale

27

modelling

28

1 Introduction

29

Coal, as a complex and heterogeneous rock, is composed of three physical components, coal

30

matrix, fractures and minerals 1. The matrix, also known as “maceral”, is fundamental for

31

combustion and methane adsorption 2. Coal is formed from the accumulation of plant materials

32

under high temperature and high pressure during long-time periods. Therefore, coal matrix is

33

divided into multiple lithotypes that originate from different plant materials

34

appear to be banded and composed of alternating bright and dull materials 5, 6. The bright band is

35

rich in vitrain, which is a brittle material with bright lustre. It is permeated with fine cracks at

36

right angles, giving a blocky appearance. In contrast, the dull band, comprised of durain, is a

37

grey to black material with a dull lustre.

3, 4

. Most coals

38

Coal fractures, also known as “cleats”, are the dominant flow pathways that determine the

39

permeability and production rate of coal seam gas (CSG) 7, 8. Cleats in different lithotypes tend to 2

ACS Paragon Plus Environment

Page 3 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

40

have different properties. For example, cleats in bright bands, named “bright cleats”, occur in

41

two orthogonal sets of sub-parallel cleats known as face and butt cleats. Face cleats form first

42

during coalification and are easier to visualise because they extend longer than butt cleats 9. Butt

43

cleats occur later due to the relaxation of the original stress field and terminate at face cleats 10, 11,

44

such that the connectivity pattern of an organised bright cleat system is mostly present as “T-

45

junctions” between face and butt cleats 12. In contrast, cleats of dull bands lack a regular pattern,

46

where face and butt cleats are hardly recognised 5, 13. Dull cleats are poorly developed and rarely

47

observed in coal samples 5, 14.

48

Minerals are the inorganic material formed during peat accumulation as well as changes to the

49

subsurface fluids 15, 16. Based on the time of formation, coal minerals are classified as syngeneic

50

or epigenetic, where syngeneic minerals formed during peat formation while the epigenetic form

51

after coal has matured 17. Epigenetic minerals are commonly deposited in cleats and are mainly

52

composed of kaolinite, illite, pyrite and calcite

53

matrix or cleats

54

morphology

55

similar to the coal cleats network. For example, “T-junction” connections can be observed in the

56

mineral phase, which resembles the connectivity pattern of face and butt cleats

57

fillings are also oriented perpendicular to the bedding plane. Thus, mineral fillings can also be

58

described by geometrical properties, such as orientation, length and mineral size. With the help

59

of X-ray micro-computed tomography (micro-CT), geometrical properties of minerals can be

60

statistically analysed by applying quantitative image analysis 24-30.

22, 23

20

15, 18-21

. Minerals can occur within either coal

. Mineral fillings distributed in cleats tend to have a well-developed

. According to Zhang et al.

22

, the morphology of mineral fillings in cleats is

22

. Mineral

61

The occurrence of minerals in cleats during diagenesis can significantly affect the conductivity

62

of coal 31. Because minerals that are lined or filled in the void cleat space 18, 19 tend to block gas 3

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 31

63

flow pathways, thus the hydraulic fracture aperture is lowered and fluid pathways are more

64

tortuous. As a result, the cleat network with mineralisation is less interconnected or even non-

65

conductive depending on the degree of mineralisation

66

to be able to prevent coal fractures from completely closing under external stress 34. Furthermore,

67

minerals can decrease the gas sorption capacity of the coal matrix by reducing the pore volume

68

and internal surface area

69

reacts with organic matrix is reduced

70

fundamental requirement for enhanced CSG recovery because the mineralisation of coal is

71

closely associated with its permeability and methane sorption capacity.

2, 20

32, 33

. However, mineral fillings are found

. Because of the presence of minerals, the contact area where gas 23

. Therefore, the study of coal minerals is a primary and

72

In this paper, a mineralised coal sample collected from Moura Mine of Bowen Basin is

73

studied. We apply micro-CT technology to obtain 3D digital images, based on which the coal

74

matrix, cleat network and mineral matter are individually analysed to obtain statistical data.

75

According to measured statistical data, an advanced digital coal model is constructed, which

76

preserves the banding information of the coal matrix, geometrical properties of cleats as well as

77

the mineralisation. To investigate the effect of mineralisation on digital coal properties, the

78

petrophysical properties of our digital coal models are calculated and compared to the original

79

micro-CT data. We find that the characterisation of mineralisation is crucial for digital coal

80

modelling since minerals have a significant impact on cleat conductivity and porosity.

4

ACS Paragon Plus Environment

Page 5 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

81

2 Methodology

82

2.1 Micro-CT Scanning and Image Processing

83

The sample studied in this work is from Moura mine of Bowen Basin, which is a medium

84

volatile bituminous coal with a vitrinite reflectance of 1.15% 35. A high-resolution, helical micro-

85

CT scanner developed at the Australian National University 36 is applied to scan the coal sample

86

to obtain its three-dimensional internal structure. The 3D micro-CT data are represented by an

87

array of voxels, where each voxel value corresponds to the X-ray attenuation coefficient of a

88

given phase

89

the investigated object

90

grey-scale values in the micro-CT data, which are also called “CT numbers”. Therefore, the main

91

coal components, including matrix, cleats and mineral matter, can be identified by the particular

92

range of CT numbers

93

by applying image segmentation methods 1. For example (Figure 1), in the grey image, cleats are

94

shown in dark black, the mineral matter is highlighted with high CT numbers (bright colours),

95

and the coal matrix is identified by medium CT numbers (light grey). For details on the image

96

segmentation process, please refer to Ramandi et al.

97

grey-scale images are partitioned into distinct phases, where each component can be extracted

98

for further analyses 42, 43.

37-40

. As the voxel values are a function of density and effective atomic number of 41

, different coal components with distinct densities will have different

24, 25, 41

, such that the grey-scale image can be partitioned to unique phases

35

. After the segmentation, the continuous

5

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 31

99 100

Figure 1: In the grey image (a), cleats are shown in dark black colour, the mineral matter are

101

highlighted with high CT numbers, and the coal matrix is identified by medium CT numbers. In

102

the segmented image (b), the continuous grey-scale image is segmented into distinct phases:

103

cleats (green), minerals (red), bright bands (grey) and dull bands (black).

104

The 3D digital micro-CT images studied in this work have a dimension of 500×500×500

105

voxels with the resolution of 16.5 µm (Figure 2a). As can be seen, bright and dull bands are

106

alternatively layered in this sample (Figure 2a). Two bright bands can be observed in this

107

domain, both of which have well-developed cleat networks (Figure 2b). Minerals (Figure 2c) are

108

widely distributed in both bright and dull bands, while minerals in bright bands have more

109

regular morphology as cleats.

110 6

ACS Paragon Plus Environment

Page 7 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

111

Figure 2: (a) The 3D micro-CT image has the dimension of 500×500×500 voxel with a

112

resolution of 16.5 µm; (b) Two bright bands are observed within the sample, both of which are

113

well-fractured; (c) Mineral matter occurs in both bright and dull bands, while minerals in bright

114

bands have more regular morphology as cleats.

115

2.2 Statistics Acquisition

116

Based on the segmented micro-CT images, all three components of coal are extracted for

117

quantitative analysis individually. For the matrix phase, a binary lithotype profile consisting of

118

bright and dull bands is obtained by applying the thresholding method

119

are measured based on the lithotype profile: (1) volume fraction of bright bands, (2) band

120

thickness and (3) band orientation. The volume fraction of bright bands is calculated by dividing

121

the number of the bright band voxels by the total number of voxels (500×500×500). To measure

122

the band thickness and orientation, an edge detector with Sobel approximation 44 is used to detect

123

the band boundaries, where points with maximum intensity gradient are highlighted as edges.

124

Then the highlighted boundaries are used for band thickness and orientation measurements,

125

where the band thickness is determined by measuring the spacing of two edges of a band . Since

126

dull cleats are poorly developed, only bright cleats are extracted for analyses by overlapping the

127

lithotype profile onto the cleat network phase. Geometrical properties, including orientation,

128

length, spacing and aperture size are measured by an automatic quantitative image analysis

129

method 29. The aperture size is defined as the opening width of a fracture. Prior to measurements,

130

a cleat-grouping algorithm

131

independent geometrical statistics since face and butt cleats have different properties

29

1, 24-26

. Three parameters

is utilised to partition face cleats and butt cleats, providing two 45

.

7

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 31

132

Measurements of geometric properties and the cleat grouping process based on micro-CT images

133

are explained in detail by Jing et al. 29.

134

Minerals occurring along the cleat surfaces are observed to have well-developed morphologies

135

that are similar to the cleat network 23. Similar with cleats, the mineral fillings are also oriented

136

perpendicular to the bedding plane and connected with “T-junctions”, which resemble the

137

morphology of face and butt cleats 22. This type of pattern indicates that minerals are deposited

138

after formation of the cleats 27. Frequency, volume fraction of minerals and length of each bright

139

band are measured separately. Frequency is the number of mineral fillings that are observed in

140

cleats. The volume fraction of minerals represents the percentage of mineral volume to the bright

141

band volume, where the mineral volume can be determined by the number of mineral-identified

142

voxels. Length refers to the extent of a mineral along the cleat surface. Herein, length is

143

measured based on the 2D slices of micro-CT images. Minerals on the micro-CT images are first

144

skeletonised with a thinning algorithm

145

Therefore, the length of a mineral can be calculated by counting the number of voxels that the

146

mineral has.

147

46

, reducing the widths of minerals to be one-voxel.

As geometric properties of cleats, such as spacing and aperture, are functions of band thickness

148

45, 47-49

149

complete set of statistical data is obtained (Figure 3). Each bright band is divided into cleats and

150

minerals, where cleats are further grouped into face cleats and butt cleats.

, individual bright bands with unique thicknesses are analysed separately. As a result, a

8

ACS Paragon Plus Environment

Page 9 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

151 152

Figure 3: The structure of mineral statistics data. Different bright bands of various band

153

thickness values are extracted for further analysing. Under individual bright bands, minerals

154

filling within different cleat groups are distinguished and statistically analysed separately.

155

2.3 Mineralised Digital Coal Model

156

A mineralised digital coal model (Figure 4c) is constructed numerically to characterise the

157

three main components of coal, including band information (Figure 4a), cleat geometrical

158

properties (Figure 4b), and mineral fillings (Figure 4c). The construction consists of three steps:

159

(1) generating layered bright and dull bands; (2) constructing DFN models in bright bands; and

160

(3) distributing minerals in the constructed DFN models.

a

b

c

2 mm

161 9

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 31

162

Figure 4: The mineralised digital coal model (c) contains band information (a) and cleat

163

geometrical properties (b). The colours indicate different phases of coal: bright bands (black),

164

dull bands (grey), cleats (green), minerals (red).

165

Based on the banding statistics of the coal matrix, a banded model is constructed where bright

166

and dull bands with stochastic band thicknesses are alternatively layered. Boundaries of each

167

band are generated according to the thickness and orientation distribution of bands. A region-

168

filling algorithm

169

procedure is repeated until the volume fraction of bright bands is equal to that of the original

170

micro-CT data with an error less than 5%. Furthermore, the bright cleat network is merged with

171

the banded model by generating discrete fracture network (DFN) models within the bright bands.

172

The DFN model is comprised of discrete cleat planes, whose geometrical properties follow

173

statistical data

174

band, DFN models are stochastically constructed within corresponding bright bands. As coal

175

cleats of bright bands have a particular network pattern, an improved DFN modelling method is

176

developed by Jing et al.

177

connectivity as the cleat connections of the original sample. Dull bands are assumed to be

178

impermeable in this work, so there are no cleat networks constructed in dull bands. This is

179

aligned with the observations of Ramandi et al.

180

dull bands is 3-4 orders of magnitude smaller than that of bright bands for the same coal sample.

50, 51

52-54

is then applied to fill the volume within boundaries. The modelling

. In this work, according to the cleat statistics of the corresponding bright

29

where face cleats and orthogonal butt cleats have “T-junction”

35

where they showed that the permeability of

181

Based on the cleat network, minerals are randomly generated along bright cleat surfaces.

182

According to the frequency statistics, a certain number of mineral fillings are initially generated,

183

where one mineral filling is made up of mineral voxels. The resulting volume fraction of

184

generated minerals is calculated and compared with that of original micro-CT data. To insure 10

ACS Paragon Plus Environment

Page 11 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

185

that the total volume of generated minerals is identical to that of the original micro-CT data, a

186

volume adjusting process is developed to vary the size of the mineral volume by changing

187

adjacent voxels to be solid (mineral phase) or void (cleat space). After each volume adjusting

188

process, the resulting volume fraction of minerals is calculated and compared to the fraction

189

value measured from the micro-CT images. The aforementioned procedure is stopped when we

190

obtain a mineralised cleat network with the desired degree of mineralisation in terms of mineral

191

frequency and mineral volume fraction. The workflow of the mineralisation process is illustrated

192

in Figure 5.

193 194

Figure 5: Flowchart of mineralisation process. Based on the cleat network, certain pieces of

195

minerals are generated. A volume adjusting process of varying the mineral size is applied until

196

the total volume of generated minerals is identical to that of original micro-CT images. Lastly,

197

mineralised cleat network with desired mineralisation degree is obtained.

11

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

198

2.4 Evaluation of Mineralised Digital Coal Models

199

2.4.1 Porosity and permeability

Page 12 of 31

200

Porosity (φ) based on a segmented binary image is determined by dividing the number of cleat-

201

identified voxels by the total number of voxels. Permeability (k) is simulated by a Navier-Stokes

202

solver in GeoDict software

203

directional permeability, a pressure difference is applied in the flow direction of the domain and

204

fluid flow is simulated through the sample. The boundary condition in the flow direction is set to

205

be periodic with a constant pressure difference, while no-flow boundary is defined in the

206

tangential direction. Additional void regions of 10-voxels are added at both inflow and outflow

207

boundaries to avoid the possibility of dead-end flow channels under periodic conditions. No-slip

208

boundary conditions are applied at other boundaries. The simulation will stop if the change of

209

permeability is less than 5% between successive iterations. Both matrix and minerals are

210

considered impermeable in the flow simulation.

211

2.4.2 Deformation

55, 56

, which uses the finite volume method

57-59

. To determine the

212

The cleat network alters when the structure is deformed under external stresses 60. In order to

213

study the impact that minerals can have on this alteration, the sample is subject to external loads

214

and simulated for permeability after structure deformation. We use the Elastodict solver of the

215

GeoDict software

216

composition, e.g. coal matrix and minerals, is assigned with its corresponding mechanical

217

parameters, such as Young’s modulus (E) and Poisson ratio (ν). Herein, the minerals are

218

assumed to be calcite, which is a common mineral type for coal cleats

219

parameters of calcite are set to be E =84 Gpa and ν = 0.32 63. For the coal matrix, we define the

61

to solve the elasticity equations

62

. For the mineralised digital model, each

18, 19

. The mechanical

12

ACS Paragon Plus Environment

Page 13 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

64

220

mechanical properties according to the work of Aziz et al.

221

load type of compression is applied on the model boundaries with specified strain (in %). Based

222

on the input information, including mechanical properties of each constituent and applied

223

boundary conditions, stresses and corresponding strains at each solid voxel are computed based

224

on Hooke’s Law 65,

225

, where E =4 Gpa and ν = 0.3. A

 = ∑ ,  ɛ , ,   {1, 2, 3}

(1)

226

where  is the stress tensor, ε is strain tensor and C is symmetric elasticity or stiffness tensor,

227

which is given by, 2μ + λ

λ

λ λ

0 0 0 μ 0 0

0 0 0 0 μ 0

0 0! 0 0 0 μ

228

  C=   

229

where λ and μ are Lamé and shear moduli, which depend on mechanical properties of specific

230

λ λ

0 0 0

λ

0 0 0

2μ + λ 0 0 0

(2)

material, including Young’s moduli (E), Poisson ratio (ν) and bulk moduli (B): %( λ'(%)

231

E=

232

ν = ((λ'%)

233

B=

234

2μ + λ

(3)

λ'% λ

(4)

λ'(%

(5)



2.4.3 Specific surface area

235

Specific surface area (SSA) is defined as the total surface area per bulk volume. Because of the

236

presence of minerals along cleat surfaces, the contact between gas and the organic matrix is

237

changed

238

digital images, boundary voxels where there is a change in the voxel values are highlighted

23

, which in turn has been shown to influence reaction rates and fluid transport

66

. In 67

.

13

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 31

239

For each boundary voxel, we identify surfaces that are exposed to the interface. Thus, SSA is

240

determined by summing the surface area of voxels that are on the solid-void boundaries 30.

241

2.4.4 Tortuosity

242

Tortuosity (τ) is a parameter that characterises the sinuosity and interconnectedness of porous 68, 69

243

media

. Herein, geometrical tortuosity is calculated, which is defined as the ratio of the

244

shortest length of the path connecting two points in the pore space to the straight-line distance.

245

So, it is a structural character of the medium, independent of any particular transport process. In

246

binary image, tortuosity is determined by the geodesic distance (*+ ) of two void voxels divided

247

by the Euclidean distance (*, ) between them 70, ./

248

τ=

249

Herein, the geodesic distance is determined by the city block distance that examines the

(3)

.0

250

absolute differences between the coordinates of two objects

251

voxels A (1 , 2 , 3 ) and B (1( , 2( , 3( ) in three dimensions is,

71

. For example, the distance of

252

*+ = |1 − 1( | + |2 − 2( | + |3 − 3( |

253

However, for 3D digital data, the geodesic and Euclidean distances between two opposite

254

surfaces rather than two void voxels are measured. The geodesic distance of every voxel in the

255

cleat space to the sample surfaces is computed

256

there are a variety of geodesic distance values where the average geodesic distance is obtained.

257

Next, the average geodesic distance values are plotted as a function of corresponding Euclidean

258

distances, which provides a linear correlation. The slope of the fitted linear relationship is the

259

tortuosity.

(4)

72

. As a result, for a certain Euclidean distance,

14

ACS Paragon Plus Environment

Page 15 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

260

3 Results and Discussion

261

3.1 Statistics

262

Table 1 summarises the measured distributions of the three components: coal matrix, cleat

263

network and mineral matter. The two bright bands (Bright Band #1 and #2) are oriented in

264

parallel and approximately 1.62mm and 3.67mm in thickness, respectively. Also, the volume

265

fraction of bright bands is 61.5%. For each bright band, cleats and minerals are individually

266

extracted for the analysis. According to the statistics, the face and butt cleats are almost equally

267

spaced, while spacing values vary within different bright bands. For example, Bright Band #1 is

268

fractured more coarsely than Bright Band #2 in terms of the development of face cleats. Because

269

of the limited size of the sample, the thickness values are measured from parts of bright bands in

270

the studied domain rather than the entire sample. Therefore, the thickness cannot represent the

271

real thickness of the bright bands. However, according to the findings of our previous work 45, 47-

272

49

273

be indicated that the whole Bright Band #1 is thicker than Bright Band #2. Furthermore, the face

274

cleats in Band #1 have an average aperture size of 0.17mm, slightly larger than those in Band #2

275

(0.16mm). Based on the linear positive correlation between aperture size and band thickness

276

provided by Close and Mavor et al. 47, Bright Band #1 is also inferred to be thicker than Bright

277

Band #2. The orientation distributions of the bands and cleat network indicate that face and butt

278

cleats are mutually orthogonal and both normal to the bedding plane. For mineralisation, Bright

279

Band #1 is more extensively mineralised compared with Bright Band #2. Bright Band #1 has

280

longer minerals fillings (the average length is 5.15mm) than Bright Band #2 (3.92mm). The

281

volume fraction of minerals of Band #1 is 3 times greater than that of Band #2. In addition, 15

, average spacing of cleats is linearly proportional to the coal band thickness. Therefore, it can

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 31

282

mineralisation is not the same for different cleat families: face cleats are filled with minerals,

283

which indicates that the minerals deposit after the creation of face cleats; but there are no

284

minerals observed in butt cleats in the studied sample, since they may be formed after water has

285

been driven from the coal, such that they are less likely to be mineral coated

286

different bright bands and cleat families tend to have a different degree of mineralisation, leading

287

to the importance of the proposed coal banding characterisation and cleat family grouping

288

process.

289

Table 1: Measured statistical data, including banding, cleat network and mineralisation, of a coal

290

sample from Moura Mine of Bowen Basin.

73-75

. Therefore,

Bright Band #1 Bright Band Thickness Ave. 1.65

Azimuth Dev. 0.08

Ave. 5.08

Dip Dev. 3.98

Ave. 92.75

Dev. 3.74

Aperture Ave. Dev. 0.17 0.03 0.04 0.02

Length Ave. Dev. 3.72 2.65 1.83 1.21

Cleat Network

Face cleats Butt cleats

Azimuth Ave. Dev. -2.41 4.61 92.26 4.78

Dip Ave. Dev. 94.41 4.59 91.23 4.67

Spacing Ave. Dev. 2.72 0.37 2.19 0.43

Mineral Volume Fraction 0.0160

Frequency Face cleats Butt cleats 3 0

Length Ave. 5.15

Dev. 3.71

Bright Band #2 Bright Band Thickness Ave. 3.67

Azimuth Dev. 0.63

Ave. 4.89

Dip Dev. 3.82

Ave. 91.45

Dev. 3.46

Aperture Ave. Dev. 0.16 0.03

Length Ave. Dev. 3.22 1.66

Cleat Network

Face cleats

Azimuth Ave. Dev. -2.52 4.54

Dip Ave. Dev. 95.42 4.65

Spacing Ave. Dev. 2.17 0.39

16

ACS Paragon Plus Environment

Page 17 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Butt cleats

91.44

4.74

90.71

4.79

2.30

0.52

0.04

0.02

1.04

0.46

Mineral Volume Fraction 0.0056

291

Frequency Face cleats Butt cleats 2 0

Length Ave. 3.92

Dev. 2.85

3.2 Study of mineral effects on digital coal models

292

The mineralised digital coal models are constructed with identical banding information, cleat

293

geometrical properties and mineral distributions as the original miro-CT images (Figure 6). It is

294

shown that the digital coal model preserves the main features of the original micro-CT images:

295

(1) the model comprises alternating bright and dull bands, where the volume fraction of bright

296

bands is 61.5%; (2) bright bands have well-developed cleat networks with “T-junction”

297

connectivity, and are constrained by adjacent dull bands; and (3) there are 5 mineral fillings with

298

the total mineral volume fraction of 0.8% within the bright cleats, which is identical to the

299

mineralisation of original micro-CT images. A total of 15 pairs of digital coal realisations are

300

stochastically constructed with the Monte Carlo method

301

one mineralised digital coal model and one without minerals. Next, porosity, permeability,

302

specific surface area and tortuosity are computed for each realisation.

76

. Each pair of realisations consists of

17

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 31

303 304

Figure 6: An example digital coal model with mineralisation (b) with a clip in y-direction. It

305

preserves the main features of original micro-CT images (a).

306

3.2.1 Petrophysical properties

307

In total 15 model pairs, i.e. with and without minerals are stochastically generated. Each model

308

pair has an identical cleat network. We then obtain the porosity and simulate flow through each

309

model pair to measure permeability. The average petrophysical properties of all model pairs are

310

plotted as a function of the number of models that are used to give the average values (Figure 7).

311

It can be seen that the average porosity and permeability of digital models are significantly

312

reduced by mineralisation. For example, the average porosity of digital models without minerals

313

reaches a plateau of 4.1%, which decreases to 3.2% after mineralisation. The effect of

314

mineralisation on permeability is more dramatic, where average permeability is reduced by

315

approximately 75%. Figure 8 provides insight into how minerals influence the flow field by

316

presenting the streamlines for a model with and without minerals. For the mineralised digital

18

ACS Paragon Plus Environment

Page 19 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

317

models, the mineral fillings block the pathways of fluid (Figure 8c), resulting in lower

318

transmissibility.

319

Additionally, the mineralised digital coal models provide permeability values that are closer to

320

the original micro-CT images than digital coal models without minerals. Specifically, the

321

average permeability of digital coal models without mineralisation stabilises at 3.96D with an

322

error of 16.8% compared with the non-mineralised original micro-CT images (4.76D).

323

Conversely, mineralised digital coal models have an average permeability of 0.94D, which is

324

near the permeability of the mineralised original micro-CT images (0.90D), with an error of only

325

4.4%.

326 327

Figure 7: The average porosity (a) and permeability (b) as a function of the number of models

328

that are used to give the average values. It can be seen that both the porosity and permeability are

329

significantly reduced by mineralisation.

19

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 31

330 331

Figure 8: Illustration of the streamlines change with the effects of mineralisation. (a) a digital

332

coal model with mineralisation; (b) and (c) shows the streamlines of models without and with

333

minerals in 3D.

334

3.2.2 Mechanical study of mineralised models

335

Since the permeability of fractured media is stress-dependent

60

, we compress the cleat

336

network of the digital coal model to study the sensitivity of permeability to external stresses. In

337

order to reduce computational time, a subsample (169×244×273) that has a well-developed cleat

338

network is cropped and investigated for an external loading test. In the subsample (Figure 9a),

339

there are two face cleats, one of which is filled with minerals, while butt cleats are not

340

mineralised. We numerically compress the subsample with a stress of 10Mpa in the X-direction

341

(normal to face cleats), while there is no pore pressure for this model. The resulting strain

342

distribution within the domain (Figure 9b) shows that the coal matrix near mineralised face cleat

343

has less strain than that of the non-mineralised face cleat. Thus, mineral fillings can mitigate the

344

geometrical deformation induced by confining stress.

345

The subsample is compressed in X- and Z-directions with external stresses ranging from 1

346

MPa to 10 MPa and the porosity of the resulting deformed models are measured (Figure 10). It is 20

ACS Paragon Plus Environment

Page 21 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

347

shown that both models with and without minerals become less porous under compression. The

348

porosity of non-mineralised subsamples presents a sharper decrease under increasing loads,

349

while porosity reduction of mineralised subsample is less significant. When applied stress is

350

greater than 5 MPa, models with mineral have greater porosity than non-mineralised models that

351

have the same cleat network. Figure 11 compares the resulting deformed geometries when the

352

subsample (Figure 11a) is subject to a confining stress of 4MPa in X- and Z-directions with and

353

without mineralisation. In Figure 11b, the non-mineralised face cleat disappears completely

354

under the compression, while the face cleat with minerals remains open. On the other hand, in

355

Figure 11c, the same face cleat is closed when there are no minerals. Therefore, mineral fillings

356

can aid in maintaining cleat aperture sizes while under confining stresses

357

less compressible than the coal matrix 22.

77

, since minerals are

358 359

Figure 9: (a) A subsample has well-developed cleat network, where minerals (in red) are

360

observed in one of the face cleats; (b) strain distribution under specified stress in the x direction

361

(as indicated by the arrows), where the colour of mineral (white) is not indication of strain.

21

ACS Paragon Plus Environment

Energy & Fuels

0.06 Models without minerals

0.05

Porosity

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 31

Models with minerals

0.04 0.03 0.02 0.01 0 1

362

3

5 7 Effective Stress, Mpa

9

363

Figure 10: The porosity variation of models with and without minerals as a function of external

364

stresses applied on surfaces. When applied stress is higher than 5 MPa, models with mineral

365

fillings inversely have greater porosity than non-mineralised one that has the identical cleat

366

structure.

367 368 369 370

Figure 11: The subsample (a) is compressed under the identical stresses (4Mpa), where the compression directions are indicated by the arrows. Resulting deformed geometries (b) and (c) are for mineralised and non-mineralised model, respectively. 22

ACS Paragon Plus Environment

Page 23 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

371

Energy & Fuels

3.2.3 Geometrical properties

372

Since the minerals in cleats tend to block fluid flow paths, the flow channels are more tortuous

373

compared with those in models without mineralisation. Figure 12 (a) plots the ratios of geodesic

374

Euclidean distances for one pair of digital models. The geodesic and Euclidean distances are

375

measured based on the surface at Z = 0 and another arbitrary parallel surface in the Z-direction. It

376

is shown that for models without minerals, the ratio values (in squares) are around 1, indicating

377

less tortuous flow channels. However, data points of the mineralised model (in circles) lie above

378

the data for models without minerals, which means that the fluid in the mineralised model travel

379

a longer pathway from inlet to outlet. Therefore, mineralisation can increase tortuosity by 21%

380

approximately. Tortuosity results of all generated models are shown in Figure 12b. Models

381

without minerals have an average tortuosity of 1.07. Conversely, after the inclusion of minerals,

382

the average tortuosity rises to 1.30 with larger variation. Compared with mineralised original

383

micro-CT images (τ = 1.32), the tortuosity estimated from mineralised models is more accurate

384

with an error of 1.5%.

385

Another geometrical property of the fluid flow channel that is influenced by mineralisation is

386

specific surface area (SSA). We find that the SSA of digital models is reduced from 30.27 cm-1

387

to 25.09cm-1 due to mineralisation. Therefore, mineralisation could have a detrimental influence

388

on the reaction rates by decreasing the contact area between adsorbed methane and coal matrix.

389

However, measured SSA of digital models is lower than that of the original micro-CT images

390

(66.00 cm-1 and 69.17 cm-1 for models with and without minerals, respectively), which is likely

391

induced by segmentation error where some regions are mislabelled as cleats. As a result, the

392

segmented micro-CT data will provide void space and SSA values that are greater than reality

393

(Figure 13b). Besides, isolated short cleats and dead-end pores of the original micro-CT images 23

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

30

Page 24 of 31

394

also result in larger internal surface area

. Since digital models are excluded of segmentation-

395

induced pores and isolated short cleats, so the estimated SSA values based on digital models are

396

more close to the area of the contact surface where fluids (e.g. methane and water) react with the

397

organic coal matrix.

398 399

Figure 12: Comparison of tortuosity between models with and without minerals. (a) the ratios of

400

geodesic distance and Euclidean distance for one pair of digital model; (b) tortuosity values of all

401

digital models and mineralised original micro-CT images.

402 24

ACS Paragon Plus Environment

Page 25 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

403

Figure 13: (a) is one slice of original grey images, where noises (in circle) have similar grey

404

value with cleats. As a result, the noises are mistakenly segmented into cleats in segmented

405

image (b).

406

4 Conclusions

407

This work develops a mineralised digital coal model that is able to characterise the core-scale

408

components of coal, including coal matrix, cleats and minerals. The developed digital model

409

comprises of alternating bright and dull bands such that the banding morphology of different

410

lithotypes are preserved. Then, a novel discrete fracture network model is integrated with the

411

banded model where geometrical properties and the particular connectivity pattern of the cleat

412

network are characterised, resulting in a digital coal model without mineralisation. Based on

413

micro-CT images of a coal sample from Moura mine of Bowen Basin, the mineralisation is

414

characterised regarding the mineral frequency, the volume fraction of mineral fillings and the

415

mineral length. Stochastic minerals are further generated within the digital coal model, with the

416

identical degree of mineralisation as the original micro-CT images. By computing the

417

petrophysical and morphological properties of the digital coal models, it is found that

418

mineralisation can significantly influence coal petrophysical properties. For instance, mineral

419

fillings act as barriers to fluid migration, resulting in less permeable coal and more tortuous flow

420

paths. The contact area between organic coal matrix and gas is reduced by 17% due to

421

mineralization, which could influence reaction rates. In addition, the less compressible minerals

422

can aid in keeping cleats open and thus reduce coal matrix deformation under confining pressure.

423

Therefore, mineralisation is crucial in the characterisation of coal since it not only retards cleat 25

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 31

424

network transmissibility but also supports coal cleat apertures sizes under external stresses.

425

Furthermore, the developed digital coal models that characterise mineralisation, are found to be

426

more representative in terms of permeability estimation. So the construction of mineralised

427

digital coal models can lead to more accurate evaluation of coal samples and thus provide more

428

reliable parameters for the prediction of coal seam gas production rates.

429

Reference

430 431 432

(1) Y. Yao, D. Liu, Y. Che, D. Tang, S. Tang, and W. Huang, Non-destructive characterization of coal samples from China using microfocus X-ray computed tomography. International Journal of Coal Geology, 2009. 80(2): p. 113-123.

433 434

(2) C.R. Ward, Analysis, origin and significance of mineral matter in coal: An updated review. International Journal of Coal Geology, 2016. 165: p. 1-27.

435

(3) R.R. Dutcher, Field Description Of Goal. Vol. 661. 1978: Astm International.

436 437 438

(4) S.R. Dindarloo, A. Bagherieh, J.C. Hower, J.H. Calder, and N.J. Wagner, Coal modeling using Markov Chain and Monte Carlo simulation: Analysis of microlithotype and lithotype succession. Sedimentary Geology, 2015. 329: p. 1-11.

439 440 441

(5) W. Riese, W.L. Pelzmann, and G.T. Snyder, New insights on the hydrocarbon system of the Fruitland Formation coal beds, northern San Juan Basin, Colorado and New Mexico, USA. Geological Society of America Special Papers, 2005. 387: p. 73-111.

442

(6) R.M. Flores, Coal and coalbed gas: fueling the future. 2013: Newnes.

443 444

(7) P.D. Gamson, B.B. Beamish, and D.P. Johnson, Coal microstructure and micropermeability and their effects on natural gas recovery. Fuel, 1993. 72(1): p. 87-99.

445 446 447

(8) C. Pattison, C. Fielding, R. McWatters, and L. Hamilton, Nature and origin of fractures in Permian coals from the Bowen Basin, Queensland, Australia. Geological Society, London, Special Publications, 1996. 109(1): p. 133-150.

448 449 450

(9) P. Gamson and B. Beamish. Coal type, microstructure and gas flow behaviour of Bowen Basin coals. in Symposium on coalbed methane research and development in Australia, Beamish, BB and Gamson, PD (eds.), James cook University of North Queensland. 1992.

451 452

(10) F. Gao, D. Stead, and H. Kang, Numerical investigation of the scale effect and anisotropy in the strength and deformability of coal. International Journal of Coal Geology, 2014. 136: p. 25-37.

26

ACS Paragon Plus Environment

Page 27 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

453 454

(11) L. Scholtès, F.-V. Donzé, and M. Khanal, Scale effects on strength of geomaterials, case study: coal. Journal of the Mechanics and Physics of Solids, 2011. 59(5): p. 1131-1146.

455 456 457

(12) J. Maillot, P. Davy, R. Le Goc, C. Darcel, J.-R. de Dreuzy, and M. Stigsson, Comparison between “Poissonian” and “mechanically-oriented” DFN models for the prediction of permeability and flow channeling, in DFNE 2014. 2014: Vancouver, Canada.

458 459

(13) G. Dawson and J. Esterle, Controls on coal cleat spacing. International Journal of Coal Geology, 2010. 82(3): p. 213-218.

460 461 462

(14) J. Esterle, Y. Kolatschek, and G. O'Brien, Relationship between in situ coal stratigraphy and particle size and composition after breakage in bituminous coals. International journal of coal geology, 2002. 49(2): p. 195-214.

463 464

(15) V.N. Misra, P. Reddy, and B. Mohapatra, Mineral Characterisation and Processing. 2004: Allied Publishers.

465 466

(16) C.R. Ward, Analysis and significance of mineral matter in coal seams. International Journal of Coal Geology, 2002. 50(1–4): p. 135-168.

467 468

(17) P.K. Singh and A. Naik, Coal microscopy as a tool to understand the beneficiation problems of Indian Gondwana coals. 2015.

469 470 471

(18) S. Dai, D. Ren, C.-L. Chou, S. Li, and Y. Jiang, Mineralogy and geochemistry of the No. 6 Coal (Pennsylvanian) in the Junger Coalfield, Ordos Basin, China. International Journal of Coal Geology, 2006. 66(4): p. 253-270.

472 473

(19) E. Daniels, S. Marshak, and S. Altaner, Use of clay-mineral alteration patterns to define syntectonic permeability of joints (cleat) in Pennsylvania anthracite coal. Tectonophysics, 1996. 263(1): p. 123-136.

474 475 476

(20) G. Dawson, S. Golding, J. Esterle, and P. Massarotto, Occurrence of minerals within fractures and matrix of selected Bowen and Ruhr Basin coals. International Journal of Coal Geology, 2012. 94: p. 150166.

477 478

(21) C.Ö. Karacan, Swelling-induced volumetric strains internal to a stressed coal associated with CO 2 sorption. International Journal of Coal Geology, 2007. 72(3): p. 209-220.

479 480 481

(22) Y. Zhang, X. Xu, M. Lebedev, M. Sarmadivaleh, A. Barifcani, and S. Iglauer, Multi-scale x-ray computed tomography analysis of coal microstructure and permeability changes as a function of effective stress. International Journal of Coal Geology, 2016. 165: p. 149-156.

482 483 484

(23) C. Deng, D. Tang, S. Liu, H. Xu, and S. Tao, Characterization of mineral composition and its influence on microstructure and sorption capacity of coal. Journal of Natural Gas Science and Engineering, 2015. 25: p. 46-57.

485 486 487

(24) F. Verhelst, P. David, W. Fermont, L. Jegers, and A. Vervoort, Correlation of 3D-computerized tomographic scans and 2D-colour image analysis of Westphalian coal by means of multivariate statistics. International Journal of Coal Geology, 1996. 29(1–3): p. 1-21.

27

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 31

488 489 490

(25) F.J. Simons, F.d.r. Verhelst, and R. Swennen, Quantitative characterization of coal by means of microfocal X-ray computed microtomography (CMT) and color image analysis (CIA). International Journal of Coal Geology, 1997. 34(1): p. 69-88.

491 492 493

(26) M. Van Geet and R. Swennen, Quantitative 3D‐fracture analysis by means of microfocus X‐Ray Computer Tomography (µCT): An example from coal. Geophysical Research Letters, 2001. 28(17): p. 3333-3336.

494 495

(27) C.O. Karacan and E. Okandan, Adsorption and gas transport in coal microstructure: investigation and evaluation by quantitative X-ray CT imaging. Fuel, 2001. 80(4): p. 509-520.

496 497 498

(28) S. Mazumder, K.-H. Wolf, K. Elewaut, and R. Ephraim, Application of X-ray computed tomography for analyzing cleat spacing and cleat aperture in coal samples. International Journal of Coal Geology, 2006. 68(3): p. 205-222.

499 500

(29) Y. Jing, R.T. Armstrong, H.L. Ramandi, and P. Mostaghimi, Coal cleat reconstruction using microcomputed tomography imaging. Fuel, 2016. 181: p. 286-299.

501 502

(30) Y. Jing, R.T. Armstrong, and P. Mostaghimi, Rough-walled discrete fracture network modelling for coal characterisation. Fuel, 2017. 191: p. 442-453.

503 504

(31) A.R. Scott, Hydrogeologic factors affecting gas content distribution in coal beds. International Journal of Coal Geology, 2002. 50(1–4): p. 363-387.

505 506

(32) S.E. Laubach, J.E. Olson, and J.F.W. Gale, Are open fractures necessarily aligned with maximum horizontal stress? Earth and Planetary Science Letters, 2004. 222(1): p. 191-195.

507 508

(33) S.E. Laubach, Practical approaches to identifying sealed and open fractures. AAPG bulletin, 2003. 87(4): p. 561-579.

509 510

(34) P. Guo and Y. Cheng, Permeability prediction in deep coal seam: a case study on the no. 3 coal seam of the Southern Qinshui Basin in China. The Scientific World Journal, 2013. 2013.

511 512 513

(35) H.L. Ramandi, P. Mostaghimi, R.T. Armstrong, M. Saadatfar, and W.V. Pinczewski, Porosity and permeability characterization of coal: a micro-computed tomography study. International Journal of Coal Geology, 2016. 154–155: p. 57-68.

514 515 516

(36) T. Varslot, A. Kingston, G. Myers, and A. Sheppard, High-resolution helical cone-beam micro-CT with theoretically-exact reconstruction from experimental data. Medical physics, 2011. 38(10): p. 54595476.

517 518 519

(37) P. Mostaghimi, R.T. Armstrong, A. Gerami, Y. Hu, Y. Jing, F. Kamali, M. Liu, Z. Liu, X. Lu, and H.L. Ramandi, Cleat-scale characterisation of coal: An overview. Journal of Natural Gas Science and Engineering, 2017. 39: p. 143-160.

520 521 522

(38) P. Iassonov, T. Gebrenegus, and M. Tuller, Segmentation of X‐ray computed tomography images of porous materials: A crucial step for characterization and quantitative analysis of pore structures. Water Resources Research, 2009. 45(9).

523 524

(39) M.J. Blunt, B. Bijeljic, H. Dong, O. Gharbi, S. Iglauer, P. Mostaghimi, A. Paluszny, and C. Pentland, Pore-scale imaging and modelling. Advances in Water Resources, 2013. 51: p. 197-216.

28

ACS Paragon Plus Environment

Page 29 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

525 526 527

(40) D. Wildenschild, C. Vaz, M. Rivers, D. Rikard, and B. Christensen, Using X-ray computed tomography in hydrology: systems, resolutions, and limitations. Journal of Hydrology, 2002. 267(3): p. 285-297.

528 529 530

(41) M. Van Geet, R. Swennen, and P. David, Quantitative coal characterisation by means of microfocus X-ray computer tomography, colour image analysis and back-scattered scanning electron microscopy. International Journal of Coal Geology, 2001. 46(1): p. 11-25.

531 532 533

(42) A.P. Sheppard, R.M. Sok, and H. Averdunk, Techniques for image enhancement and segmentation of tomographic images of porous materials. Physica A: Statistical mechanics and its applications, 2004. 339(1): p. 145-151.

534 535

(43) S. Schlüter, A. Sheppard, K. Brown, and D. Wildenschild, Image processing of multiphase images obtained via X‐ray microtomography: a review. Water Resources Research, 2014. 50(4): p. 3615-3639.

536 537

(44) I. Sobel, An isotropic 3× 3 image gradient operator. Machine Vision for three-demensional Sciences, 1990.

538 539

(45) S. Laubach, R. Marrett, J. Olson, and A. Scott, Characteristics and origins of coal cleat: a review. International Journal of Coal Geology, 1998. 35(1): p. 175-207.

540 541

(46) T.-C. Lee, R.L. Kashyap, and C.-N. Chu, Building skeleton models via 3-D medial surface/axis thinning algorithms. CVGIP: Graph. Models Image Process., 1994. 56(6): p. 462-478.

542 543 544

(47) J.C. Close and M.J. Mavor, Influence of Coal Composition and Rank on Fracture Development in Fruitland Coal Gas Reservoirs of San Juan Basin, in Assoc. Geol. Field Conf. 1991: Coalbed Methane of Western North America, Rk. Mt. p. 109-121.

545 546 547

(48) B. Law. The relationship between coal rank and cleat spacing: implications for the prediction of permeability in coal. in Proceedings of the 1993 International Coalbed Methane Symposium. 1993. Twscaloosa, ALabma: University of ALabama.

548 549 550

(49) C.M. Tremain, S.E. Laubach, and N.H. Whitehead III, Coal fracture (cleat) patterns in Upper Cretaceous Fruitland formation, San Juan Basin, Colorado and New Mexico-implications for coalbed methane exploration and development. 1991.

551 552

(50) Y. Jing, R.T. Armstrong, and P. Mostaghimi, Digital coal: Generation of fractured cores with microscale features. Fuel, 2017. 207: p. 93-101.

553 554

(51) Z. Zhao, B. Li, and Y. Jiang, Effects of Fracture Surface Roughness on Macroscopic Fluid Flow and Solute Transport in Fracture Networks. Rock Mechanics and Rock Engineering, 2014. 47(6): p. 2279.

555 556

(52) G. Baecher, N. Lanney, and H. Einstein, Statistical description of rock properties and sampling, in The 18th US Symposium on Rock Mechanics (USRMS). 1977, American Rock Mechanics Association.

557 558

(53) A. Robertson, The interpretation of geological factors for use in slope theory, in Planning Open Pit Mines, Proceedings, Johannesburg. 1970. p. 55-71.

559 560 561

(54) M. Cacas, E. Ledoux, G.d. Marsily, B. Tillie, A. Barbreau, E. Durand, B. Feuga, and P. Peaudecerf, Modeling fracture flow with a stochastic discrete fracture network: Calibration and validation: 1. The flow model. Water Resources Research, 1990. 26(3): p. 479-489.

29

ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 31

562 563 564

(55) S. Rief, D. Kehrwald, K. SCHMDT, and A. Wiegmann, Fraunhofer Software Tools GeoDict/FilterDict for the Simulation of Diesel Particulate Filters. NAFEM: Reliable Use of Numerical Methods in Upfront Simulations. Wiesbaden, Germany, 2007.

565 566

(56) A. Wiegmann, Computation of the permeability of porous materials from their microstructure by FFF-Stokes. 2007: Fraunhofer-Institut für Techno-und Wirtschaftsmathematik, Fraunhofer (ITWM).

567 568

(57) D. McBride, J. Gebhardt, and M. Cross, A comprehensive gold oxide heap leach model: Development and validation. Hydrometallurgy, 2012. 113: p. 98-108.

569 570 571

(58) P. Mostaghimi, B.S. Tollit, S.J. Neethling, G.J. Gorman, and C.C. Pain, A control volume finite element method for adaptive mesh simulation of flow in heap leaching. Journal of Engineering Mathematics, 2014. 87(1): p. 111-121.

572 573

(59) M. Liu and P. Mostaghimi. Coupling lattice Boltzmann and finite volume methods for simulation of reactive transport. in International Symposium of the Society of Core Analysts, Colorado, USA. 2016.

574 575 576

(60) K.-B. Min, J. Rutqvist, C.-F. Tsang, and L. Jing, Stress-dependent permeability of fractured rock masses: a numerical study. International Journal of Rock Mechanics and Mining Sciences, 2004. 41(7): p. 1191-1210.

577 578

(61) V. Rutka, H. Andrä, and A. Wiegmann, EJIIM for calculation of effective elastic moduli in 3D linear elasticity. Vol. 1. 2006: ITWM.

579 580

(62) M. Liu and P. Mostaghimi, Pore-scale simulation of dissolution-induced variations in rock mechanical properties. International Journal of Heat and Mass Transfer, 2017. 111: p. 842-851.

581 582 583

(63) A. Wiegmann, S. Rief, and A. Latz, Geodict and filterdict: Software for the virtual material design of new filter media. Proc. New Developments in Filtration Technology (Loughborough, Angleterre), 2006.

584 585

(64) N. Aziz, B. Kininmonth, J. Nemcik, T. Ren, and J. Hoelle, Proceedings of the 2013 Coal Operators' Conference. 2013.

586 587

(65) H. Neuber, Theory of stress concentration for shear-strained prismatical bodies with arbitrary nonlinear stress-strain law. Journal of Applied Mechanics, 1961. 28(4): p. 544-550.

588 589

(66) J. Thomas and H.H. Damberger, Internal surface area, moisture content, and porosity of Illinois coals: Variations with coal rank. Circular no. 493, 1976.

590

(67) J.R. Parker, Algorithms for image processing and computer vision. 2010: John Wiley & Sons.

591 592

(68) M.B. Clennell, Tortuosity: a guide through the maze. Geological Society, London, Special Publications, 1997. 122(1): p. 299-344.

593 594

(69) E. Fjar, R.M. Holt, A. Raaen, R. Risnes, and P. Horsrud, Petroleum related rock mechanics. Vol. 53. 2008: Elsevier.

595 596 597

(70) C.J. Gommes, A.J. Bons, S. Blacher, J.H. Dunsmuir, and A.H. Tsou, Practical methods for measuring the tortuosity of porous materials from binary or gray‐tone tomographic reconstructions. AIChE Journal, 2009. 55(8): p. 2000-2012.

30

ACS Paragon Plus Environment

Page 31 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

598 599

(71) E.F. Krause, Taxicab geometry: An adventure in non-Euclidean geometry. 2012: Courier Corporation.

600 601

(72) P. Soille, Morphological image analysis: principles and applications. 2013: Springer Science & Business Media.

602 603

(73) B. Ryan and N.V. Branch, Cleat development in some British Columbia coals. Geological fieldwork, 2002: p. 2003-1.

604 605 606

(74) B.S. Faraj, C.R. Fielding, and I.D. Mackinnon, Cleat mineralization of upper permian Baralaba/Rangal coal measures, Bowen Basin, Australia. Geological Society, London, Special Publications, 1996. 109(1): p. 151-164.

607 608 609

(75) J.C. Hower, D.A. Williams, C.F. Eble, T. Sakulpitakphon, and D.P. Moecher, Brecciated and mineralized coals in Union County, Western Kentucky coal field. International Journal of Coal Geology, 2001. 47(3): p. 223-234.

610 611 612

(76) S.H. Chen, X.M. Feng, and S. Isam, Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method. International journal for numerical and analytical methods in geomechanics, 2008. 32(12): p. 1459-1477.

613 614 615

(77) H. Kumar, E. Lester, S. Kingman, R. Bourne, C. Avila, A. Jones, J. Robinson, P.M. Halleck, and J.P. Mathews, Inducing fractures and increasing cleat apertures in a bituminous coal under isotropic stress via application of microwave energy. International Journal of Coal Geology, 2011. 88(1): p. 75-82.

616

31

ACS Paragon Plus Environment