Sedimentological Influence on Physical Properties of a Tight Sandstone Reservoir: The Cretaceous Nenjiang Formation, Southern Songliao Basin, Northeast China

Abstract: The Nenjiang Formation, south of Songliao Basin, has many hydrocarbon bearing units, but currently the understanding of the distribution of viable reservoir sandstones is too limited to support a development strategy. Therefore, a detailed study on the sedimentary microfacies and reservoir properties has been completed in order to reduce uncertainty and improve subsurface predictions. Nine lithofacies and five lithofacies associations were identified supporting the development of a sedimentary model of a river-dominated delta front setting, which could be divided into four sedimentary environments: subaqueous distributary channel-fill, mouth bar, sand sheet, and interdistributary bay. The distribution sandbodies extend to the south in a tongue-like form, and they thin and pinch out. Finally, the influence of sedimentary process on properties was assessed by establishing the correlation between microfacies and reservoir physical parameters, such as porosity, permeability, pore radius, throat radius, and clay minerals. It is revealed that the correspondence between reservoir physical properties and microfacies types is strong; the physical properties of the subaqueous distributary channel and mouth bar are the best.


INTRODUCTION
Lithofacies are identified by mineral composition, structure, color, and sedimentary structure of rocks, which is directly related to the sedimentary environment and facies of silicic clastic rocks. Early lithofacies type analysis and application were based on fluvial facies sandbodies (Miall, 1985). The changes of lithofacies characteristics reflect the changes of various conditions during the formation of sandbodies of different scales, such as hydrodynamic conditions, sediment transport mode, and the degree of external interference (Qi et al., 2009;Wu D. et al., 2019). The combination of these single lithofacies forms lithofacies associations, which is the result of multistage sedimentation and superposition (Eyles et al., 2010;Ding et al., 2014). Lithofacies associations and the arrangement of architectural elements can be used to interpret the climatic and environmental conditions during deposition and the paleohydrodynamic conditions (Jin et al., 2019).
According to the literature, many researchers have carried out conventional, but not detailed, sedimentary facies research in the Songliao Basin. These studies were usually at a large scale and the fourth member of the Cretaceous Nenjiang Formation was presumed to be a single prograding sequence (Liu et al., 1993;Sun et al., 2006). Most authors interpreted that tectonic movement controls the sedimentary facies types in the southern Songliao Basin (Li et al., 2009;Huang et al., 2013;Wang J. et al., 2014); the sedimentary facies types of the Nenjiang Formation in the Songliao Basin are gravity flows, alluvial fans, deltas, and rivers Li et al., 2015;Lv et al., 2016;Wang D.-d. et al., 2016). In this sedimentary period, there are two main river deltas, one comes from the northeast, which plays an important role in the sedimentary process of the first member of the Nenjiang Formation, the other comes from the east Wang W. et al., 2016). Other interpret that the Nenjiang Formation in the Songliao Basin is a mud-prone delta with poor continuity, which has minor sand components and a lower angle, compared with the traditional delta. In addition to delta sedimentary environments, other sedimentary environments of the Nenjiang Formation include deep and shallow lakes (Gao and Wang, 2010;Wang L. et al., 2014) and even deep-water turbidite fans (Pan et al., 2017;Wang et al., 2018). Recently, some researchers established a new ultrafine-grained sedimentary model in the Songliao Basin: the channel fan, which includes straight and sinuous channels extending from north to south. This system is dominated by fine sediment and contains a large amount of sandstone, fluvial sedimentary structure, and internal erosion surfaces. In addition, the system also contains abundant terrestrial organic debris, showing bedload and suspended load transport (Tong et al., 2018;Mo et al., 2019). A series of diagenetic processes will occur after sediment deposition, which will directly affect the size and sorting of reservoir pores (Guo and Mao, 2019). Compaction and cementation greatly reduce the porosity of the reservoir, and at the same time the pore throats are blocked and the rock permeability is also greatly reduced, which is the main factor behind the formation of tight reservoir (Shi et al., 2015;Yang et al., 2017). However, in the middle and late diagenesis, some organic acid solutions will dissolve mineral particles and cements in pores, thus expanding the pores of rocks and improving its permeability. The main minerals of dissolution are potash feldspar and calcite (Wang et al., 2020). However, the location and degree of dissolution are quite different, which will lead to a different degree of dissolution of rock particles and cements, resulting in less strength in the correlation between porosity and permeability, and being prone to high porosity, low permeability, or high permeability and low porosity .
Many research results have made the structural evolution characteristics and sedimentary environment of Changling fault depression in the south of Songliao Basin clear, which provide a wealth of sources for this study (Meng et al., 2016;Xu et al., 2019;Cai et al., 2020). However, the characterization of sedimentary facies and evaluation of the fourth member of Nenjiang Formation in the south of Songliao Basin are not sufficient and need detailed research Zhang and Wang, 2019). Here, we refine the characterization of microfacies in the areas with high well density to establish corresponding lithofacies models (either single or composite lithofacies model) and investigate the influence of sedimentation and diagenesis on reservoir physical properties. Through the study of fine reservoir description, we establish a new integrated sedimentary model and lithofacies model, discuss the factors controlling the sandstone physical properties, and play a positive guiding role in the effective development of this reservoir.

General Situation of the Study Area
Songliao Basin is the largest inland sedimentary basin fill of Mesozoic Cenozoic in eastern China, which is divided into two parts, the south and the north, with Songhua River as the boundary. The southern Songliao Basin comprises four firstorder structural units: western slope, central depression, southeast uplift, and southwest uplift, among which the central depression and southeast uplift are important geological units for hydrocarbon exploration ( Figures 1A,B) (Li et al., 2007;Zhang and Huang, 2010). The study area (H168 block) is located in Qian'an County, in Songyuan City of the Jilin province northeast of China, and its regional tectonic location is in the south of the Changling depression of the southern Songliao Basin. The structural characteristics of the Changling sag are relatively gentle. Its tectonic evolution has gone through three stages: early fault subsidence, middle subsidence, and tectonic quiescence. It mainly includes the Qian'an subsag and Heidimiao subsag, and the Hei168 block is located in the Heidimiao subsag, its east part is Huazijing terrace, and its northwest part is Daan-Honggang terrace. At present, 122 horizontal wells have been drilled and 117 wells have been put into production in the study area, which is a typical thin-layer reservoir developed by horizontal wells ( Figures 1C,D).

Stratigraphic and Tectonic Evolution in the Study Area
Overall, the Songliao Basin experienced warm, humid climate during the Cretaceous period (Liu et al., 2003). The tectonic process includes three stages: early fault depression, middle depression, and late shrinkage. During the basin expansion stage (the early fault depression and middle depression stages), the Denglouku Formation, Quantou Formation, Qingshankou Formation, and Yaojia Formation were successively deposited in the study area. In the low subsidence stage (the late shrinkage stage) of the basin evolution, deltaic and lacustrine strata are Nenjiang Formation (Figure 2A), Sifangtai Formation, and Mingshui Formation. The Yanshan V episode tectonic movement in the late Mingshui Formation resulted in strong inversion and uplift, forming the present structural form. The period of the Nenjiang Formation is the last stage . The Nenjiang Formation comprises a large delta that built out into a lacustrine water body during the second flooding period Frontiers in Earth Science | www.frontiersin.org February 2021 | Volume 9 | Article 606701 (Feng and Zhang, 2012;Huang et al., 2013). The relative lake level decreased continuously during the Nenjiang Formation period except for a short rise in the middle of the fourth member, and the main sedimentary types are semideep lake, delta, and shallow lake ( Figure 2B).
Using these rock samples, a large number of micro characterization experiments have been carried out by means of optical or electronic microscope, such as thin section analysis, scanning electron microscope, and X-ray diffraction analysis experiment. On the basis of these experimental results, the micro components, pore structure, and diagenesis of reservoir sandstone have been characterized in detail, and the single and combined lithofacies models are established. Then, the characteristics of rock lithology, sandstone type, grain size, and sorting are studied in detail, the plan form and vertical distribution of microfacies are characterized, and the Frontiers in Earth Science | www.frontiersin.org February 2021 | Volume 9 | Article 606701 3 sedimentary model is established. Finally, the sandstone diagenesis in the study area is characterized by using the micro pore structure experimental parameters, and the variation law of reservoir pore and throat, as well as their internal filler type, is clarified.

Lithologic Characteristics
The lithology of the fourth member of the Nenjiang Formation is mainly mudstone, whose content is as high as 90% ( Figure 3A).  The lithology of sandstone is mainly gray medium-and finegrained sandstone, siltstone, and argillaceous siltstone, showing normal and inverse grading. The thickness range of a single sand body is 1-5 m, with an average of 2.2 m ( Figure 3B). According to Folk's three terminal sandstone classification method (Folk, 1968), the relative content of quartz, feldspar, and rock debris shows that the main types of sandstone are lithic feldspar sandstone ( Figure 3C), of which the average content of quartz is 32%; feldspar content is 31% (potash feldspar), and lithics is 37%, which includes igneous rocks, metamorphic rocks, and mica debris. All these suggest that the sediment source area is far away from the provenance area and belongs to the terminal position of the river (Wu J. et al., 2019), and the source area lithology is complex, mainly metamorphic rock and intermediate magmatic rock ( Figure 3D).

Sedimentary Structure
Through core observation, it can be found that the sedimentary structures shown in the rocks in the study area mainly include cross bedding, parallel bedding, wavy bedding, and sand ripple bedding, which are consistent with an interpretation of the sedimentary environment of river-dominated delta front ( Figure 4).
The sedimentary structures show weak hydrodynamic characteristics. The strong scour structures at the bottom of the subaqueous distributary channel are rare; therefore, the trough cross bedding in the core is not obvious. The main sedimentary structure includes planer cross bedding, parallel bedding, ripple lamination, and other sedimentary structures. The single sand body of the mouth bar is relatively thick, mainly composed of fine sand and silty sand; the sedimentary structures include current sand ripple bedding, parallel bedding, low angle cross bedding, and wave formed sand ripple bedding. The sand sheet is well sorted and is mostly fine sandstone and siltstone. There is no large-scale cross bedding. Parallel bedding and low angle cross bedding are mainly developed. There are also sand ripple bedding and wavy bedding. The biogenetic structure is not developed.
The grain size parameters of sandstone in the study area were obtained by using the casting slice observation method. The grain size of the reservoir rock in the study area is generally small, and the cumulative curve of grain size probability is mostly two-stage (composed of jump population and suspension population). The particle size of the sediment is small, but the distribution range is large. The cut-off point between the saltating population and the suspended population is between 2.5Φ and 4.0Φ. The suspended population content is mostly between 10 and 40%, and the jumping population content is generally more than 60% ( Figures 5A-C). In the C-M diagram, the sample points are concentrated in the QR segment and a few in the RS segment. The QR segment is parallel to the baseline of C M, which represents the deposition mechanism of gradual suspension. In the fluid, the particle size of suspended materials gradually becomes finer from the base to the top, which is often caused by the development of turbulence, and its C and M increase in proportion. RS section is a uniform suspension, which is a complete suspension with no change of particle size and density with depth, reflecting a weak and calm sedimentary environment (Figures 5D).

Lithofacies Characteristics
The grain size range of sediments is relatively small, the types of sedimentary rocks are relatively unified, and the sedimentary sequences formed are relatively similar. It is difficult to distinguish their microfacies types from rock types alone and to understand the relationship between different microfacies types. Therefore, this study identifies sedimentary microfacies based on the combination of lithofacies characteristics. On the basis of identified lithofacies types, a model of single lithofacies and composite lithofacies is established.
(1) Single Lithofacies Nine lithofacies types were identified in the delta facies according to the characteristics of rock fabric, bedding, and grain order of the study area Kallepalli and Richardson, 2017); the morphology and characteristics of each lithofacies are shown in Figure 6.
These lithofacies include (1) the massive fine sandstone facies (FIM), which was formed under rapid deposition; (2) trough cross bedding sandstone facies (FIT), which was formed by migration of 3D dune-scale bedforms at the base of subaqueous distributary channel; (3) oblique bedding fine sandstone facies (FIO), which was usually formed in the deposition environment of mouth bar or sand sheet; (4) lowangle cross-bedding fine sandstone facies (FIL), which was usually formed by deposition in subaqueous distributary channels with weak current action; (5) inverse graded sandstone facies (SI), which was usually formed by waxing hyperpycnal flows that deposit on the mouth bar; (6) ripple cross lamination fine sandstone facies (PR), which was formed in a relatively quiet water environment and mainly composed of argillaceous siltstone; (7) ripple cross lamination and lenticular bedding siltstone facies (SiC), which was usually formed by the interdistributary bay setting; (8) ripple cross lamination siltstone facies (SiR), which was formed in the front of the river controlled delta with weak hydrodynamic force (Gugliotta et al., 2015;Xu et al., 2020); and (9) the horizontal bedding mudstone facies (MH), which was formed by suspension deposition in the interdistributary bay-fill.

Sedimentary Facies and Microfacies Characteristics
Based on the analysis of lithological characteristics, sedimentary structure, and other sedimentary facies markers, combined with the sedimentary background and environment of the fourth member of the Nenjiang Formation of the Songliao Basin, the types of sedimentary facies are determined as a transgressive delta. By analyzing the lithology, rock color, lithofacies type, and other facies indicators, it can be concluded that the study area is located in the delta front. The delta front can be divided into four different parts in the study area: subaqueous distributary channelfill, mouth bar, sand sheet, and interdistributary bay-fill. The log response characteristics, curve shape, sand body shape, and lithology of each sedimentary microfacies sand body are very different, which are their main distinguishing features ( Table 1).

(1) Subaqueous Distributary Channel-Fill
The distributary channel is the subaqueous extension part of the onshore branch channel. In the process of extension, the channel widens, and the bifurcation increases, the flow velocity decreases, and the deposition rate increases. The lithology is mainly grayish white fine to medium sandstone, with trough cross bedding, parallel bedding, and small trough cross bedding which can be seen from the bottom to the top; SP log curve is box shaped; sand body vertically presents a positive rhythm of coarse and fine from the bottom to the top. The subaqueous distributary channel-fill thickness in the study area is generally 2-4 m. The base of these sand bodies is an erosional overlain by trough cross bedding and muddy gravel. The middle is oblique bedded sandstone and the upper is laminated sandstone. The top is sand ripple bedding with an increased mud content ( Figure 8A).
(2) Mouth Bar The mouth bar formed after the branch channel entered the lake where the distributive channel bifurcates. Most of the mouth bar in the study area coarsens upward, increasing sand content and better sorting. SP log curve of the mouth bar sand body is funnel-shaped with abrupt change at the top and gradual change at the base. The mouth bar of the study area is mainly composed of fine sand and silty sandstone; its thickness is generally more than 2.5 m. The base of the sand body has more gradual contact with the underlying rock layer and the top contact is abrupt with the overlying mudstone layer. The bedding structures include small cross bedding, parallel bedding, and low angle oblique bedding ( Figure 8B).

(3) Sand Sheet
The sand sheet is distributed in the front and flank of the mouth bar. It is characterized by a wide distribution area, but the particle size becomes smaller and the sandbodies become thinner. It is mostly distributed in the front of the delta front, where the hydrodynamic conditions are weak. The shape of SP log curve is usually fingerand cone-like, and the amplitude is generally high and medium. The sand sheet sandstone is well sorted; its lithology is mostly argillaceous siltstone and siltstone, with a thickness of 1-2 m. Small cross bedding, parallel bedding, and flow sand ripple bedding are usually found in the sand sheet, in addition to some biogenic structures ( Figure 8C).

(4) Interdistributary Bay-Fill
Because of the river overflow and the change of water depth, an interdistributary bay forms. Longitudinally, it is usually composed of multiple sand bodies. There are mudstone interbeds in the middle of the sand bodies, and the thickness is uneven. In terms of internal structure, there are obvious interbeds between each single sand body and other sand bodies in the sequence. Log curves are similar through the sand bodies but have correlated peaks, which are generally tooth-shaped, conical, or elliptical, with moderate or low amplitudes, at the mudstone interbeds ( Figure 8D).

Distribution of Sedimentary Microfacies
According to the distribution law of sedimentary microfacies, the Sanzhao and Gulong depressions in the northern part of the study area are the main source areas. The progradational delta distributed sediment southward along the long axis of the basin. The delta front facies sand bodies extend to the inner part of the study area, forming a reservoir mainly composed of sheet sand. The catchment area of Songliao Basin is in the south, basinward of the delta front subfacies, and the sedimentary microfacies are dominated by the sand sheet. Most of the thin layers in the fourth member of Nenjiang Formation are poor in physical properties, and there is no effective oil-bearing layer ( Figure 9B). The k 2 n 4 6 sandbody thickness is 3-5 m, and the two microfacies, distributary channel-fill and mouth bar, account for a relatively high proportion. Therefore, its sand body has good physical properties, which is the main oil-bearing layer in the fourth member of Nenjiang Formation and the target member of our study ( Figure 9C).
Analysis of previous research results shows that the sediments of the study area are derived from the northeast (Jia et al., 2007;Huang et al., 2013;Bao and Niu 2017). Affected by this, the thickness of sandstone in the northeast is larger than that in the southwest. The thickness of sand body decreases gradually from southwest to northeast until it pinches out. However, there is no boundary of sandstone disappearance in the study area. The distribution area of sandstone here is very large, and its thickness is between 1 and 4 m, with an average thickness of 3.5 m.

Sedimentary Model
According to the structural evolution characteristics of Songliao Basin and the sedimentary facies, a sedimentary model of river delta in the study area is established (Figure 10).
The sedimentary model is derived from the interrelationship between the sedimentary environment, sedimentation, and the results (sedimentary facies). The sedimentary model not only shows the most typical characteristics of sedimentary facies, but also shows the formation mechanism and process of sedimentary facies (Ye, 2005;Feng and Zhang, 2012;Gao et al., 2017). The sedimentary system of Songliao Basin was dominated by a process of progradation; during the first member and the second member of Nenjiang Formation, the relative lake level rose and the water body deepened continuously. From the late stage of the second member of the Nenjiang Formation, the relative lake level began to decrease until the third member ended.
In the early stage of the fourth member of Nenjiang Formation, the basement tectonic activity in the southern Songliao Basin FIGURE 10 | Sedimentary model of river delta (the study area is located in the delta front area).
Frontiers in Earth Science | www.frontiersin.org February 2021 | Volume 9 | Article 606701 10 gradually weakened, which entered the active period of depression. At the end of the fourth member of Nenjiang Formation, the relative lake-level fallen and the delta continued to prograde southward, forming the deposition dominated by mouth bar and sand sheet. The fifth member of Nenjiang Formation is a period of continuous rising of relative lake level.

Analysis of the Influence Degree of Sedimentation on Reservoir Properties
(1) Because the target area is located at the end of the delta front, distal to the source area, the deposits are predominantly finegrained, and mudstone accounts for a relatively high proportion, sandstone are very fine grained, and the deposits are well sorted. Under the influence of these conditions, the pores of sandstone are relatively small, and the smaller intergranular pores are dominant (46%). In addition, due to the influence of parent rock lithology in the source area, the content of feldspar and lithics in sandstone is higher, which increases the probability of particle dissolution. Therefore, the content of dissolution pores and mold pores related to dissolution is 29 and 17% respectively. Dissolution depends largely on the cement in the pores, and the content of early carbonate cement in the pores of the rock is very high, with an average of 18%, which is easily dissoluted by the organic acid. But the clay mineral content is low (15%) mainly composed of illite/ montmorillonite (I/M) mixed layer, which has little effect on the physical properties of the reservoir and only plays a role of interstitial filling, respectively. In addition, other types of pores in the rock are relatively small ( Figure 11).
The probability distribution of sandstone pore and throat radius is obtained by using the constant velocity mercury injection experiment. The correlation between these data and FIGURE 11 | The different pore types of the sandstone (① intergranular pores; ② dissolution pore of feldspar; and ③ mold pore of the debris dissolution).
Frontiers in Earth Science | www.frontiersin.org February 2021 | Volume 9 | Article 606701 11 sedimentary microfacies shows that the physical conditions of distributary channel sandbodies are the best, which have good storage and seepage performance. The average of pore radius is 56.6 um, and the average of throat radius is 8.6 um. The physical conditions of the mouth bar sandbodies are the second, the average of pore radius is 33.8 um, and the throat radius is 5.7 um. The physical properties of the sand sheet body at the delta front are very poor, with an average pore radius of 21.5 um and an average throat radius of 3.1 um. The physical properties of the interdistributary sandbodies are the worst, with an average pore radius of 15.3 um and an average throat radius of 2.5 um, almost without storage and seepage performance. This shows that the size and distribution of pores and throats are greatly influenced by sedimentary microfacies (Figure 12).
(2) Due to the influence of sedimentary environment, the sandstone in the study area has finer grain size, more argillaceous sediments, and high clay mineral content. The sandstone interstitial material includes mud, carbonate, and clay minerals; the mud is mostly recrystallized and locally distributed in a muddy strip Li et al., 2019;Zhang et al., 2020). Carbonate cementation (calcite and dolomite) is an important cementation, which is one of the main factors to destroy the reservoir conditions in the study FIGURE 12 | Statistical chart of relationship between sandstone pore and throat size and microfacies. Frontiers in Earth Science | www.frontiersin.org February 2021 | Volume 9 | Article 606701 12 area. According to the precipitation sequence and composition changes, they can be divided into early and late stages.
The carbonate minerals in the study area are mainly early cements, which are mainly formed in the late stage of early diagenesis. The crystal morphology of calcite is generally poor, while that of dolomite is better (Figures 13A,B). The results of whole rock X-ray diffraction show that the mineral content of carbonate rock is between 5 and 25%, with an average of 18% ( Figure 13G). The early carbonate cement fills the primary pores, reduces the porosity and permeability of rocks, and destroys the reservoir performance. The clay mineral content in the fourth member of Nenjiang Formation is relatively high, and clay cementation is also the most important cementation in the study area. Clay cements include the I/M mixed layer ( Figure 13C), illite ( Figure 13D), chlorite ( Figure 13E), and kaolinite ( Figure 13F) of which the content of the I/M mixed layer is the highest, reaching 55% ( Figure 13H). This shows that the hydrodynamic conditions are weak, but later diagenesis is strong.
(3) Different microfacies lead to different types of sandbodies, which leads to great differences in porosity and permeability of rocks, so the distribution range of sandstone porosity and permeability of each microfacies can directly reflect the matching relationship between different sedimentary microfacies and reservoir physical properties. For the reservoir in the study area, the physical properties of the subaqueous distributary channel-fill and mouth bar are the best, with an average porosity of more than 13% and a permeability of more than 17 × 10 −3 μm 2 ; the sand body of the sand sheet is the second, with a porosity of 10.48% and a permeability of 9.36 × 10 −3 μm 2 ; the sandstone physical properties of the interdistributary bay-fill are the worst, with an average porosity of only 8.33% and a permeability of only 3.55 × 10 −3 μm 2 . The correspondence between the physical properties of the reservoir and the types of sedimentary microfacies is very good, which reflects the control of the depositional processes and detrital mineralogy on the reservoir performance, as have been interpreted from deep-water settings (Porten et al., 2016).

CONCLUSION
(1) The sandstone content of the K 2 n 4 strata in southern Songliao Basin is merely 10%, which indicates long transport distance for the sediments. The sedimentary facies support interpretation of a river-dominated delta front which includes four microfacies: subaqueous distributary channel, mouth bar, sand sheet, and interdistributary bay-fill. The lake basin in the study area was gradually shrinking during the deposition period of Nenjiang Formation, and the river delta continued to prograde southward, which causes the sand body to become tongue-like and gradually thinned to the south.
(2) There are 9 single lithofacies types and five lithofacies associations in the target layer of the study area. The single lithofacies mainly consist of SIT, SR and PR, and the main combination of lithofacies is MH-FIT-FIO-FIT-FIO-MH, which show the vertical superposition of positive rhythm channel sand bodies. This is consistent with the characteristics of microfacies in the study area.
(3) Influenced by sedimentary environment and distance from source area, the grain size of sediment is very fine, which determines that the primary porosity and permeability of sandstone are low. In addition, the compaction and cementation in the diagenesis stage are strong, which reduces reservoir pore space. Therefore, the pores of the sandstone are mainly formed by lithic and feldspar dissolution in the study area. Dissolution depends largely on the cement in the pores, and the content of early carbonate cement in the pores of the rock is very high, with an average of 18%, but the clay mineral content of is low (15%) which is mainly composed of illite/ montmorillonite-(I/M-) mixed layer. The absolute content of clay mineral has little effect on the physical properties of the reservoir and only plays a role of interstitial filling.

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 author.

AUTHOR CONTRIBUTIONS
JW was mainly responsible for the conception of the article. He put forward a good idea and completed the research content related to sedimentary facies in the article. KZ was in charge of sedimentary facies research and carried out rock analysis experiments. JZ gave some good ideas to the lithofacies model establishing. Xie has been engaged in the research on the reservoir property influence factors.

FUNDING
Projects ZR2020MD035 was supported by the Natural Science Foundation of Shandong Province. Projects 51504143 and 51674156 were supported by the National Natural Science Foundation of China.