Determination of Paleocurrent Directions Based on Well Logging Technology Aiming at the Lower Third Member of the Shahejie Forma

water

Article Determination of Paleocurrent Directions Based on Well Logging Technology Aiming at the Lower Third Member of the Shahejie Formation in the Chezhen Depression and Its Implications

Yangjun Gao 1,2, Furong Li 2,3, Shilong Shi 2 and Ye Chen 1,4,*

1 School of Earth Sciences and Resources, China University of Geosciences, Beijing 100083, China; [email protected] 2 Shengli Oilﬁeld Branch Company, SINOPEC, Dongying 257001, China; [email protected] (F.L.); [email protected] (S.S.) 3 Faculty of Land and Resources Engineering, Kunming University of Science and Technology, Kunming 650093, China 4 School of Water Resources and Environment, China University of Geosciences, Beijing 100083, China * Correspondence: [email protected]

Abstract: The Bohai Bay basin, mainly formed in the Cenozoic, is an important storehouse of groundwater in the North China Plain. The sedimentary deposits transported by paleocurrents often provided favorable conditions for the enrichment of modern liquid reservoirs. However, due to limited seismic and well logging data, studies focused on the macroscopic directions of paleocurrents L are scarce. In this study, we obtained a series of well logging data for the sedimentary layers of Es3 Formation in the Chezhen depression. The results indicate the sources of paleocurrents from the northeast, northwest, and west to a center of subsidence in the northern Chezhen depression at that Citation: Gao, Y.; Li, F.; Shi, S.; time. Based on the well testing data, the physical properties of the layers from Es L Formation in Chen, Y. Determination of 3 Paleocurrent Directions Based on this region were generally poor, but two abnormal overpressure zones were found at 3700–3800 m Well Logging Technology Aiming at and 4100–4300 m deep intervals, suggesting potential high-quality underground liquid reservoirs. the Lower Third Member of the By combining with other geological evidence, we suggest that the Paciﬁc Plate was retreating and Shahejie Formation in the Chezhen changing its direction from NE–SE to W–E and the Bohai–Luxi block was suffering an extrusion from Depression and Its Implications. NE induced by the Lan–Liao and Tan–Lu strike-slip faults in the early Paleogene. Water 2021, 13, 408. https://doi.org/ 10.3390/w13040408 Keywords: well logging; paleocurrent; provenance; Bohai Bay basin; Paciﬁc Plate; Tan–Lu strike-slip fault Academic Editor: Constantinos V. Chrysikopoulos Received: 16 December 2020 Accepted: 1 February 2021 1. Introduction Published: 4 February 2021 The North China Plain has one of the largest groundwater stores in the world. The

Publisher’s Note: MDPI stays neutral groundwater reservoirs are widely distributed in the sedimentary interlayers of the post- with regard to jurisdictional claims in Cenozoic age in the Bohai Bay basin (Figure1a), laying in the north of the North China Plain. published maps and institutional affil- The Chezhen depression (Figure1b) is a representative secondary sub-basin tectonic unit in iations. the Bohai Bay basin [1–3]. The study of this depression is of great help to understand the distribution of groundwater reservoirs and the tectonic genesis of the Bohai Bay basin. In previously published papers, many studies have focused on the stratigraphy, sedimentology, microtectonics, oil–gas exploration, and thermal evolution history of the eastern Bohai Bay basin [4–9]. However, due to limited seismic and well logging data, studies focused on the Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. macroscopic direction of paleocurrents are scarce, and the coupling relationship between This article is an open access article paleocurrents and the tectonic background has not received enough attention [10,11]. distributed under the terms and The migration of paleocurrents controls the sedimentary filling process and the compo- conditions of the Creative Commons sition of sediments in the basin, and it further influences the scale and physical properties Attribution (CC BY) license (https:// of underground liquid reservoirs [12–14]. In particular, the paleowater flow was always creativecommons.org/licenses/by/ accompanied by mineral transportation and sedimentary screening, leading to the depo- 4.0/). sition of mature sediments in the downstream area to benefit the modern preservation

Water 2021, 13, 408. https://doi.org/10.3390/w13040408 https://www.mdpi.com/journal/water Water 2021, 13, 408 2 of 18

of groundwater. Therefore, the direction of the provenance of sedimentary deposits is a reliable indicator of the macroscopic direction of paleocurrents [15–17]. Furthermore, the state of paleocurrents reflects the paleogeomorphologic characteristics. In general, paleocurrents always flow from higher to lower altitudes, so the paleogeomorphologic features can be restored by using the results of paleocurrents to illustrate the dynamic changes of the surrounding tectonics [18–20]. Traditional methods for provenance analysis mainly include heavy mineral analysis, clastic rock analysis, sedimentary analysis, and fission track methods [21–26]. However, the experimental results of traditional methods are easily limited by the geographical locations of the sampling sites, the number of samples, and the erosion and chemical weathering of outcrops, leading to inconclusive results. In well logging technology, data such as gamma- ray and photoelectric absorption data are the sum of the specific information obtained from the target formations buried underground. Therefore, the required data for provenance analysis can be obtained directly from the original rocks [27,28], which greatly improves the accuracy. In recent years, well logging methods have been used to determine the directions of paleocurrents and the sedimentary provenance by quantitative and semi-quantitative calculation, and good results have been achieved [29,30]. In this paper, we obtained the seismic and well logging data from seven wells in the Chezhen depression focusing on the sandstone layers in the lower part of the third member L of the Shahejie Formation (Es3 ) in the early Paleogene. The directions of paleocurrents and the sedimentary provenance in this area were determined, and then favorable locations of underground liquid reservoirs were predicted. By combining the results with the tectonic background of the Bohai Bay basin in the early Paleogene, we interpreted the geotectonic dynamic of the Chezhen depression and the block it belongs to.

Figure 1. (a) A simplified tectonic map of the Bohai Bay basin, modified from Li et al. (2011) [31]. (b) A simplified geological map of the Chezhen depression, modified from Lao et al. (2011) [3]. It should be noted that the logging data from the seven wells (yellow points) were used in most results; however, the reconstruction of the paleogeomorphologic map was performed based on all 40 wells in this region. The locations of the other 33 wells are hidden in this map following privacy policies. Water 2021, 13, 408 3 of 18

2. Geological Information Located in the east of the North China block, the Bohai Bay basin is the central region where the Paleo-Asian Ocean, Tethys Ocean, and Pacific Ocean interacted [32–35]. It is a large-scale basin where the Mesozoic-Cenozoic rocks are superimposed on the pre- Paleozoic basement of the North China craton [36–39]. The tectonic evolution of the Bohai Bay basin can be divided into three scenes: the first syn-rift and post-rift scene in the Jurassic (201.3–145.0 Ma); the second syn-rift and post-rift scene in the Cretaceous (145.0–65.0 Ma); and the third syn-rift and post-rift scene from the Paleogene (~65.0 Ma) to the present [40]. Our research was carried out in the Chezhen depression, laying at the southeast of the Bohai Bay basin (Figure1b) [ 1]. It is one portion of a series of tectonic combinations of uplifts and depressions arranged in a northwest direction; it is adjacent to the Bohai bay and is bounded by the Tan–Lu strike-slip fault in the east and the Lan–Liao strike- slip fault in the west [41–43]. The Chezhen depression is a Cenozoic continental faulted basin developed based on the long-term reconstruction of the cratonic basement [1,3]. It is a synsedimentary half-graben controlled by the Chengnan normal fault that connects with the Chengning uplift and a partially exposed Precambrian crystalline basement to the north, and the Cenozoic sedimentary overlapped cover extends to the Yihezhuang uplift to the south [44]. This pattern of W–E trending fault in the north and stratigraphic overlapping in the south may suggest a crustal stretching. The Chezhen depression can be further divided into three sub-tectonic regions: the steep south-dipping slope formed by the fault in the north, the gentle slope belt extending to the south, and the depression area formed inside [1]. The interior of the Chezhen depression has stable sedimentation and weak faulting, and only four weak secondary shear faults have developed in it. Thus, the Chezhen basin is suitable for drilling and well logging research in the depression. The Chezhen depression contains the ancient basement in the lower part and the strata of sedimentary cover generated during the three rift cycles in the upper part [1,3]. The basement is composed of a Precambrian metamorphic basement and metamorphic or unmetamorphosed Paleozoic sedimentary rocks, including Cambrian–Ordovician marine carbonates and Carboniferous–Permian clastic sediments and limestones. The cover consists of Jurassic to Cretaceous volcanic rocks and terrestrial sediments, as well as lacustrine and fluvial sediments in the Cenozoic (Figure2). The thick Paleogene strata are mainly composed of thick Shahejie Formation (Es) sediments, which are further divided into four members. Es4 is composed of a 200–1000 m thick red bed, mudstones with interbedded gypsum, and salt rocks. Es3 consists of 300 m to 1200 m gray-black sandstones and conglomerates. Es2 includes 100–600 m variegated mudstones and inter-layered siltstones. Es1 contains 50 m to 400 m sandstones, mudstones, shales, and biogenic limestones [3,45,46]. Among them, the lower sub-member of Es3, the target formation in our study, is considered to have the largest high-quality reservoirs in this region, making it of great exploration value. Water 2021, 13, 408 4 of 18

Figure 2. Stratigraphic sequence of the Chezhen depression in the Cenozoic. Modiﬁed from Zhang et al. (2019) [2].

3. Methods Well logging is one of the most commonly used geophysical wireline techniques applied to hydrocarbon exploration and reservoir description [47]. It is widely used for applications including lithotype discrimination [48], estimation of shale content [49,50], lithological correlations from well to well in sediments [51], interpretation of the paleo- Water 2021, 13, 408 5 of 18

climate and/or paleo-environment [52–54], sedimentary facies characterization, and reservoir analysis [55]. In this study, the well logging technique was applied as an indicator of paleocurrents and sedimentary provenance directions because the variations of index components from sediments are sensitive to erosion during migration. Where lithological variations and deformations are gentle, the changes in the distribution and maturity of sedimentary deposits and abundance of particular active materials can thus be deemed a corresponding kinestate of the flowage of paleocurrents. Gamma-ray (GR) logging data are measured from the “total” gamma activity surrounding a borehole. This activity is a sum of radioactive contributions from the decay series of U238 (mainly from the intermediate daughter Bi214) and Th232 (mainly from Tl208), and from K40 [56]. Th mainly exists in heavy/resistate minerals (e.g., thorite, monazite, and rutile) [57] and is generally presumed to be contained by clay minerals in clastic sediments [58,59]. Because Th is mainly transported in a state of solid particles [58,59], the concentration of Th is decreased by hydraulic sorting and recycling of sedimentary environments along the paleocurrents [60]. U is present in uraninite (only in the case of sediments deposited before oxy-atmospheric conversion), allanite, zircon, and monazite in sedimentary rocks; in various rock-forming minerals; and in association with various colloidal and carbonaceous substances in sediments [57]. In the migration of sediments, U may be extracted by chemical weathering, recycling [61], and the pH and Eh of ancient fluvial systems and transport to the reservoirs along paleocurrents [57]. K is contained mainly in K-feldspars and clay minerals. The survival of K-feldspars in a sedimentary environment is a function of the relief and climatic factors that prevailed in the provenance [62]. The process of sedimentary transport from high to low positions decreases the concentration of K in sediments. Because the K-bearing minerals are less stable compared to Th-bearing minerals, K decays faster than Th in a natural sedimentary environment [63]. Additionally, the effective elimination of K-feldspars renders sediments more quartz-rich and mature. Therefore, along the direction of paleocurrents, the Th/U ratio decreases and the Th/K ratio increases. The photoelectric absorption (Pe) index of the sedimentary rocks can be measured by lithology-density logging. The Pe index can effectively distinguish the typical lithology, and it suffers little influence from the properties and other contents of fluids in the source rocks [64,65]. It is especially sensitive to the total organic carbon content in mud- stone [66]. With paleocurrent transport, the argillaceous components decrease and the Pe index decreases. Therefore, the Pe index decreases following the fluxion of paleocurrents. Through dip logging, the tendency and dip of the sedimentary layers around the wells can be measured. The variations of tendency and dip corrected by tectonic dip angle are classified into five patterns [30,67]. Among them, the blue pattern presents a similar tendency and decreasing of dip angle related to the increase of depth, suggesting a result of the flow of paleocurrents. The paleocurrent direction is determined by the main tendency. Others are named as red, green, yellow, and white patterns, presenting faults, stable sedimentation, overturned strata, and meaninglessness, respectively. The permeability, porosity, and thickness of reservoirs were derived from seven wells L for the sedimentary layers of the Es3 Formation. The well testing process was performed to determine the variations of porosity and permeability with depth.

4. Results L According to the W–E seismic interpretation profile, the deformation of Es3 For- mation tends to be gentle in this region, and the fault scarp with a large dip only exists near the eastern boundary (Figure3a). Combining this with the stratigraphic correlation L profile derived from the logging data, the Es3 Formation can be divided into six sandsets (Figure3b) . Because Sandset 6 was only identified in the easternmost well drilling and there were no comparable ones in other wells, we excluded it and performed the analyses focused on the data from Sandsets 1–5 in this region. Water 2021, 13, 408 6 of 18

Figure 3. (a) W–E trending seismic profile interpretations of the northern Chezhen depression. (b) W–E trending well logging profile. L The variations of the Th/U, Th/K, and Pe values of Sandsets 1–5 in the Es3 Formation obtained by GR logging and lithology-density logging from seven wells in this region are shown in Figure4 and Table1. There is a strong similarity in the statistics for the different sandsets from each well. The standard deviations of data from each well are relatively small, and the scope of AVG ± SD covers most of the data (Figure4). These results indicate L the stability of sedimentary layers in the Es3 Formation. Thus, we use the averaged Th/U, Th/K, and Pe values derived from five sandsets to get further results.

Figure 4. The variations of Th/U, Th/K, and Pe in each well, corresponding to Table1. Green shades are intervals of ± standard deviation (SD). Water 2021, 13, 408 7 of 18

L Table 1. Th/U, Th/K, and Pe derived from seven wells from the Es3 Formation in the Chezhen depression.

Type/ Wells Sandset No. C01 C8 C61 C0 C6 C10 C2 Th/U S1 4.77 3.25 3.06 4.18 2.16 5.05 2.38 S2 4.22 3.11 2.76 4.34 1.78 4.93 2.31 S3 4.79 3.15 2.94 4.64 1.79 4.40 2.64 S4 4.32 3.94 2.97 4.43 1.98 5.25 2.33 S5 4.42 3.05 2.88 4.62 2.04 4.38 2.26 AVG a 4.50 3.30 2.92 4.44 1.95 4.80 2.38 SD 0.23 0.33 0.10 0.17 0.15 0.35 0.13 Th/K S1 4.98 4.79 5.86 4.99 6.38 4.83 4.60 S2 5.40 5.55 5.28 5.26 5.94 4.88 5.00 S3 5.59 5.59 6.05 5.59 6.36 4.53 4.92 S4 5.42 5.28 5.66 5.30 6.82 4.74 4.58 S5 5.11 5.41 5.85 5.68 6.46 4.23 4.51 AVG b 5.30 5.32 5.74 5.36 6.39 4.64 4.72 SD 0.22 0.29 0.26 0.25 0.28 0.24 0.20 Pe S1 5.62 7.83 4.41 7.00 3.72 7.68 5.66 S2 5.75 8.38 4.62 7.13 4.21 7.53 5.57 S3 5.85 7.31 4.45 7.68 4.37 8.04 4.84 S4 5.72 7.75 4.38 7.49 3.38 7.69 5.22 S5 6.11 9.23 4.71 7.37 3.90 8.69 5.18 AVG c 5.81 8.10 4.51 7.33 3.92 7.93 5.29 SD 0.17 0.66 0.13 0.24 0.35 0.42 0.29 Note: S1–5 = Sandsets 1–5. AVG = Averaged data. SD = Standard deviation. a–c = used in Figures4 and5.

The variation of the averaged Th/U ratio of the five sandsets indicates a decline from Wells C0 and C10 in the north to C6 in the south, and from Well C8 in the west to C2 in the east (Figure5a). The averaged Th/K ratio indicates an increase from Wells C10 and C0 in the northwest to Well C6 in the south, from Well C8 in the west to Well C6, and from Well C2 in the northeast to Well C6 (Figure5b). The variation of the averaged Pe index is similar to that of Th/U, which indicates a decline from Wells C10 and C0 in the northeast to Well C6 in the south, from Well C8 in the west to Well C6, and from Well C2 in the northeast to Well C6 (Figure5c). The dip logging data were mainly obtained from the cross beddings, and the patterns of the stratigraphic dip variations were identified for each sandset from Wells C10, C01, and C61 (Table2). The dip angles are generally slight in this region, and the main tendencies are concentrated. In Wells C10, C01, and C61, the main tendencies are directed southeast, southeast/southwest, and southeast/southwest, respectively. Although secondary green and red patterns were recognized in some sandsets, the dominating dip pattern in these wells is the blue one, suggesting a strong correlation of the scouring of paleocurrents in this region. The porosity and reservoir thickness for Sandsets 1–5 were obtained from each well (Table3). Based on these data, we made the averaged reservoir thickness and averaged porosity variation maps for this region (Figure6). The averaged reservoir thickness map and revised well logging profile show a thickening of sedimentary deposits from Wells C8 and C6 in the southwest to Well C2 in the northeast (Figure6a). The averaged porosity map and the revised well logging profile show that the porosity for the sandstone of the L Es3 Formation is lower in Wells C8, C10, and C2 at the edges and gradually increases to Wells C6 and C61 in the south-central area (Figure6b). Water 2021, 13, 408 8 of 18

Figure 5. Schematic maps for the directions of paleocurrents based on (a) Th/U, (b) Th/K, and (c) Pe. Blue bars indicate the data for each well, corresponding to Table1; however, 1, 4, and 3 were subtracted from averaged Th/U, Th/K, and Pe data, respectively, to clarify the differences. Black lines are paleoelevations. Red arrows indicate the directions of paleocurrents.

L Table 2. Dip logging data derived from seven wells from the Es3 Formation in the Chezhen depression.

Well/Sandset No. Color Patterns Dip Angle (◦) Azimuth Scope (◦) Main Azimuth (◦) C10 S1 Blue 17–27 91.1–224.1 126.3 S2 Blue, green 17–29 88.2–211.6 138.8 S3 Blue, green 14–26 101.1–243.2 146.4 S4 Blue 15–26 108.2–179.5 155.3 S5 Blue, red 10–27 114.7–273.6 160.2 C01 S1 Blue, green 20–32 98.4–241.4 103.9/240.2 S2 Blue 15–28 95.3–259.1 108.1/235.5 S3 Blue, green 17–31 103.7–246.8 110.7/233.9 S4 Blue, green 15–29 121.3–179.5 98.6/179.5 S5 Blue 16–29 115.7–252.6 116.3/236.5 C61 S1 Blue, green 11–19 68.7–187.6 93.6/208.8 S2 Blue 11–20 82.8–213.2 101.2/206.3 S3 Blue, green 9–21 90.1–196.2 118.0/202.6 S4 Blue 8–18 123.8–179.5 126.7/179.5 S5 Blue 8–19 102.7–210.5 123.5/203.4 Water 2021, 13, 408 9 of 18

L Table 3. The porosity, permeability, and reservoir thickness derived from seven wells from the Es3 Formation in the Chezhen depression.

Well/Sandset No. Porosity (%) Averaged Porosity (%) Permeability (mD) Reservoir Thickness (m) Total Thickness (m) C10 S1 3.41 0.35 30.4 S2 3.17 0.36 20.8 S3 2.17 0.12 39.0 S4 2.95 0.19 28.0 S5 2.51 2.84 a 0.12 12.0 130.2 g C01 S1 4.52 0.33 19.5 S2 5.74 0.58 19.0 S3 5.86 0.76 8.0 S4 7.14 1.65 35.5 S5 6.05 5.86 b 0.73 29.5 111.5 h C2 S1 2.60 0.10 27.4 S2 2.44 0.10 12.4 S3 1.46 0.10 10.2 S4 3.90 0.13 64.5 S5 3.15 0.10 42.6 S6 3.20 2.80 c 0.10 53.2 210.3 i C6 S1 - - 0.0 S2 5.38 0.40 6.0 S3 2.35 0.10 7.2 S4 - - 0.0 S5 4.70 4.14 d 0.77 4.5 17.7 j C61 S1 4.83 0.20 1.0 S2 4.85 0.30 26.8 S3 8.01 2.37 18.8 S4 6.40 1.20 28.1 S5 8.34 6.49 e 2.25 29.5 104.3 k C0 S1 3.52 0.10 5.7 S2 3.86 0.26 36.6 S3 4.67 1.66 37.8 S4 5.29 1.95 50.9 S5 4.97 4.46 f 2.35 18.3 149.3 l Note: “S1–5” = Sandset 1–5. “-” = no data. a–f = data used in Figure6b, g–l = data used in Figure6a.

Figure 6. (a) Averaged reservoir thickness contour map. (b) Averaged porosity contour map. The data correspond to Table3.

The well testing data from a depth interval of 3356.2 m to 4552 m show that the porosity ranges from 0.1% to 10.6% and the permeability from 0.1 to 5.535 milli-darcy (abbreviated to mD) (Figure7; Table4), suggesting generally poor physical properties L of the layers from the Es3 Formation in this region. However, two positive anomalies were found from 3700–3800 m and 4100–4300 m intervals (Figure7), suggesting improved physical properties. Water 2021, 13, 408 10 of 18

Figure 7. The variations of porosity (a) and permeability (b) with depth, corresponding to Table4. Thick black lines indicate the ﬁtted variational curves. Red shades indicate the abnormal overpressure zones. WAVG = Weighted-average data.

L Table 4. Well testing data of the Es3 Formation.

Porosity (%) Permeability (mD) Thickness/Plies Ply Thickness (m) Deep (m) Min Max WAVG Min Max WAVG 36.3/1 36.3 3356.2 1.44 5.23 3.88 0.10 0.10 0.10 2.0/1 2.0 3425.4 1.55 6.79 5.12 0.10 0.39 0.22 161.0/1 161.0 4462.0 0.42 10.10 6.09 0.10 3.96 0.24 14.8/3 1.9 4233.0 5.46 10.90 8.91 0.27 5.54 2.93 7.6 4243.0 3.98 7.94 5.56 0.10 1.38 0.45 5.3 4253.8 2.92 4.46 4.12 0.10 0.35 0.18 4.0/1 4.0 3759.6 2.53 10.62 7.58 0.10 4.94 2.13 12.6/5 2.0 3883.7 2.05 6.95 4.06 0.10 0.76 0.22 3.0 4066.5 2.03 3.87 2.16 0.10 0.10 0.10 2.5 4077.3 0.46 5.04 1.57 0.10 0.19 0.13 2.4 4091.6 1.75 6.20 2.30 0.10 0.46 0.23 2.7 4117.7 3.98 6.60 5.23 0.10 0.61 0.27 23.8/6 6.0 4234.5 3.45 3.63 3.50 0.14 0.36 0.29 1.7 4239.7 0.99 5.45 3.31 0.10 1.05 0.22 2.4 4243.4 1.59 2.78 2.54 0.10 0.10 0.10 10.1 4254.9 0.42 4.92 1.80 0.10 0.67 0.13 1.9 4261.7 0.75 5.15 2.40 0.10 0.82 0.22 1.7 4264.5 1.58 4.58 3.15 0.10 0.49 0.16 16.0/7 3.0 4527.0 0.10 3.81 2.88 0.10 0.10 0.10 6.0 4536.0 2.14 4.10 3.14 0.10 0.10 0.10 1.0 4543.0 0.11 2.12 1.98 0.10 0.10 0.10 1.0 4546.0 0.80 2.66 2.10 0.10 0.10 0.10 2.0 4549.0 1.13 2.99 1.89 0.10 0.10 0.10 2.0 4552.0 1.19 4.02 3.08 0.10 0.10 0.10 1.0 4557.0 0.10 0.75 0.35 0.10 0.10 0.10 Note: WAVG = Weighted-average data. Water 2021, 13, 408 11 of 18

5. Discussion 5.1. Analysis of Provenance Directions The concentrations of these heavy metals are determined by the abundance of the minerals that host them. Besides the paleocurrents, the diagenesis may also influence them. In the Chezhen depression, the sedimentary deposits experienced gentle compaction, cementation, and dissolution in the burial diagenetic evolution [68]. However, the study region was located in the center of the Chezhen depression and was shown to be a single and stable lacustrine environment [3], suggesting that the influence of diagenesis on the variations of indicators was homogeneous in different wells. The original distribution of heavy minerals was basically finished during sedimentation and syndiagenesis, and was then changed by the sustaining paleocurrent flows. Some regional studies suggested L that the sedimentary layers of the Es3 Formation were mainly in the middle diagenetic stage A [68–70]. The alteration to the mother rocks was of not very significant. Analyses of lab data and rock slices suggested that the primary pores were effectively preserved, and the secondary pores affected amounted to 3–8% [69], which did not critically affect the sedimentary deposits. The paleogeotemperature was estimated to be less than 140 ◦C. It was at a low thermal evolution stage of the clay minerals and the tectonic zones [71]. The cementation and dissolution were not sufficient to result in an extensive migration L of heavy minerals. The minerals in the deposits of the Es3 Formation in this region were mainly disordered [71]. Additionally, the assemblages of heavy minerals were valid in the provenance analysis even in regions with stronger diagenesis [72]. The tectonic L activities were slight and infrequent after the age of the Es3 Formation, and the depocenters L almost unchanged [1]. The sedimentary layers of the Es3 Formation have been kept in good condition since its formation, also suggesting a lower effect of epidiagenesis and other exogenous processes. Therefore, we suggest that the influence of diagenesis on the variations of indicators can be neglected in this study, and the heavy minerals were still effective in the provenance analysis. Because the particles of heavy minerals are small, they are often wrapped in muds, and the concentration is enriched in different areas during the transportation of paleocurrents [58,59]. The heavy mineral analysis is also a well-established method used for the determination of palaeocurrent directions and provenance analysis [21,22]. Different regions of provenance have specific characteristics of parent rocks, which are manifested as different assemblages of heavy minerals in the sedimentary deposits transported by paleocurrents [21,22]. The gamma logging technology measures the value of the total concentration of elements from all minerals around the underground wellbores, which is a better indicator than the ones derived from individual samples [27,28]. The effect of dissolution on the results is also greatly limited. Therefore, we suggest that results obtained by GR logging and lithology-density logging can effectively represent the variation of ionic concentration after the redistribution of paleocurrents in the study area. The results of GR logging and lithology-density logging from seven wells in this area were plotted in Figure5. The directions of paleocurrent sources were recognized as northwest and west; northeast, northwest and west; and northeast, northwest and west, L based on the variations of the Th/U, Th/K, and Pe values of Sandsets 1–5 in the Es3 Formation, respectively. The results of dip logging are in good agreement with those of GR logging and lithology-density logging. Because the sedimentation and regional deformation were gentle after the formation of the Chezhen depression [1,3] and the blue pattern of dip variation was predominant in this region, the paleocurrent directions can be easily recognized by the main tendencies of the strata. Therefore, the paleocurrents from the northeast, northwest and west were identified from the results in Wells C10, C01, and C61 (Figure8). Water 2021, 13, 408 12 of 18

Figure 8. (a) Cumulative seismic amplitude map. Blue to red indicates the velocity of seismic waves from fast to slow. Red arrows indicate the directions of paleocurrents. A and B indicate the predicted high-quality liquid reservoirs (presented as buckets). (b–d) The paleocurrents implied by the results of stratigraphic dip logging of Wells C10, C01, and C61 are presented in rose charts.

The variations of reservoir thickness and porosity are reliable auxiliary indexes in a region with a stable sedimentary environment and strong paleocurrents. The reservoir thickness has a high correlation with paleogeomorphology and depositional sequences. From the statistical results, the variation of reservoir thickness proves the existence of a source of sediments and paleocurrents in the further northeast (Figure6a). As the sediments were transported by paleocurrents, the muds and clays decreased. The porosity of the sedimentary deposits should increase from the provenance to the depocenter. The variation of porosity in this region agrees well with the inferences from well loggings, suggesting the provenance of paleocurrents from three directions (northeast, northwest and west) (Figure6b). Two positive abnormal zones of porosity and permeability were found at the depth intervals of 3700–3800 m and 4100–4300 m (Figure7 and Table4). Because the controlling action of overpressure on the reservoirs is stronger than in other reasons for the Paleogene strata in the Chezhen Depression [70], the anomalies of porosity and permeability may be caused by the overpressure of liquid reservoirs and the gas produced by organic materials. Combining these results and the provenance analysis, we suggest that high-quality liquid reservoirs should be located at the intersection of multiple paleocurrents, namely, between Wells C8 and C61 and between Wells C2 and C61, with depths of 3700–3800 m and 4100–4300 m (A and B in Figure8). Based on the logging data from all 40 wells (33 of them are hidden following privacy policies) in the Chezhen depression, we restored the paleogeomorphology associated with L the deposition of the Es3 Formation in this region (Figure9). The paleogeomorphologic restoration map shows that the central to northeastern part of the region was the subsidence center of the Chezhen depression with lower altitude, collecting the sediments from the terrains with higher altitude in the northeast, northwest, and west. Although there is an Water 2021, 13, 408 13 of 18

area of middle-high paleoelevation in the southern part of this region, we consider it may be the stratigraphic overlapping of sediments deposited along the southern margin of the half-graben based on our results regarding paleocurrents (Figure9b).

L Figure 9. (a) Reconstruction of the paleogeomorphology for the Es3 Formation in the Chezhen area. (b) Three-dimensional L reconstruction of the paleogeomorphology for the Es3 Formation in the Chezhen area. Blue to red indicates the paleoelevation from low to high. Red arrows indicate the directions of paleocurrents.

5.2. Interpretation of the Tectonic Dynamics The macroscopical paleocurrents and sedimentary transport in basins are often controlled by the surrounding tectonic activities. In a faulted basin, the macroscopic directions of paleocurrents and provenance usually indicate the sources of tectonic stress in the region [61,73,74]. As a typical sub-basin of the Bohai Bay basin, the Chezhen depression, belonging to the Bohai–Luxi block, is adjacent to the N–S-trending Tan–Lu strike-slip fault, the Izanagi Plate, and the Pacific Plate in the east; is adjacent to the NE–SW-trending Lan–Liao strike-slip fault in the west; and it is remotely affected by the compression-stretch of the Yinshan–Yanshan fold belt and the collision of the North China block and South China block in the N–S direction (Figure 10)[41–43,75]. Therefore, the multidirectional tectonic genesis makes it of great research value. There are a series of NW–SE-ranged Cenozoic secondary basin and uplift tectonic assemblages, including the Chezhen depression, in the Bohai Bay basin, in the east of the North China block, which are referred to as the “wave basining” phenomenon [10,76]. Dur- ing the Cenozoic, the mantle thermal equilibrium in the Bohai Bay area was destroyed due to the subduction of the Pacific Plate from the Eurasian Plate and the collision of the Indian Plate with the South China block. After the dissolution of the old, cold lithosphere (earlier than ca. 450 Ma) in the middle and late Mesozoic, the thermal convection in the mantle of this area was disturbed and resulted in regional uplifts, while the continental crusts above the new underplating mantle stretched and resulted in rifting depressions [77,78]. On the other hand, the plate tectonic activities and interactions were still ongoing at the same time. The Pacific Plate began subducting to Eurasia before the Jurassic [42,79–81]. The velocity of subduction reached its peak at ~10 cm/y at the end of the Cretaceous and then declined rapidly after the Cenozoic (ca. 50 Ma) [42,82,83]. The rapid attenuation of the subduction velocity led to the extension of the crusts on both sides of the subduction zone, resulting in a large number of collapse and pull-apart depressions [42,84,85]. The distribution of volcanic rocks decreased and migrated from the western inland to the eastern coastal area during the Paleogene, which also suggests the retracement of the Pacific Plate [86–88]. Despite the subduction zone of the Pacific Plate retreating on the eastern margin of the North China block, the mid-ocean ridge was still expanding [82,83]. The alternations of subductions and rollbacks of the Pacific Plate in the eastern margin of the North China block and the re-activation of the ancient mantle under the North China craton might have Water 2021, 13, 408 14 of 18

led to the “wave basining” phenomenon in the Bohai Bay basin. The analyses of provenance and paleocurrents indicated that the Chezhen depression received a large amount L of sediments transported by paleocurrents from northwest in the age of Es3 , suggesting that a southeastern retracement of the Paciﬁc Plate was happening in the early Paleogene (ca. 50.5–41.7 Ma) (Figure 10). The observation of absolute moving velocities of the Paciﬁc Plate relative to the lower mantle suggested its dynamic changed from nearly NW–SE to W–E [82,89], which resulted in the western provenance for the Chezhen depression.

Figure 10. A simplified tectonic map of eastern China. Modified from Li et al. (2010) [89]. BL = Bohai–Luxi block; ES = Erlian–Songliao block; JLJ = Jiaobei–Liaodong–Jinan block; JO = Jinji–Ordos block; muY = middle–upper Yangtze block; YH = lower Yangtze–southern Huanghai block; YYOB = Yinshan–Yanshan orogenic belt; QDOB = Qinling–Dabie orogenic belt; F1 = Tan–Lu strike-slip fault; F2 = Lan–Liao strike-slip fault. Red arrows indicate the directions of tectonic stress.

In the eastern margin of North China, N–NE-trending Mesozoic and Cenozoic strike- slip structures are widely developed in the Bohai Bay basin, including the Tan–Lu fault, Lan–Liao fault, Dagang fault, and Shulu fault [42,76,77]. These large-scale NNE strike-slip fault systems played an important role in the tectonic formation of the Bohai Bay basin. These strike-slip faults have very obvious piecewise evolutionary characteristics: the initial formation time of these fault groups is vague, possibly during the early Indochina period; then they were basically formed in the early Indochina period. In the late Yanshan period, most of them were left-lateral strike-slip faults; in the Cenozoic, they mostly manifested as right-lateral extensional strike-slip faults [42,76,77,89]. Among them, the Lan–Liao Water 2021, 13, 408 15 of 18

fault showed stable right-lateral strike-slipping during the entire Cenozoic, while the Tan–Lu fault showed left-lateral strike-slipping during the Paleocene to early Eocene and changed to right-lateral until the middle Eocene [31,77,89]. During the Cenozoic, the subduction front of the Paciﬁc Plate reached the Bohai Bay and even further north [86,89], which intensiﬁed the strike-slip movement of the fault zones. The provenance from northeast may suggest that the Chezhen depression (Bohai–Luxi block) was subjected to the southwestward escape tectonic movement induced by the Lan–Liao and Tan–Lu strike-slip faults in the early Paleogene (Figure 10).

6. Conclusions In this study, we obtained the Th/U and Th/K ratios, Pe index, porosity, permeability, stratigraphic dips, and the thickness of reservoirs by using well logging technology on the L sedimentary layers of the Es3 Formation in the Chezhen depression; and we reconstructed the paleogeomorphologic maps. The combined results indicate the existence of paleocurrents from the NE, W, and NW and convergence in the northern Chezhen depression at that time. Based on the well testing data, the physical properties of the layers from the L Es3 Formation in this region were generally poor; however, two abnormal overpressure zones were found at 3700–3800 m and 4100–4300 m deep intervals, suggesting potential high-quality underground liquid reservoirs. By combining with other geological evidence, we suggest that the Paciﬁc Plate was retreating and changing its direction from NE–SW to W–E and the Bohai–Luxi block was suffering an extrusion from NE induced by the Lan–Liao and Tan–Lu strike-slip faults in the early Paleogene.

Author Contributions: Conceptualization, Y.G. and S.S.; writing—original draft, Y.G.; software, supervision, F.L.; writing—review & editing, Y.C. All authors have read and agreed to the published version of the manuscript. Funding: This research was supported by the National Natural Science Foundation of China (8519768). Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: All data necessary to carry out the work in this paper are included in the figures, tables or are available in the cited references. Acknowledgments: We thank Lanshuang Gao and Shuling Yang for their financial assistance; Changtian Guan and Min Mu for the collection of partial data; and the members of Zhan–Che institute for their assistance in the fieldwork. Conflicts of Interest: The authors declare no conflict of interest.

References 1. Su, J.; Zhu, W.; Wei, J.; Xu, L.; Yang, Y.; Wang, Z.; Zhang, Z. Fault Growth and Linkage: Implications for Tectonosedimentary Evolution in the Chezhen Basin of Bohai Bay, Eastern China. AAPG Bull. 2011, 95, 1–26. [CrossRef] 2. Zhang, J.; He, S.; Wang, Y.; Wang, Y.; Hao, X.; Luo, S.; Li, P.; Dang, X.; Yang, R. Main Mechanism for Generating Overpressure in the Paleogene Source Rock Series of the Chezhen Depression, Bohai Bay Basin, China. J. Earth Sci. 2019, 30, 775–787. [CrossRef] 3. Lao, H.; Wang, Y.; Shan, Y.; Hao, X.; Li, Q. Hydrocarbon Downward Accumulation from an Upper Oil Source to the Oil Reservoir Below in an Extensional Basin: A Case Study of Chezhen Depression in the Bohai Bay Basin. Mar. Pet. Geol. 2019, 103, 516–525. [CrossRef] 4. Hu, S.; O’Sullivan, P.B.; Raza, A.; Kohn, B.P. Thermal History and Tectonic Subsidence of the Bohai Basin, Northern China: A Cenozoic Rifted and Local Pull-Apart Basin. Phys. Earth Planet. Inter. 2001, 126, 221–235. [CrossRef] 5. Yang, Y.; Xu, T. Hydrocarbon Habitat of the Offshore Bohai Basin, China. Mar. Pet. Geol. 2004, 21, 691–708. [CrossRef] 6. Hao, F.; Zhou, X.; Zhu, Y.; Bao, X.; Yang, Y. Charging of the Neogene Penglai 19-3 Field, Bohai Bay Basin, China: Oil Accumulation in a Young Trap in an Active Fault Zone. AAPG Bull. 2009, 93, 155–179. [CrossRef] 7. Zhu, H.; Yang, X.; Liu, K.; Zhou, X. Seismic-Based Sediment Provenance Analysis in Continental Lacustrine Rift Basins: An Example from the Bohai Bay Basin, Chinaseismic-Based Sediment Provenance Analysis. AAPG Bull. 2014, 98, 1995–2018. [CrossRef] 8. Sun, Z.; Zhu, H.; Xu, C.; Yang, X.; Du, X. Reconstructing Provenance Interaction of Multiple Sediment Sources in Continental Down- Warped Lacustrine Basins: An Example from the Bodong Area, Bohai Bay Basin, China. Mar. Pet. Geol. 2020, 113, 104142. [CrossRef] Water 2021, 13, 408 16 of 18

9. Li, S.; Zhu, H.; Xu, C.; Zeng, H.; Liu, Q.; Yang, X. Seismic-Based Identification and Stage Analysis of Overlapped Compound Sedimentary Units in Rifted Lacustrine Basins: An Example from the Bozhong Sag, Bohai Bay Basin, China. AAPG Bull. 2019, 103, 2521–2543. [CrossRef] 10. Zhao, L.; Zheng, T. Seismic Structure of the Bohai Bay Basin, Northern China: Implications for Basin Evolution. Earth Planet. Sci. Lett. 2005, 231, 9–22. [CrossRef] 11. Liu, S.; Lin, G.; Liu, Y.; Zhou, Y.; Gong, F.; Yan, Y. Geochemistry of Middle Oligocene-Pliocene Sandstones from the Nanpu Sag, Bohai Bay Basin (Eastern China): Implications for Provenance, Weathering, and Tectonic Setting. Geochem. J. 2007, 41, 359–378. [CrossRef] 12. Pang, X.; Wang, Q.; Du, X.; Dai, L.; Li, H.; Wang, M. Matter Provenance Evolution and Its Influences on Paleogene Reservoirs in the Northwestern Margin of Bozhong Sag. Pet. Geol. Oilfield Dev. Daqing 2016, 35, 34–41. 13. Du, X.; Wang, Q.; Zhao, M.; Liu, X. Sedimentary Characteristics and Mechanism Analysis of a Large-Scale Fan Delta System in the Paleocene Kongdian Formation, Southwestern Bohai Sea, China. Interpretation 2020, 8, SF81–SF94. [CrossRef] 14. Li, H.; Du, X.; Wang, Q.; Yang, X.; Zhu, H.; Wang, F. Formation of Abnormally High Porosity/Permeability in Deltaic Sandstones (Oligocene), Bozhong Depression, Offshore Bohai Bay Basin, China. Mar. Pet. Geol. 2020, 121, 104616. [CrossRef] 15. Miao, X.; Lu, H.; Li, Z.; Cao, G. Paleocurrent and Fabric Analyses of the Imbricated Fluvial Gravel Deposits in Huangshui Valley, the Northeastern Tibetan Plateau, China. Geomorphology 2008.[CrossRef] 16. Amajor, L.C. Paleocurrent, Petrography and Provenance Analyses of the Ajali Sandstone (Upper Cretaceous), Southeastern Benue Trough, Nigeria. Sediment. Geol. 1987, 54, 47–60. [CrossRef] 17. Das, R.; Pandya, K.L. Palaeocurrent Pattern and Provenance of a Part of Gondwana Succession, Talchir Basin, Orissa. J. Geol. Soc. India 1997, 50, 425–433. 18. Condie, C.K. Chemical Composition and Evolution of the Upper Continental Crust: Contrasting Results from Surface Samples and Shales. Chem. Geol. 1993, 104, 1–37. [CrossRef] 19. Naqvi, S.M.; Sawkar, R.H.; Rao, D.S.; Govil, P.K.; Rao, T.G. Geology, Geochemistry and Tectonic Setting of Archaean Greywackes from Karnataka Nucleus, India. Precambrian Res. 1988, 39, 193–216. [CrossRef] 20. Raza, M.; Casshyap, S.M.; Khan, A. Geochemistry of Mesoproterozoic Lower Vindhyan Shales from Chittaurgarh, Southeastern Rajasthan and Its Bearing on Source Rock Composition, Palaeoweathering Conditions and Tectonosedimentary Environments. J. Geol. Soc. India 2002, 60, 505–518. 21. Morton, C.A. Geochemical Studies of Detrital Heavy Minerals and Their Application to Provenance Research. Geol. Soc. Lond. Spec. Publ. 1991, 57, 31–45. [CrossRef] 22. Morton, A.C.; Hallsworth, C.R. Processes Controlling the Composition of Heavy Mineral Assemblages in Sandstones. Sediment. Geol. 1999, 124, 3–29. [CrossRef] 23. Morton, A.C.; Whitham, A.G.; Fanning, C.M. Provenance of Late Cretaceous to Paleocene Submarine Fan Sandstones in the Norwegian Sea: Integration of Heavy Mineral, Mineral Chemical and Zircon Age Data. Sediment. Geol. 2005, 182, 3–28. [CrossRef] 24. Morton, A.; Hurst, A. Correlation of Sandstones Using Heavy Minerals: An Example from the Statfjord Formation of the Snorre Field, Northern North Sea. Geol. Soc. Lond. Spec. Publ. 1995, 89, 3–22. [CrossRef] 25. Verma, S.P.; Armstrong-Altrin, J.S. New Multi-Dimensional Diagrams for Tectonic Discrimination of Siliciclastic Sediments and Their Application to Precambrian Basins. Chem. Geol. 2013, 355, 117–133. [CrossRef] 26. Verma, S.P.; Armstrong-Altrin, J.S. Geochemical Discrimination of Siliciclastic Sediments from Active and Passive Margin Settings. Sediment. Geol. 2016, 332, 1–12. [CrossRef] 27. Miroslav, K.; Stanislav, M.; Paillet, F. Geophysical Well Logging. In Hydrogeophysics; Springer: Berlin/Heidelberg, Germany, 2005. 28. Working Group on Nuclear Techniques in, Hydrology. Nuclear Well Logging in Hydrology; International Atomic Energy Agency: Vienna, Austria, 1971. 29. Hu, C.; Hu, G.; He, Y.; Li, J.; Li, H.; Xu, M. Application of Radioactive Thorium Logging in Provenance Analysis. Geol. Sci. Technol. Inf. 2019, 38, 54–63. 30. Zhao, R.; Chen, S.; Wang, H.; Yan, D.; Cao, H.; Gong, Y.; He, J.; Wu, Z. Paleogene Sedimentation Changes in Lenghu Area, Qaidam Basin in Response to the India–Eurasia Collision. Int. J. Earth Sci. 2019, 108, 27–48. [CrossRef] 31. Li, S.; Zhao, G.; Dai, L.; Zhou, L.; Liu, X.; Suo, Y.; Santosh, M. Cenozoic Faulting of the Bohai Bay Basin and Its Bearing on the Destruction of the Eastern North China Craton. J. Asian Earth Sci. 2012, 47, 80–93. [CrossRef] 32. Tapponnier, P.; Peltzer, G.L.D.A.Y.; Le Dain, A.Y.; Armijo, R.; Cobbold, P. Propagating Extrusion Tectonics in Asia: New Insights from Simple Experiments with Plasticine. Geology 1982, 10, 611–616. [CrossRef] 33. Liu, Z.; Huang, C.; Algeo, T.J.; Liu, H.; Hao, Y.; Du, X.; Lu, Y.; Chen, P.; Guo, L.; Peng, L. High-Resolution Astrochrono- logical Record for the Paleocene-Oligocene (66–23 Ma) from the Rapidly Subsiding Bohai Bay Basin, Northeastern China. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2018, 510, 78–92. [CrossRef] 34. Liu, M.; Cui, X.; Liu, F. Cenozoic Rifting and Volcanism in Eastern China: A Mantle Dynamic Link to the Indo–Asian Collision? Tectonophysics 2004, 393, 29–42. [CrossRef] 35. Chen, S.; Liu, X.; Cui, Y.; Zhao, X.; Zhang, J. Palaeogene Structural Evolution of Dongying Depression, Bohai Bay Basin, Ne China. Int. Geol. Rev. 2017, 59, 259–273. [CrossRef] 36. Allen, M.B.; Macdonald, D.I.M.; Xun, Z.; Vincent, S.J.; Brouet-Menzies, C. Early Cenozoic Two-Phase Extension and Late Cenozoic Thermal Subsidence and Inversion of the Bohai Basin, Northern China. Mar. Pet. Geol. 1997, 14, 951–972. [CrossRef] Water 2021, 13, 408 17 of 18

37. Huang, H.; Pearson, M.J. Source Rock Palaeoenvironments and Controls on the Distribution of Dibenzothiophenes in Lacustrine Crude Oils, Bohai Bay Basin, Eastern China. Org. Geochem. 1999, 30, 1455–1470. [CrossRef] 38. Huang, L.; Liu, C.; Wang, Y.; Zhao, J.; Mountney, N.P. Neogene–Quaternary Postrift Tectonic Reactivation of the Bohai Bay Basin, Eastern China. AAPG Bull. 2014, 98, 1377–1400. [CrossRef] 39. Dong, Y.; Xiao, L.; Zhou, H.; Wang, C.; Zheng, J.; Zhang, N.; Xia, W.; Ma, Q.; Du, J.; Zhao, Z. The Tertiary Evolution of the Prolific Nanpu Sag of Bohai Bay Basin, China: Constraints from Volcanic Records and Tectono-Stratigraphic Sequences. Bulletin 2010, 122, 609–626. [CrossRef] 40. Huang, C.; Wang, H.; Wu, Y.; Wang, J.; Chen, S.; Ren, P.; Liao, Y.; Xia, C. Genetic Types and Sequence Stratigraphy Models of Palaeogene Slope Break Belts in Qikou Sag, Huanghua Depression, Bohai Bay Basin, Eastern China. Sediment. Geol. 2012, 261, 65–75. [CrossRef] 41. Li, S.; Li, S.; Wang, Y.; Li, X.; Suo, Y.; Wang, P.; Wang, Q.; Zhang, J.; Somerville, I. Cenozoic Faulting Response of Eastern North China to Subduction of the Pacific Plate: A Case of Study of the Luxi Block. Geol. J. 2017, 52, 70–80. [CrossRef] 42. Li, S.; Suo, Y.; Li, X.; Zhou, J.; Santosh, M.; Wang, P.; Wang, G.; Guo, L.; Yu, S.; Lan, H. Mesozoic Tectono-Magmatic Response in the East Asian Ocean-Continent Connection Zone to Subduction of the Paleo-Pacific Plate. Earth-Sci. Rev. 2019, 192, 91–137. [CrossRef] 43. Cao, X.; Flament, N.; Müller, D.; Li, S. The Dynamic Topography of Eastern China since the Latest Jurassic Period. Tectonics 2018, 37, 1274–1291. [CrossRef] 44. Min, H.; Yingchang, C.; Yanzhong, W. Sedimentary Facies of Sha4 to Lower Sha3 Member in Paleogene in Chezhen Depression. Pet. Geol. Recovery Effic. 2009, 16, 16–18. 45. Hongwen, D.; Yi, X.; Lixiang, M.; Zhenglong, J. Genetic Type, Distribution Patterns and Controlling Factors of Beach and Bars in the Second Member of the Shahejie Formation in the Dawangbei Sag, Bohai Bay, China. Geol. J. 2011, 46, 380–389. [CrossRef] 46. Su, S.; Jiang, Z.; Gao, Z.; Ning, C.; Wang, Z.; Li, Z.; Zhu, R. A New Method for Continental Shale Oil Enrichment Evaluation. Interpretation 2017, 5, T209–T217. [CrossRef] 47. Wahl, J.S. Gamma-Ray Logging. Geophysics 1983, 48, 1536–1550. [CrossRef] 48. Galbraith, J.H.; Saunders, D.F. Rock Classification by Characteristics of Aerial Gamma-Ray Measurements. J. Geochem. Explor. 1983, 18, 49–73. [CrossRef] 49. Poupon, A.; Gaymard, R. The Evaluation of Clay Content from Logs. In Proceedings of the SPWLA 11th Annual Logging Symposium, Los Angeles, CA, USA, 3–6 May 1970. 50. Wonik, T. Gamma-Ray Measurements in the Kirchrode I and Ii Boreholes. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2001, 174, 97–105. [CrossRef] 51. Thibal, J.; Etchecopar, A.; Pozzi, J.P.; Barthès, V.; Pocachard, J. Comparison of Magnetic and Gamma Ray Logging for Correlations in Chronology and Lithology: Example from the Aquitanian Basin (France). Geophys. J. Int. 1999, 137, 839–846. [CrossRef] 52. Herbert, D.T.; Mayer, L.A. Long Climatic Time Series from Sediment Physical Property Measurements. J. Sediment. Res. 1991, 61, 1089–1108. 53. Liu, Z.C.; Chen, Y.; Yuan, L.W.; Zhou, C.L.; Wang, Y.J.; Li, J.Q.; Yang, P.; Xi, P. The Application of Natural Gamma Logging Curve to Recovering Palaeoclimatic Changes from 2.85 Ma Bp. Sci. China (Ser. D) 2000, 30, 609–618. 54. Zhong-Hong, C.; Zha, M.; Jin, Q. Application of Natural Gamma Ray Logging and Natural Gamma Spectrometry Logging to Recovering Paleoenvironment of Sedimentary Basins. Chin. J. Geophys. 2004, 47, 1286–1290. 55. Song, Z.Q.; Li, W.-f.; Tang, C.-j.; Li, W.-f.; Pang, Z.-y.; Wang, Y. Dividing Sedimentary Facies and Reservoir Distributions by Using Natural Potential and Natural Gamma Ray Logging Curves. Prog. Geophys 2009, 24, 651–656. 56. Bücker, C.; Rybach, L. A Simple Method to Determine Heat Production from Gamma-Ray Logs. Mar. Pet. Geol. 1996, 13, 373–375. [CrossRef] 57. Boyle, R.W. Geochemical Prospecting for Thorium and Uranium Deposits; Elsevier: Amsterdam, The Netherlands, 2013. 58. Taylor, S.R.; McLennan, S.M. The Continental Crust: Its Composition and Evolution; U.S. Department of Energy: Washington, DC, USA, 1985. 59. Condie, C.K. Another Look at Rare Earth Elements in Shales. Geochim. Cosmochim. Acta 1991, 55, 2527–2531. [CrossRef] 60. Bauluz, B.; Mayayo, M.J.; Fernandez-Nieto, C.; Lopez, J.M.G. Geochemistry of Precambrian and Paleozoic Siliciclastic Rocks from the Iberian Range (Ne Spain): Implications for Source-Area Weathering, Sorting, Provenance, and Tectonic Setting. Chem. Geol. 2000, 168, 135–150. [CrossRef] 61. McLennan, S.M.; Hemming, S.; McDaniel, D.K.; Hanson, G.N. Geochemical Approaches to Sedimentation, Provenance, and Tectonics. Spec. Pap. Geol. Soc. Am. 1993, 21. [CrossRef] 62. Pettijohn, F.J. Sedimentary Rocks; Harper and Brothers: New York, NY, USA, 1957; 718p. 63. Basu, H.; Kumar, K.M.; Paneerselvam, S.; Chaki, A. Study of Provenance Characteristics and Depositional History on the Basis of U, Th and K Abundances in the Gulcheru Formation, Cuddapah Basin in Tummalapalle-Somalollapalle Areas, Cuddapah- Anantapur Districts, Andhra Pradesh. J. Geol. Soc. India 2009, 74, 318. [CrossRef] 64. Pan, H.P.; Niu, Y.X.; Wang, W.X. Ccsd Well Logging Engineering Program. J. China Univ. Geosci. 2002, 13, 91–94. 65. Niu, Y.X.; Pan, H.P.; Wang, W.X.; Zhu, L.F.; Xu, D.H. Geophysical Well Logging in Main Hole (0~2000 M) of Chinese Continental Scientific Drilling. Acta Petrol. Sin. 2004, 20, 165–178. 66. He, J.; Ding, W.; Zhang, J.; Li, A.; Zhao, W.; Dai, P. Logging Identification and Characteristic Analysis of Marine–Continental Transitional Organic-Rich Shale in the Carboniferous-Permian Strata, Bohai Bay Basin. Mar. Pet. Geol. 2016, 70, 273–293. [CrossRef] Water 2021, 13, 408 18 of 18

67. Wen, J.; He, Y.-b.; Wu, C.-x.; Luo, W. Geologic Application of Diplogging. J. Oil Gas. Technol. 2008, 1, 263–265. 68. Jin, J.; Guan, D. Anlyse the Characteristic of Diagenesis and Porosity in Dawangbei Sub-Sag, Chezhen Sag. Sichuan Nonferrous Met. 2016, 3, 14–17. 69. Wang, L. The Characteristics of Reservoirs and Hydrocarbon Accumulation in the Middle Section of the Steep Slope Zone in Chezhen Depression. Master’s Thesis, China University of Petroleum, Beijing, China, 2016. 70. Wang, Y.; Cao, Y. Lower Property Limit and Controls on Deep Effective Clastic Reservoirs of Paleogene in Chezhen Depression. Acta Sedimentol. Sin. 2010, 28, 752–761. 71. Ni, J.L.; Zhang, H.; Tang, X.L.; Han, S. Diagenetic and Paleogeothermal Evolution of the Clay Minerals in the Center Uplift Belt and Its Adjacent Region of Huimin Sag, Jiyang Depression. Nat. Gas Geosci. 2016, 27, 1778–1788. 72. Fu, H.; Hu, X.; Liang, W.; Eduardo, G. Heavy Mineral Assemblages of Jurassic-Paleogene Sandstones in Southern Tibet: Implications for Provenance Interpretations of Magmatic Arc and Continental Block. Geol. J. China Univ. 2020, 26, 530–539. 73. Mondal, M.E.A.; Wani, H.; Mondal, B. Geochemical Signature of Provenance, Tectonics and Chemical Weathering in the Quaternary Flood Plain Sediments of the Hindon River, Gangetic Plain, India. Tectonophysics 2012, 566, 87–94. [CrossRef] 74. Davis, W.D. U–Pb Geochronology of Archean Metasedimentary Rocks in the Pontiac and Abitibi Subprovinces, Quebec, Constraints on Timing, Provenance and Regional Tectonics. Precambrian Res. 2002, 115, 97–117. [CrossRef] 75. LiQing, X.; SanZhong, L.; LingLi, G.; YanHui, S.; XianZhi, C.; LiMing, D.; PengCheng, W.; GeGe, H. Impaction of the Tan-Lu Fault Zone on Uplift of the Luxi Rise: Constraints from Apatite Fission Track Thermochronology. Acta Petrol. Sin. 2016, 32, 1153–1170. 76. Li, S.Z.; Suo, Y.H.; Zhou, L.H.; Dai, L.M.; Zhou, J.T.; Zhao, F.M. Pull-Apart Basins within the North China Craton: Structural Pattern and Evolution of Huanghua Depression in Bohai Bay Basin. J. Jilin Univ. (Earth Sci. Ed.) 2011, 41, 1362–1379. 77. Li, S.Z.; Suo, Y.H.; Santosh, M.; Dai, L.M.; Liu, X.; Yu, S.; Zhao, S.J.; Jin, C. Mesozoic to Cenozoic Intracontinental Deformation and Dynamics of the North China Craton. Geol. J. 2013, 48, 543–560. [CrossRef] 78. Jiang, F.; Pang, X.; Yu, S.; Hu, T.; Bai, J.; Han, G.; Li, B. Charging History of Paleogene Deep Gas in the Qibei Sag, Bohai Bay Basin, China. Mar. Pet. Geol. 2015, 67, 617–634. [CrossRef] 79. Zhou, J.-B.; Li, L. The Mesozoic Accretionary Complex in Northeast China: Evidence for the Accretion History of Paleo-Pacific Subduction. J. Asian Earth Sci. 2017, 145, 91–100. [CrossRef] 80. Zhou, J.B.; Cao, J.L.; Wilde, S.A.; Zhao, G.C.; Zhang, J.J.; Wang, B. Paleo-Pacific Subduction-Accretion: Evidence from Geochemical and U-Pb Zircon Dating of the Nadanhada Accretionary Complex, Ne China. Tectonics 2014, 33, 2444–2466. [CrossRef] 81. Liu, L.; Xu, X.; Zou, H. Episodic Eruptions of the Late Mesozoic Volcanic Sequences in Southeastern Zhejiang, Se China: Petrogenesis and Implications for the Geodynamics of Paleo-Pacific Subduction. Lithos 2012, 154, 166–180. [CrossRef] 82. Müller, R.D.; Seton, M.; Zahirovic, S.; Williams, S.E.; Matthews, K.J.; Wright, N.M.; Shephard, G.E.; Maloney, K.T.; Barnett-Moore, N.; Hosseinpour, M. Ocean Basin Evolution and Global-Scale Plate Reorganization Events since Pangea Breakup. Annu. Rev. Earth Planet. Sci. 2016, 44, 107–138. [CrossRef] 83. Zahirovic, S.; Matthews, K.J.; Flament, N.; Müller, R.D.; Hill, K.C.; Seton, M.; Gurnis, M. Tectonic Evolution and Deep Mantle Structure of the Eastern Tethys since the Latest Jurassic. Earth-Sci. Rev. 2016, 162, 293–337. [CrossRef] 84. Liu, L.; Xu, X.; Xia, Y. Asynchronizing Paleo-Pacific Slab Rollback beneath Se China: Insights from the Episodic Late Mesozoic Volcanism. Gondwana Res. 2016, 37, 397–407. [CrossRef] 85. Qiu, L.; Yan, D.P.; Xu, H.; Shi, H.; Dong, W.; Sun, S. Late Cretaceous Mud Volcanism in the Southwestern Songliao Basin Records Slab Rollback of the Subducted Paleo-Pacific Plate Underneath Ne China. J. Asian Earth Sci. X 2020, 3, 100028. [CrossRef] 86. Tang, J.; Xu, W.; Wang, F.; Ge, W. Subduction History of the Paleo-Pacific Slab beneath Eurasian Continent: Mesozoic-Paleogene Magmatic Records in Northeast Asia. Sci. China Earth Sci. 2018, 61, 527–559. [CrossRef] 87. Tang, Z.Y.; Sun, D.Y.; Mao, A.Q.; Yang, D.G.; Deng, C.Z. Timing and Evolution of Mesozoic Volcanism in the Central Great Xing’an Range, Northeastern China. Geol. J. 2019, 54, 3737–3754. [CrossRef] 88. Tang, J.; Zhang, Z.; Ji, Z.; Chen, Y.; Li, K.; Jinfu, Y.; Liu, P. Early Cretaceous Volcanic and Sub-Volcanic Rocks in the Erlian Basin and Adjacent Areas, Northeast China: New Geochemistry, Geochronology and Zircon Hf Isotope Constraints on Petrogenesis and Tectonic Setting. Int. Geol. Rev. 2019, 61, 1479–1503. [CrossRef] 89. Li, S.Z.; Suo, Y.H.; Dai, L.M.; Liu, L.P.; Jiao, Q. Development of the Bohai Bay Basin and Destruction of the North China Craton. Earth Sci. Front. 2010, 17, 64–89.