# Applications of Vitrinite Anisotropy in the Tectonics: A Case Study of Huaibei Coalfield, Southern North China

Authors: Yudong Wu; Quanlin Hou; Yiwen Ju; Daiyong Cao; Junjia Fan; Wei Wei
Journal: Journal of Geological Research
 (2012)

---

## Abstract

29 oriented and 10 nonoriented coal samples are collected in the study from three different regions of the Huaibei coalfield, eastern China, and their vitrinite reflectance indicating surface (RIS) parameters are systematically calculated and analyzed. Using the available methods, Kilby’s transformations and RIS triaxial orientations are obtained. The magnitudes and orientations of the RIS axes of the three regions were respectively projected on the horizontal planes and vertical sections. The results show that the samples in high deformed region have significant anisotropy magnitudes (higherBw/Rmax and Ram values) with a biaxial negative style, whereas the samples in the slightly deformed area have unimpressive anisotropy magnitudes with a biaxial negative style. Thermal metamorphism superposed might enhance the complication and variation of RIS style. RIS projection analysis deduced that the RIS orientation is mainly controlled by regional tectonic stress, and likely influenced by deformation mechanisms of coal.


---

## Body

## 1. Introduction

During the coalification, tectonic-thermal events could influence the vitrinite reflectance of coal [1, 2]. Effects of temperature and stress on modifying vitrinite reflectance anisotropy (VRA) were proved by high temperature-pressure experiments on bituminous anthracite [3–13]. Patterns of VRA [14, 15] in coal are usually represented by a three-dimensional graphical form resembling an ellipsoid and termed as RIS (reflectance indicating surface).The investigation of the RIS for vitrinite in coals is performed by means of sections cut from oriented samples, and the magnitudes and orientations of the ellipsoid parameters are obtained by calculating the data measured from three mutually perpendicular surfaces. This method has been reported [14] and verified by high temperature-pressure experiments [16] and tectonic stress field analysis [11, 15, 17–20]. The anisotropy of RIS in this method is quantified by a function of rank [15], including in terms of bireflectance ratio [21], anisotropy ratio [17] and bireflectance [22].The RIS styles may also be obtained by nonoriented sections from crushed coal samples. The method requires measurement of the maximum and minimum apparent reflectance from a series of randomly oriented vitrinite particles, and the magnitudes of RIS principal reflectance axes were estimated by using the reflectance crossplots [23, 24]. A series of parameters named as Kilby’s transformations in the method are proposed more effectively to estimate the anisotropy magnitude and shape of RIS ellipsoid, including RIS style (Rst), anisotropy magnitude (Ram), and reflectance of equivalent volume isotropic RIS (Rev). Compared with Levine’s method, basing on the application of vitrinite maximum reflectance versus vitrinite random reflectance [25], less time consuming is found during Kilby’s data processing [13, 23, 26–30]. In the nonoriented sections, however, the orientations of three principal reflectance axes are missing.Previous works suggested that the anisotropic strains developed by lithostatic pressure, especially by tectonic strains, including oriented pressures and shears. High rank coals and anthracites show evidence of the presence of biaxial negative materials [15, 19, 23, 31, 32]. As the result of tectonic dislocations, in many cases associated with fractures in the plane normal to the coal bedding, it is possible to find microtextures responsible for the biaxial positive character [29, 30]. It is found that minimum reflectance develops incrementally parallel to the direction of maximum compressive stress (σ1) during coalification, and maximum reflectance develops incrementally parallel to the direction of minimum compressive stress (σ3) [1, 14]. These viewpoints are verified by simple shear experiments in high temperature-pressure environment [28, 33, 34]. On this basis, the RIS may approximate to finite strain ellipsoid [35], and the orientations and magnitudes of RIS principal reflectance axes may correspondingly have further usages in the tectonic stress field analysis.The Huaibei Coalfield in Anhui Province of China has abundant coal resources. Thermal model research showed that the coal experienced low temperature (<230°C) after being buried during the coalification process [36] and, therefore, thermal function was limited in the coalification history. Although the coalification process was influenced by other geological agencies such as tectonic stress, the RIS parameters may provide semiquantitative information about the magnitudes and orientations of tectonic deformation.In this paper, 29 oriented coal samples from the working face of the underground mines and 10 nonoriented coal samples from the core drilling were collected. In order to estimate the magnitudes and orientations of the three principal reflectance axes, the RIS was built on the oriented coal sample blocks. The orientation of RIS was obtained by calculating the data measured from three mutually perpendicular surfaces. Meanwhile, the anisotropy magnitude and shape were analyzed with the parameters such as bireflectance ratio, Flinn’s diagram, and Kilby’s transformations. The effects of tectonic stress and thermal evolution on coal are therefore expected to be evaluated.


## 2. Geology of the Study Area

The Huaibei Coalfield in the southeast area of North China Craton was influenced by the Tancheng–Lujiang sinistral strike-slip fault in Mesozoic [39–42] and a related arcuate thrust system [37]. Many researchers believed that the arcuate thrust system was a west-thrust imbricate fan with the arc-top near to Huaibei City [37, 43]. A syncline formed in the margin of the thrust system had a certain influence on the burial depth of coal seams and the level of coal rank [44]. Coal-bearing strata mainly occur in Permo-Carboniferous age.The Huaibei Coalfield is divided into three regions representing different stress-thermal environments. Vitrinite reflectance of coal samples collected from these regions was measured and sampling locations are shown in Figure1(a).(1)
Region A, including the coal mines of Taoyuan (label 3), Qinan (label 4), Qidong (label 5), and Xutuan (label 9). The strata were likely influenced less by magmatism and tectonic agency and more by burial metamorphism.(2)
Region B, including the coal mines of Haizi (label 6), Linhuan (label 7), Tongting (label 8) and Baishan (label 10). Both Mesozoic deformation and magmatism affected deeply the strata and likely impacted on the samples collected in the region.(3)
Region C, including the coal mines of Luling (label 1), Zhuxianzhuang (label 2), Shitai (label 11) and Yangzhuang (label 12). The strata in the region strongly deformed, and samples are probably highly influenced by effects of tectonic deformation and thermal-shearing.(4)
The lithostratigraphic sequence of the study area is shown in Figure1(b). The major coal seams occur in the lower Permian system (P1) and minor coal seams locally occur in the upper Permian System (P2). Southward marine regression during early Permian time led to deposition in an island-lagoon environment and followed by terrestrial delta-plain deposits that contained extensive and minable coal seams [44]. The deltaic environment became unstable in early Late Permian, and coal seams became thinned, discontinuous and unminable. The alluvial deposits in Late Permian marked the end of coal accumulation in the Huaibei coalfield [44]. All the oriented samples collected from the working face in coal seams of 3rd, 7th, 8th and 10th and the nonoriented samples collected from the well-drills belonging to the Permian strata.Simplified geologic map with late Paleozoic stratigraphic column of the Huaibei Coalfield. Coal mine label: (1) Luling mine, (2) Zhuxianzhuang mine, (3) Taoyuan mine, (4) Qinan mine, (5) Qidong mine, (6) Haizi mine, (7) Linhuan mine, (8) Tongting mine, (9) Xutuan mine, (10) Baishan mine, (11) Yangzhuang mine, (12) Shitai mine, The location is shown in the inset of Figure1(a): Areas in gray where Upper Permian coal bearing strata are preserved in the Huaibei coalfield. Oriented coal samples were gathered from the Luling (LL), Taoyuan (TY), Haizi (HZ), Linhuan (LH), Xutuan (XT), Shitai (ST), and Yangzhuang (YZ) mines. The nonoriented coal samples were collected from bore holes in TY, XT, LH and Qidong (QD) coal mines. As shown in Figure 1(b), the Mesozoic and lower Cenozoic strata were suffered from erosion and missed. The Permian strata mainly consist of sandy mudstone and fine sandstone with 10 major coal seams (Nos. 1~10). Geological data are compiled based on [37, 38] the stratigraphic column is modified from China National Administration of coal Geology (1996).
(a)
(b)


## 3. Experimental Section

Coal macerals of samples were determined by spot checks before the reflectance measurement, and the reflectance was measured by means of the method described in GB/T 8899-1998 (equated with ISO7404/3 standard). At least 500 points for each sample were measured on the same polished specimen. The contents of vitrinite, inertinite, liptinite, and mineral matter were obtained and are shown in Table1. Coal lithotype in the Huaibei coalfield is dominated by clarain, with secondary durian and vitrain, and coal maceral is dominated by vitrinite (40.6~92.5%), which is followed by inertinite and liptinite (Table 1).Table 1
Maceral composition of coal samples studied.
Sample number
V (%)
I (%)
E (%)
MM (%)
TYM03
40.6
22.5
36.7
0.2
TY02
66.4
18.4
14.4
0.8
XTM08
89.2
9.2
0.8
0.8
LHM03
67.4
31.7
x
0.9
LHM02
81.9
17.4
x
0.7
LHM10
76.8
22.8
x
0.4
HZM03
92.5
7.0
x
0.5
LLM05
58.2
24.2
15.2
2.4
LLM08
65.0
17.9
16.3
0.8
LLM01
67.8
18.3
12.4
1.5
YZM12
67.3
32.6
x
0.1
STM05
67.7
30.8
x
1.5
See Figure1 for sample location. Abbreviations: V→vitrinite content (vol.%); I→inertinite content (vol.%); E→Exinite content (vol.%); MM→mineral matter content (vol.%). x = content < 0.1%.


### 3.1. Sample Preparation

The 29 oriented samples were prepared following the rules described in GB/T 16773-2008 (equated with ISO7404/2 standard) and cut into cubes with approximately 3-4 cm in dimensions (Figure2). One pair of planes are paralleling to the horizontal plane and the other two mutually perpendicular pairs normal to the horizontal one. Each sample, therefore, has three pairs of mutually perpendicular surfaces (a/a′, b/b′, and c/c′ as shown in Figure 2) under the process of polishing surface.Figure 2
Diagram of three sections cut from oriented blocks. Arrowhead on the top surface of the block shows the spatial orientation. Surfaces marked witha/a′, b/b′, c/c′ were mutually perpendicular with the surfaces a and a′ horizontal and c and c′ vertical and parallel to the arrowhead.The 11 nonoriented samples were prepared by a modified procedure of GB/T 16773-2008 (equated with ISO 7404/02 standard). These samples were also used in our previous work for thermal simulation of maturation for coal [36].


### 3.2. Optical Microscopic Measurements

The measurements were performed for each oriented sample using an MPV-3 Combi (Leitz) microscope, reflected light, and oil immersion objectives. The maximum and minimum apparent reflectance (Rmax and Rmin) on each polished surface were measured in polarized light (λ=546μ) and by rotating the microscope stage through 360° and following the analytical method described in GB/T 6948-2008 (equated with ISO7404/5 standard). at least 20 points on each surface of each sample were measured to obtain the mean statistical values and azimuths of Rmax and Rmin.As for nonoriented samples, at least 20 points were measured with the same analytical method, and the mean statistical values ofRmax and Rmin were obtained.


### 3.3. Data Processing

On the basis of measured mean values, the values and occurrences of maximum reflectance (Rmax), intermediate reflectance (Rint), and minimum reflectance (Rmin) axes may be obtained by solving ellipsoid equation in analytic geometry of space [35], which is the magnitudes and orientations of RIS. The computing method is in Figure 3 as follows:
(1)λx′x2+λy′y2+λz′z2-2γxy′xy-2γyz′yz-2γxz′zx=1,
where
(2)λx′=1λx,λy′=1λy,λz′=1λz,γ′=γλ=γλ′,λx′=cos2θ¯x(R¯max)xy2+sin2θ¯x(R¯min)xy2,λy′=sin2θ¯x(R¯max)xy2+cos2θ¯x(R¯min)xy2,γxy′=[1(R¯max)xy2-1(R¯min)xy2]*sinθ¯x*cosθ¯x,λy′=cos2θ¯y(R¯max)yz2+sin2θ¯y(R¯min)yz2,λz′=sin2θ¯y(R¯max)yz2+cos2θ¯z(R¯min)yz2,γyz′=[1(R¯max)yz2-1(R¯min)yz2]*sinθ¯y*cosθ¯y,λz′=cos2θ¯z(R¯max)zx2+sin2θ¯z(R¯min)zx2,λx′=sin2θ¯z(R¯max)zx2+cos2θ¯z(R¯min)zx2,γzx′=[1(R¯max)zx2-1(R¯min)zx2]*sinθ¯z*cosθ¯z.Coordinates and coal brick VR ellipses [47]. A coal sample brick was demonstrated in Figure 3(a), the measuring result were quantitatively represented in the YXZ coordinate system. The value and orientation of Rmax/Rmin on each polished surface were shown on three mutually perpendicular cube face, represented by the long/minor axis of ellipse. For example, as shown in Figure 3(b), the measured value (RXOY·max and RXOY·min) and their orientations (θx) are expressed as their coordinates and the RIS ellipsoid equation can be built.
(a)
(b)This ellipsoid equation has three eigenvalues, corresponding with the values of three RIS ellipsoid axes (Rmax>Rint>Rmin):
(3)1+ei=λi;i=1,2,3.The occurrence of three RIS ellipsoid axis is calculated according to the formula(4)Ai=arctanMiLi,Bi=arccosNi,i=1,2,3,
where Ai and Bi represent the direction and dip angle of the axis in the coordinate system, respectively. 								(5)Li=cosRix=Xiλi;Mi=cosRiy=Yiλi,Ni=cosRiz=Ziλi,i=1,2,3,Xi+Yi+Zi=λi,Xiλy′λz′-λi′λy′-λi′λz′+λi'2-γyz'2=-Yi-γxy′λz′+γxy′λi′-γzx′γyz′=-Ziγxy′γyz′+γzx′λy′-γzx′λi′.The geological occurrence of RIS ellipsoid axis (Diri∠Dipi), however, needs a data conversion as follows:
(6)ifBi≥0,Diri=Ci-Ai±360°,Dipi=Bi,ifBi<0,Diri=Ci-Ai±180°,Dipi=-Bi,
where 0<Diri<360°, 0<Dipi<90°; where Ci represents the geological direction of arrowhead in Figure 2.Therefore, the reflectance anisotropy (bireflectanceBw and Ram coefficient values) and RIS shapes (such as Rev, Rst coefficient value, Flinn’s parameters, etc.) can be determined. The details of the parameters are given as follows.(1)
Bw coefficient value [22], representing bireflectance value, is calculated according to the formula. 											(7)Bw=Rmax-Rmin.Bw/Rmax [21] represents the bireflectance ratio.(2)
Flinn’s parameter was firstly proposed for examining the geometry of three-dimensional homogeneous strain [35, 45]. Formula of Ramsay was chosen here
(8)K=ε1-ε2ε2-ε3,where(9)ε1=ln(Rmax),ε2=ln(Rint),ε3=ln(Rmin).(3)
Ram coefficient value [23, 24], representing the RIS anisotropy magnitude, is calculated according to the formula
(10)Ram=(x2-y2)1/2,where(11)x=RintRmax+Rint+Rmin+[y+13cos30°]tan30°-12,y=[Rmax(Rmax+Rint+Rmin)-13]cos30°.This value is the distance between the plotted position of RIS and the position of an isotropic RIS, a single measure of anisotropy for all styles of RIS [23].Isotropic constituent of coal is characterized by Ram=0. The higher coefficient value is, the stronger optical anisotropy will be. Typical Ram values for raw vitrinites are from 0.03 to 0.05 and the value of 0.1 represents strong anisotropy of a constituent.(4)
R
st coefficient value [23, 24], representing RIS style, is calculated according to the formula
(12)Rst=30-arctan(xy),
wherex and y are determined in the same way as shown above. Rst values range from (−30) to (+30) and the RIS shape is described as
R
st
=
(
-
30
): uniaxial negative; Rst=(+30): uniaxial positive; -30<Rst<0: biaxial negative; 0<Rst<+30: biaxial positive; Rst=0: biaxial neutral.
(5)
R
ev coefficient value [23, 24], representing the reflectance of equivalent volume of isotropic RIS, is calculated according to the formula
(13)Rev=(Rmax×Rint×Rmin)1/3.


## 3.1. Sample Preparation

The 29 oriented samples were prepared following the rules described in GB/T 16773-2008 (equated with ISO7404/2 standard) and cut into cubes with approximately 3-4 cm in dimensions (Figure2). One pair of planes are paralleling to the horizontal plane and the other two mutually perpendicular pairs normal to the horizontal one. Each sample, therefore, has three pairs of mutually perpendicular surfaces (a/a′, b/b′, and c/c′ as shown in Figure 2) under the process of polishing surface.Figure 2
Diagram of three sections cut from oriented blocks. Arrowhead on the top surface of the block shows the spatial orientation. Surfaces marked witha/a′, b/b′, c/c′ were mutually perpendicular with the surfaces a and a′ horizontal and c and c′ vertical and parallel to the arrowhead.The 11 nonoriented samples were prepared by a modified procedure of GB/T 16773-2008 (equated with ISO 7404/02 standard). These samples were also used in our previous work for thermal simulation of maturation for coal [36].


## 3.2. Optical Microscopic Measurements

The measurements were performed for each oriented sample using an MPV-3 Combi (Leitz) microscope, reflected light, and oil immersion objectives. The maximum and minimum apparent reflectance (Rmax and Rmin) on each polished surface were measured in polarized light (λ=546μ) and by rotating the microscope stage through 360° and following the analytical method described in GB/T 6948-2008 (equated with ISO7404/5 standard). at least 20 points on each surface of each sample were measured to obtain the mean statistical values and azimuths of Rmax and Rmin.As for nonoriented samples, at least 20 points were measured with the same analytical method, and the mean statistical values ofRmax and Rmin were obtained.


## 3.3. Data Processing

On the basis of measured mean values, the values and occurrences of maximum reflectance (Rmax), intermediate reflectance (Rint), and minimum reflectance (Rmin) axes may be obtained by solving ellipsoid equation in analytic geometry of space [35], which is the magnitudes and orientations of RIS. The computing method is in Figure 3 as follows:
(1)λx′x2+λy′y2+λz′z2-2γxy′xy-2γyz′yz-2γxz′zx=1,
where
(2)λx′=1λx,λy′=1λy,λz′=1λz,γ′=γλ=γλ′,λx′=cos2θ¯x(R¯max)xy2+sin2θ¯x(R¯min)xy2,λy′=sin2θ¯x(R¯max)xy2+cos2θ¯x(R¯min)xy2,γxy′=[1(R¯max)xy2-1(R¯min)xy2]*sinθ¯x*cosθ¯x,λy′=cos2θ¯y(R¯max)yz2+sin2θ¯y(R¯min)yz2,λz′=sin2θ¯y(R¯max)yz2+cos2θ¯z(R¯min)yz2,γyz′=[1(R¯max)yz2-1(R¯min)yz2]*sinθ¯y*cosθ¯y,λz′=cos2θ¯z(R¯max)zx2+sin2θ¯z(R¯min)zx2,λx′=sin2θ¯z(R¯max)zx2+cos2θ¯z(R¯min)zx2,γzx′=[1(R¯max)zx2-1(R¯min)zx2]*sinθ¯z*cosθ¯z.Coordinates and coal brick VR ellipses [47]. A coal sample brick was demonstrated in Figure 3(a), the measuring result were quantitatively represented in the YXZ coordinate system. The value and orientation of Rmax/Rmin on each polished surface were shown on three mutually perpendicular cube face, represented by the long/minor axis of ellipse. For example, as shown in Figure 3(b), the measured value (RXOY·max and RXOY·min) and their orientations (θx) are expressed as their coordinates and the RIS ellipsoid equation can be built.
(a)
(b)This ellipsoid equation has three eigenvalues, corresponding with the values of three RIS ellipsoid axes (Rmax>Rint>Rmin):
(3)1+ei=λi;i=1,2,3.The occurrence of three RIS ellipsoid axis is calculated according to the formula(4)Ai=arctanMiLi,Bi=arccosNi,i=1,2,3,
where Ai and Bi represent the direction and dip angle of the axis in the coordinate system, respectively. 								(5)Li=cosRix=Xiλi;Mi=cosRiy=Yiλi,Ni=cosRiz=Ziλi,i=1,2,3,Xi+Yi+Zi=λi,Xiλy′λz′-λi′λy′-λi′λz′+λi'2-γyz'2=-Yi-γxy′λz′+γxy′λi′-γzx′γyz′=-Ziγxy′γyz′+γzx′λy′-γzx′λi′.The geological occurrence of RIS ellipsoid axis (Diri∠Dipi), however, needs a data conversion as follows:
(6)ifBi≥0,Diri=Ci-Ai±360°,Dipi=Bi,ifBi<0,Diri=Ci-Ai±180°,Dipi=-Bi,
where 0<Diri<360°, 0<Dipi<90°; where Ci represents the geological direction of arrowhead in Figure 2.Therefore, the reflectance anisotropy (bireflectanceBw and Ram coefficient values) and RIS shapes (such as Rev, Rst coefficient value, Flinn’s parameters, etc.) can be determined. The details of the parameters are given as follows.(1)
Bw coefficient value [22], representing bireflectance value, is calculated according to the formula. 											(7)Bw=Rmax-Rmin.Bw/Rmax [21] represents the bireflectance ratio.(2)
Flinn’s parameter was firstly proposed for examining the geometry of three-dimensional homogeneous strain [35, 45]. Formula of Ramsay was chosen here
(8)K=ε1-ε2ε2-ε3,where(9)ε1=ln(Rmax),ε2=ln(Rint),ε3=ln(Rmin).(3)
Ram coefficient value [23, 24], representing the RIS anisotropy magnitude, is calculated according to the formula
(10)Ram=(x2-y2)1/2,where(11)x=RintRmax+Rint+Rmin+[y+13cos30°]tan30°-12,y=[Rmax(Rmax+Rint+Rmin)-13]cos30°.This value is the distance between the plotted position of RIS and the position of an isotropic RIS, a single measure of anisotropy for all styles of RIS [23].Isotropic constituent of coal is characterized by Ram=0. The higher coefficient value is, the stronger optical anisotropy will be. Typical Ram values for raw vitrinites are from 0.03 to 0.05 and the value of 0.1 represents strong anisotropy of a constituent.(4)
R
st coefficient value [23, 24], representing RIS style, is calculated according to the formula
(12)Rst=30-arctan(xy),
wherex and y are determined in the same way as shown above. Rst values range from (−30) to (+30) and the RIS shape is described as
R
st
=
(
-
30
): uniaxial negative; Rst=(+30): uniaxial positive; -30<Rst<0: biaxial negative; 0<Rst<+30: biaxial positive; Rst=0: biaxial neutral.
(5)
R
ev coefficient value [23, 24], representing the reflectance of equivalent volume of isotropic RIS, is calculated according to the formula
(13)Rev=(Rmax×Rint×Rmin)1/3.


## 4. Results

### 4.1. Anisotropic Properties of Nonoriented Samples

R
max and Rmin of 11 non-oriented samples were measured and the anisotropic properties of the material such as Bw and Bw/Rmax were calculated as shown in Table 2.Table 2
The vitrinite reflectance anisotropy of non-oriented samples, Huaibei coalfield.
Coal mine
Sample number
Coal seam
R
max
R
min
B
w
B
w/Rmax
R
max/Rmin
Linhuan
LH10
9
1.37
1.24
0.13
0.09
1.10
Shitai
ST04
3
1.44
1.29
0.15
0.10
1.12
Taoyuan
TY01
7
1.02
0.9
0.12
0.12
1.13
Taoyuan
TY02
10
1.00
0.91
0.09
0.09
1.10
Taoyuan
TY03
10
1.14
1.02
0.12
0.11
1.12
Taoyuan
TY04
7
1.06
0.95
0.11
0.10
1.12
Taoyuan
TY05
3
0.87
0.79
0.08
0.09
1.10
Qidong
QD01
9
1.09
0.97
0.12
0.11
1.12
Qidong
QD02
7
0.94
0.87
0.07
0.07
1.08
Qidong
QD03
9
0.89
0.77
0.12
0.13
1.16
Qidong
QD04
10
0.98
0.88
0.10
0.10
1.11
Abbreviations:Rmax: maximum reflectance; Rint: intermediate reflectance; Rmin: minimum reflectance; Bw: bireflectance, Rmax−Rmin.


### 4.2. Magnitude and Orientation of RIS

The magnitudes and orientations of RIS for 29 oriented samples were calculated by measuring the values and azimuths of theRmax and Rmin. Many of the coal samples suffered from tectonic deformation during coalification. The deformation types were identified by the proposed tectonic coals standard [38] and divided into three series: brittle-ductile transition, ductile shear, and brittle crush as shown in Table 3.Table 3
Magnitude and orientation of RIS, Huaibei coalfield.
Region
Coal number
Coal seam
DT
R
max axes
R
int axes
R
min axes
Value
Strike
Dip
Value
Strike
Dip
Value
Strike
Dip
Region A
TYM03
8
BD
1.01
41°
81°
0.95
177°
21°
0.77
266°
6°
TYM04
8
BC
1.00
323°
69°
0.90
153°
55°
0.81
37°
44°
TYM06
10
BC
0.87
234°
89°
0.85
28°
49°
0.74
160°
16°
TYM07
10
BC
0.79
164°
4°
0.75
291°
72°
0.72
82°
89°
XTM01
7
BC
1.11
151°
69°
0.99
239°
36°
0.89
37°
63°
XTM02
7
BD
1.11
106°
87°
1.01
157°
1°
0.97
248°
77°
XTM06
7
DS
1.13
248°
85°
1.07
107°
40°
0.98
25°
63°
XTM07
7
BC
1.11
242°
64°
1.05
126°
4°
0.98
60°
45°
Region B
LHM02
7
DS
1.31
339°
87°
1.26
250°
67°
1.14
139°
86°
LHM03
7
DS
1.55
300°
88°
1.25
108°
6°
1.17
196°
54°
LHM04
9
BD
1.46
256°
54°
1.40
344°
22°
1.29
82°
67°
LHM05
9
BD
3.09
56°
74°
2.69
310°
54°
2.36
222°
49°
LHM09
9
BD
1.44
295°
56°
1.31
174°
36°
1.15
25°
82°
LHM12
10
BD
1.24
62°
85°
1.15
196°
5°
1.10
272°
82°
LHM14
10
BD
1.50
344°
15°
1.36
205°
72°
1.18
73°
55°
LHC01*
—
—
1.37
337°
41°
1.15
246°
1°
1.13
156°
41°
HZM03
10
BC
1.92
121°
32°
1.70
47°
89°
1.63
222°
35°
HZM04
10
BC
2.01
98°
27°
1.98
333°
82°
1.82
186°
89°
HZM10
8
DS
2.66
332°
33°
2.54
233°
16°
2.51
151°
84°
HZC02*
—
—
1.46
93°
26°
1.32
345°
31°
1.17
214°
46°
HZC03*
—
—
1.45
3°
36°
1.26
101°
11°
1.13
205°
51°
BSC04*
—
—
2.61
24°
1°
2.16
289°
87°
2.08
114°
3°
BSC05*
—
—
2.54
249°
8°
2.21
343°
26°
1.96
149°
62°
Region C
YZM11
3
BC
1.70
30°
77°
1.59
160°
89°
1.34
215°
63°
YZM14
3
BC
1.69
132°
88°
1.61
248°
29°
1.38
335°
88°
YZM16
3
BC
1.82
164°
24°
1.60
3°
65°
1.30
254°
74°
STM01
3
BD
1.52
128°
84°
1.47
4°
11°
1.21
275°
75°
STM02
3
BC
1.41
215°
33°
1.36
17°
83°
1.10
97°
42°
STM03
3
BD
1.72
94°
2°
1.25
209°
64°
1.18
16°
2°
STC06*
—
—
1.52
186°
22°
1.37
277°
3°
1.19
15°
67°
LLM03
8
BD
1.06
47°
50°
0.99
261°
38°
0.88
171°
77°
LLM05
8
DS
0.84
322°
74°
0.76
247°
17°
0.59
118°
63°
LLM07
10
BC
0.95
156°
45°
0.86
345°
47°
0.74
256°
57°
LLM08
10
BC
0.83
146°
18°
0.80
357°
79°
0.75
235°
83°
LLM10
8
DS
1.11
307°
67°
0.92
35°
42°
0.83
168°
68°
*Data complied from [18].
Abbreviations: DT: deformation type; BD: brittle-ductile transition; DS: ductile shear; BC: brittle crush;Rmax: maximum reflectance; Rint: intermediate reflectance; Rmin: minimum reflectance.


### 4.3. Anisotropic Property of Oriented Samples

Based on magnitudes and orientations of RIS, Levine’s rank and Kilby’s transformations are obtained as shown in Table4.Table 4
Calculated results of RIS parameters, Huaibei coalfield.
Region
Coal number
Coal seam
DT
R
st
R
am
R
ev
B
w
B
w/Rmax
Axis ratio
K
R
max/Rint
R
int/Rmin
R
max/Rmin
Region A
TYM03
8
BD
−14.43
0.05
0.90
0.24
0.24
1.07
1.23
1.31
0.33
TYM04
8
BC
1.06
0.04
0.90
0.19
0.19
1.11
1.11
1.23
0.96
TYM06
10
BC
−20.43
0.03
0.82
0.13
0.15
1.03
1.15
1.18
0.20
TYM07
10
BC
−0.05
0.02
0.75
0.07
0.09
1.05
1.05
1.10
0.95
XTM01
7
BC
2.52
0.04
0.99
0.23
0.20
1.12
1.12
1.26
1.04
XTM02
7
BD
14.66
0.03
1.03
0.14
0.12
1.10
1.04
1.14
2.49
XTM06
7
DS
−6.05
0.03
1.06
0.16
0.14
1.06
1.10
1.16
0.64
XTM07
7
BC
−4.73
0.02
1.05
0.12
0.11
1.05
1.07
1.13
0.71
Region B
LHM02
7
DS
−12.79
0.03
1.23
0.18
0.13
1.04
1.11
1.16
0.41
LHM03
7
DS
17.42
0.06
1.31
0.38
0.24
1.23
1.07
1.32
2.93
LHM04
9
BD
−9.73
0.02
1.38
0.17
0.12
1.04
1.08
1.13
0.51
LHM05
9
BD
2.89
0.05
2.70
0.72
0.23
1.15
1.14
1.31
1.04
LHM09
9
BD
−3.89
0.04
1.29
0.30
0.21
1.10
1.15
1.26
0.70
LHM12
10
BD
12.26
0.02
1.16
0.14
0.11
1.08
1.04
1.13
2.08
LHM14
10
BD
−4.52
0.03
1.34
0.33
0.22
1.10
1.16
1.28
0.67
LHC01*
—
—
25.69
0.04
1.21
0.24
0.18
1.19
1.02
1.21
9.98
HZM03
10
BC
15.63
0.03
1.75
0.30
0.15
1.13
1.05
1.18
2.65
HZM04
10
BC
−21.15
0.02
1.94
0.19
0.09
1.02
1.09
1.10
0.19
HZM10
8
DS
18.44
0.01
2.57
0.15
0.06
1.05
1.01
1.06
3.63
HZC02*
—
—
−1.14
0.04
1.31
0.29
0.20
1.11
1.13
1.25
0.84
HZC03*
—
—
6.18
0.05
1.27
0.32
0.22
1.15
1.12
1.28
1.29
BSC04*
—
—
21.95
0.05
2.27
0.53
0.20
1.21
1.04
1.25
5.01
BSC05*
—
—
4.55
0.05
2.22
0.58
0.23
1.15
1.13
1.30
1.16
Region C
LLM03
8
BC
−7.01
0.04
0.98
0.18
0.17
1.07
1.12
1.21
0.59
LLM05
8
BC
−11.51
0.07
0.72
0.24
0.29
1.10
1.28
1.41
0.40
LLM07
10
BC
−6.37
0.05
0.85
0.21
0.22
1.10
1.17
1.28
0.60
LLM08
10
BD
−8.05
0.02
0.79
0.08
0.10
1.04
1.07
1.11
0.58
LLM10
8
BC
11.71
0.06
0.94
0.27
0.25
1.20
1.11
1.33
1.84
YZM11
3
BC
−12.87
0.05
1.54
0.37
0.21
1.07
1.19
1.27
0.38
YZM14
3
BC
−15.35
0.04
1.55
0.32
0.19
1.05
1.17
1.23
0.32
YZM16
3
BC
−4.68
0.06
1.56
0.53
0.29
1.14
1.23
1.41
0.63
STM01
3
BD
−21.59
0.05
1.39
0.32
0.21
1.03
1.22
1.26
0.17
STM02
3
BC
−20.95
0.05
1.28
0.31
0.22
1.04
1.24
1.29
0.18
STM03
3
BD
23.20
0.08
1.37
0.54
0.31
1.38
1.06
1.46
5.57
STC06*
—
—
−3.00
0.05
1.35
0.33
0.22
1.11
1.15
1.28
0.74
*Data complied from [18] except for the coal seam and deformation type of these four samples. Abbreviations: BD: brittle-ductile transition; DS: ductile shear; BC: brittle crush; DT: deformation type; Rmax: maximum reflectance; Rint: intermediate reflectance; Rmin: minimum reflectance; Bw: bireflectance, Rmax−Rmin; K: Flinn’s parameter, (ε1-ε2)/(ε2-ε3), where ε1= ln (Rmax), ε2 = ln (Rint), ε3 = ln (Rmin); Ram: RIS anisotropy magnitude, (x2-y2)1/2; Rst → RIS style, 30-arctan (x/y); Rev: reflectance of equivalent volume isotropic RIS, (Rmax×Rint×Rmin)1/3; where x = Rint/(Rmax+Rint+Rmin) + [y + 1/3 cos 30°] × tan 30° − 1/2, y = [Rmax/(Rmax+Rint+Rmin) − 1/3] × cos 30°.


## 4.1. Anisotropic Properties of Nonoriented Samples

R
max and Rmin of 11 non-oriented samples were measured and the anisotropic properties of the material such as Bw and Bw/Rmax were calculated as shown in Table 2.Table 2
The vitrinite reflectance anisotropy of non-oriented samples, Huaibei coalfield.
Coal mine
Sample number
Coal seam
R
max
R
min
B
w
B
w/Rmax
R
max/Rmin
Linhuan
LH10
9
1.37
1.24
0.13
0.09
1.10
Shitai
ST04
3
1.44
1.29
0.15
0.10
1.12
Taoyuan
TY01
7
1.02
0.9
0.12
0.12
1.13
Taoyuan
TY02
10
1.00
0.91
0.09
0.09
1.10
Taoyuan
TY03
10
1.14
1.02
0.12
0.11
1.12
Taoyuan
TY04
7
1.06
0.95
0.11
0.10
1.12
Taoyuan
TY05
3
0.87
0.79
0.08
0.09
1.10
Qidong
QD01
9
1.09
0.97
0.12
0.11
1.12
Qidong
QD02
7
0.94
0.87
0.07
0.07
1.08
Qidong
QD03
9
0.89
0.77
0.12
0.13
1.16
Qidong
QD04
10
0.98
0.88
0.10
0.10
1.11
Abbreviations:Rmax: maximum reflectance; Rint: intermediate reflectance; Rmin: minimum reflectance; Bw: bireflectance, Rmax−Rmin.


## 4.2. Magnitude and Orientation of RIS

The magnitudes and orientations of RIS for 29 oriented samples were calculated by measuring the values and azimuths of theRmax and Rmin. Many of the coal samples suffered from tectonic deformation during coalification. The deformation types were identified by the proposed tectonic coals standard [38] and divided into three series: brittle-ductile transition, ductile shear, and brittle crush as shown in Table 3.Table 3
Magnitude and orientation of RIS, Huaibei coalfield.
Region
Coal number
Coal seam
DT
R
max axes
R
int axes
R
min axes
Value
Strike
Dip
Value
Strike
Dip
Value
Strike
Dip
Region A
TYM03
8
BD
1.01
41°
81°
0.95
177°
21°
0.77
266°
6°
TYM04
8
BC
1.00
323°
69°
0.90
153°
55°
0.81
37°
44°
TYM06
10
BC
0.87
234°
89°
0.85
28°
49°
0.74
160°
16°
TYM07
10
BC
0.79
164°
4°
0.75
291°
72°
0.72
82°
89°
XTM01
7
BC
1.11
151°
69°
0.99
239°
36°
0.89
37°
63°
XTM02
7
BD
1.11
106°
87°
1.01
157°
1°
0.97
248°
77°
XTM06
7
DS
1.13
248°
85°
1.07
107°
40°
0.98
25°
63°
XTM07
7
BC
1.11
242°
64°
1.05
126°
4°
0.98
60°
45°
Region B
LHM02
7
DS
1.31
339°
87°
1.26
250°
67°
1.14
139°
86°
LHM03
7
DS
1.55
300°
88°
1.25
108°
6°
1.17
196°
54°
LHM04
9
BD
1.46
256°
54°
1.40
344°
22°
1.29
82°
67°
LHM05
9
BD
3.09
56°
74°
2.69
310°
54°
2.36
222°
49°
LHM09
9
BD
1.44
295°
56°
1.31
174°
36°
1.15
25°
82°
LHM12
10
BD
1.24
62°
85°
1.15
196°
5°
1.10
272°
82°
LHM14
10
BD
1.50
344°
15°
1.36
205°
72°
1.18
73°
55°
LHC01*
—
—
1.37
337°
41°
1.15
246°
1°
1.13
156°
41°
HZM03
10
BC
1.92
121°
32°
1.70
47°
89°
1.63
222°
35°
HZM04
10
BC
2.01
98°
27°
1.98
333°
82°
1.82
186°
89°
HZM10
8
DS
2.66
332°
33°
2.54
233°
16°
2.51
151°
84°
HZC02*
—
—
1.46
93°
26°
1.32
345°
31°
1.17
214°
46°
HZC03*
—
—
1.45
3°
36°
1.26
101°
11°
1.13
205°
51°
BSC04*
—
—
2.61
24°
1°
2.16
289°
87°
2.08
114°
3°
BSC05*
—
—
2.54
249°
8°
2.21
343°
26°
1.96
149°
62°
Region C
YZM11
3
BC
1.70
30°
77°
1.59
160°
89°
1.34
215°
63°
YZM14
3
BC
1.69
132°
88°
1.61
248°
29°
1.38
335°
88°
YZM16
3
BC
1.82
164°
24°
1.60
3°
65°
1.30
254°
74°
STM01
3
BD
1.52
128°
84°
1.47
4°
11°
1.21
275°
75°
STM02
3
BC
1.41
215°
33°
1.36
17°
83°
1.10
97°
42°
STM03
3
BD
1.72
94°
2°
1.25
209°
64°
1.18
16°
2°
STC06*
—
—
1.52
186°
22°
1.37
277°
3°
1.19
15°
67°
LLM03
8
BD
1.06
47°
50°
0.99
261°
38°
0.88
171°
77°
LLM05
8
DS
0.84
322°
74°
0.76
247°
17°
0.59
118°
63°
LLM07
10
BC
0.95
156°
45°
0.86
345°
47°
0.74
256°
57°
LLM08
10
BC
0.83
146°
18°
0.80
357°
79°
0.75
235°
83°
LLM10
8
DS
1.11
307°
67°
0.92
35°
42°
0.83
168°
68°
*Data complied from [18].
Abbreviations: DT: deformation type; BD: brittle-ductile transition; DS: ductile shear; BC: brittle crush;Rmax: maximum reflectance; Rint: intermediate reflectance; Rmin: minimum reflectance.


## 4.3. Anisotropic Property of Oriented Samples

Based on magnitudes and orientations of RIS, Levine’s rank and Kilby’s transformations are obtained as shown in Table4.Table 4
Calculated results of RIS parameters, Huaibei coalfield.
Region
Coal number
Coal seam
DT
R
st
R
am
R
ev
B
w
B
w/Rmax
Axis ratio
K
R
max/Rint
R
int/Rmin
R
max/Rmin
Region A
TYM03
8
BD
−14.43
0.05
0.90
0.24
0.24
1.07
1.23
1.31
0.33
TYM04
8
BC
1.06
0.04
0.90
0.19
0.19
1.11
1.11
1.23
0.96
TYM06
10
BC
−20.43
0.03
0.82
0.13
0.15
1.03
1.15
1.18
0.20
TYM07
10
BC
−0.05
0.02
0.75
0.07
0.09
1.05
1.05
1.10
0.95
XTM01
7
BC
2.52
0.04
0.99
0.23
0.20
1.12
1.12
1.26
1.04
XTM02
7
BD
14.66
0.03
1.03
0.14
0.12
1.10
1.04
1.14
2.49
XTM06
7
DS
−6.05
0.03
1.06
0.16
0.14
1.06
1.10
1.16
0.64
XTM07
7
BC
−4.73
0.02
1.05
0.12
0.11
1.05
1.07
1.13
0.71
Region B
LHM02
7
DS
−12.79
0.03
1.23
0.18
0.13
1.04
1.11
1.16
0.41
LHM03
7
DS
17.42
0.06
1.31
0.38
0.24
1.23
1.07
1.32
2.93
LHM04
9
BD
−9.73
0.02
1.38
0.17
0.12
1.04
1.08
1.13
0.51
LHM05
9
BD
2.89
0.05
2.70
0.72
0.23
1.15
1.14
1.31
1.04
LHM09
9
BD
−3.89
0.04
1.29
0.30
0.21
1.10
1.15
1.26
0.70
LHM12
10
BD
12.26
0.02
1.16
0.14
0.11
1.08
1.04
1.13
2.08
LHM14
10
BD
−4.52
0.03
1.34
0.33
0.22
1.10
1.16
1.28
0.67
LHC01*
—
—
25.69
0.04
1.21
0.24
0.18
1.19
1.02
1.21
9.98
HZM03
10
BC
15.63
0.03
1.75
0.30
0.15
1.13
1.05
1.18
2.65
HZM04
10
BC
−21.15
0.02
1.94
0.19
0.09
1.02
1.09
1.10
0.19
HZM10
8
DS
18.44
0.01
2.57
0.15
0.06
1.05
1.01
1.06
3.63
HZC02*
—
—
−1.14
0.04
1.31
0.29
0.20
1.11
1.13
1.25
0.84
HZC03*
—
—
6.18
0.05
1.27
0.32
0.22
1.15
1.12
1.28
1.29
BSC04*
—
—
21.95
0.05
2.27
0.53
0.20
1.21
1.04
1.25
5.01
BSC05*
—
—
4.55
0.05
2.22
0.58
0.23
1.15
1.13
1.30
1.16
Region C
LLM03
8
BC
−7.01
0.04
0.98
0.18
0.17
1.07
1.12
1.21
0.59
LLM05
8
BC
−11.51
0.07
0.72
0.24
0.29
1.10
1.28
1.41
0.40
LLM07
10
BC
−6.37
0.05
0.85
0.21
0.22
1.10
1.17
1.28
0.60
LLM08
10
BD
−8.05
0.02
0.79
0.08
0.10
1.04
1.07
1.11
0.58
LLM10
8
BC
11.71
0.06
0.94
0.27
0.25
1.20
1.11
1.33
1.84
YZM11
3
BC
−12.87
0.05
1.54
0.37
0.21
1.07
1.19
1.27
0.38
YZM14
3
BC
−15.35
0.04
1.55
0.32
0.19
1.05
1.17
1.23
0.32
YZM16
3
BC
−4.68
0.06
1.56
0.53
0.29
1.14
1.23
1.41
0.63
STM01
3
BD
−21.59
0.05
1.39
0.32
0.21
1.03
1.22
1.26
0.17
STM02
3
BC
−20.95
0.05
1.28
0.31
0.22
1.04
1.24
1.29
0.18
STM03
3
BD
23.20
0.08
1.37
0.54
0.31
1.38
1.06
1.46
5.57
STC06*
—
—
−3.00
0.05
1.35
0.33
0.22
1.11
1.15
1.28
0.74
*Data complied from [18] except for the coal seam and deformation type of these four samples. Abbreviations: BD: brittle-ductile transition; DS: ductile shear; BC: brittle crush; DT: deformation type; Rmax: maximum reflectance; Rint: intermediate reflectance; Rmin: minimum reflectance; Bw: bireflectance, Rmax−Rmin; K: Flinn’s parameter, (ε1-ε2)/(ε2-ε3), where ε1= ln (Rmax), ε2 = ln (Rint), ε3 = ln (Rmin); Ram: RIS anisotropy magnitude, (x2-y2)1/2; Rst → RIS style, 30-arctan (x/y); Rev: reflectance of equivalent volume isotropic RIS, (Rmax×Rint×Rmin)1/3; where x = Rint/(Rmax+Rint+Rmin) + [y + 1/3 cos 30°] × tan 30° − 1/2, y = [Rmax/(Rmax+Rint+Rmin) − 1/3] × cos 30°.


## 5. Discussion

### 5.1. Major Facts Affecting VRA

Both temperature and stress have influences on modifying VRA. It is therefore necessary to estimate roughly the thermal effects during coalification before analyzing RIS.In previous works, the mean random vitrinite reflectance (Rm %) shows a strong correlation (r2=0.7, n>600) with maximum burial temperature (Tmax in °C). These data are modeled by the linear regression equation [46]:
(14)ln(Rm)=0.0078Tmax-1.2.The conversion betweenRmax and Rm is following the analytical method described in GB/T 6948-2008 (equated with ISO7404/5 standard). On the basis of correlation given above, the maximum temperatures sustained during burial are estimated.As shown in Table5, the measured reflectance values range from 0.74% to 3.01, indicating conversion of medium-volatile bituminous coal to anthracite coal. Most of the studied samples experienced a maximum temperature between 120~20°C, less than 290°C The available high temperature-pressure experiments on bituminous anthracite show that there is little change in vitrinite reflectance anisotropy under the temperature less than 400°C [47, 48]. It is reasonable to believe that the vitrinite reflectance anisotropy here was mainly caused by the tectonic differential stress.Table 5
Maximum temperatures sustained during burial, Upper Permian strata, Huaibei coalfield.
Region A
Region B
Region C
Number
R
max
T
max (°C)
Number
R
max
T
max (°C)
Number
R
max
T
max (°C)
TYM03
1.01
147
LHM02
1.31
180
STM01
1.52
200
TYM04
0.94
138
LHM03
1.52
200
STM02
1.41
190
TYM06
0.88
129
LHM04
1.41
190
STM03
1.72
215
TYM07
0.81
119
LHM05
3.13
283
STC06*
1.52
200
TY01
1.02
148
LHM09
1.34
183
YZM11
1.70
214
TY02
1.00
146
LHM12
1.21
170
YZM14
1.69
213
TY03
1.14
163
LHM14
1.37
186
YZM16
1.82
223
TY04
1.06
153
LHC01*
1.37
186
ST01
1.44
193
TY05
0.87
128
LH01
1.37
186
LLM03
1.11
159
QD01
1.09
157
HZM03
1.82
223
LLM05
0.83
122
QD02
0.94
138
HZM04
1.99
234
LLM07
0.96
141
QD03
0.89
131
HZM10
2.58
261
LLM08
0.88
129
QD04
0.98
143
HZC02*
1.46
194
LLM10
1.04
151
XTM01
1.13
162
HZC03*
1.45
193
XTM02
1.12
160
BSC04*
2.61
263
XTM06
1.13
162
BSC05*
2.54
260
XTM07
1.21
170
*Data complied from [18]. During the calculation, Rmax is converted into random vitrinite reflectance at first, and then the maximum buried temperature are obtained by the formula given above. Abbreviations: Rmax: maximum reflectance; Tmax: maximum buried temperature.


### 5.2. Magnitudes of RIS

The three sampling regions represent different stress-thermal environments and, therefore, the corresponding coal samples experienced different effects of deformation and metamorphism. The parameters of reflectance anisotropy (bireflectanceBw and Ram coefficient values) and RIS shapes (such as Rev and Rst coefficient values, Flinn’s parameters) are analyzed as follows.


#### 5.2.1.  Bireflectance

The bireflectance increases progressively with upgrading of coal rank under the heat treatment without differential stress involved. In that case, the bireflectance for the coal ranking from medium-volatile bituminous to anthracite coal should be at the same level and theBw/Rmax coefficient, representing bireflectance ratio, should be suitable for the bireflectance analysis among various coal ranks.Relationships betweenRmax and Bw/Rmax are shown in Figure 4. The range of Rmax value is 0.81~1.21% in Region A (Taoyuan mine, Qidong mine, and Xutuan mine) and 1.21~3.13% in Region B (Linhuan mine and Haizi mine). This obvious difference above indicates that samples from Region B was dominated by the magmatic thermal activity. However, compared with Region A (Luling mine, Shitai mine, and Yangzhuang mine), samples in Region C where experienced extra effects of tectonic deformation and shearing heat, are manifested by the various Rmax ranges of 0.83~1.72%.Correlations betweenRmax and Bw/Rmax for the anthracites. Rmax: maximum reflectance; Bw/Rmax: bireflectance ratio.
(a)
(b)In contrast to Region B, theBw/Rmax value in Region C and Region A simultaneously increased with the Rmax and indicates that the reflectance in the two regions is dominated by the tectonic deformation and shearing heat.The 95% confidence intervals ofBw/Rmax value are 0.13±0.0195 in Region A and 0.23±0.0333 in Region C. The difference between two regions represents the effects of tectonic stress. The same indictor in Region B is 0.17±0.0326. The coal samples in the region were partly influenced by Mesozoic magmatism, and represent polarization of Ro value. Meanwhile, some strata underwent stronger deformation with increasing Bw/Rmax values.


#### 5.2.2. Reflectance of Equivalent Volume Isotropic RIS

RIS parameterRev, the reflectance of the isotropic RIS of equivalent volume, as a result of the chemical structure ordering during heating, is suggested to characterize the basic structure unites (BSU) [49].As shown in Figure5, positive linear regression relationships between Rmax and Rev(r2) for Regions A, B, and C are 0.956, 0.980, and 0.965, respectively. The correlations between Rev values and the hydrogen contents (Hdaf) have been verified that the Rev depends on the chemical composition of anthracites rather than the three-dimensional BSU [30]. Since the chemical composition and structure mainly evolved under the effect of temperature and coincided with the results, the linear- regression slope of Region B seems to be mainly influenced by the thermal effect.Correlation betweenRmax and Rev for the anthracites. Region A: y=(0.866*x)+0.057 (coef. of determination, r2=0.956); Region B: y=0.888x+0.031 (coef. of determination, r2=0.980); Region C: y=(0.838*x)+0.063 (coef. of determination, r2=0.965); Rmax: maximum reflectance; Rev: Reflectance of equivalent volume isotropic RIS, (Rmax×Rint×Rmin)1/3.
(a)
(b)


#### 5.2.3. RIS Style and Anisotropy Magnitude

The coefficient values ofRst and Ram construct polar scatter plots proposed by [23], and the three-dimensional style and anisotropy magnitude of RIS can be analyzed together. As shown in Figure 6(a), samples from Region A and C are characterized by substantial RIS style of biaxial negative (Rst ranges between −30~15). The anisotropy magnitude in Region C (Ram ranges between 0.04~0.08), however, is higher than that in Region A (Ram ranges between 0.02~0.04), which coincides with the distribution of bireflectance ratio. Conversely, as shown in Figure 6(b), the samples from Region B are characterized by significant biaxial positive character or unconspicuous biaxial negative (Rst ranges between −15~30) and the anisotropy magnitude (Ram) is various, ranging from 0.02 to 0.06.R
st-Ram polar plots of anthracites. Ram-RIS anisotropy magnitude, (x2-y2)1/2; Rst-RIS style, 30-arctan(x/y)where x=Rint/(Rmax+Rint+Rmin)+[y+1/3cos30°]×tan30°-1/2,y=[Rmax/(Rmax+Rint+Rmin)-1/3]×cos30°; RIS shape described as -30<Rst<0-biaxial negative; 0<Rst<+30-biaxial positive; Rst=(−30) uniaxial negative; Rst=(+30) uniaxial positive, Rst=0-biaxial neutral.
(a)
(b)


#### 5.2.4. RIS-Logarithmic Flinn Diagrams

For a more complete characterization of the optical properties of both raw and carbonized anthracites, the modified Flinn’sK parameters [35], related to the anisotropy as well as the optical character, are also calculated and shown in. The samples from Region C and A are characterized by evolvement of RIS style from constriction types to flattening types (0<k<1 shown in Figure 7(a)). However, RIS styles of the samples from Region B are complicated and various (as shown in Figure 7(b)), corresponding to plane strain in finite strain analysis.RIS-logarithmic Flinn diagrams. The arrowhead directs the evolutionary path of RIS style;K-Flinn’S parameter [35], (ε1-ε2)/(ε2-ε3), where ε1ln(Rmax), ε2=ln(Rint), ε3=ln(Rmin); K=∞-uniaxial prolate types; K=0-uniaxial oblate types; K=1-plane strain types;1<K<∞-constriction types; 0<K<1-flattening types.
(a)
(b)


#### 5.2.5. RIS Anisotropy Evolution Stages

The deformation path of finite strain analyses proposed by Ramsay suggested that the cleavage is developed by the route of “sphere types → uniaxial oblate types → uniaxial prolate types →uniaxial oblate types → flattening types”. The anisotropic coal samples may have the same evolution process. A related evolutionary path of RIS style during coalification was reported by Levine [15] and the stages are given as follows.The first stage was after the deposition. The peat was subjected to mild geothermal process due to broad regional subsidence and burial of overlying rocks. During this period, the ambient geologic stresses were nontectonic and only due to vertical static pressure loading by the overlying strata. Owing to the vertical downward lithostatic pressure, RIS would show the style of uniaxial negative. At the second stage, the anisotropic strain developed by tectonic differential stress (in a lateral direction) and lithostatic pressure (in the vertical direction) and the coal seams would suffer from extrusion in two directions and the anisotropic RIS represents the style of biaxial positive. At the third stage, with the increase of buried depth (more than 1000 m), the pressure on the coal seams from the overlying strata may approximately represent as isotropic hydrostatic pressure (Heim’s hypothesis). The tectonic differential stress is upgraded with the enhancement of tectonic agency. All of these factors would make the anisotropic RIS style as biaxial negative.As the suggested route given above, the style of uniaxial negative might take place either in the third stage or in the transition between the first and the second stages. More evidence is needed to constrain the actual process of physical deformation and metamorphism. During coalification, the coal samples collected from three regions likely belong to different evolution stages.(1)
Samples of Region C experienced more effects of tectonic deformation and shearing heat. With significant anisotropy magnitude (higherBw/Rmax and Ram value) and biaxial negative style, they have more likely reached into the third stage.(2)
Samples of Region A experienced more effect of burial metamorphism and less influenced by magmatism and tectonic agency. With unimpressive anisotropy magnitude (lowerBw/Rmax and Ram value) and biaxial negative style, they more likely belong to the transition between the first to the second stage.(3)
Samples of Region B influenced by the Mesozoic magmatism, representing polarized anisotropy magnitude (diversified value ofRo, Bw/Rmax and Ram). With complicated and various RIS style, they probably belong to the transition between the second and the third stage. Scanning electron microscopic (SEM) photos of coals in the study area indicated that some coal seams have strong deformation (Figure 8). The local stress field was presumably influenced by the complicated volcanic activity such as granite pegmatite and diorite porphyrite during the Yanshanian orogeny [50].SEM photographs showing different types of deformation. (a) YZM16, strongly folded deformation; (b) LHM05, ductile shear deformation; (c) HZM03, slight deformation, brittle fractures; (d) HZM10, strong deformation, locally folded.
(a)
(b)
(c)
(d)


### 5.3. Implications of RIS-Orientations

As mentioned in Section1, previous works suggest that minimum reflectance develops incrementally parallel to σ1 and maximum reflectance develops incrementally parallel to σ3 during coalification [1, 14]. These viewpoints were verified by simple shear experiments under high temperature-pressure environment [28, 33, 34]. On this basis, the vitrinite reflectance may provide important information about tectonic stress in each period during the coalification, and the σ1-direction could be indicated by the minimum reflectance axis.Eastern North China Craton experienced an important tectonic inversion during Mesozoic. The EW-trending tectonic grain was transformed to NE-NNE-trending and the contractional regime to an extensional regime during Jurassic-Middle Cretaceous [51]. The Huaibei Coalfield lies in the southeastern margin of North China Craton, and influenced by multistage and complicated tectonic events, which could be represented by the triaxial rotation of RIS.The orientations of RIS axes are respectively projected on the horizontal plane, (Figure9) and vertical section (Figure 10), and analyzed as follows:Figure 9
RIS-triaxial horizontal Projection. Coal mine label: 1-Luling mine, 2-Yangzhuang mine, 3-Shitai mine, 4-Taoyuan mine, 5-Xutuan mine, 6-Haizi mine, 7-Linhuan mine, 8-Baishan mine. Structure label: (1) Luling syncline; (2) Xisipo fault; (3) Sunan syncline; (4) Sunan anticline; (5) Nanping syncline; (6) Tongting anticline; (7) Wugou syncline.The radius of base circle represents the value of Rint; long (middle and short) arrows representRmax(Rint and Rmin).Figure 10
RIS-vertical projection. Generally speaking, stronger tectonic agency played the primary role and brought about coal seams to be RIS-orientation basically alike.


#### 5.3.1. Horizontal Projection

The horizontal projections are based on the stereographic polar method. The radius of base circle represents the value of Rint. Based on the feature that minimum reflectance develops incrementally parallel toσ1 [1, 14], the cumulative effect of multistage tectonism on each sample is indicated by the orientation of Rmin axis.The samples collected from different working faces have different RIS orientation, even though they are in the same coal mine. The samples of Region C are fan-shaped distributing on the thrust arc, and their Rmin axial directions are mainly pointing to the center of the circle such as, from north to south in order, NW-SW or NNW-SSE direction in the Luling mine), NEE-SWW direction in Yangzhuang mine, SSE-NNW or NNE-SSW direction in the Shitai mine. In other words, theRmin axial direction was parallel to the thrust- direction.The samples of Region A and B are distributing in the front of the thrust arc and a series of NNW-NE-trending folds. TheirRmin axial directions are mainly perpendicular to the corresponding fold-hinge as shown in Figure 9. These features suggest that the RIS-orientation in these areas was primarily controlled by the thrust-fold system.


#### 5.3.2. Vertical Projections

The cross-section (A-B-C) is about perpendicular to the hinges of the major folds and its position is shown in Figure9. The vertical projections on the cross-section are based on RIS-analyses of the samples from five coal mines. Since most of the Rmin axial directions are located in the cross-section, the projection provides information on RIS-dip angles.Samples of the Luling mine belong to Region C, and their RIS orientations were following the direction of thrust. Number 8 coal seam has similar RIS dips as Number 10. However, samples of the Taoyuan mine, Xutuan mine, Linhuan mine and Haizi mine seem to be more influenced by the regional tectonic stress, the distinctions between samples of different coal seam are obvious.Combining with the deformation type of coal samples (as listed in Table3), No. 7, No. 8, and No. 9 coal seams are mainly characterized as ductile-type or brittle-ductile-type, with lower Rmin axial dip angles. In contrast, N0. 10 coal seam is mainly characterized as brittle-type or brittle-ductile type, and their RIS Rmin axial dip angles are higher. The contrastive fact indicates that the coal samples in different deformation types have remarkably different RIS orientation, although they belong to the same structural unit.


## 5.1. Major Facts Affecting VRA

Both temperature and stress have influences on modifying VRA. It is therefore necessary to estimate roughly the thermal effects during coalification before analyzing RIS.In previous works, the mean random vitrinite reflectance (Rm %) shows a strong correlation (r2=0.7, n>600) with maximum burial temperature (Tmax in °C). These data are modeled by the linear regression equation [46]:
(14)ln(Rm)=0.0078Tmax-1.2.The conversion betweenRmax and Rm is following the analytical method described in GB/T 6948-2008 (equated with ISO7404/5 standard). On the basis of correlation given above, the maximum temperatures sustained during burial are estimated.As shown in Table5, the measured reflectance values range from 0.74% to 3.01, indicating conversion of medium-volatile bituminous coal to anthracite coal. Most of the studied samples experienced a maximum temperature between 120~20°C, less than 290°C The available high temperature-pressure experiments on bituminous anthracite show that there is little change in vitrinite reflectance anisotropy under the temperature less than 400°C [47, 48]. It is reasonable to believe that the vitrinite reflectance anisotropy here was mainly caused by the tectonic differential stress.Table 5
Maximum temperatures sustained during burial, Upper Permian strata, Huaibei coalfield.
Region A
Region B
Region C
Number
R
max
T
max (°C)
Number
R
max
T
max (°C)
Number
R
max
T
max (°C)
TYM03
1.01
147
LHM02
1.31
180
STM01
1.52
200
TYM04
0.94
138
LHM03
1.52
200
STM02
1.41
190
TYM06
0.88
129
LHM04
1.41
190
STM03
1.72
215
TYM07
0.81
119
LHM05
3.13
283
STC06*
1.52
200
TY01
1.02
148
LHM09
1.34
183
YZM11
1.70
214
TY02
1.00
146
LHM12
1.21
170
YZM14
1.69
213
TY03
1.14
163
LHM14
1.37
186
YZM16
1.82
223
TY04
1.06
153
LHC01*
1.37
186
ST01
1.44
193
TY05
0.87
128
LH01
1.37
186
LLM03
1.11
159
QD01
1.09
157
HZM03
1.82
223
LLM05
0.83
122
QD02
0.94
138
HZM04
1.99
234
LLM07
0.96
141
QD03
0.89
131
HZM10
2.58
261
LLM08
0.88
129
QD04
0.98
143
HZC02*
1.46
194
LLM10
1.04
151
XTM01
1.13
162
HZC03*
1.45
193
XTM02
1.12
160
BSC04*
2.61
263
XTM06
1.13
162
BSC05*
2.54
260
XTM07
1.21
170
*Data complied from [18]. During the calculation, Rmax is converted into random vitrinite reflectance at first, and then the maximum buried temperature are obtained by the formula given above. Abbreviations: Rmax: maximum reflectance; Tmax: maximum buried temperature.


## 5.2. Magnitudes of RIS

The three sampling regions represent different stress-thermal environments and, therefore, the corresponding coal samples experienced different effects of deformation and metamorphism. The parameters of reflectance anisotropy (bireflectanceBw and Ram coefficient values) and RIS shapes (such as Rev and Rst coefficient values, Flinn’s parameters) are analyzed as follows.


### 5.2.1.  Bireflectance

The bireflectance increases progressively with upgrading of coal rank under the heat treatment without differential stress involved. In that case, the bireflectance for the coal ranking from medium-volatile bituminous to anthracite coal should be at the same level and theBw/Rmax coefficient, representing bireflectance ratio, should be suitable for the bireflectance analysis among various coal ranks.Relationships betweenRmax and Bw/Rmax are shown in Figure 4. The range of Rmax value is 0.81~1.21% in Region A (Taoyuan mine, Qidong mine, and Xutuan mine) and 1.21~3.13% in Region B (Linhuan mine and Haizi mine). This obvious difference above indicates that samples from Region B was dominated by the magmatic thermal activity. However, compared with Region A (Luling mine, Shitai mine, and Yangzhuang mine), samples in Region C where experienced extra effects of tectonic deformation and shearing heat, are manifested by the various Rmax ranges of 0.83~1.72%.Correlations betweenRmax and Bw/Rmax for the anthracites. Rmax: maximum reflectance; Bw/Rmax: bireflectance ratio.
(a)
(b)In contrast to Region B, theBw/Rmax value in Region C and Region A simultaneously increased with the Rmax and indicates that the reflectance in the two regions is dominated by the tectonic deformation and shearing heat.The 95% confidence intervals ofBw/Rmax value are 0.13±0.0195 in Region A and 0.23±0.0333 in Region C. The difference between two regions represents the effects of tectonic stress. The same indictor in Region B is 0.17±0.0326. The coal samples in the region were partly influenced by Mesozoic magmatism, and represent polarization of Ro value. Meanwhile, some strata underwent stronger deformation with increasing Bw/Rmax values.


### 5.2.2. Reflectance of Equivalent Volume Isotropic RIS

RIS parameterRev, the reflectance of the isotropic RIS of equivalent volume, as a result of the chemical structure ordering during heating, is suggested to characterize the basic structure unites (BSU) [49].As shown in Figure5, positive linear regression relationships between Rmax and Rev(r2) for Regions A, B, and C are 0.956, 0.980, and 0.965, respectively. The correlations between Rev values and the hydrogen contents (Hdaf) have been verified that the Rev depends on the chemical composition of anthracites rather than the three-dimensional BSU [30]. Since the chemical composition and structure mainly evolved under the effect of temperature and coincided with the results, the linear- regression slope of Region B seems to be mainly influenced by the thermal effect.Correlation betweenRmax and Rev for the anthracites. Region A: y=(0.866*x)+0.057 (coef. of determination, r2=0.956); Region B: y=0.888x+0.031 (coef. of determination, r2=0.980); Region C: y=(0.838*x)+0.063 (coef. of determination, r2=0.965); Rmax: maximum reflectance; Rev: Reflectance of equivalent volume isotropic RIS, (Rmax×Rint×Rmin)1/3.
(a)
(b)


### 5.2.3. RIS Style and Anisotropy Magnitude

The coefficient values ofRst and Ram construct polar scatter plots proposed by [23], and the three-dimensional style and anisotropy magnitude of RIS can be analyzed together. As shown in Figure 6(a), samples from Region A and C are characterized by substantial RIS style of biaxial negative (Rst ranges between −30~15). The anisotropy magnitude in Region C (Ram ranges between 0.04~0.08), however, is higher than that in Region A (Ram ranges between 0.02~0.04), which coincides with the distribution of bireflectance ratio. Conversely, as shown in Figure 6(b), the samples from Region B are characterized by significant biaxial positive character or unconspicuous biaxial negative (Rst ranges between −15~30) and the anisotropy magnitude (Ram) is various, ranging from 0.02 to 0.06.R
st-Ram polar plots of anthracites. Ram-RIS anisotropy magnitude, (x2-y2)1/2; Rst-RIS style, 30-arctan(x/y)where x=Rint/(Rmax+Rint+Rmin)+[y+1/3cos30°]×tan30°-1/2,y=[Rmax/(Rmax+Rint+Rmin)-1/3]×cos30°; RIS shape described as -30<Rst<0-biaxial negative; 0<Rst<+30-biaxial positive; Rst=(−30) uniaxial negative; Rst=(+30) uniaxial positive, Rst=0-biaxial neutral.
(a)
(b)


### 5.2.4. RIS-Logarithmic Flinn Diagrams

For a more complete characterization of the optical properties of both raw and carbonized anthracites, the modified Flinn’sK parameters [35], related to the anisotropy as well as the optical character, are also calculated and shown in. The samples from Region C and A are characterized by evolvement of RIS style from constriction types to flattening types (0<k<1 shown in Figure 7(a)). However, RIS styles of the samples from Region B are complicated and various (as shown in Figure 7(b)), corresponding to plane strain in finite strain analysis.RIS-logarithmic Flinn diagrams. The arrowhead directs the evolutionary path of RIS style;K-Flinn’S parameter [35], (ε1-ε2)/(ε2-ε3), where ε1ln(Rmax), ε2=ln(Rint), ε3=ln(Rmin); K=∞-uniaxial prolate types; K=0-uniaxial oblate types; K=1-plane strain types;1<K<∞-constriction types; 0<K<1-flattening types.
(a)
(b)


### 5.2.5. RIS Anisotropy Evolution Stages

The deformation path of finite strain analyses proposed by Ramsay suggested that the cleavage is developed by the route of “sphere types → uniaxial oblate types → uniaxial prolate types →uniaxial oblate types → flattening types”. The anisotropic coal samples may have the same evolution process. A related evolutionary path of RIS style during coalification was reported by Levine [15] and the stages are given as follows.The first stage was after the deposition. The peat was subjected to mild geothermal process due to broad regional subsidence and burial of overlying rocks. During this period, the ambient geologic stresses were nontectonic and only due to vertical static pressure loading by the overlying strata. Owing to the vertical downward lithostatic pressure, RIS would show the style of uniaxial negative. At the second stage, the anisotropic strain developed by tectonic differential stress (in a lateral direction) and lithostatic pressure (in the vertical direction) and the coal seams would suffer from extrusion in two directions and the anisotropic RIS represents the style of biaxial positive. At the third stage, with the increase of buried depth (more than 1000 m), the pressure on the coal seams from the overlying strata may approximately represent as isotropic hydrostatic pressure (Heim’s hypothesis). The tectonic differential stress is upgraded with the enhancement of tectonic agency. All of these factors would make the anisotropic RIS style as biaxial negative.As the suggested route given above, the style of uniaxial negative might take place either in the third stage or in the transition between the first and the second stages. More evidence is needed to constrain the actual process of physical deformation and metamorphism. During coalification, the coal samples collected from three regions likely belong to different evolution stages.(1)
Samples of Region C experienced more effects of tectonic deformation and shearing heat. With significant anisotropy magnitude (higherBw/Rmax and Ram value) and biaxial negative style, they have more likely reached into the third stage.(2)
Samples of Region A experienced more effect of burial metamorphism and less influenced by magmatism and tectonic agency. With unimpressive anisotropy magnitude (lowerBw/Rmax and Ram value) and biaxial negative style, they more likely belong to the transition between the first to the second stage.(3)
Samples of Region B influenced by the Mesozoic magmatism, representing polarized anisotropy magnitude (diversified value ofRo, Bw/Rmax and Ram). With complicated and various RIS style, they probably belong to the transition between the second and the third stage. Scanning electron microscopic (SEM) photos of coals in the study area indicated that some coal seams have strong deformation (Figure 8). The local stress field was presumably influenced by the complicated volcanic activity such as granite pegmatite and diorite porphyrite during the Yanshanian orogeny [50].SEM photographs showing different types of deformation. (a) YZM16, strongly folded deformation; (b) LHM05, ductile shear deformation; (c) HZM03, slight deformation, brittle fractures; (d) HZM10, strong deformation, locally folded.
(a)
(b)
(c)
(d)


## 5.2.1.  Bireflectance

The bireflectance increases progressively with upgrading of coal rank under the heat treatment without differential stress involved. In that case, the bireflectance for the coal ranking from medium-volatile bituminous to anthracite coal should be at the same level and theBw/Rmax coefficient, representing bireflectance ratio, should be suitable for the bireflectance analysis among various coal ranks.Relationships betweenRmax and Bw/Rmax are shown in Figure 4. The range of Rmax value is 0.81~1.21% in Region A (Taoyuan mine, Qidong mine, and Xutuan mine) and 1.21~3.13% in Region B (Linhuan mine and Haizi mine). This obvious difference above indicates that samples from Region B was dominated by the magmatic thermal activity. However, compared with Region A (Luling mine, Shitai mine, and Yangzhuang mine), samples in Region C where experienced extra effects of tectonic deformation and shearing heat, are manifested by the various Rmax ranges of 0.83~1.72%.Correlations betweenRmax and Bw/Rmax for the anthracites. Rmax: maximum reflectance; Bw/Rmax: bireflectance ratio.
(a)
(b)In contrast to Region B, theBw/Rmax value in Region C and Region A simultaneously increased with the Rmax and indicates that the reflectance in the two regions is dominated by the tectonic deformation and shearing heat.The 95% confidence intervals ofBw/Rmax value are 0.13±0.0195 in Region A and 0.23±0.0333 in Region C. The difference between two regions represents the effects of tectonic stress. The same indictor in Region B is 0.17±0.0326. The coal samples in the region were partly influenced by Mesozoic magmatism, and represent polarization of Ro value. Meanwhile, some strata underwent stronger deformation with increasing Bw/Rmax values.


## 5.2.2. Reflectance of Equivalent Volume Isotropic RIS

RIS parameterRev, the reflectance of the isotropic RIS of equivalent volume, as a result of the chemical structure ordering during heating, is suggested to characterize the basic structure unites (BSU) [49].As shown in Figure5, positive linear regression relationships between Rmax and Rev(r2) for Regions A, B, and C are 0.956, 0.980, and 0.965, respectively. The correlations between Rev values and the hydrogen contents (Hdaf) have been verified that the Rev depends on the chemical composition of anthracites rather than the three-dimensional BSU [30]. Since the chemical composition and structure mainly evolved under the effect of temperature and coincided with the results, the linear- regression slope of Region B seems to be mainly influenced by the thermal effect.Correlation betweenRmax and Rev for the anthracites. Region A: y=(0.866*x)+0.057 (coef. of determination, r2=0.956); Region B: y=0.888x+0.031 (coef. of determination, r2=0.980); Region C: y=(0.838*x)+0.063 (coef. of determination, r2=0.965); Rmax: maximum reflectance; Rev: Reflectance of equivalent volume isotropic RIS, (Rmax×Rint×Rmin)1/3.
(a)
(b)


## 5.2.3. RIS Style and Anisotropy Magnitude

The coefficient values ofRst and Ram construct polar scatter plots proposed by [23], and the three-dimensional style and anisotropy magnitude of RIS can be analyzed together. As shown in Figure 6(a), samples from Region A and C are characterized by substantial RIS style of biaxial negative (Rst ranges between −30~15). The anisotropy magnitude in Region C (Ram ranges between 0.04~0.08), however, is higher than that in Region A (Ram ranges between 0.02~0.04), which coincides with the distribution of bireflectance ratio. Conversely, as shown in Figure 6(b), the samples from Region B are characterized by significant biaxial positive character or unconspicuous biaxial negative (Rst ranges between −15~30) and the anisotropy magnitude (Ram) is various, ranging from 0.02 to 0.06.R
st-Ram polar plots of anthracites. Ram-RIS anisotropy magnitude, (x2-y2)1/2; Rst-RIS style, 30-arctan(x/y)where x=Rint/(Rmax+Rint+Rmin)+[y+1/3cos30°]×tan30°-1/2,y=[Rmax/(Rmax+Rint+Rmin)-1/3]×cos30°; RIS shape described as -30<Rst<0-biaxial negative; 0<Rst<+30-biaxial positive; Rst=(−30) uniaxial negative; Rst=(+30) uniaxial positive, Rst=0-biaxial neutral.
(a)
(b)


## 5.2.4. RIS-Logarithmic Flinn Diagrams

For a more complete characterization of the optical properties of both raw and carbonized anthracites, the modified Flinn’sK parameters [35], related to the anisotropy as well as the optical character, are also calculated and shown in. The samples from Region C and A are characterized by evolvement of RIS style from constriction types to flattening types (0<k<1 shown in Figure 7(a)). However, RIS styles of the samples from Region B are complicated and various (as shown in Figure 7(b)), corresponding to plane strain in finite strain analysis.RIS-logarithmic Flinn diagrams. The arrowhead directs the evolutionary path of RIS style;K-Flinn’S parameter [35], (ε1-ε2)/(ε2-ε3), where ε1ln(Rmax), ε2=ln(Rint), ε3=ln(Rmin); K=∞-uniaxial prolate types; K=0-uniaxial oblate types; K=1-plane strain types;1<K<∞-constriction types; 0<K<1-flattening types.
(a)
(b)


## 5.2.5. RIS Anisotropy Evolution Stages

The deformation path of finite strain analyses proposed by Ramsay suggested that the cleavage is developed by the route of “sphere types → uniaxial oblate types → uniaxial prolate types →uniaxial oblate types → flattening types”. The anisotropic coal samples may have the same evolution process. A related evolutionary path of RIS style during coalification was reported by Levine [15] and the stages are given as follows.The first stage was after the deposition. The peat was subjected to mild geothermal process due to broad regional subsidence and burial of overlying rocks. During this period, the ambient geologic stresses were nontectonic and only due to vertical static pressure loading by the overlying strata. Owing to the vertical downward lithostatic pressure, RIS would show the style of uniaxial negative. At the second stage, the anisotropic strain developed by tectonic differential stress (in a lateral direction) and lithostatic pressure (in the vertical direction) and the coal seams would suffer from extrusion in two directions and the anisotropic RIS represents the style of biaxial positive. At the third stage, with the increase of buried depth (more than 1000 m), the pressure on the coal seams from the overlying strata may approximately represent as isotropic hydrostatic pressure (Heim’s hypothesis). The tectonic differential stress is upgraded with the enhancement of tectonic agency. All of these factors would make the anisotropic RIS style as biaxial negative.As the suggested route given above, the style of uniaxial negative might take place either in the third stage or in the transition between the first and the second stages. More evidence is needed to constrain the actual process of physical deformation and metamorphism. During coalification, the coal samples collected from three regions likely belong to different evolution stages.(1)
Samples of Region C experienced more effects of tectonic deformation and shearing heat. With significant anisotropy magnitude (higherBw/Rmax and Ram value) and biaxial negative style, they have more likely reached into the third stage.(2)
Samples of Region A experienced more effect of burial metamorphism and less influenced by magmatism and tectonic agency. With unimpressive anisotropy magnitude (lowerBw/Rmax and Ram value) and biaxial negative style, they more likely belong to the transition between the first to the second stage.(3)
Samples of Region B influenced by the Mesozoic magmatism, representing polarized anisotropy magnitude (diversified value ofRo, Bw/Rmax and Ram). With complicated and various RIS style, they probably belong to the transition between the second and the third stage. Scanning electron microscopic (SEM) photos of coals in the study area indicated that some coal seams have strong deformation (Figure 8). The local stress field was presumably influenced by the complicated volcanic activity such as granite pegmatite and diorite porphyrite during the Yanshanian orogeny [50].SEM photographs showing different types of deformation. (a) YZM16, strongly folded deformation; (b) LHM05, ductile shear deformation; (c) HZM03, slight deformation, brittle fractures; (d) HZM10, strong deformation, locally folded.
(a)
(b)
(c)
(d)


## 5.3. Implications of RIS-Orientations

As mentioned in Section1, previous works suggest that minimum reflectance develops incrementally parallel to σ1 and maximum reflectance develops incrementally parallel to σ3 during coalification [1, 14]. These viewpoints were verified by simple shear experiments under high temperature-pressure environment [28, 33, 34]. On this basis, the vitrinite reflectance may provide important information about tectonic stress in each period during the coalification, and the σ1-direction could be indicated by the minimum reflectance axis.Eastern North China Craton experienced an important tectonic inversion during Mesozoic. The EW-trending tectonic grain was transformed to NE-NNE-trending and the contractional regime to an extensional regime during Jurassic-Middle Cretaceous [51]. The Huaibei Coalfield lies in the southeastern margin of North China Craton, and influenced by multistage and complicated tectonic events, which could be represented by the triaxial rotation of RIS.The orientations of RIS axes are respectively projected on the horizontal plane, (Figure9) and vertical section (Figure 10), and analyzed as follows:Figure 9
RIS-triaxial horizontal Projection. Coal mine label: 1-Luling mine, 2-Yangzhuang mine, 3-Shitai mine, 4-Taoyuan mine, 5-Xutuan mine, 6-Haizi mine, 7-Linhuan mine, 8-Baishan mine. Structure label: (1) Luling syncline; (2) Xisipo fault; (3) Sunan syncline; (4) Sunan anticline; (5) Nanping syncline; (6) Tongting anticline; (7) Wugou syncline.The radius of base circle represents the value of Rint; long (middle and short) arrows representRmax(Rint and Rmin).Figure 10
RIS-vertical projection. Generally speaking, stronger tectonic agency played the primary role and brought about coal seams to be RIS-orientation basically alike.


### 5.3.1. Horizontal Projection

The horizontal projections are based on the stereographic polar method. The radius of base circle represents the value of Rint. Based on the feature that minimum reflectance develops incrementally parallel toσ1 [1, 14], the cumulative effect of multistage tectonism on each sample is indicated by the orientation of Rmin axis.The samples collected from different working faces have different RIS orientation, even though they are in the same coal mine. The samples of Region C are fan-shaped distributing on the thrust arc, and their Rmin axial directions are mainly pointing to the center of the circle such as, from north to south in order, NW-SW or NNW-SSE direction in the Luling mine), NEE-SWW direction in Yangzhuang mine, SSE-NNW or NNE-SSW direction in the Shitai mine. In other words, theRmin axial direction was parallel to the thrust- direction.The samples of Region A and B are distributing in the front of the thrust arc and a series of NNW-NE-trending folds. TheirRmin axial directions are mainly perpendicular to the corresponding fold-hinge as shown in Figure 9. These features suggest that the RIS-orientation in these areas was primarily controlled by the thrust-fold system.


### 5.3.2. Vertical Projections

The cross-section (A-B-C) is about perpendicular to the hinges of the major folds and its position is shown in Figure9. The vertical projections on the cross-section are based on RIS-analyses of the samples from five coal mines. Since most of the Rmin axial directions are located in the cross-section, the projection provides information on RIS-dip angles.Samples of the Luling mine belong to Region C, and their RIS orientations were following the direction of thrust. Number 8 coal seam has similar RIS dips as Number 10. However, samples of the Taoyuan mine, Xutuan mine, Linhuan mine and Haizi mine seem to be more influenced by the regional tectonic stress, the distinctions between samples of different coal seam are obvious.Combining with the deformation type of coal samples (as listed in Table3), No. 7, No. 8, and No. 9 coal seams are mainly characterized as ductile-type or brittle-ductile-type, with lower Rmin axial dip angles. In contrast, N0. 10 coal seam is mainly characterized as brittle-type or brittle-ductile type, and their RIS Rmin axial dip angles are higher. The contrastive fact indicates that the coal samples in different deformation types have remarkably different RIS orientation, although they belong to the same structural unit.


## 5.3.1. Horizontal Projection

The horizontal projections are based on the stereographic polar method. The radius of base circle represents the value of Rint. Based on the feature that minimum reflectance develops incrementally parallel toσ1 [1, 14], the cumulative effect of multistage tectonism on each sample is indicated by the orientation of Rmin axis.The samples collected from different working faces have different RIS orientation, even though they are in the same coal mine. The samples of Region C are fan-shaped distributing on the thrust arc, and their Rmin axial directions are mainly pointing to the center of the circle such as, from north to south in order, NW-SW or NNW-SSE direction in the Luling mine), NEE-SWW direction in Yangzhuang mine, SSE-NNW or NNE-SSW direction in the Shitai mine. In other words, theRmin axial direction was parallel to the thrust- direction.The samples of Region A and B are distributing in the front of the thrust arc and a series of NNW-NE-trending folds. TheirRmin axial directions are mainly perpendicular to the corresponding fold-hinge as shown in Figure 9. These features suggest that the RIS-orientation in these areas was primarily controlled by the thrust-fold system.


## 5.3.2. Vertical Projections

The cross-section (A-B-C) is about perpendicular to the hinges of the major folds and its position is shown in Figure9. The vertical projections on the cross-section are based on RIS-analyses of the samples from five coal mines. Since most of the Rmin axial directions are located in the cross-section, the projection provides information on RIS-dip angles.Samples of the Luling mine belong to Region C, and their RIS orientations were following the direction of thrust. Number 8 coal seam has similar RIS dips as Number 10. However, samples of the Taoyuan mine, Xutuan mine, Linhuan mine and Haizi mine seem to be more influenced by the regional tectonic stress, the distinctions between samples of different coal seam are obvious.Combining with the deformation type of coal samples (as listed in Table3), No. 7, No. 8, and No. 9 coal seams are mainly characterized as ductile-type or brittle-ductile-type, with lower Rmin axial dip angles. In contrast, N0. 10 coal seam is mainly characterized as brittle-type or brittle-ductile type, and their RIS Rmin axial dip angles are higher. The contrastive fact indicates that the coal samples in different deformation types have remarkably different RIS orientation, although they belong to the same structural unit.


## 6. Conclusions

Our results indicate that:(1)
there are obvious relationships between the anisotropic parameters and tectonic stress action. Based on the tectonic situation, the Huaibei Coalfield is divided into three regions. Samples of each region have special RIS-shape and magnitude. Samples of Region C have significant anisotropy magnitude (higherBw/Rmax and Ram value) and biaxial negative style. By contrast, samples of the Region A have unimpressive anisotropy magnitude and biaxial negative style. Samples of Region B represent polarized anisotropy magnitude (diversified value of Ro, Bw/Rmax and Ram) with complicated and various RIS styles,(2)
orientation of RIS is mainly controlled by regional tectonic stress: either parallel to the direction of thrusting or perpendicular to the hinge of folds. It is likely influenced by deformation mechanisms of coal.



---

*Source: 526016-2012-08-23.xml*