Identification of Favorable Zones of Gas Accumulation via Fault Distribution and Sedimentary Facies: Insights From Hangjinqi Area, Northern Ordos Basin

The Hangjinqi area was explored for natural gas around 40 years ago, but the efficient consideration in this area was started around a decade ago for pure gas exploration. Many wells have been drilled, yet the Hangjinqi area remains an exploration area, and the potential zones are still unclear. The Lower Shihezi Formation is a proven reservoir in the northern Ordos Basin. This study focuses on the second and third members of the Lower Shihezi Formation to understand the controlling factors of faults and sedimentary facies distribution, aimed to identify the favorable zones of gas accumulation within the Hangjinqi area. The research is conducted on a regional level by incorporating the 3D seismic grid of about 2500 km2, 62 well logs, and several cores using seismic stratigraphy, geological modeling, seismic attribute analysis, and well logging for the delineation of gas accumulation zones. The integrated results of structural maps, thickness maps, sand-ratio maps, and root mean square map showed that the northwestern region was uplifted compared to the southern part. The natural gas accumulated in southern zones was migrated through Porjianghaizi fault toward the northern region. Well J45 from the north zone and J77 from the south zone were chosen to compare the favorable zones of pure gas accumulation, proving that J45 lies in the pure gas zone compared to J77. Based on the faults and sedimentary facies distribution research, we suggest that the favorable zones of gas accumulation lie toward the northern region within the Hangjinqi area.


INTRODUCTION
The Ordos Basin (OB), located in central China, is one of the most productive and petroliferous craton basins with the highest annual gas output in China (Yang et al., 2005;Xu et al., 2018;Wang et al., 2020). The gas exploration in the OB started in the 1980s. Since then, numerous Permian and Ordovician gas fields have been discovered in the central and north regions of the OB. These incorporate Sulige, Yulin, Changqing, Jinbian, and Daniudi gas fields (Duan et al., 2008;Yang et al., 2015a). Our study is focused on the Hangjinqi area, which is situated in the Yimeng Uplift and northern Yishan slope, covering an area of 9825 km 2 within northern Ordos Basin   (Figure 1A). In recent times, exploratory wells within the Hangjinqi area found a considerable volume of gas in the sandstones of the Upper Paleozoic era (Xu et al., 2018). Hangjinqi area offers essential geologic conditions for making a large-scale gas reservoir, and thus is a potential area for becoming a vast gas exploration zone . The Lower Shihezi Formation (LSF) of the Middle Permian Period is considered the most significant reservoir of the OB (Sun et al., 2009). We are focusing on understanding the controls of faults and sedimentary facies distribution with the help of geological modeling, seismic stratigraphy, seismic attribute analysis, and petrophysical analysis on the second and third members of LSF in the Hangjinqi area.
The published literature on the sedimentary facies and faults distribution of the second and third members of the LSF within the Hangjinqi exploration area is still scarce. In addition, the growth history of fault distribution of central fault systems in the OB is still poorly understood (Xu et al., 2018). The main issues in identifying and understanding the controlling factors of faults and sedimentary facies distribution are as follows: (1) seismic amplitude resolution is limited (Liu et al., 2007), (2) presence of coal layers present beneath LSF , (3) uplifting during Paleozoic Era causes various lateral variations which were further triggered during Yanshan structural movements (Sun et al., 2009). Numerous faults were established, which made the study area more complicated. Former studies conducted within the study area were focused on the structural evolution (Yang et al., 2015a;Yang et al., 2015b), tectonic controls , channel identification Anees et al., 2022), evaluation of genetic types and source of Upper Paleozoic tight gas , hydrothermal mineralization of uranium deposits , and evaluation of clay mineral content, and type of reservoir properties .
The current study aims to understand the controlling factors of faults and sedimentary facies distribution for the identification of favorable zones for gas accumulation of the second member (T9e) and third member (T9f) of the Middle Permian sediments of the Lower Shihezi Formation in the Hangjinqi area, northern Ordos Basin, China. We aimed to utilize the vast post-stacked 3D seismic and well log data along with cores to better understand facies distribution through seismic attributes, seismic stratigraphy, petrophysical analysis, and geological modeling. We also aimed to put forward new facts about the sedimentation of narrow Permian sections. This study will also demarcate the favorable gas zones and channel recognition through structure maps, thickness maps, sand-ratio maps, and depositional facies map. We anticipate that our research will offer more insights within the second and third members of the Middle Permian sediments of LSF within the Hangjinqi exploration area.

Structural Characteristics
The Ordos Basin (OB) was established in the west margin of the North China Block (NCB) of the Paleozoic Era on the Archean-Proterozoic basement . During the Late Cretaceous, the OB was uplifted, resulting in a plateau . However, due to large-scale strike-slip events during the Cenozoic Era, the OB was separated from the NCB, and numerous graben structures were formed around the basin (Zhang et al., 1998;Liu et al., 2020). The basin consists of six tectonic regions and five secondary structural regions characterized by fluctuations in the structure at the top surface of the basement . The six tectonic units include Weibei uplift, Yimeng uplift, Western edge overthrust belt, Jinxi flexural belt, Yishan ramp, and Tianhuan depression (Hao et al., 2007;Liu et al., 2020). The study area of Hangjinqi is situated at the junction of the Yishan ramp and Yimeng uplift. The five structural units incorporate the Gongkahan uplift, the Hangjinqi fault, the Wulangeer uplift, the Yishan slope, and the Tianhuan depression. The Hangjinqi area is considered a promising region for hydrocarbon migration because of the paleohigh in northern OB. From the west to east region of the study area, three main faults are disseminated named as Wulanjilinmiao fault, Sanyanjing fault, and Porjianghaizi fault (PF) (also referred to as Boerjianghaizi fault in the literature) . These three faults have W-E to SW-NE dipping orientations . The PF is the main and the largest fault that separated the Hangjinqi area into north and south regions. The PF experienced various evolution stages, including formation stage, thrusting stage, and tearing stage during Caledonian Yanshanian, Indosinian Yanshanian, and Late Yanshanian (Li et al., 2015).

Stratigraphic Characteristics
The Upper Paleozoic rocks are comprised of Carboniferous Taiyuan Formation (C 3 t), Lower Permian Shanxi Formation (P 1 s), Middle Permian Lower Shihezi Formation (P 2 x), Middle Permian Upper Shihezi Formation (P 2 s), and Upper Permian Shiqianfeng Formation (P 3 s) (Ritts et al., 2004). The coal and mudstones of Taiyuan and Shanxi Formations are the primary source rocks . The Shanxi Formation (P 1 s) comprises of Lower Shanxi member (P 1 s 1 ) and Upper Shanxi member (P 1 s 2 ). The Lower Shihezi Formation (LSF) (P 2 x) consists of three members. These members are categorized as a first member (P 2 x 1 ), second member (P 2 x 2 ), and third member (P 2 x 3 )  ( Figure 1B). The tight sandstones of LSF are acting as the main reservoir within the study area, incorporating sandy fluvial conglomerate and coarse-to-finegrained sandstones (Zhang et al., 2009;Anees et al., 2019). These tight sandstones have low porosity (less than 12%) and low permeability (less than 1 mD) (Dai et al., 2012). However, these values of tight reservoir sands may differ due to lithological variations (Yang et al., 2005). The thick mudstones of Upper Shihezi and Shiqianfeng Formations are main seal rocks and consist of thick lacustrine layers of mudstone and sandstone intercalations .

MATERIALS AND METHODS
Materials used in our study to understand the faults and facies distribution for identifying favorable gas zones of second and third members of LSF include cores, post-stacked seismic grid, and well log data. A total of 62 wells were utilized to evaluate the most favorable zones of gas accumulation of the Hangjinqi area. Several geophysical logs such as gamma-ray (GR), spontaneous potential (SP), deep resistivity (LLD), shallow resistivity (LLS), caliper (CAL), and sonic (DT) logs were incorporated (Ali et al., 2020;Ali et al., 2021). The regional seismic grid of about 2500 km 2 is included to study the distribution characteristics and favorable gas accumulation zones. A total of 2280 inlines and 1880 crosslines are used. The intervals between the inlines and crosslines are 25 m. The cores samples from wells J45 and J77, including fine sandstone, coarse sandstone, and mudstone, were utilized in our study. The wells J45 and J77 were chosen based on their location nearby braided channels.
The methodology presented in our study incorporates the following steps: Initially, seismic and well data were employed to establish synthetic seismogram, and targeted horizons were marked. The two interpreted horizons are interpolated and utilized to study the faults distribution, topography, and thickness of tight reservoir sand with the help of geological modeling. Seismic attribute analysis interpreted braided channels' lateral distribution and geometry. Seismic sections were utilized to study the behavior of sediment distribution and channel distribution. Finally, we establish the depositional facies distribution map of targeted T9e-T9f horizons to interpret the favorable zones of sand deposition by marking the favorable gas accumulation zones. Geophysical logging and their geological cores through wells J45 and J77 were used to evaluate gas and water volume. These wells were finally employed to identify the favorable zones of gas accumulation.

Geological Modeling
Geological models are an effective and advanced tool to image the subsurface facies. Geoscientists are continuously developing geological modeling techniques for better prediction and simulation of subsurface systems (Radwan, 2022). The geological model is the numerical equivalent of a 3D geological map that gives the targeted subsurface geobody. Structural maps, thickness maps, and sand-ratio maps of Frontiers in Earth Science | www.frontiersin.org February 2022 | Volume 9 | Article 822670 targeted T9e and T9f horizons are created to study the distribution characteristics of faults and sedimentary facies to demarcate the favorable zones.

Seismic Attribute Analysis
Sedimentary facies classification through seismic attribute analysis has been proficiently employed in geosciences for decades (Wang et al., 2017). Seismic attributes provide valuable, reliable information regarding bedding thickness, continuity, fractures, porosity, depositional environment, lithofacies, and sequence boundaries (Chopra and Marfurt, 2005;Abbas et al., 2019). Various horizon attributes were computed between targeted T9e and T9f horizons to identify the sedimentary facies such as fluvial channels. The Root Mean Square (RMS) attribute and Continuous Wavelet Transform (CWT) method of spectral decomposition provided the best results and were incorporated in our study.

Petrophysical Analysis
For a tight sand reservoir region of our study area, the fluid information within the reservoir can be ambiguous because the logging curve is disturbed by the borehole environment and formation (Jiang et al., 2017). The relationship between log data and test data is not accurate, and the corrected logging curve data depends on the test data for standard inspection (Xiang et al., 2021). Consequently, consistency amid log data and test data is based on the preliminary logging interpretation, gas logging, and production data. Integrating the test results with LLD, porosity, and other gas logging data of the targeted interval, the high consistency intervals are chosen for the standard sample with enhanced test quality. In the end, we executed the standard inspection to interpret the gas-water and pure gas layers within the targeted horizons of the LSF.

Depositional Facies Map and its Indication
Facies maps are drawn at an isochronous surface. The depositional facies map is an effective tool to interpret the distribution of sedimentary facies. The aim of incorporating the depositional facies map in our study was to calibrate the favorable sand distribution and accumulation zones. The depositional facies map of the targeted T9e-T9f horizons was made by integrating the resulted information of the RMS attribute map, structural maps, thickness distribution map, and sand-ratio maps of T9e and T9f.

Seismic-To-Well Tie and Horizon Interpretation
We have interpreted the horizons based on the stratigraphic column model. Sedimentary deposits are represented as stratigraphic sequences, each constituted by a vertical succession of parallel layers. Seismic and well log data were jointly used to interpret and identify the different horizons with more clarity. Synthetic seismogram was calibrated via seismic-to-well tie (Ashraf et al., 2016;Anees et al., 2017), which is a process of acquiring the association amid seismic reflections and geological features. The attained synthetic traces were matched with the corresponding seismic traces to build the connection among time and depth (Ashraf et al., 2020b). Two reflectors were marked named as T9e (P 2 x 2 ) and T9f (P 2 x 3 ), which denote the second member and third member of the LSF (Figure 2). The interpolations of the interpreted horizons show that the area is tectonically unstable. The study area is separated into northern and southern parts by a central Porjianghaizi Fault (PF) (Figures 3A,B). The targeted horizons T9f (P 2 x 3 ) and T9e (P 2 x 2 ) are also shown on the uninterpreted and interpreted seismic sections near well J77 ( Figures 4A,B).

Distribution Characteristics of Fault in the Hangjinqi Area
The Hangjinqi Fault Zone (HFZ) consists of two different wide relay ramps (RRs) in the fault trace in the middle of the data coverage. The present-day fault trace is discontinuous for over 200 km (about 150 km within the Hangjinqi area) (Yang et al., 2013). Three main segments of HFZ, east (PF), central (WF), and west (SF), coincide with left steps in the fault geometry (Yang et al., 2013). HFZ consists of three main types of faults; there are also three levels and three periods of these faults. The three stages of faults prominently developed in the east Hangjinqi area are as follows: (1) (3) late Yanshan early fault (weak extension and positive inversion) developed in the PF zone, the Lijiaqu fault zone, and the Shiguhao and Azhen areas (Yang et al., 2015b). These three levels were characterized based on the extension of the faults developed toward the east of the Hangjinqi area. The tectonic activity in the study area is controlled and triggered by the main fault of PF, which is extending laterally within the northern OB. The faults distribution across horizons T9e and T9f show evidence of more tectonic activity at the north of the PF than the south region of the study area ( Figures 5A,B). In addition to significant regional faults, there are also local minor faults. There are substantial differences in the development of structures on both sides of PF. There are many minor faults in the north of the PF, while the faults in the south are sparse, and minor faults are mainly locally developed, mainly in the Lower Carboniferous-Permian strata.

Structural Mapping
A structure map shows the topography of the unit. The contours within the maps indicate the relief of the specific strata, usually beneath the surface. The structural styles such as syncline or anticline, faults and other geologic features in space are displayed. The structure map shows the topography of targeted T9e, and T9f sedimentary layers lie toward the east of the Hangjinqi area.
The structural maps of Middle Permian layers of T9e and T9f show a deepening trend toward the southern region, having depths of about 3000 m. The structural topography is divided into uplifted northern and deepened southern regions through the PF. The northeastern region shows stratigraphic thinning   with depths of about 1800-1900 m. The regions of the north of the structural maps show more deformation due to frequent lowto-medium rank faults that offset the sedimentary layers within the northern area. However, there are fewer faults in the southern region. The northeastern regions of T9e and T9f are uplifted compared to other maps' areas ( Figures 6A,B). The overall results indicate that the preliminary shortening deformation and uplifting took place during this period.

Thickness Mapping
Thickness mapping is an essential tool in geosciences to identify the thickness variations and thickness trends of different layers. T9e was marked at the bottom of the second member of LSF and T9f at the top of the third member of the Lower Shihezi Formation, so a thickness map between T9e-T9f (Lower Shihezi 2 + 3) was generated. The resulting thickness map shows significant thickness variations across the Hangjinqi area, large in the northern area compared to the southern region. The thickness is lesser in the northwestern (NW) area than in the south region, which indicates erosion in the NW region due to uplifting. The uplifting at the NW region suggests the origin of the braided channels. These results indicate a gentle slope from the north toward the south (Figure 7). Corollary, the results show that the braided channels originate from NW and flow toward the southern region.

Sand-Ratio Mapping
The sand-ratio provides the distribution of the sand to mud ratio. Sand-ratio distribution maps of the second member (T9e) and third member (T9f) were plotted with the help of structure and thickness maps to study the favorable zones of tight sand accumulation. The light blue color shows that more sand is deposited in these regions and the regions with yellow color indicate that more mud is deposited in these areas. The results of sand-ratio maps of targeted members show that maximum sand was deposited toward the northwestern region, which flows toward the south direction ( Figure 8). These maximum values of sand-ratio indicate braided channels. In addition, the extreme  eastern region also shows maximum sand-ratio values, indicating the favorable zone for braided channel deposition (Figures 8  and 9). The result of the sand-ratio map of T9e and T9f almost shows similar results, and their response agrees with each other. The corollary, the maximum sand-ratio values indicate that the braided channel originates from the north and flows toward the southern direction.

Sediment Distribution via Seismic Stratigraphy
The OB is saturated with numerous sedimentary layers, ranging from Paleozoic to Cenozoic Era. Six major unconformities (T3, T4, T5, T7, T8, and T9) were marked on the seismic data to show the behavior of sediment distribution characteristics ( Figures  10A,B). The Cambrian and Devonian rocks are missing within the Hangjinqi area. In addition, Lower Jurassic, Upper Jurrasic, Upper Cretaceous, Paleogene, and Neogene strata are also missing ( Figure 1B). These sediments of different ages are missing because of regional uplifting and hiatus in the study area (Yang et al., 2015b). The strata between T4 and T3 consists of Middle Jurassic Zhiluo (J 2 z) and An'ding (J 2 a) Formations. The J 2 z is composed of sandstone, while J 2 a is composed of mudstone facies. The sediments amid T5 and T4 are also composed of Middle Jurassic Yan'an (J2y) Formation consisting of sand and gravel facies. The Middle Jurassic sediments were deposited in the fluvial and lacustrine environments (Yang et al., 2015b). The rocks among T7 and T5 are comprised of Upper Triassic Yachang (T 3 y) Formation that varies between sandstone and mudstone facies. The sediments between T8 and T7 are composed of Lower-to-Middle Triassic Liujiagou (T 1 l), Heshangou (T 1 h), and Ermaying (T 2 e) Formations. The lithologies of T 1 l and T 1 h sediments vary from muddy siltstone interbedded with sandstone facies, whereas T 2 e consists of mudstone facies. The Permian to Triassic sediments show an fluvial, deltaic, and lacustrine depositional environment (Yang et al., 2013). The sediments between T9 and T8 are composed of Upper Permian (Shiqianfeng (P 3 s) Formation) and Upper Shihezi (P 2 s) Formation of Middle Permian age. P 3 s are composed of thick beds of mudstone interbedded with fine-sand facies, whereas P 2 s are composed of thick lacustrine beds of mudstone with thin intercalations of sandstone facies . The sediments above T9 consist of Pennsylvanian Taiyuan (C 3 t) Formation, Lower Permian Shanxi (P 1 s) Formation, and Middle Permian Lower Shihezi (P 2 x) Formation (LSF). Several coal beds are interbedded with mud and sand facies within the C 3 t and P 1 s. The LSF consists of three members (P 2 x 1 , P 2 x 2 , P 2 x 3 ), whereas the area of interest of our study is the second member and third member marked as T9e (P 2 x 2 ) and T9f (P 2 x 3 ) within

Distribution Characteristics of Channels Through Seismic Attribute Analysis
Seismic attribute analysis was performed to identify the braided channels between the targeted T9e-T9f horizons. Several amplitude and instantaneous attributes were applied to get the optimized results. The Root Mean Square (RMS) horizon attribute and Continuous Wavelet Transform (CWT) method of spectral decomposition provided the best results, and hence are included in the results. The horizon RMS attribute map confirms the presence of numerous braided channels within the study area. The sinuous zones having relatively high RMS values showing narrow parallel lines were interpreted as channels (highlighted by red arrows in Figure 11A). The RMS values are higher above the major PF (highlighted by the green arrow) than the lower region.  (Li et al., 2017). Hence, a 45-Hz slice was included in our results. The dark-colored narrow sinuous features running parallel on the map were interpreted as channels which were marked by red arrows on the map ( Figure 11B). We have also interpreted the channels on the corresponding seismic sections at various locations for a reliable interpretation. We have chosen two locations, each from northern (J45) and southern (J77) region separated by PF (highlighted by light red squares). These wells were chosen because of their optimal location nearby the channels on the RMS attribute map. The corresponding location of J45 on the seismic section also confirms the convex bottom, which suggests undercutting close to the targeted T9e-T9f horizons ( Figure 11C). In addition, the corresponding location of J77 also shows convex bottom seismic features of the channels ( Figure 11D). The associated vertical sections of J45 and J77 confirm the presence of channels which show the reliability of the study.

Comparison of Favorable Zones via Formation Evaluation
There is one perforation interval within the well J45 ( Figure 12). The depth of the perforated layer ranges from 2483.5 to 2487.5 m. The total thickness of the perforated interval is 4 m. The average porosity value of the interval is 14.87%. The average permeability value of the interval is 2.14 mD. The perforated layer is less affected by mud intrusion and diameter expansion. The log curves of the reservoir section show medium resistivity, low GR, low SP, medium density, relatively low AC, and high porositypermeability values. The test gas production volume is 4375 m 3 /day, and the water volume is 0 m 3 /day. The gas logging curve of total hydrocarbon value (TG) is maximum close to reservoir depth. These log interpretation results suggest that the reservoir is a pure gas layer. Henceforth, the reservoir within J45 of the Lower Shihezi Formation belongs to a good gas reservoir. There are three perforation intervals (marked as 1, 2, and 3) within the well J77 ( Figure 13). The depths of the perforated layers range from 2609 to 2611.5 m for interval-1, 2614-2616 m for interval-2, and 2619-2621 m for interval-3. The thicknesses of these three perforated intervals are 2.5 m, 2 m, and 2 m. The average porosity values of these intervals are 10.05%, 11.1%, and 9.71%. The average permeability values of the three intervals are 1.40 mD, 1.43 mD, and 1.21 mD. The perforated layers are less affected by mud intrusion and diameter expansion. The log curves of the reservoir section show high resistivity, low GR, low density, relatively high AC, and low porosity-permeability values. The test gas production volume is 3547 m 3 /day, and the water volume is 2.74 m 3 /day. The gas logging curve of the total hydrocarbon value is close to 10 at the reservoir level. These log interpretation results suggest that the reservoir is not a pure gas layer but a gas-water layer. Henceforth, the well J77 within the LSF belongs to a tight sandstone reservoir.
The results of the geological cores of wells J45 and J77 from the second and third members of the LSF are also incorporated to evaluate the distribution of sedimentary facies. The primary lithologies in the study area are sandstone and mudstone ( Figures 14A-H). The cores of J77 at depths 2618.89-2619.08 m shows the presence of light gray gravel-bearing coarse sandstone (i.e., quartz gravel and lithic gravel) ( Figures 14A,B). The gravel content is increasing from top to bottom. The facies consists of light gray gravelly coarse sandstone and gray mudstone at a depth of 2619.37-2619.53 m (Figures 14C,D).
The cores of well J45 at depths 2488.35-2488.56 m show that the structural features of the upper light gray medium sandstone are in gradual contact with the middle and lower light gray fine sandstone ( Figure 14E). The cores at depth 2483.18-2483.36 m show the presence of light gray gravelly medium sandstone ( Figure 14F). The cores at depths 2488.76-2488.92 m show that the gray mudstone is 2 cm away from the top, while light gray gravelly medium sandstone is at the lower part ( Figure 14G), while depths 2486.97-2487.09 m show the presence of gray mudstone ( Figure 14H).
The sand has moderately sorted sub-angular grains with compact cementation. It shows horizontal-bedding to trough Frontiers in Earth Science | www.frontiersin.org February 2022 | Volume 9 | Article 822670 11 cross-bedding of normal grading sequences, which shows the characteristics of braided channels (Lin et al., 2017). The gradual contact between sandstone and mudstone facies indicates changes in facies from the fluvial channel to the floodplain. The presence of different colors in mudstone facies shows different depositional environments (Jiang et al., 2011). The brown color suggests that the mudstone was deposited in the oxidizing depositional environment, while gray color in mudstone indicates a reducing environment.

Distribution Characteristics of Sedimentary Facies of Lower Shihezi Formation
The northern region of OB was uplifted for an extended time due to the presence of the Yinshan fold belt, which affected the sedimentary units of the Late Carboniferous Taiyuan-Permian Shihezi Period (Yang et al., 2015a). Due to strong tectonic activity, many regional faults were developed that controlled the depositional pattern of the Taiyuan Formation Period. During early Permian times, the delta-plain depositional environment was developed from the northern to the southern region. The sedimentary filling configuration got thinner in the north area and thicker in the south. The pattern goes on, and the Lower Shihezi-Upper Shihezi Formation revealed an overall filling-up effect on the uplifted depressions due to solid tectonic activities.
During the Late-Carboniferous-Early Permian Period of the Taiyuan and Lower Shihezi Formation, the depositional environment of the basin completely changed to continental facies from marine facies. During the early stage of the LSF, the tectonic activity was intense toward the northern part of OB. Due to this, the supply of the provenances was adequate, rivers were developed, and the grain size of the sandstones was mainly coarser, i.e., gravel-bearing sandstone. During the late stage of the LSF, the tectonic activity was less intense; sedimentary slope decreased, river construction was weakened which gave rise to the distribution of mediumto fine-grained sediments. Henceforth, the depositional pattern of the LSF consisted of braided river deposits.
The structural maps showed that the paleogeomorphology at the NW direction was uplifted, and a gentle slope was formed toward the SE direction ( Figures 6A,B). The distribution of the sedimentary facies is mainly governed by paleogeomorphology. The integrated results of structural, thickness, and sand-ratio mapping show that the sedimentary facies of LSF were distributed in braided river deposits. The main sedimentary subfacies were channels and floodplains. The sand bodies were deposited in the river channels and the sub-aqueous braided channels. At the same time, the mudstone facies were deposited in the floodplains. The braided channels originate from the NW direction and flow in the SE direction ( Figure 15).

Identification of Favorable Zones Within Hangjinqi Exploration Area
The current study is focused on identifying the favorable zones of gas accumulation that can be targeted for future studies. The sand-ratio, RMS attribute, and thickness maps indicated the presence of uplifted paleohigh regions toward the northern section of the Hangjinqi area. Due to this uplifting, the accumulated gas in the southern region was affected and migrated laterally toward the north of the region due to the presence of open areas in the PF . In our study, we have chosen two wells (J77 from the southern area and J45 from the northern region) to identify the favorable zones in the exploration area. We studied the formation evaluation of these two wells to compare the test results for promising zones (Figures 12 and 13). The test results of J77 showed that the reservoir region is composed of a mixed gas-water layer, whereas the test results of well J45 revealed that the reservoir zone is Frontiers in Earth Science | www.frontiersin.org February 2022 | Volume 9 | Article 822670 13 composed of pure gas. Hence, we proposed that the northern region of the PF is more favorable for gas exploration ( Figure 15).

CONCLUSION
The research has been conducted on the Middle Permian sediments of the second and third members of the Lower Shihezi Formation interpreted as T9e (P 2 x 2 ) and T9f (P 2 x 3 ) within the Hangjinqi area, northern Ordos Basin, China. The fault distribution characteristics of T9e-T9f showed that the significant Porjianghaizi Fault mainly controls the study area, and many minor faults developed in the lower Carboniferous-Permian strata. The structural maps revealed a deepening trend of about 3000 m in the southern region while thinning trend of about 1900 m in the northern region. The thickness map also suggested that the thickness is lesser in the north of the region compared to southern areas, which revealed that the NW region is more uplifted. These results were supported by sand-ratio mapping and RMS attribute map, which also revealed the presence of maximum sand deposition toward the NW region.
Cores of wells J77 and J45 show the trough crossbedding of gravel-bearing coarse sandstone and mudstone, which indicated the presence of braided channels confirmed by the RMS attribute map and spectral decomposition of a 45-Hz slice. The locations of J45 and J77 wells also confirmed the presence of a convex bottom and undercutting on the seismic sections close to the targeted T9e-T9f horizons. The integrated results of drilling cores, attribute analysis, petrophysical analysis, and geological modeling confirmed the presence of braided channels in the Lower Shihezi Formation. The test results of wells J45 (NW region) and J77 (southern region) showed that the NW region is more favorable for pure gas exploration than the south region. Corollary, we proposed that more attention should be given to the NW region for pure gas exploration.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

AUTHOR CONTRIBUTIONS
AA, WS, and UA conceived this study. WS supervised this project. RW and KL undertook the responsibility of arranging the data for this project. AA and UA wrote the manuscript. HM and RJ helped in formal analysis and software. XZ and QL discussed the data. ST and HZ reviewed the manuscript and provided the necessary fundings. All authors contributed to the article and approved the submitted version.