Formation Characterization and Type Prediction Based on Geophysical Well Log Data in Horizontal Well: A Case Study of Triassic Chang 8 Formation in Shunning Region, Central Ordos Basin

: The role of the horizontal well in developing unconventional oil and gas reservoirs is particularly significant. Different from vertical wells, horizontal wells are greatly affected by many factors, e.g., well track, surrounding mudstone, resistivity, and pore structure heterogeneity in horizontal and vertical directions. These make it difficult to evaluate reservoir parameters and determine optimized test layers. In order to improve formation evaluation in horizontal wells, it is necessary to carry out the research of analyzing formation anisotropy, predict physical property parameters, and classify formation to determine high-quality intervals. In this study, taking Triassic Chang 8 Formation in Shunning Region, Central Ordos Basin as an example, 40 core samples were drilled and applied for nuclear magnetic resonance (NMR), mercury injection capillary pressure (MICP), and resistivity experiments. The porosity, permeability, resistivity, and pore structure anisotropy are analyzed. Results illustrate that the physical properties and pore structure in horizontal direction are superior to vertical direction. Meanwhile, NMR log loses its role in pore structure characterization in horizontal wells. Afterward, methods of correcting interval transit time (DT) and resistivity anisotropy are raised, and models of predicting formation physical property parameters, such as porosity, permeability, and water saturation, are established. A parameter, named as the formation type indicator ( FTI ), is proposed to reflect reservoir oil-bearing properties and pore structure. Finally, our target horizontal intervals are classified into four types, and the highest-quality “sweet spot” is determined.


Introduction
A drilling horizontal well is considered the primary method to explore unconventional oil and gas resources and increase hydrocarbon production. Compared with the vertical well, the horizontal well can traverse thousands of meters of intended intervals. Hence, it can reduce production costs and significantly improve oil and gas production [1]. The horizontal well technique has become an essential method for efficiently developing low permeability oil and gas reservoirs in the Ordos Basin [2]. In recent years, techniques and methods that focused on horizontal well log interpretation and formation evaluation were rare; models that were raised in the vertical well were still used [3]. However, log responses in horizontal intervals are absolutely different from that in vertical intervals

Geological Setting
Ordos Basin is the second-largest sedimentary basin in China. It used to be a part of the Greater North China continental sedimentary basin [19]. However, it was separated and transformed into a large inland basin in the late Triassic [20]. The Ordos Basin is characterized by a large asymmetric syncline with a gentle eastern flank and a steep western flank [21]. The Zhijing-Ansai oilfield is located in the west and central part of the Yinshan slope. It contains six regions, which are Jingbian in the north, Yongning in the south, Ansai in the east, Wuqi in the west, and Zhidan and Shunning in the central sections. It covers about 1.0 × 10 4 km 2 [22] (Figure 1). The Shunning Region is a part of the Zhijing-Ansai oilfield. It is located in the Yinshan slope structural belt and belonged to delta front deposition, close to Chang 7 high-quality source rock. It is a typical forming pattern of "upper generation and lower storage" [23,24]. Triassic Chang 8 Formation is the main oil-production zone, and formation depth ranges from 1500 to 2100 m. It is subdivided into Chang 81 and Chang 82 Formation. The formation is tight and has little natural productivity. Drilling horizontal well in the Chang 81 Formation is the main route that increases hydrocarbon production.
- Figure 1. Location map of our target Shunning Region in Central Ordos Basin (a); and the formation member, thickness and lithology (b) [23].

Lithologic Features
Based on the experimental data of the thin section that are acquired from 145 core samples, the lithology is mainly feldspar sandstone and lithic feldspar sandstone; the total content of quartz and felspar reaches to 68.88% (Figure 2a). The Chang 8 sandstones, which are close to the Chang 7 source rock, have a high content of feldspar and lithic dissolution and grain collapse from the acidic fluid injection. This also indirectly leads to the deterioration of reservoir physical properties [25]. The detrital composition is mainly metamorphic rocks and igneous rock (Figure 2b), and the interstitial matter is mainly composed of ferrocalcite and chlorite (Figure 2c).

Physical Characteristics
In Chang 8 Formation of the Shunning Region, 1988 core samples were recovered from 33 straight boreholes and applied for routine physical properties measurements. Figure 3a,b show the statistical histogram of core-derived porosity and permeability, separately. These figures illustrate that the porosities range from 4% to 12%, and the average porosity is 7.42% (Figure 3a). Permeabilities are predominantly distributed from 0.01 mD to 2 mD, and the average permeability is only 0.062 mD (Figure 3b). The relationship between the porosity and permeability is poor (Figure 3c). Figure 3a-c indicate that our target Chang 8 Member belongs to typical tight sandstone reservoirs and contains strong heterogeneity.

Anisotropy Characteristics
In Chang 8 Formation, we drilled 40 samples along a 90-degree angle to separately obtain 20 horizontal and 20 vertical core samples, following with the method displayed in Figure 4. These core samples were applied for nuclear magnetic resonance (NMR), mercury injection capillary pressure (MICP), and resistivity experimental measurements in the laboratory, respectively. The NMR Response Figure 5a,b exhibit the experimented NMR T2 distribution for 20 core samples separately. T2 spectra displayed in Figure 5a are acquired from core samples that were drilled from the horizontal direction, whereas Figure 5b shows the NMR T2 spectra that were acquired from vertical core samples. A comparison of these two figures illustrates that the shapes of NMR T2 distributions are similar. Irreducible water saturation (Swi) and T2 logarithmic mean (T2lm) that extracted from experimented NMR data for the horizontal and vertical core samples also exhibit approximate values ( Figure 6). These mean that NMR T2 distributions are invalid in indicating formation anisotropy.  The failure of NMR T2 distribution in reflecting the reservoir anisotropy can be explained as follows: for core samples with full water saturation, the NMR T2 transverse relaxation time is dominated by surface relaxation, while the bulk relaxation and diffusion relaxation can be ignored [26]. Surface relaxation reflects the energy exchange between hydrogen nuclei and pore surfaces. Rocks with big pores always need a long relaxation time, and vice versa. The energy exchange process is independent of the angle and direction but only affected by pore size (Figure 7) [27].

The MICP Curve
We compare the shapes of MICP curves for two sets of core samples that are extracted from horizontal and vertical core samples separately ( Figure 8). The location of MICP curves that are extracted from horizontal core samples is on the bottom, and the corresponding threshold pressures are low. This means that the pore structure in the horizontal direction is greater than that in the vertical direction. In Figure 9, we compare the pore structure characterization parameters (the average pore-throat radius Rm and median pore-throat radius R50) that are extracted from two different directions. A similar conclusion can be obtained. Through Figures 5-9, we can conclude that the pore structure anisotropy is strong in our target Chang 8 Formation in Shunning Region. The NMR log loses its role in pore structure characterization in the horizontal well, whereas with the capillary pressure, the pore-throat radius distribution will be of great importance.  Comparisons of experimental the average pore-throat radius (a) and median pore-throat radius (b) for 20 horizontal and vertical core samples. These two figures illustrate that the pore structure characterization parameters extracted from horizontal core samples are superior to those in vertical direction.

Resistivity Anisotropy
Forty core samples are applied for resistivity experiments by a gas-drive method in the laboratory. Rock resistivities under full water and different oil saturation are acquired. The cementation exponent and saturation exponent are obtained by using cross-plots of porosity and formation factor as well as water saturation versus resistivity index, separately. In Figure 10, we compare rock resistivity under full water saturation that is measured from horizontal and vertical core samples. It can be observed that the resistivities of horizontal core samples are lower than that of vertical core samples. This is caused by the better pore structure of vertical core samples. Except for rock resistivity, the cementation exponent m and saturation exponent n are also heavily affected by formation anisotropy (Figure 11a,b). Figure 11 indicates that the values of m and n are obviously lower in the horizontal direction than in the vertical direction. This coincided with the pore structure heterogeneity. Based on the negative correlation between pore structure and rock electrical conductivity, the pore structure of horizontal core sample is superior to vertical core samples, and the corresponding values of m and n should be lower [28].

Anisotropy Correction in Horizontal Well
Above experimental results and theoretical analysis have illustrated that the log responses of horizontal well are absolutely different from that of the vertical well in the same layer due to formation anisotropy. This leads to the established models in vertical intervals, which cannot be directly used in horizontal intervals. To predict precise formation parameters in horizontal intervals, the measured conventional log curves should be first corrected. In our target region, the measured conventional log curves contain natural gamma (GR), spontaneous potential (SP), caliper (CAL), interval travel time (DT), and deep and shallow lateral resistivity (RLLD and RLLS). Density and neutron logs are not acquired. GR, SP, and CAL have a negligible effect by the direction. However, DT, RT, and RXO are heavily affected. In this study, we correct DT and RT before they are used in horizontal interval characterization.

Anisotropy Correction of Interval Travel Time (DT or Δt)
In our target formation, no core samples are applied for acoustic experiments in laboratory, and the number of key wells is fewer. These means traditional interval travel time heterogeneity correction methods based on statistical regression cannot be used [9,10]. In this study, we raise a method to correct the effect of formation anisotropy to Δt, and this method covers several procedures.
First, selecting several vertical wells that are located around the target horizontal wells as control wells, collecting the Δt data of Chang 8 Formation, making a statistical histogram, collecting Δt value that represents the Chang 8 Formation based on the highest frequency, and naming it as Δth (Figure 12a).
Second, collecting Δt data of horizontal wells in Chang 8 Formation, making a statistical histogram, recording the Δt value that corresponds to the highest frequency, and defining it as Δtv (Figure 12b).
Third, calculating the correlation value of Δtv and Δth by using Equation (1).
where Δtd is the correlation value of interval travel time between horizontal and vertical directions, separately. Δth is the interval travel time that is acquired from horizontal direction. Δtv is the interval travel time that is acquired from vertical direction. The unit of them is μs/m. Combined with Figure 12a,b, the value of Δtd is determined as +10.0 μs/m. Finally, correcting the interval travel time of horizontal wells by adding Δtd to the Δt curve (habitually marked as DT) and making a statistical histogram of corrected interval travel time (Figure 12c). The correction procedures will stop once the shapes of the statistical histograms of interval travel time before and after the correction are similar.

Anisotropy Correction of Deep Lateral Resistivity
In this study, after the formation is assumed to be composed of a bedded sequence of parallel and discrete sandstone and shale (or thin beds), and sandstone anisotropy is considered small enough to be negligible, a method is raised to correct the anisotropy of deep lateral resistivity (RLLD). This correction method covers several procedures.
where GR is the measured gamma ray in the horizontal well, GRmin is the gamma ray of clean sandstone, and GRmax is the gamma ray of pure shale; the unit of them is API. SH is the relative shale content, Vsh is the shale volume content, and Vss is the sandstone volume content; the unit of them is v/v. φ is the porosity of the horizontal well in v/v. GCUR is a constant that is associated with the geologic period, and its value is 3.7. Second, assigning an initial value of the horizontal resistivity of shale as Rshho, and calculating formation anisotropy coefficient by using Equation (5).
where λ is the formation anisotropy coefficient. Rss is the sandstone resistivity in Ω·m.
Rshho is the initial value of shale horizontal resistivity in Ω·m. Third, substituting measured apparent resistivity Ra, borehole deviation angle α, and anisotropy coefficient into Equation (6) where Rh is the corrected horizontal resistivity and Ra is the measured apparent resistivity in the horizontal well; the unit of them is Ω·m. α is the borehole deviation angle in ( o ). Fourth, calculating vertical resistivity Rv and horizontal resistivity of shale Rshh by using Equations (7) and (8).
where Rv is the vertical resistivity with anisotropy correction, Rshh is the shale resistivity in horizontal direction, and the unit of them is Ω·m. Finally, the difference between Rshh and Rshho is compared to determine if the heterogeneity correction is finished. Once the difference between Rshh and Rshho meets the accuracy requirement (lower than 0.01 Ω·m), anisotropy correction is finished and the values of Rh and Rv are considered as true formation resistivity in horizontal and vertical directions.
Based on continuous iterative calculation, the optimal values of Rss and Rshh are defined as 110 Ω·m and 25 Ω·m, separately. By using the above-mentioned method, horizontal wells in the Shunning Region are processed, and the corrected Rh and Rv are consecutively determined. Figure 13 displays a comparison of corrected Rh and Rv with measured apparent deep lateral resistivity (RLLD). It can be observed that RLLD is absolutely different from Rh and Rv. Rh is lower than RLLD and Rv, which coincides with the experimental resistivity and can be explained by pore structure.

Porosity Calculation
We establish the porosity estimation model based on DT data because DEN and CNL are not acquired in the Chang 8 Formation. In the adjacent vertical wells, 19 core samples were recovered for routine experimental measurement. Core-derived porosity and permeability are extracted. After analyzing the core-derived and log data, we find that Raymer's equation is optimal in expressing the relationship between core-derived porosity and DT ( Figure 14). Hence, porosity estimation model is established as follows [29]: where φ is the porosity in %. Δtc is the interval transit time with anisotropy correction in μs/m. Δtma is the interval transit time of rock matrix in μs/m. Generally, its value is 182 μs/m. x is the undetermined coefficient. Its value needs to be calibrated by using the experimented data of core samples. In Chang 8 Formation of Shunning Region, its value is −0.612.

Permeability Prediction
Considering the poor relationship between core-derived porosity and permeability (Figure 3c), common permeability prediction method based on porosity cannot be used. Generally, permeability is heavily associated with the pore structure, and many models of calculating the permeability based on pore structure had been raised [30][31][32]. These models were proposed based on the MICP experimented data of core samples, such as the Swanson-based model, the R10 (pore-throat radius that corresponds to 10% mercury injection saturation)-based model, and the Rmax (maximal pore-throat radius)-based model [30,31]. These models are verified to be effective in tight sandstone reservoirs in Ordos Basin [32]. In this study, we analyze the MICP experimental results of 20 core samples that were drilled from horizontal direction and find that permeability is heavily associated with porosity and R10. Hence, a permeability prediction model based on R10 is raised and displayed in Figure 15. The model of calculating permeability from R10 and porosity is expressed as Equation (10).
where K is the permeability in horizontal direction in mD; φ is the porosity in horizontal direction in %; R10 is the pore-throat radius that corresponds to 10.0% mercury injection saturation in μm.  Figure 15. Relationship between core-derived porosity, permeability, and R10.

Water Saturation Calculation
Water saturation calculation models are abundant after the Archie's equation was first raised in 1942 [33]. These models contain the Simandoux equation, the Indonesian equation, the Waxman-Smits model, the dual-water model, etc. [34][35][36][37][38][39]. These models are all raised to solve the water saturation calculation in shaly sands. However, no specific models have been proposed in tight sandstone reservoirs, and the Archie's equation is still used. The derivative Archie's equation is expressed: where Rw is the resistivity of formation water in Ω·m; Rt is the resistivity of rock under hydrocarbon saturation in Ω·m; Sw is the water saturation in v/v; a is the pore tortuosity factor; m is the cementation exponent; b is a constant that associated with lithology; and n is the saturation exponent. a, b, m, and n are defined as the rock-resistivity parameters, and their values are determined by the resistivity experiments of rocks.
In the Chang 8 Formation of the Shunning Region, we have chosen a total of 12 core samples (six from the horizontal direction and six from the vertical direction) to apply for resistivity experiments in laboratory. The relationships between porosity and formation factor, and water saturation versus resistivity index, are established in horizontal and vertical directions and displayed in Figures 16 and 17 separately. Comparing these two figures, some regularity can be concluded: (1) The resistivity properties in horizontal and vertical directions are absolutely different; this diversity can be observed by the difference of the values of a, b, m, and n. The values of m and n in the horizontal direction are lower than that of in vertical direction. (2) The relationships between water saturation versus resistivity index are all perfect. This means that the pore structure slightly affects the saturation exponent, and they can be well-used in water saturation calculation (Figures 16b  and 17b). (3) Although the correlation coefficients between porosity and formation factor are high, the regressed trendline cannot well express the relation, especially for the part with porosity lower than 8% (Figures 16a and 17a). This can be demonstrated by the fact that the trend line does not cross the boundary point (1,1). This regularity is similar to that extracted from tight gas sands in the Xujiahe Formation in the Central Sichuan Basin [28].  To improve the accuracy of the water saturation calculation by using the Archie's equation, the cementation exponent should be optimized. Based on the theoretical analysis,  concluded a relation between porosity and the cementation exponent [28]. This theoretical relation can be well-used to express the variation between pore structure and cementation exponent, especially for tight rocks with porosities lower than 8%. The theoretical relation raised by  has been verified to be widely usable in different types of tight reservoirs [40]. In this study, using the resistivity experimental data that are displayed in Figures 16 and 17, we also analyze the effect of pore structure to cementation exponent in Chang 8 Formation of Shunning Region and determine that the cementation exponent is heavily affected by porosity. The cementation exponent decreases significantly with the increase of porosity. Based on this analysis, we optimize the relationships between porosity and formation factor and deduce two formulae to calculate various cementation exponents from porosity in horizontal and vertical directions separately (Figure 18a,b). These formulae can be expressed by Equations (12) and (13).
Horizontal direction: m = 0.35 × log(φ) + 2.23 (12) Vertical direction: m = 0.52 × log(φ) + 2.54 (13) where m is the cementation exponent; φ is porosity in v/v. It should be noted that the value of a is defined as 1 in optimized relationships. Comparing Figures 16 and 17 with Figure 18, it can be easily observed that the trend lines between porosity and formation factor are deflexed once the porosities of core samples decrease. The optimized cementation exponent can well express this bendability. However, unreliable water saturation should be calculated if we directly define fixed cementation exponent in the Archie's equation. Combining Equation (11) and Figure 16 to Figure 18, we separately conclude optimal water saturation calculation models in horizontal and vertical directions based on the Archie's equation, and these two models are expressed as:  (15) where Sw_h and Sw_v are the calculated water saturation in horizontal and vertical directions, separately. The unit of them is v/v.
By substituting the anisotropy corrected resistivity Rh and Rv, porosity into Equations (14) and (15), Sw_h and Sw_v can be calculated. To verify the reliability of resistivity correction in the horizontal well, we also directly use the measured apparent deep lateral resistivity to calculate water saturation (Sw). Comparison of calculated Sw_h and Sw_v with Sw in a horizontal well is displayed in Figure 19. It can be observed that Sw_h and Sw_v are almost equal, whereas the value of Sw is about 5% to 8% higher. This means that the resistivity anisotropy correction is essential before it is applied for formation evaluation, or else the hydrocarbon volume will be underestimated. Either Sw_h or Sw_v can represent true formation water saturation.

Formation Classification Based on Pore Structure
The ultimate purpose of formation evaluation is to determine the highest-quality intervals, and offer the optimal development zones. Hence, establishing formation classification criteria is indispensable. Generally, formation quality is heavily controlled by pore structure, oil volume, and permeability, and they significantly affect the movable fluid saturation. The greater the change in permeability, the more significant the change in movable fluid saturation. Pore structure is the main factor that contributes to the obviously different occurrence characteristics of movable fluid [41]. Formation pore structure can be characterized by the shape of the pore-throat radius. Generally, rocks with good pore structure always contain wide pore-throat radius distribution, and it exhibits bimodality or multimodality. However, narrow unimodal pore-throat radius always corresponds to relatively simple pore structure. Hence, many parameters that are extracted from pore-throat radius distributions, e.g., the average pore-throat radius (Rm), Rmax, the median pore-throat radius (R50), are used to represent the pore structure. In the Chang 8 Formation of Shunning Region, the pore-throat radius distributions are acquired from MICP experimental data of 20 core samples in the vertical direction, and the relationship between formation quality and pore-throat radius distribution is compared. Figure 20a,c,e display the pore-throat radius spectra of three typical core samples of No.1, No.5, and No.13, and Figure 20b,d,f display the experimented NMR T2 distributions of these three core samples. Comparisons clearly illustrate that the pore-throat radius distribution match well with NMR. Hence, the pore-throat radius distributions can replace NMR T2 spectra to characterize formation pore structure once the NMR data had not been consecutively acquired. Figure 20. Comparison of pore-throat radius spectra and the corresponding NMR T2 distributions of rocks with good pore structure (a,b), medium pore structure (c,d) and simple pore structure (e,f) in Chang 8 Formation of Shunning Region. Li et al. (2023) pointed out that the pore-throat radius distribution can be reconstructed from conventional well logs based on the deep-learning method [42]. This method has been verified to be widely usable in the Ordos Basin [43]. Hence, the raised method by Li et al. (2023) is directly used in our target Chang 8 Formation in the horizontal well. The pore-throat radius distributions are consecutively constructed. Afterward, Rm is consecutively calculated to reflect formation pore structure characteristic.
In addition, hydrocarbon volume is also an important parameter that controls formation quality. Generally, high-quality reservoirs always contain high hydrocarbon volume and, thus, good production. On the contrary, once the hydrocarbon volume is low, the corresponding reservoir quality is also poor [44]. Hence, hydrocarbon (oil in the Chang 8 Formation) volume is also considered as an input parameter in establishing formation classification criteria in the horizontal well. Based on the definition, oil volume is calculated as: where fo is the oil volume in %; So is the oil saturation in v/v.
Based on the above analysis, we chose two parameters of Rm and fo, which respectively reflect the pore structure and formation quality to establish the formation type indicator parameter. This parameter is calculated as follows: where FTI is the formation type indicator parameter in the horizontal well. By using Equation (17), we process the Chang 8 Formation in 17 horizontal wells, which contain drill stem test (DST) data, in our target region, and FTIs are calculated. Relationships among FTI and formation physical property parameters (porosity, permeability, and DT) are analyzed. Finally, we establish formation classification criteria ( Figure  21). Based on the criteria, our target Chang 8 Formation is classified into four types. The classification criteria are listed in Table 1. The first type of formation contains the best pore structure, the highest fo, and, thus, the highest-quality and oil production. The pore structure and fo of the second type of formation are relatively weaker than that of the first type of formation, and the quality is also poorer. The third type of formation is poorer oil-bearing formation due to lower FTI and physical property parameter. The last type of formation is the dry layer; no development program should be made for this type of formation.

Case Studies
The proposed models and formation classification criteria are applied to our target Chang 8 Formation in two horizontal wells B and C. After the anisotropy of DT and RT are corrected, formation porosity, water saturation, and oil volume are calculated. Afterwards, FTI is calculated to classify horizontal intervals and determine the highest-quality "sweet spot" area. Figures 22 and 23 display the formation classification results based on the criteria that are listed in Table 1. The evaluation results illustrate that the first and second types of formations are dominant in the horizontal interval in Well B; the total amount reaches to 76.4%. Especially for the first type of formation, it contributes the most to oil production capacity. The fourth type of formation can be ignored. This means that the reservoir quality is greater. However, in the horizontal interval in Well C, the proportion of the first type reservoir is only 3.4%, but the amount of the third reservoir reaches to 52.6%.
Comparison of formation classification in these two wells illustrates that Well B has great0er reservoir quality despite the longer horizontal interval of Well C. Hence, we predict that Well B will contain much more oil production capacity. This prediction is verified by the drill stem test data. In Well B, approximately 74.16 bbl of oil is produced per day. However, in Chang 8 Formation of Well C, only 26.28 bbl of oil is produced per day; the daily oil production of Well C is only one-third as much as that of Well B.  To quantize the relationship between formation quality and oil production, we extract the total proportion of first and second types of formation (these two types of formation have a dominant contribution to oil production capacity in Chang 8 Formation) and find that it is heavily positively associated with daily oil production ( Figure 24). With the increase of the total proportion of the first and second types of formation, daily oil production per meter increases. Figure 24. Cross-plot of the total proportion of 1st and 2nd types of formation versus daily oil production. Positive relationship means the dominant contribution of formation quality to oil production capacity, especially the high-quality formation.

Discussions
In this study, we raise methods and models of correcting the anisotropy of interval transit time and deep lateral resistivity, calculating formation parameter and classifying formation type in horizontal wells, and well results are acquired in our target Chang Formation in Shunning Region, Central Ordos Basin. However, these methods and models are all established based on experimented results of typical core samples and drill stem test data. This means they are only valuable in the Chang 8 Formation. In other regions, these methods and models, especially the formation classification criteria, may not be available unless they contain similar formation properties. If we want to apply the established methods and models in any other formation, a certain amount of core samples should be first drilled for petrophysical experiments to calibrate the involved parameters in our raised models. In addition, the corresponding drill stem test data should also be acquired to adjust the established criteria in Table 1.

Conclusions
The Triassic Chang 8 Formation in Shunning Region belongs to typical tight sandstone reservoirs, and the heterogeneity of interval travel time, resistivity, and pore structure are strong. It is difficult to accurately calculate formation physical property parameters and characterize the pore structure in horizontal wells before the anisotropy of interval travel time and resistivity is first corrected.
After the anisotropy of apparent resistivity of the horizontal well is corrected, two resistivities, the horizontal and vertical resistivities, are acquired. The Archie's equation is respectively optimized, and the water saturations in horizontal and vertical directions are calculated. Consistency of these two water saturations verifies the reliability of resistivity anisotropy correction. Meanwhile, interval travel time anisotropy is also corrected to predict formation porosity in horizontal intervals.
Relationships among hydrocarbon volume, pore structure, and formation quality are analyzed, and we raise an assembly parameter of FTI to identify high-quality formation in Chang 8 Member. Combined with FTI and drill stem test data, our target horizontal intervals of Chang 8 Formation are classified into four types. The first type of formation contains the highest quality, and it has the greatest contribution to oil production capacity. The pore structure and FTI of second and third types of formation are poorer than that of the first type of formation, the formation quality and production capacity are also poorer. The fourth type of formation is dry layer; no hydrocarbon can be produced.
The total proportion of first and second types of formation is positively associated with daily oil production. Hence, once the pore structure is first characterized, the proportion of four types of formation can be identified. Afterward, the total proportion of first and second types of formation can be calculated to predict oil production capacity in our target horizontal intervals in Chang 8 Formation of the Shunning Region.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to restrictions of privacy.

Conflicts of Interest:
The authors declare no conflict of interest.