Journal of Geological ResearchVolume 2012 (2012), Article ID 526016, 16 pagesdoi:10.1155/2012/526016

Research Article

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

Yudong Wu,1,2 Quanlin Hou,² Yiwen Ju,² Daiyong Cao,³ Junjia Fan,⁴ and Wei Wei²

¹MLR Key Laboratory of Metallogeny and Mineral Assessment, Institute of Mineral Resources, CAGS, Beijing 100037, China²Key Lab of Computational Geodynamics of Chinese Academy of Sciences and College of Earth Science, Graduate University of the Chinese Academy of Sciences, Beijing 100049, China³Key Laboratory of Coal Resources, China University of Mining and Technology, Beijing 100083, China⁴Key laboratory of Basin Structure and Petroleum Accumulation, CNPC and PetroChina Research Institute of Petroleum Exploration and Development, Beijing 100083, China

Received 28 February 2012; Accepted 6 April 2012

Academic Editor: Hongyuan Zhang

Copyright © 2012 Yudong Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 (higher and 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.

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 (), anisotropy magnitude (), and reflectance of equivalent volume isotropic RIS (). 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 during coalification, and maximum reflectance develops incrementally parallel to the direction of minimum compressive stress [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 Figure 1(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 Figure 1(b). The major coal seams occur in the lower Permian system () and minor coal seams locally occur in the upper Permian System (). 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.

(a)

(b)

Figure 1: 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 Figure 1(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).

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 Table 1. 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 Figure 1 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 (Figure 2). 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 (, , and 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 with , , were mutually perpendicular with the surfaces and horizontal and and 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 ( and ) on each polished surface were measured in polarized light () and by rotating the microscope stage through 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 and .

As for nonoriented samples, at least 20 points were measured with the same analytical method, and the mean statistical values of and were obtained.

3.3. Data Processing

On the basis of measured mean values, the values and occurrences of maximum reflectance (), intermediate reflectance (), and minimum reflectance () 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: where

(a)

(b)

Figure 3: 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 coordinate system. The value and orientation of 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 ( and ) and their orientations () are expressed as their coordinates and the RIS ellipsoid equation can be built.

This ellipsoid equation has three eigenvalues, corresponding with the values of three RIS ellipsoid axes ():

The occurrence of three RIS ellipsoid axis is calculated according to the formula where and represent the direction and dip angle of the axis in the coordinate system, respectively.

The geological occurrence of RIS ellipsoid axis (), however, needs a data conversion as follows: where , ; where represents the geological direction of arrowhead in Figure 2.

Therefore, the reflectance anisotropy (bireflectance and coefficient values) and RIS shapes (such as , coefficient value, Flinn’s parameters, etc.) can be determined. The details of the parameters are given as follows.(1) coefficient value [22], representing bireflectance value, is calculated according to the formula.    [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  where (3) coefficient value [23, 24], representing the RIS anisotropy magnitude, is calculated according to the formula  where  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 . The higher coefficient value is, the stronger optical anisotropy will be. Typical values for raw vitrinites are from 0.03 to 0.05 and the value of 0.1 represents strong anisotropy of a constituent.

(4) coefficient value [23, 24], representing RIS style, is calculated according to the formula  where and are determined in the same way as shown above. values range from (−30) to () and the RIS shape is described as : uniaxial negative; : uniaxial positive; : biaxial negative; : biaxial positive; : biaxial neutral.(5) coefficient value [23, 24], representing the reflectance of equivalent volume of isotropic RIS, is calculated according to the formula

4. Results

4.1. Anisotropic Properties of Nonoriented Samples

and of 11 non-oriented samples were measured and the anisotropic properties of the material such as and / 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

/
/

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: : maximum reflectance; : intermediate reflectance; : minimum reflectance; : bireflectance, −.

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 the and . 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
axes
axes
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; : maximum reflectance; : intermediate reflectance; : 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 Table 4.

Table 4: Calculated results of RIS parameters, Huaibei coalfield.

Region
Coal number
Coal seam
DT
R st
R am
R ev

/

Axis ratio

/
/
/

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; : maximum reflectance; : intermediate reflectance; : minimum reflectance; : bireflectance, −; K: Flinn’s parameter, (-)/(-), where = ln (), = ln (), = ln (); : RIS anisotropy magnitude, ; → style, 30-arctan (/); : reflectance of equivalent volume isotropic RIS, ; where = /() + [ + 1/3 cos 30°] × tan 30° − 1/2, = [/() − 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 ( %) shows a strong correlation (, ) with maximum burial temperature ( in °C). These data are modeled by the linear regression equation [46]:

The conversion between and 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 Table 5, 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

(°C)
Number

(°C)
Number

(°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, is converted into random vitrinite reflectance at first, and then the maximum buried temperature are obtained by the formula given above. Abbreviations: : maximum reflectance; : 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 (bireflectance and coefficient values) and RIS shapes (such as and 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 the coefficient, representing bireflectance ratio, should be suitable for the bireflectance analysis among various coal ranks.

Relationships between and / are shown in Figure 4. The range of 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 ranges of 0.83~1.72%.

(a)

(b)

Figure 4: Correlations between and for the anthracites. : maximum reflectance; : bireflectance ratio.

In contrast to Region B, the value in Region C and Region A simultaneously increased with the and indicates that the reflectance in the two regions is dominated by the tectonic deformation and shearing heat.

The 95% confidence intervals of value are in Region A and in Region C. The difference between two regions represents the effects of tectonic stress. The same indictor in Region B is . The coal samples in the region were partly influenced by Mesozoic magmatism, and represent polarization of value. Meanwhile, some strata underwent stronger deformation with increasing values.

5.2.2. Reflectance of Equivalent Volume Isotropic RIS

RIS parameter , 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 Figure 5, positive linear regression relationships between and for Regions A, B, and C are 0.956, 0.980, and 0.965, respectively. The correlations between values and the hydrogen contents () have been verified that the 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.

(a)

(b)

Figure 5: Correlation between and for the anthracites. Region A: (coef. of determination, ); Region B: (coef. of determination, ); Region C: (coef. of determination, ); : maximum reflectance; : Reflectance of equivalent volume isotropic RIS, .

5.2.3. RIS Style and Anisotropy Magnitude

The coefficient values of and 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 ( ranges between −30~15). The anisotropy magnitude in Region C ( ranges between 0.04~0.08), however, is higher than that in Region A ( 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 ( ranges between −15~30) and the anisotropy magnitude () is various, ranging from 0.02 to 0.06.

(a)

(b)

Figure 6: - polar plots of anthracites. -RIS anisotropy magnitude, ; -RIS style, 30-where ,; RIS shape described as -biaxial negative; -biaxial positive; (−30) uniaxial negative; uniaxial positive, 0-biaxial neutral.

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’s 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 ( 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.

(a)

(b)

Figure 7: RIS-logarithmic Flinn diagrams. The arrowhead directs the evolutionary path of RIS style; -Flinn’S parameter [35], ()/(), where , , ; -uniaxial prolate types; 0-uniaxial oblate types; -plane strain types;-constriction types; -flattening types.

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 (higher and 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 (lower and 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 of , and ). 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].

(a)

(b)

(c)

(d)

Figure 8: 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.

5.3. Implications of RIS-Orientations

As mentioned in Section 1, previous works suggest that minimum reflectance develops incrementally parallel to and maximum reflectance develops incrementally parallel to 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 -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, (Figure 9) 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 represent ( and ).

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, 14], the cumulative effect of multistage tectonism on each sample is indicated by the orientation of 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, the 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. Their 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 Figure 9. The vertical projections on the cross-section are based on RIS-analyses of the samples from five coal mines. Since most of the 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 Table 3), No. 7, No. 8, and No. 9 coal seams are mainly characterized as ductile-type or brittle-ductile-type, with lower axial dip angles. In contrast, N0. 10 coal seam is mainly characterized as brittle-type or brittle-ductile type, and their RIS 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 (higher and 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 , and ) 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.

Acknowledgments

These research results are part of a key project carried out in 2006–2011 and financially supported by The national Natural Science Foundation of China (Grant no. 41030422; 40940014; 40772135) and National Basic Research Program of China (Grant no. 2009CB219601). The authors are grateful to Dai Jiming and Ai Tianjie at China University of Mining and Technology (Beijing) for their assistance in the experiments. Many thanks to Zheng Yadong, Wang Guiliang, Jin Weijun and Liu Qinfu, Liu Dameng, and three reviewers for their help with improving the paper.

References

I. Stone and A. Cook, “The influence of some tectonic structures upon vitrinite reflectance,” The Journal of Geology, vol. 87, pp. 497–508, 1979.

E. Stach, M.-T. Mackowsky, M. Techmuller, G. H. Taylor, D. Chandra, and R. Techmuller, Stach's Textbook of Coal Petrology, 3rd edition, 1982.

J. V. Ross, R. M. Bustin, and J. N. Rouzaud, “Graphitization of high rank coals—the role of shear strain: experimental considerations,” Organic Geochemistry, vol. 17, no. 5, pp. 585–596, 1991.

R. M. Bustin, J. V. Ross, and J. N. Rouzaud, “Mechanisms of graphite formation from kerogen: experimental evidence,” International Journal of Coal Geology, vol. 28, no. 1, pp. 1–36, 1995.

J. Komorek and R. Morga, “Relationship between the maximum and the random reflectance of vitrinite for coal from the Upper Silesian Coal Basin (Poland),” Fuel, vol. 81, no. 7, pp. 969–971, 2002.

S. Pusz, S. Duber, and B. K. Kwiecińska, “The study of textural and structural transformations of carbonized anthracites,” Fuel Processing Technology, vol. 77-78, pp. 173–180, 2002.

J. Komorek and R. Morga, “Vitrinite reflectance property change during heating under inert conditions,” International Journal of Coal Geology, vol. 54, no. 1-2, pp. 125–136, 2003.

S. Pusz, B. K. Kwiecińska, and S. Duber, “Textural transformation of thermally treated anthracites,” International Journal of Coal Geology, vol. 54, no. 1-2, pp. 115–123, 2003.

M. Krzesińska, S. Pusz, and A. Koszorek, “Elastic and optical anisotropy of the single-coal monolithic high-temperature (HT) carbonization products obtained on a laboratory scale,” Energy and Fuels, vol. 19, no. 5, pp. 1962–1970, 2005.

D. Y. Cao, X. M. Li, and S. R. Zhang, “Influence of tectonic stress on coalification: stress degradation mechanism and stress polycondensation mechanism,” Science in China, Series D, vol. 50, no. 1, pp. 43–54, 2007.

J. P. Jiang, G. Y. Gao, and J. W. Kang, “Tests on vitrinite reflectance of coal and analysis of tectonic stress field,” Chinese Journal of Geophysics (Acta Geophysica Sinica), vol. 50, no. 1, pp. 138–145, 2007.

J. Komorek and R. Morga, “Evolution of optical properties of vitrinite, sporinite and semifusinite in response to heating under inert conditions,” International Journal of Coal Geology, vol. 71, no. 4, pp. 389–404, 2007.

I. Suárez-Ruiz and A. B. García, “Optical parameters as a tool to study the microstructural evolution of carbonized anthracites during high-temperature treatment,” Energy and Fuels, vol. 21, no. 5, pp. 2935–2941, 2007.

J. R. Levine and A. Davis, “Optical anisotropy of coals as an indicator of tectonic deformation, Broad Top Coal Field, Pennsylvania ( USA),” Geological Society of America Bulletin, vol. 95, no. 1, pp. 100–108, 1984.

J. R. Levine and A. Davis, “Reflectance anisotropy of Upper Carboniferous coals in the Appalachian foreland basin, Pennsylvania, U.S.A,” International Journal of Coal Geology, vol. 13, no. 1-4, pp. 341–373, 1989.

B. Jiang and Y. Qin, “Experimental research on deformation of optical fabric of coal vitrinite reflectance experimental research on deformation of optical fabric of coal vitrinite reflectance,” Coal Geology & Exploration, vol. 25, pp. 11–15, 1997 (Chinese).

J. C. Hower and A. Davis, “Vitrinite reflectance anisotropy as a tectonic fabric element,” Geology, vol. 9, pp. 165–168, 1981.

D. Y. Cao, “The vitrinite reflectance anisotropy in the nappe structure in the Huaibei coalfield, Anhni province,” Geology Review, vol. 36, pp. 333–340, 1990 (Chinese).

W. Langenberg and W. Kalkreuth, “Reflectance anisotropy and syn-deformational coalification of the Jewel seam in the Cadomin area, Alberta, Canada,” International Journal of Coal Geology, vol. 19, no. 1-4, pp. 303–317, 1991.

B. Jiang, F. Y. Xu, Y. Liu, and F. L. Jin, “Vitrinite optical fabric and stress-strain analysis of northern margin of Chaidamu Basin,” Journal of China University of Mining & Technology, vol. 31, pp. 561–564, 2002.

J. Jones and S. Creaney, “Optical character of thermally metamorphosed coals of northern England,” Journal of Microscopy, vol. 109, pp. 105–118, 1977.

G. K. Khorasani, D. G. Murchison, and A. C. Raymond, “Molecular disordering in natural cokes approaching dyke and sill contacts,” Fuel, vol. 69, no. 8, pp. 1037–1046, 1990.

W. E. Kilby, “Recognition of vitrinite with non-uniaxial negative reflectance characteristics,” International Journal of Coal Geology, vol. 9, no. 3, pp. 267–285, 1988.

W. E. Kilby, “Vitrinite reflectance measurement—some technique enhancements and relationships,” International Journal of Coal Geology, vol. 19, no. 1-4, pp. 201–218, 1991.

J. C. Hower, R. F. Rathbone, G. D. Wild, and A. Davis, “Observations on the use of vitrinite maximum reflectance versus vitrinite random reflectance for high volatile bituminous coals,” Journal of Coal Quality, vol. 13, pp. 71–76, 1994.

K. R. Wilks, M. Mastalerz, R. M. Bustin, and J. V. Ross, “The role of shear strain in the graphitization of a high-volatile bituminous and an anthracitic coal,” International Journal of Coal Geology, vol. 22, no. 3-4, pp. 247–277, 1993.

D. W. Houseknecht and C. M. B. Weesner, “Rotational reflectance of dispersed vitrinite from the Arkoma basin,” Organic Geochemistry, vol. 26, no. 3-4, pp. 191–206, 1997.

J. V. Ross and R. M. Bustin, “Vitrinite anisotropy resulting from simple shear experiments at high temperature and high confining pressure,” International Journal of Coal Geology, vol. 33, no. 2, pp. 153–168, 1997.

S. Duber and J. N. Rouzaud, “Calculation of relectance values for two models of texture of carbon materials,” International Journal of Coal Geology, vol. 38, no. 3-4, pp. 333–348, 1999.

S. Duber, S. Pusz, B. K. Kwiecińska, and J. N. Rouzaud, “On the optically biaxial character and heterogeneity of anthracites,” International Journal of Coal Geology, vol. 44, no. 3-4, pp. 227–250, 2000.

A. C. Cook, D. G. Murchison, and E. Scott, “Optically biaxial anthracitic vitrinites,” Fuel, vol. 51, no. 3, pp. 180–184, 1972.

F. Ting, “Uniaxial and biaxial vitrinite reflectance models and their relationship to paleotectonics,” in Organic Maturation Studies and Fossil Fuel Exploration, J. Brooks, Ed., pp. 379–392, Academic Press, London, UK, 1981.

R. M. Bustin, J. V. Ross, and I. Moffat, “Vitrinite anisotropy under differential stress and high confining pressure and temperature: preliminary observations,” International Journal of Coal Geology, vol. 6, no. 4, pp. 343–351, 1986.

B. Jiang, F. Jin, Q. Zhou, and W. Wang, “Experimental research on deformation of optical fabric of coal vitrinite reflectance,” Coal Geology & Exploration, vol. 25, pp. 11–15, 1997 (Chinese).

J. G. Ramsay, Folding and Fracturing of Rocks, Mc Graw-Hill, New York, NY, USA, 1967.

Y. D. Wu, Y. W. Ju, Q. L. Hou, S. Hu, S. Q. Ni, and J. J. Fan, “Characteristics of tectono-thermal modeling and restriction on coalbed-gas generation in Sulin mining area, Huaibei coalfield,” Progress in Natural Science, vol. 19, pp. 1134–1141, 2009 (Chinese).

W. Guiliang, J. Bo, C. Daiyong, Z. Hai, and J. Weijun, “On the Xuzhou-Suzhou arcuate duplex-imbricate fan thrust system,” Acta Geologica Sinica, vol. 72, no. 3, pp. 235–236, 1998.

Y. W. Ju, Q. L. Hou, B. Jiang, G. L. Wang, and A. M. Fang, “Tectonic coals: structure and physical properties of reservoirs,” in Proceedings of the Proceedings of the 6th International Workshop on CBM/CMM in China, Beijing, China, 2006.

G. Zhu, C. Song, D. Wang, G. Liu, and J. Xu, “Studies on ⁴⁰Ar/ ³⁹Ar thermochronology of strike-slip time of the Tan-Lu fault zone and their tectonic implications,” Science in China, Series D, vol. 44, no. 11, pp. 1002–1009, 2001.

G. Zhu, Y. Wang, G. Liu, M. Niu, C. Xie, and C. Li, “⁴⁰Ar/ ³⁹Ar dating of strike-slip motion on the Tan-Lu fault zone, East China,” Journal of Structural Geology, vol. 27, no. 8, pp. 1379–1398, 2005.

G. Zhu, C. L. Xie, Y. S. Wang, M. L. Niu, and G. S. Liu, “Characteristics of the Tan-Lu high-pressure strike-slip ductile shear zone and its ⁴⁰Ar/³⁹Ar dating,” Acta Petrologica Sinica, vol. 21, no. 6, pp. 1687–1702, 2005.

A. Yin and S. Nie, “An indentation model for the north and south China collision and development of the Tan-Lu and Honam fault systems, Eastern Asia,” Tectonics, vol. 12, pp. 810–813, 1993.

G. Zhu, Y. S. Wang, N. M. Lan, G. S. Liu, and C. L. Xie, “Synorogenic movement of the Tan-Lu fault zone,” Earth Science Frontiers, vol. 11, pp. 169–182, 2004 (Chinese).

D. Liu, Y. Yao, D. Tang, S. Tang, Y. Che, and W. Huang, “Coal reservoir characteristics and coalbed methane resource assessment in Huainan and Huaibei coalfields, Southern North China,” International Journal of Coal Geology, vol. 79, no. 3, pp. 97–112, 2009.

D. Flinn, “On folding during three-dimensional progressive deformation,” Quarterly Journal of the Geological Society, vol. 118, pp. 385–428, 1962.

C. Barker and M. Pawlewicz, “The correlation of vitrinite reflectance with maximum temperature in humic organic matter,” Lecture Notes in Earth Science, vol. 5, pp. 79–81, 1986.

F. Goodarzi and D. G. Murchison, “Optical properties of carbonized vitrinites,” Fuel, vol. 51, no. 4, pp. 322–328, 1972.

D. G. Murchison, “Petrographic aspects of coal structure: reactivity of macerals in laboratory and natural environments,” Fuel, vol. 70, no. 3, pp. 296–315, 1991.

S. Duber, J. N. Rouzaud, C. Clinard, and S. Pusz, “Microporosity and optical properties of some activated chars,” Fuel Processing Technology, vol. 77-78, pp. 221–227, 2002.

G. L. Wang, B. Jiang, D. Y. Cao, H. Zou, and W. J. Jin, “On the Xuzhou-Suzhou arcuate duplex-imbricate fan thrust system.,” Acta Geologica Sinica (in Chinese with English summary), vol. 72, pp. 228–236, 1998.

M. Zhai, R. Zhu, J. Liu et al., “Time range of Mesozoic tectonic regime inversion in eastern North China Block,” Science in China, Series D, vol. 47, no. 2, pp. 151–159, 2004.