Multiple Influences of Paleogeography, Sequence Stratigraphy, and Fold Structure on Coalbed Methane Accumulation in the Lopingian (Late Permian) of Junlian Coalfield, SW China

Junlian coalfield is one of the main targets for coalbed methane (CBM) exploration and development in the southwest China. Based on field geological survey, core observation, gas content statistics, coal maceral composition, vitrinite reflectance (
 
 
 
 R
 
 
 O
 
 
 
 ), proximate analysis and trace element test, lithological types, lithofacies, sedimentary environment, and structural analysis, this research established the sequence stratigraphy frame, revealed the plane distribution characteristics of sedimentary facies, and defined CBM accumulation mode. The results show that six rock types were identified and further subdivided into twenty lithofacies types. Four types of sedimentary systems such as alluvial plains, delta, lagoon-tidal flat, and carbonate platform were summarized according to their combination characteristics. Additionally, 12-14 fourth-order sequences and three third-order sequences CSI, CSII, and CSIII were divided, and a sequence stratigraphic framework of the Lopingian coal-bearing series was established. Among them, the features of third-order sequence CSIII paleogeography from west to east are alluvial plains, deltas, lagoon-tidal flats, and limited carbonate platforms. Thick coal seams are mainly developed in the sedimentary environment of tidal flats, delta plains, and floodplains behind banks. Closely related to coal seam thickness, gas contents of Lopingian coal seams are generally higher than 8 m3/t, except the low level in northwest and partial denudation areas. CBM accumulation is significantly controlled by the fold structure, and the hydraulic plugging effect makes the syncline core favorable for CBM accumulation. Furthermore, favorable geostress conditions enable the secondary anticline to become a favorable area for CBM accumulation when the sealing conditions are better. This research will provide a theoretical guide for the exploration and development of CBM in the study area.


Introduction
During the past 20 years, CBM has emerged as an important energy resource and is expected to be a significant component in the future world energy portfolio. CBM is considered a clean fuel because its combustion releases no toxins and ash and less carbon dioxide per unit of energy than the com-bustion of coal, oil, or wood [1,2]. These advantages, as well as rapidly growing worldwide energy demands, are prompting many countries having large coal resources to evaluate the CBM potential of coal basins [3,4]. China is rich in CBM resources, and the economic production of CBM has been successfully achieved from shallow coal seams [5,6]. However, the CBM development activities in China have an extremely uneven geographical distribution concentrated in North China, especially in the southern Qinshui Basin and eastern Ordos Basin [7,8]. Southern Sichuan is an important coal-producing region in South China. With a high degree of coal thermal evolution and a wide development area, the CBM resources are abundant in the Lopingian coal-bearing rock system [9]. Junlian coalfield is an important target area for CBM exploration and development in this area [10]. Nevertheless, the CBM exploration and development have not reached commercial breakthroughs. One of the main reasons is the lack of understanding of the CBM accumulation under the complex geological conditions in the study area.
Many scholars have conducted in-depth research on CBM geology in southern Sichuan [11,12]. They have realized that the resource potential of CBM is affected by the thickness of the coalbed, buried depth, gas content, coal and rock composition, roof and floor lithology and hydrodynamic conditions, and other factors. The depositional environment of coal-bearing rock series controls the development of coal seams, which in turn affects the generation, storage, and preservation conditions of CBM, and directly affects the level of CBM content [13]. Structure is an important controlling factor for the evolution of coal reservoirs and the adjustment of the occurrence of CBM, and it is also a key factor affecting the current gas content of coal reservoirs [14]. At present, some scholars tend to agree that the CBM enrichment areas are mainly located in coalaccumulating structural units dominated by synclines [15,16]. Therefore, it is very necessary to clarify the depositional environment of coal accumulation period and current structural features of coal reservoirs.
Based on field geological survey, core observation, gas content statistics, coal maceral composition, vitrinite reflectance (R O ), proximate analysis and trace element test, lithological types, lithofacies, sedimentary environment, and structural analysis, this study is aimed at building the sequence stratigraphy frame, clarifying the sedimentary environment of coal formation period, revealing distribution characteristics of gas content in coal reservoirs, deepening CBM accumulation mode, and providing a theoretical guide for the exploration and development of CBM in the study area.

Geological Background
The Junlian coalfield, located in the south of Sichuan Province, is an important part of the southern Sichuan coalfield and one of the main targets of CBM exploration and development in the southwest China (Figure 1(a)) [10]. In the Late Permian, it was located at the southwest margin of the Upper Yangzi Craton Basin tectonically, which is part of the intact South China Plate, and the Lopingian coal-bearing series was developed on a stable sedimentary basement [17,18]. After the Lopingian deposition, Junlian area experienced superposition of multiple phases of tectonic movements, forming a pattern of alternate distribution of multiple obliques and back slopes (Figure 1(b)) [19]. The coal seams are mainly located inside these obliques or on the two flanks of the back slope, with the dip angle of the coal seams larger at the edge of the oblique and gentler at the core of the oblique. The coal seams are mostly exposed or stripped at the edge of the oblique and back slope, and the burial depth of the main coal seams is generally less than 2000 m. Coal types range from fat, coking, lean, poor to anthracite and are mainly high-rank bituminous coals and low-rank anthracite coals [10], with the maximum vitrinite reflectance between 0.67% and 3.49% and the volatile fraction yield generally less than 10%. The coal rock types dominate semidark to semibright with strong luster. Mirror coal and bright coal show diamond luster, in black, gray-black, or lead gray color and mostly with linear or fine to medium stripes and few flake structure [20].
The Longtan Formation and Changxing Formation are the main coal-bearing series in the Lopingian. The Longtan Formation was mainly formed in a sea-land transitional environment with a stratigraphic thickness of 50-130 m mainly consisting of conglomerate, sandstone, mudstone, coal, and limestone. The bottom is usually unconformably overlain on the basal Emeishan basalt or the tuff of the Maokou Formation, and the top was in integrated contact with the Changxing Formation [21,22]. There are about 1~30 coal layers developed in the Longtan Formation, and the main mineable coal seams are C 2 , C 3 , C 7 , and C 8 . The Changxing Formation dominantly comprises marine carbonate rocks and sea-land transitional deposits, including limestone, mudstone, coal seam, and fine sandstone, with stratigraphic thickness of 40-80 m and No. 0-10 coal seams [23] as well as abundant paleontological fossils including plants, brachiopods, petrous gills, foraminifera, Fusulinida, moss worms, and coral species, specifically Gigantopteris dictyophylloides Gu and Zhi, Lepidodendron acutangulum, G. nicotianaefolia, Lobatannularia heianensis, Plagiozamites oblongifolius, Phricodothyris asiatica, Oldhamina squamosa, and Aviculopecten sp., etc. [24]. The top layer was in integrated contact with the Lower Triassic Feixianguan Formation.

Methods
The enrichment of CBM is mainly dependent on the development characteristics of coal seams, which are inextricably linked to the depositional environment [25,26]. Petrographic paleogeographic analysis can reconstruct the geological landscape and environment based on the petrographic features [27], and the petrographic parameters can independently reflect the characteristics of the depositional environment in a particular geological period [17].
This research established a database for analyzing coal accumulation and CBM enrichment law. The data consists of 39 boreholes and outcrop profiles covering the entire study area. Paleogeographic parameters such as stratum thickness, sandstone thickness, mudstone thickness, limestone percentage, coal seam percentage, and the ratio of sand to mud were analyzed. The lithofacies paleogeography was reconstructed by the single-factor analysis and multiple-factor comparison method. Additionally, coal seam thickness, coal maceral content, ash yield, and sulfur content, as well as 30 gas content 2 Geofluids data of C 2 + C 3 coal seam in Changxing Formation and 39 gas content of C 7 + C 8 coal seam in Longtan Formation, were calculated. R O measurement and coal maceral composition identification were both implemented with a Leitz MPV-3 photometer microscope. The R O measurement followed the China National Standards GB/T 6948-1998. Coal maceral composition identification followed the China National Standards GB/T 8899-1998. The proximate analysis, which was used to obtain ash content, moisture, and volatile content of the coal, was performed following the China National Standards GB/T212-2008. Trace elements were analyzed by plasma mass spectrometer ICP-MS, with a test temperature of 23.6°C and a relative humidity of 41.9%. In the analysis, the sample was dissolved in a closed sample dissolver with hydrofluoric acid and nitric acid, and the hydrofluoric acid was evaporated on an electric heating plate; then the samples were sealed and dissolved with nitric acid. After dilution, the content of elements was determined by ICP-MS external standard method. Finally, the calculation and correction are carried out, and the error is less than 5%. The sequential process to study the influence of coal forming environment on methane accumulation in the Lopingian coal seams of Late Permian in Junlian area is as follows: (1) clarify the controlling effects of sedimentary environment on the thickness distribution of coal seam, coal maceral content changes, and ash and sulfur content as well       (Table 1).
Conglomerates: the conglomerates of the Lopingian in the southern Sichuan area are mainly fine and medium conglomerates, with occasional coarse conglomerates. Medium conglomerates and fine conglomerates are mainly distributed in the lower member of the Longtan Formation, with scour-fill tectonics developed at the bottom and inside and two-way interstratification visible. The thinning sequence upward corresponds to fluvial retention deposition, divergent fluvial retention deposition, or fluvial flood deposition, and the upwardly thickening sequence corresponds to estuarine dam deposition, while the two-way interstratification corresponds to tide-controlled divergent fluvial or tidal channel deposition.
Sandstone: the sandstone in south Sichuan is generally light gray to dark gray in color, and the debris composition is mainly basaltic rock debris. The sandstone filler is mainly  5 Geofluids clay mineral miscellaneous base, including chlorite, kaolinite, montmorillonite, and illite, followed by a small amount of calcareous, siliceous, and rhodochrosite cement. The debris particles are generally rounded, mostly subrounded, and subangular, with average to good sorting.
Siltstone: the siltstone is widely distributed in south Sichuan, with similar clastic composition to sandstone and deeper clastic alteration which is mostly in the form of chloritic clasts. The filler is relatively more than that in sandstone, dominated by cryptocrystalline chlorite and kaolinite fragments, and transi-tional to mudstone. According to the sedimentary structure, it can be divided into horizontal laminated siltstone interbedded with mudstone, tidally laminated siltstone interbedded with mudstone, soft sedimentary deformed siltstone interbedded with sandstone and mudstone, and rhodolitic siltstone interbedded with mudstone.
Mudstones: Lopingian coal-bearing mudstones are developed in various environments and can be classified as carbonaceous mudstones, dark mudstones, light-colored mudstones, rhodochrosite mudstones, bauxite and bauxite  6 Geofluids mudstones, root mudstones, and volcanic ash altered claystones according to the differences in composition, color, and genesis of the mudstones. The clay mineral composition is mainly chlorite, kaolinite, illite, montmorillonite and illitemontmorillonite mixed-layer minerals, containing quartz, rhodochrosite, pyrite, and other minerals, with horizontal lamination and plant or marine animal fossils. When the rhodolite mudstone and siltstone are interstratified, it can form a very characteristic "row of bones layer." Coal: the main coals of the Lopingian in southern Sichuan are mainly formed in the peat swamp environments in the delta plain divergent interdistributary bays and lagoon-tide flat environment. It is mainly anthracite, and the macroscopic coal rock type is mainly semibright to semidark coal, with strong luster, adamantine luster, and metallic luster. These coals are mostly developed striated structure or linear structure and few with flake structure. Pyrite, calcite, and quartz are mostly filled in coal fissures or distributed along the coal seam strips.
Limestone: the Lopingian limestone was mainly developed in the Changxing Formation and dominated by tuffs. It is generally dark gray-light gray in color, mainly medium-thick bedded, locally thin bedded, or lenticular, and occasionally has siliceous nodules, oolitic structures, sutures, and other sedimentary structures. The muddy tuffs are mainly composed of illite and eumonite, with occasional chlorite, ranging from 10% to 50%, mostly in thin laminations or bands, interspersed with bioclastic tuffs or dark mudstones, containing marine animals. Bioplastics are abundant and mainly are algae, followed by foraminifera, Fusulinida, mesomorphs, brachiopods, bivalves, sea lily stems, chrysocolla, etc.

Sedimentary Environment.
Based on the combination of lithofacies and their spatial distribution, the environment of the Lopingian in the study area was identified as a sea-land transitional craton basin environment consisting of river flood plains, shallow-water deltas, and tidal flat-lagoon sedimentary systems (Table 2).
(1) Facies Association A: Fluvial Plain. The fluvial plain was mainly developed in the western part of Junlian and Yanjin areas in the middle and upper member of Xuanwei Formation. It can be divided into fluvial channel, overbank, and flood plain deposit. According to the characteristics of channel sand, it can be divided into braided river and meandering river system, corresponding to the evolution of from early river life to late.   7 Geofluids ash and sulfur content in the coal bed. The lower delta plain has a relatively low degree of superposition of sand bodies in tidal-controlled distributary channels, with rich sea greenstone and more developed tidal bedding and bidirectional cross-bedding. Abundant marine animal fossils are visible in the distributary bays; There is a transition zone between the upper and lower delta plains, and the delta plains of the transition zone are affected by the rivers and tides, forming relatively stable peat swamp environment and the horizontally stable coal seams with large thickness.
The delta front includes mouth bars and distal bars. The shape of the mouth bar sand is mostly lenticular in the section and elliptical in the plane parallel to the direction of the tidal current, which is the result of the progradation of the mouth bar. The distal bar has a finer grain size than the mouth bar and is farther away from the delta plain.
(3) Facies Association C: Lagoon-Tidal Flat. The barrier unit corresponding to lagoon with brackish water in the transitional environment of Lopingian coal-bearing series is the transitional barrier sand dam (beach). In addition, the reefs in the basin edge and the tidal sand in epeiric sea shelf can also act as barriers. Due to the extremely gentle slope of the epeiric sea shelf, the direct distance between the lagoon and barrier unit may be so long; thus, the vertical sedimentary sequence of them might not be in direct contact, and the barrier sand bar can be sandwiched between carbonate platform deposits or between carbonate platform deposits and tidal flats deposits, and lagoon sediments can also be sandwiched in tidal flats.
Tidal flat deposits are mainly distributed on the outer edge of shallow-water delta facies belts. Sandstone, siltstone, mudstone, and other deposits from shallow-water deltas are transformed by tidal action to form subtidal tidal channels and intertidal tidal flats. Tidal flats mostly change into fining sedimentary sequence upward comprising sand flat-mixed flat-mud flat. The supratidal zone is generally a salt marsh environment, and peat swamps can develop.

Geochemical Analysis.
The Sr/Ba ratio value is a geochemical indicator to distinguish terrestrial and marine sedimentary environments [28]. The geochemical index V/ ðV + NiÞ shows that Lopingian coal-bearing series was formed in an oscillatory anoxic depositional environment, and the paleosalinity index Sr/Ba shows that the Lopingian was a semisaline-brackish water environment, where Mn elements could exist in seawater in a relatively stable manner as Mn 2+ , and Fe was easily oxidized to Fe 3+ . Most of the Fe element transported to the ocean that appeared in suspension. Mn element content and Mn/Fe ratio indicate that the sea level changed frequently during the Lopingian period ( Figure 3).

Sequence Stratigraphic Framework.
In this paper, the sedimentary sequence classification scheme proposed by Van Wagoner (1987) was adopted for the establishment of the sea-land transitional deposits of the inner craton basin in Junlian mining area [29,30]. The boundaries of sedimentary sequence are the terrestrial unconformity at the basin margin and the corresponding integration at the basin center [31]. Based on the sedimentary characteristics of the

Geofluids
Lopingian coal-bearing series in Junlian area, the key stratigraphic interfaces reflecting sea-level changes include regional exposed unconformity, river cut valley erosion basement, paleosols, sedimentary facies transitional surfaces, and coal seams, and four third-order stratigraphic interfaces can be identified. According to the superimposition style of the fourthorder sequence (progradation, aggradation, and retrogradation) or the shoreline migration characteristics corresponding to the fourth-order largest ocean flooding surface, the third-order composite stratigraphic sequence can be further divided into Lowstand Systems Tracts (LST), Transgressive Systems Tract (TST), and Highstand Systems Tract (HST). These three third-order composite stratigraphic sequences CSI, CSII, and CSIII are divided ascendingly based on the four third-order stratigraphic interfaces, among which CSIII includes seven fourth-order stratigraphic sequences, CSII includes five fourth-order stratigraphic sequences, while due to the restricted development and diachronism of the bottom surface of the Lopingian in some terrestrial sedimentary areas, CSI include up to four fourth-order sequences in marine sedimentary areas (Figure 4).

Lithofacies Paleogeography.
In this paper, the petrographic paleogeographic analysis is conducted only for CSIII, which is based on the sandstone/mudstone contour map with reference to the single-factor parameter contours of stratigraphic thickness, sandstone thickness, mudstone thickness, limestone thickness, etc. The sedimentary zones of CSIII are spreading nearly east-west, with the river alluvial plain distributed to the west in Gaoxian-Junlian area and the shallow-water delta in Junlian-Gongxian area, respectively. The river-controlled deltaic plain was distributed in the area of Junlian-Gongxian, and the tide-controlled deltaic plain was spread in the eastern part of the deltaic plain, and the deltaic foreland-tide flat zone was located in the area of Weixin, and the tide flat-restricted platform deposit was developed in eastern part of Xingwen-Xuyong ( Figure 5).

The Characteristics of Coal Seam Distribution, Coal
Macerals, and Coal Quality 4.4.1. Coal Seam Distribution. Coal seams are developed when the growth rate of basin accommodative space is slightly greater than the peat accumulation rate. The Lopingian sedimentary environments in Junlian area mainly comprise river alluvial plains, shallow-water deltas, and lagoontidal flat environments. Junlian area mainly develops C 7 and C 8 and C 2 and C 3 coal seams, of which the thickness of C 7 + C 8 and C 2 + C 3 coal is 3-10.5 m ( Figure 6) and 1-4.5 m (Figure 7 YSL151  YSL2-1  YSL44  YSL41  YSL50   0  150  1  3  0  150  1  3  0  150  1  3  0  150  1  3  0  150  1

Coal Maceral and Coal
Quality. The organic maceral components not only affect the pore structure and pore volume of coal reservoirs but also affect the adsorption capacity and desorption rate of coal for methane [32,33]. Ten exploration reports and drilling core test results show that the organic maceral components in the Junlian area are mainly vitrinite group, followed by inertinite group, and liptinite group content is tiny. The average content of the vitrinite group in C 7 coal seam is 71.11%, and the average content of the inertinite group is 21.51%. The total amount of organic matter accounts for 92.62% of the microscopic coal composition. The vitrinite content of the C 8 coal seam averages 70.13%; the average content of the inertinite group is 23.50%. The total amount of organic matter accounts for 93.63% of the microscopic coal composition. Compared with Junlian coals, vitrinite contents in coals of Xuyong area are low, but other characteristics of coals are basically similar ( Table 3).
The R O is an effective indicator of the coal coalification degree. According to the analysis results of coal exploration reports, the metamorphism level of the Lopingian coal seams in the Junlian area is generally as high as the lean coal-anthracite. The R O values of C 7 and C 8 coal seams are between 2.86% and 3.84%, corresponding to anthracite (Figure 8). On the whole, R O values of the main mineral coal seams in this region show a decreasing trend from west to east. Among these values, the borehole YSL6 has a minimum value of 2.86% and the borehole YSL11 has a maximum value of 3.84%.
The R O values of main coal seams C 2 and C 3 in the study area are between 2.77% and 3.74%, equal to anthracite (Figure 9). It shows a decreasing trend of main coal seams C 2 and C 3 R O values from west to east in this region on the whole, among which the borehole YSL6 has a minimum value of 2.77% and the borehole YSL11 well has a maximum value of 3.74%.

Spatial Distribution of Gas Content.
Based on the data from 39 CBM wells and borehole data, the gas-bearing characteristics of the main coal seams C 7 and C 8 and C 2 and C 3 were statistically analyzed, respectively. The R O values of both C 7 and C 8 coal seams in the southern Sichuan area are greater than 2.86% representing a high coal rank, and the gas content is greater than 8 m 3 /t in most areas, with low level only in the stripping areas and the areas south of borehole YSL4 (Figure 10). Among them, the gas content of C 7 coal seam in the Junlian area ranges from 6.96 to 18.94 m 3 /t, with an average of 14.19 m 3 /t, and the gas content of coal seam C 7 in the core area of Muai is more than 16 m 3 /t. The gas content of the C 8 coal seam is similar to that of the C 7 coal seam, and the gas content of the C 8 coal seam is between 4.9 and 20.45 m 3 /t, with an average of 13.22 m 3 /t. The gas content of the C 8 coal seam in the Muai area also exceeds 16 m 3 /t. In summary, the C 7 and C 8 coal seams can be considered as the same gas-bearing system.
The R O values of both C 2 and C 3 coal seams in the study area are greater than 2.77%, indicating a high rank, and the gas content in most areas is also greater than 8 m 3 /t, and those in the stripping area and the area northwest of the borehole YSL25 are less than this level (Figure 11). 10 Geofluids framework and the paleogeographic depositional environment. In different paleogeographic backgrounds of different stratigraphic frameworks, R P (rate of peat accumulation) and R A (increasing rate of accommodation) have different performances, and the spatial distribution of coal seam thickness in the vertical and planar directions has different characteristics [18]. The R A -R P balance in the stratigraphic frameworks is most significantly influenced by sea-level change, differential basin subsidence, and sedimentary supply in the east-west direction. The thickness of coal seams generally has a positive influence on the gas content of coal seams [34], and there is an overall positive correlation between the gas content and coal thickness in the study area. Coal seam thickness not only controls the gas potential of  coal reservoirs but also has an impact on CBM preservation. In the case of unsatisfactory storage conditions for the roof and floor of the coal seam, the thickness of the coal seam has a more significant influence on the gas content. Under the condition of offshore coal accumulation on the margin of the craton epeiric sea, most of river alluvial plains and river-controlled upper delta plains near the source area were in a compensation or overcompensation state, due to the relatively sufficient land source supply, that is, R P > R A . When the base level rise rate and the basin structure subsidence rate were fast, the accommodation space increased rapidly, and R P and R A could maintain a stable equilibrium state for a long time. Therefore, thick coal seams generally developed near the MFS (maximum flooding surface) or in the HST of CSIII and CSII. The coal accumulation center in the delta system with relatively low accommodation space was generally located in the lower delta plain or the interdistributary bays of the upper and lower delta plain transition zone.

Discussion
According to the analysis results of the sequence stratigraphic framework in the east-west direction, it shows that, from CSI and CSII upward to CSIII, coal accumulation grad-ually increased and also the thickness of coal seams. During the CSI stage, coal seams were mainly developed in TST and HST, mainly distributed in the Xuyong area on the north and south sides of the north Guizhou-south Sichuan uplift belt, where the subsidence rate was relatively high. The thickness of the coal seams ranged from 0.2 to 5.69 m and gradually pinched westward. The coal-accumulating environment was mainly lagoon-tidal flat environment, with low ash content and high sulfur content.
In the CSII stage, the coal accumulation effect was stronger than that in the CSI stage and developed in both TST and HST. The thickness of the coal seam ranges from 0.1 to 7.67 m and gradually decreases from east to west also with decreased number of coal seams. The coal-accumulating environment was tidal-controlled lower delta plain and the interdistributary bay of the upper and lower delta plain transition zone, with high ash content and medium sulfur content.
During the CSIII stage, coal accumulation was the strongest, and coal accumulation was generally developed in the whole area. The thickness of the coal seam ranges from 1 to 13 m and increases first and then decreases from east to   13 Geofluids seams, with characteristics of large coal seam thickness and good lateral continuity. The maximum thickness of the coal seam can reach 6 m, with an average of about 2.5 m. In stage of HST, it mainly develops C 2 , C 3 , and C 7 coal seams, among which the C 7 coal seam has the largest thickness and also the best horizontal continuity as well as the maximum coal thickness reaching 13 m with the average about 5.6 m. These coals were formed in the transition zone of the upper and lower delta plains, with medium ash content and medium sulfur content. The coal seams in the western part of Junlian area were mainly developed in the floodplains of the river alluvial plains of HST and the peat swamps in the interdistributary bays of the upper delta plains, with high ash content and low sulfur content.

5.2.
Influence of Fold Structure on CBM Accumulation. Fold structures have a significant control effect on the migration and accumulation of CBM [35,36]. On the one hand, coalbed in the fold structure rises and the static pressure decreases accordingly, which is conducive to the desorption of CBM. On the other hand, fissures developed in the anticline and syncline axis with stronger ground stresses, coal reservoir porosity and permeability will increase, and reservoir pressure decreases, which is also conducive to the desorption of CBM. The barrier-type fold structures with open syncline and closed anticline were mostly developed in southern Sichuan. The two flanks and the axis above the neutral plane of the syncline structure are compressive stress fields, and below the neutral surface are tensile stress fields, and because of the deep buried depth of coal seams, only a few open fissures are produced and part of the stress is released to form a relatively low-pressure area. Therefore, the two flanks and the axis above neutral surface of the syncline are conducive to the storage and accumulation of CBM. In particular, the axis of the syncline is often an area with high gas content. The secondary anticlines without weathering and denudation are also conducive to the enrichment and accumulation of CBM. The secondary anticlines in the Junlian area are located on the two flanks of a large, wide, and gentle compound syncline. The anticline is small in amplitude, the flanks are gentle, and the fissures are not developed. Gas formed in coal seams migrates upwards, forming free gas reservoirs under favorable sealing conditions. Above the neutral surface of the anticline axis, there is a tensile stress field, which is a low-pressure zone, with developed fissures and increased coal seam permeability. The large accommodation space is conducive to the accumulation of CBM.

Enrichment Model of CBM in the Junlian
Coalfield. Sedimentary facies belts have a significant control effect on the development of the gas-bearing system combination. From the river alluvial plain to the tidal flat-restricted platform environment, the CBM gas-bearing system combination is becoming more and more complicated [37]. The Late Permian coal-bearing strata in the study area are dominated by delta and tidal flat deposits. The distribution of the depositional system has obvious east-west zoning characteristics, which are river alluvial plain, delta, tidal flat-lagoon, and tidal flat-carbonate platform from west to east.
The coal-bearing strata receive atmospheric precipitation in the outcropping zone of the syncline flank to form groundwater, which flows from the flank to the axis. The hydraulic effect is gradually weakened, and eventually, a 14 Geofluids groundwater retention zone is formed in the syncline core. Groundwater runoff produced a hydraulic plugging effect on the escape of CBM, and the syncline core area is conducive to the accumulation of CBM [38]. The Muai core exploration and development area in Junlian is mainly composed of a wide and gentle syncline. The groundwater flow channels provided by the two flanks of the syncline are strong runoff areas, and the coalbed water gradually converges downwards to the retention area at the syncline axis. During the flow of groundwater, the direction of runoff is opposite to the upward escape of CBM, forming a hydrodynamic blockage. The groundwater in the retention zone gradually increases, and the pressure of the coal reservoir increases correspondingly to form a confined groundwater zone. Affected by the lithology of the clastic rock and the small stratum dip angle, the groundwater flow in the coal-bearing strata is very limited and weak, forming a medium-weak runoff zone. CBM is blocked by water and caprocks, and it is easy to form a confined water-blocked CBM reservoir. Large-scale anticline structures are not conducive to the accumulation of CBM, and the axis of secondary anticlines can work as free gas accumulation areas. Some secondary anticlines developed in the Junlian mining area; within them, the coal seams are buried at a depth of about 300 m with mostly mudstone roof, which exhibits high plasticity. The tensile stress will only increase the plastic deformation and not produce open cracks; thus, the roof of coal still maintains good sealing performance. Under the action of compression stress in the anticline flanks, the coal reservoir pressure increases. Affected by the tectonic stress, the CBM in the anticline flanks migrates to the upper part, forming a CBM enrichment area in the secondary anticline axis. At this time, CBM is still enriched above the neutral surface, and free gas accounts for a considerable proportion. The core of the secondary anticline has a high degree of groundwater salinity (>10000 mg/L) with slow groundwater flow, forming a hydraulic seal on the CBM, which is conducive to the accumulation of CBM.

Conclusion
(1) Six rock types and twenty lithofacies types were identified of the Lopingian coal-bearing series in southern Sichuan. Four types of sedimentary systems were summarized, and five key sequence interfaces were identified. Combined with the spatial distribution of regional marine markers, fourteen to sixteen fourth-order sequences and three third-order sequences were divided. The stratigraphic framework of the Lopingian coal-bearing series was established (2) The third-order sequence CSIII paleogeography features from west to east are alluvial plains, deltas, lagoon-tidal flats, and limited carbonate platforms. Thick coal seams are mainly developed in the sedimentary environment of tidal flats, delta plains, and floodplains behind banks. The vitrinite content in the continental coal-accumulating environment is significantly higher than that in the transitional coal-accumulating environment, while the inertinite content is the opposite (3) The gas content of Lopingian coal seams in the Junlian area, which is generally higher than 8 m 3 /t, is closely related to coal seam thickness. Relative sealevel changes and sedimentary paleogeography affect the spatial distribution of coal-accumulating sedimentary environment and the degree of coal seam development, which in turn control gas generation potential and preservation conditions of CBM (4) CBM accumulation is significantly controlled by the fold structure. The hydraulic plugging effect makes the syncline core favorable for CBM accumulation.
In addition, favorable geostress conditions enable the secondary anticline to become a favorable area for CBM accumulation when the sealing conditions are better

Data Availability
This research established a database consisting of 39 boreholes and outcrop profiles for coal accumulation and CBM enrichment rule analysis of the Lopingian coal-bearing series.

Conflicts of Interest
The authors declare that they have no conflicts of interest.