A New Evaluation Method of Shale Oil Sweet Spots in Chinese Lacustrine Basin and Its Application

: Compared with the marine shale from North America, Chinese lacustrine basins have more complex geological and reservoir conditions, which makes the selection of sweet spot intervals in the shale oil reservoir particularly difﬁcult. It is thus crucial to accurately predict the shale oil sweet spots for effective exploration and development of shale oil in a lacustrine basin. In this paper, we propose an innovative evaluation method of shale oil sweet spots, which considers ﬁve primary parameters (i.e., total oil content, movable oil ratio, reservoir pressure coefﬁcient, permeability, and rock brittleness index) to construct a comprehensive weighting factor, which is used to quantitatively identify a favorable shale oil reservoir. This method ﬁrstly constructs an evaluation function for each of the parameters, and then calculates a comprehensive weighting factor to determine the shale oil sweet spot. Statistical results show that the oil production of formation testing intervals have a good positive correlation with the average value of the comprehensive weighting factor of the intervals, which veriﬁes the feasibility of the method. Based on this method, one of the key exploratory wells, Qiang21 in the Raoyang Sag of Bohai Bay basin, was selected as a case study and was determined to be a sweet spot interval of the shale oil reservoir in the upper third member of the Shahejie Formation. This study provides a new way to obtain a favorable exploration interval of shale oil reservoirs and serves shale oil development.


Introduction
As a promising alternative to conventional petroleum resources, shale oil has received a lot of attention in recent years [1][2][3][4]. However, the exploration and development of shale oil is still at a relatively early stage due to the extremely great challenge of shale oil recovery [5]. Shale oil refers to the petroleum that occurs in the shale-dominated strata in a free, dissolved, or adsorbed state, of which free and adsorbed oil are the two main constituents [6][7][8]. Shale oil is the accumulation of residual hydrocarbons in organic-rich shales or juxtaposed (overlying, underlying, or embedded) organic-lean lithofacies after crude oil is expelled from source rocks, which contributes to conventional oil reserves [9,10]. Usually, shale oil underwent a primary migration over a short distance, or even underwent no migration, within the shale reservoir. Compared with the shale oil in North America, which is mainly derived from marine strata [9,11], China's shale oil mainly comes from the mudstone-shale strata in lacustrine basins, such as the Cretaceous Qingshankou Formation

Geological Settings
The Raoyang Sag is a secondary structural unit in the central Jizhong Depression of the Bohai Bay Basin (Figure 1), and covers an area of about 5280 km 2 [23,24]. The Sag is one of the most promising shale oil accumulation and development areas in the Jizhong Depression. There are two sets of main source rocks in the Lower first member of the Shahejie Formation (Es 1 L ) and the upper third member of the Shahejie Formation (Es 3 U ) in the Sag, which are main sources for the formation of shale oil. After more than 40 years of exploration and development in the Raoyang Sag, oil shows have been commonly encountered in shale formations. Among them, more than 200 wells in the shale strata in the Es 1 L and Es 3 U have presented oil shows. In this study, a key exploratory well, Qiang21, was primarily studied to determine favorable intervals of shale oil sweet spots. A total of 48 samples located at the Es 3 U were obtained from the well. Based on the basic measurements of organic geochemistry, they show that: (a) the total organic carbon (TOC) content ranges from 0.112 to 5.536% (mean 2.320%); (2) S 1 content is the pyrolysis hydrocarbon content obtained from Rock-Eval measurement at 573.15K and varies from 0.01 to 1.35 mg/g (mean 0.48 mg/g); and (3) the production index (S 1 /(S 1 + S 2 ), where S 2 is Rock-Eval measured kerogen yield) has a range of 0.17-41.13 mg/g (mean 15.12 mg/g). Overall, studied shales belong to good source rocks and are favorable to shale oil enrichment ( Figure 2). Inorganic mineral components of shales are composed of clay minerals, quartz, and calcite, with weight percentages In this study, a key exploratory well, Qiang21, was primarily studied to determine favorable intervals of shale oil sweet spots. A total of 48 samples located at the Es3 U were obtained from the well. Based on the basic measurements of organic geochemistry, they show that: (a) the total organic carbon (TOC) content ranges from 0.112 to 5.536% (mean 2.320%); (2) S1 content is the pyrolysis hydrocarbon content obtained from Rock-Eval measurement at 573.15K and varies from 0.01 to 1.35 mg/g (mean 0.48 mg/g); and (3) the production index (S1/(S1 + S2), where S2 is Rock-Eval measured kerogen yield) has a range of 0.17-41.13 mg/g (mean 15.12 mg/g). Overall, studied shales belong to good source rocks and are favorable to shale oil enrichment ( Figure 2). Inorganic mineral components of shales are composed of clay minerals, quartz, and calcite, with weight percentages of 37.9 wt.%, 18.7 wt.%, and 17.9 wt.%, respectively. Shales contain a small amount of aragonite

Determination of Primary Parameters
A comprehensive sweet spot of shale oil is a superposition of geological and engineering sweet spots. However, due to the inconsistency between various types of sweet spots, it is difficult to quantitatively determine the comprehensive sweet spot for shale oil recovery. It is necessary to select primary parameters to characterize the various sweet spots. Geological sweet spots mainly include good oiliness, strong mobility of oil, and high permeability, which can be expressed by the four parameters of oil content, movable oil ratio, pressure coefficient of reservoir fluid, and shale permeability. The engineering sweet spot mainly focuses on the evaluation of fracturing performance of rock, which can be reflected by a brittleness index.

Shale Oil Content
Total oil content in shale can directly characterize the degree of shale oil enrichment, and is a direct manifestation of the abundance of shale oil resources. There are two parameters that can be used to reflect shale oil content, which are pyrolysis hydrocarbon (S1) and chloroform bitumen "A". The former is obtained by the rock pyrolysis analyzer (Rock-Eval); the latter can be obtained by the chloroform extraction method [25][26][27]. It should be noted that the pyrolysis hydrocarbon (S1) is obtained when the heating temperature does not exceed 573.15K in a Rock-Eval experiment and is the result of a loss of light and heavy hydrocarbons. Thus, it needs to be restored to be able to evaluate the resource abundance of shale oil. However, chloroform bitumen "A" includes pyrolysis hydrocarbon S1 and some heavy components. According to the relationship of chloroform bitumen "A" with S1 (Figure 4a), the recovery coefficient (K1) of heavy hydrocarbons can be obtained. This was determined to be about 2.1 in this study. Thus, the corrected S1 value ( 1 S ' ) after the recovery of the heavy components can be calculated as: where 1 S ' is the corrected S1 value with a heavy hydrocarbon recovery, mg/g, and K1 is the recovery coefficient of heavy hydrocarbons, dimensionless.

Determination of Primary Parameters
A comprehensive sweet spot of shale oil is a superposition of geological and engineering sweet spots. However, due to the inconsistency between various types of sweet spots, it is difficult to quantitatively determine the comprehensive sweet spot for shale oil recovery. It is necessary to select primary parameters to characterize the various sweet spots. Geological sweet spots mainly include good oiliness, strong mobility of oil, and high permeability, which can be expressed by the four parameters of oil content, movable oil ratio, pressure coefficient of reservoir fluid, and shale permeability. The engineering sweet spot mainly focuses on the evaluation of fracturing performance of rock, which can be reflected by a brittleness index.

Shale Oil Content
Total oil content in shale can directly characterize the degree of shale oil enrichment, and is a direct manifestation of the abundance of shale oil resources. There are two parameters that can be used to reflect shale oil content, which are pyrolysis hydrocarbon (S 1 ) and chloroform bitumen "A". The former is obtained by the rock pyrolysis analyzer (Rock-Eval); the latter can be obtained by the chloroform extraction method [25][26][27]. It should be noted that the pyrolysis hydrocarbon (S 1 ) is obtained when the heating temperature does not exceed 573.15K in a Rock-Eval experiment and is the result of a loss of light and heavy hydrocarbons. Thus, it needs to be restored to be able to evaluate the resource abundance of shale oil. However, chloroform bitumen "A" includes pyrolysis hydrocarbon S 1 and some heavy components. According to the relationship of chloroform bitumen "A" with S 1 (Figure 4a), the recovery coefficient (K 1 ) of heavy hydrocarbons can be obtained. This was determined to be about 2.1 in this study. Thus, the corrected S 1 value (S 1 ) after the recovery of the heavy components can be calculated as: where S 1 is the corrected S 1 value with a heavy hydrocarbon recovery, mg/g, and K 1 is the recovery coefficient of heavy hydrocarbons, dimensionless. Based on the recovery of heavy hydrocarbons, the relationship between the ratio of C 5-14 to C 14+ and the thermal maturity (R o ) of shale is determined according to the experimental data from the thermal simulation analysis of hydrocarbon generation. Further, according to the relationship between burial depth and R o of shale, the recovery coefficient (1 + K 2 ) of light hydrocarbon can be determined at different burial depths (Figure 4b), and the total oil content after the recovery of light hydrocarbons can be obtained, which is expressed as: where S 1 is the total oil content, mg/g, and 1 + K 2 is the recovery coefficient of light hydrocarbons, dimensionless.
hydrocarbons, dimensionless. According to the above analysis, the burial depth of the shale oil reservoir shoul an important parameter. However, we did not choose burial depth as one of the m geological parameters, mainly because burial depth is closely related to thermal matu (Ro) and influences the shale oil content. Thus, the burial depth has been indirectly sidered in the proposed evaluation method. Due to the constraints of analysis cost, and especially the limited number of s samples that can be provided, it is difficult to continuously describe the organic heter neity of shales with measured data. The characteristics of organic matter and oiliness h obvious effects on many logging responses (acoustic wave, resistivity, neutron, den gamma, etc.) [20,28], and the logging data are characterized by continuity and high r lution. This makes it possible to use logging data to evaluate the organic heterogenei the shale reservoir. The most commonly used prediction method in the industry is ΔLogR method, which can continuously predict the two important parameters of organic carbon (TOC) and S1 in shales based on logging curves of acoustic wave and sistivity [28,29]. In this study, the S1 and TOC values of shales in the Well Qiang21 w determined as linear functions of ΔLogR. Calculated results are presented in Figu From the figure, it can be seen that there is a good match between calculated and m ured results. Further, total oil content can be obtained. It can be seen that oil conte relatively high at the depth of 2926-2955 m, where a shale oil-rich interval exists. According to the above analysis, the burial depth of the shale oil reservoir should be an important parameter. However, we did not choose burial depth as one of the main geological parameters, mainly because burial depth is closely related to thermal maturity (R o ) and influences the shale oil content. Thus, the burial depth has been indirectly considered in the proposed evaluation method.
Due to the constraints of analysis cost, and especially the limited number of shale samples that can be provided, it is difficult to continuously describe the organic heterogeneity of shales with measured data. The characteristics of organic matter and oiliness have obvious effects on many logging responses (acoustic wave, resistivity, neutron, density, gamma, etc.) [20,28], and the logging data are characterized by continuity and high resolution. This makes it possible to use logging data to evaluate the organic heterogeneity of the shale reservoir. The most commonly used prediction method in the industry is the ∆LogR method, which can continuously predict the two important parameters of total organic carbon (TOC) and S 1 in shales based on logging curves of acoustic wave and resistivity [28,29]. In this study, the S 1 and TOC values of shales in the Well Qiang21 were determined as linear functions of ∆LogR. Calculated results are presented in Figure 5. From the figure, it can be seen that there is a good match between calculated and measured results. Further, total oil content can be obtained. It can be seen that oil content is relatively high at the depth of 2926-2955 m, where a shale oil-rich interval exists.

Movable Oil Ratio and Reservoir Pressure Coefficient
In addition to considering the degree of oil resource enrichment, the mobility of shale oil is also closely related to productivity [30]. The mobility of shale oil mainly includes the two aspects of flow potential and flow capacity of oil, which is similar to the mobility of water [31]. Flow potential is characterized by the movable oil ratio. The higher the movable ratio of shale oil, the more favorable it may be for shale oil extraction. Flow capacity is characterized by the pressure coefficient of reservoir fluid. Abnormally high pressure (i.e., overpressure) can provide sufficient elastic energy for shale oil seepage through the reservoir [32].

Movable Oil Ratio
Theoretically, total oil content includes two parts: movable oil and bound oil, of which the movable oil is composed of the free hydrocarbons in the measured S 1 and the recovered light hydrocarbons, and the bound oil includes the adsorbed hydrocarbons in the measured S 1 and the recovered heavy hydrocarbons. The key is to determine the adsorbed and free oil contents in S 1 . Jarvie showed that the ratio of S 1 to TOC can be used to characterize the lower limit of shale adsorbed hydrocarbons [9], so it can be used to evaluate the adsorbed and free amounts of shale oil. Determination of these values was based on the understanding that oil begins expulsing from organic matter as the shale reservoir submerges down to the hydrocarbon expulsion threshold depth during burial history; that is, oil achieves a saturated sorption on the organic matter and begins migrating into various pores at the hydrocarbon expulsion threshold. In this study, several relationships of burial depth with S 1 /TOC, (S 1 + S 2 )/TOC, S 1 , and hydrogen index (HI) were obtained, respectively, as shown in Figure 6. It can be seen that the depth of the hydrocarbon expulsion threshold occurs at around 2850 m, which means that hydrocarbons begin to expel and the shale has the maximum adsorbed oil content at this depth. It can be determined that if the saturated adsorption oil content (Q c ) of shales is equal to about 50 mg/g, then the adsorbed oil amount in S 1 is equal to the product of TOC multiplied by Q c . Therefore, the amounts of bound and movable oil and the ratio of movable oil are obtained as: where Q b is the bound oil content, mg/g; Q c is the saturated adsorption oil content, mg/g; Q m is the movable oil content, mg/g; and R m is the movable oil ratio, %.

Movable Oil Ratio and Reservoir Pressure Coefficient
In addition to considering the degree of oil resource enrichment, the mobility of shale oil is also closely related to productivity [30]. The mobility of shale oil mainly includes the two aspects of flow potential and flow capacity of oil, which is similar to the mobility of water [31]. Flow potential is characterized by the movable oil ratio. The higher the mova- where Qb is the bound oil content, mg/g; Qc is the saturated adsorption oil content, mg/g; Qm is the movable oil content, mg/g; and Rm is the movable oil ratio, %. Calculated results were presented in Figure 5. It can be seen that shale oil has a high mobility and movable ratio, with an average movable ratio of about 35% in the depth section of 2926-2955 m, where a shale oil-rich interval is also found.

Pressure Coefficient of Reservoir Fluid
There are many logging methods predicting reservoir fluid pressure, mainly including: (a) the equivalent depth method and Eaton method, based on undercompaction theory, (b) the effective stress method and Bowers method, based on rock mechanics theory, and (c) the empirical statistical method [33][34][35][36][37][38]. Among them, the Eaton method is an effective method of reservoir pressure prediction based on acoustic logging data, rock density logging data, and vertical effective stress of the rock. In addition, the method is based on the actual measured reservoir pressure [33,34]. Calculated results were presented in Figure 5. It can be seen that shale oil has a high mobility and movable ratio, with an average movable ratio of about 35% in the depth section of 2926-2955 m, where a shale oil-rich interval is also found.

Pressure Coefficient of Reservoir Fluid
There are many logging methods predicting reservoir fluid pressure, mainly including: (a) the equivalent depth method and Eaton method, based on undercompaction theory, (b) the effective stress method and Bowers method, based on rock mechanics theory, and (c) the empirical statistical method [33][34][35][36][37][38]. Among them, the Eaton method is an effective method of reservoir pressure prediction based on acoustic logging data, rock density logging data, and vertical effective stress of the rock. In addition, the method is based on the actual measured reservoir pressure [33,34].
In the study area, reservoir fluid pressure was divided into four classes according to the pressure coefficient: atmospheric pressure (0.9-1.1), weak overpressure (1.1-1.2), medium overpressure (1.2-1.5), and strong overpressure (>1.5). The reservoir fluid pressure at a certain depth is predicted by an Eaton equation (Equation (6)) combined with a normal compaction trend line. The Eaton index in the equation is obtained from the measured reservoir pressure data and is expressed as Equation (7). After determining the Eaton index and normal compaction trend line, reservoir fluid pressure can be calculated by Equation (6). Several main equations are as follows: where P p is the reservoir fluid pressure, MPa; P 0 is the overburden pressure, MPa; P n is the hydrostatic column pressure, MPa; AC n is the sonic differential time on the normal compaction trend line, µs/ft; AC is the logging data of sonic differential time, µs/ft; and c is the Eaton index. According to the result of reservoir pressure, the pressure coefficient (σ p ) can be obtained as: Calculated results were presented in Figure 5. Overall, the reservoir pressure rank belongs to atmospheric pressure with a reservoir coefficient of 0.9-1.1. Moreover, reservoir pressure has an increasing trend with the increase in burial depth.

Shale Permeability
Previous studies found that the production of shale oil is positively correlated with the physical properties of the reservoir [39,40]. The better the physical properties, the higher the production. One of the key parameters that characterizes the physical properties of the reservoir is permeability. Permeability is one of the important parameters for reservoir evaluation. It is not only related to porosity, but also closely related to the pore structure of the rock. The view that the higher the porosity is, the better the permeability is has been guiding the production of medium-high porosity/permeability reservoirs [5,41]. However, as the porosity decreases, the relationship between porosity and permeability becomes more complicated. When the reservoir has low porosity and low permeability, the difference in permeability reaches nearly 4 orders of magnitude under similar porosity situations [42]. Usually, reservoir porosity is relatively easily obtained by logging evaluation. However, due to the complexity of shale oil reservoir permeability, the relationship of porosity with permeability as established by traditional methods is not effective. Hearn et al. [43] proposed the concept of hydraulic flow unit (HFU), that is, different hydraulic flow units have varied petrophysical characteristics. The proposal of hydraulic flow unit effectively solves the complexity of shale oil reservoir porosity and permeability. Porosity and permeability of the reservoir have a good correlation for the same HFU, and the relationships of porosity with permeability in different HFUs are established in this study. Based on the porositypermeability relationships, shale permeability can be predicted [44][45][46]. The calculation method is expressed as follows: where FZI is the flow zone index, 10 −6 µm 2 ; RQI is the reservoir quality index; K is the permeability, mD; φ is the porosity, %; and PMR is the normalized porosity, dimensionless. According to the above method, this study divided the shale reservoirs of Raoyang Sag into 4 types of HFU based on the cumulative density distribution of the reservoir flow zone index (FZI) (Figure 7a). For each of the HFUs, there is a good correlation between porosity and permeability (Figure 7b). Based on the BP neural network method, the relationships of logging data with FZI and porosity are established to predict shale oil reservoir permeability (Figure 7c,d). Calculated results are shown in Figure 8. It can be seen that there is a low value interval located at 2870-2960 m, where shale has low porosity and permeability.

Brittleness Index
The brittleness of the shale reservoir is one of the key parameters of shale oil development and determines the fracturability of shale oil reservoirs. The brittleness index of rock is closely related to shale's hydraulic fracturing, which directly affects the seepage channels of shale oil. The existence of a large number of brittle minerals is the reason why rock can obtain high-yield oil flow through fracturing fractures [17,47]. An area with a high brittleness index is called an engineering sweet spot. Sag into 4 types of HFU based on the cumulative density distribution of the reservoir flow zone index (FZI) (Figure 7a). For each of the HFUs, there is a good correlation between porosity and permeability (Figure 7b). Based on the BP neural network method, the relationships of logging data with FZI and porosity are established to predict shale oil reservoir permeability (Figure 7c,d). Calculated results are shown in Figure 8. It can be seen that there is a low value interval located at 2870-2960 m, where shale has low porosity and permeability.

Brittleness Index
The brittleness of the shale reservoir is one of the key parameters of shale oil development and determines the fracturability of shale oil reservoirs. The brittleness index of rock is closely related to shale's hydraulic fracturing, which directly affects the seepage channels of shale oil. The existence of a large number of brittle minerals is the reason why rock can obtain high-yield oil flow through fracturing fractures [17,47]. An area with a high brittleness index is called an engineering sweet spot.
Evaluation of rock brittleness is mainly performed by calculating the brittleness index based on mineral compositions, or by calculating Young's modulus and Poisson's ratio based on rock mechanics data. In this paper, the brittleness index is used to characterize the rock brittleness of the shale reservoir in the study area. The inorganic mineral com- Evaluation of rock brittleness is mainly performed by calculating the brittleness index based on mineral compositions, or by calculating Young's modulus and Poisson's ratio based on rock mechanics data. In this paper, the brittleness index is used to characterize the rock brittleness of the shale reservoir in the study area. The inorganic mineral composition of shale is simplified into three parts: clay minerals, felsic minerals (quartz and feldspar), and calcareous minerals (calcite and dolomite). Based on the measured and logging data, the BP neural network method was used to establish prediction models for clay minerals, felsic minerals, and calcareous minerals, respectively. It should be noted that most studies in the past have neglected the influence of organic matter content on rock brittleness, however, as a rock component, the mechanical properties of organic matter are quite different from those of inorganic minerals [22]. Therefore, this research carried out the calculation of brittleness index considering the TOC content, and the expression was written as: where B is the rock brittleness index, dimensionless; C 1 is the felsic mineral content, %; C 2 is the calcareous mineral content, %; C TOC is the TOC content, %; and C 3 is the content of clay minerals, %. As shown in Figure 9, continuous inorganic minerals in the shale interval were predicted by the BP neural network method. There is a good match between the calculated and the measured result. Based on organic and inorganic compositions, the brittleness index of shales and their distribution in the well Qiang21 were calculated. Overall, there is a good brittleness in the studied shales. The brittleness index of shales varies from 0.147 to 0.829, with an average value of 0.514.

Calculation of Comprehensive Weighting Factor
Taking into account the entire study interval/area, the five primary parameters are normalized. The oil content and brittleness index are linearly normalized with the maximum value, the permeability is normalized by the logarithmic range, and evaluation values of the movable oil ratio and pressure coefficient can be used directly. They can be mathematically expressed as:

Calculation of Comprehensive Weighting Factor
Taking into account the entire study interval/area, the five primary parameters are normalized. The oil content and brittleness index are linearly normalized with the maximum value, the permeability is normalized by the logarithmic range, and evaluation values of the movable oil ratio and pressure coefficient can be used directly. They can be mathematically expressed as: where U 1 is the weighting factor of oil content, dimensionless; S 1 is the shale oil content after the recovery of light and heavy hydrocarbons, mg/g; S 1 max is the maximum shale oil content in the study area, mg/g; U 2 is the weighting factor of the movable oil ratio, dimensionless; S m is the movable content of shale oil, mg/g; U 3 is the weighting factor of the pressure coefficient of reservoir fluid, dimensionless; f is the pressure coefficient of reservoir fluid, dimensionless; U 4 is the weighting factor of permeability, dimensionless; lgk is the logarithm of shale permeability, mD; lgk max is the maximum logarithm of the shale permeability in the study area, mD; lgk min is the minimum logarithm of the shale permeability in the study area, mD; U 5 is the weighting factor of the shale brittleness index, dimensionless; B is the shale brittleness index, dimensionless; and B max is the maximum value of the shale brittleness index in the study area, dimensionless.
In order to quantitatively determine the comprehensive shale oil sweet spot, the weighting factors of five primary parameters are multiplied together, and then the order of magnitude is enhanced by 1000 times to obtain the final comprehensive weighting factor U, which is expressed as: According to the vertical distribution of comprehensive weighting factors, the shale oil comprehensive sweet spot can be found.

Comprehensive Evaluation Results
According to statistical results from some exploration wells in the study area, it can be seen that initial production of shale oil will increase with an increasing comprehensive weighting factor, especially when the comprehensive weight factor is larger than about 20 ( Figure 10). Correspondingly, there is a decrease in the water production with an increasing comprehensive weighting factor. Based on the results calculated by Equation (18), the sweet spot of shale oil was determined. The average comprehensive weighting factor of the interval 2926-2955 m is 65.7, which is a favorable interval for shale oil development ( Figure 11). The shale reservoir, within this interval, has a relatively high oil content, a high mobility ratio, and a high brittleness index, which are favorable factors for shale oil extraction; however, the shale has a normal formation pressure coefficient and relatively poor physical properties, which may increase the difficulty in fracturing simulation of the shale reservoir and mining the shale oil. It is important to note that we consider the shale strata as a whole, and whether it is an interlayer or a pure shale interval, all are regarded as part of the shale strata.
Structurally, Well Qiang21 is located at the edge of the Yangwuzhai Sag (see Figure 1), with a shallow burial depth interval of 2732.3-3035.0 m for the upper third number of the Shahejie Formation. This causes poor physical properties of the shale oil. As shown in Figure 12, as the depth increases, the viscosity and density of shale oil changes significantly. When the depth is greater than 3500 m, the viscosity and density decrease sharply. Relatively, the interval from 2732.3 to 3035.0 m in the Qing21 well has a high viscosity and density, which is not favorable for shale oil seepage. In addition, shale oil contains a high colloidal constituent and low saturated and aromatic hydrocarbons at the studied interval ( Figure 13). Usually, saturated and aromatic hydrocarbons have strong flowability, while colloidal is composed by heavier hydrocarbons and has poor flowability. Therefore, shale oil extraction will present a great challenge for the Qiang21 well.
( Figure 11). The shale reservoir, within this interval, has a relatively high oil content, a high mobility ratio, and a high brittleness index, which are favorable factors for shale oil extraction; however, the shale has a normal formation pressure coefficient and relatively poor physical properties, which may increase the difficulty in fracturing simulation of the shale reservoir and mining the shale oil. It is important to note that we consider the shale strata as a whole, and whether it is an interlayer or a pure shale interval, all are regarded as part of the shale strata.  Structurally, Well Qiang21 is located at the edge of the Yangwuzhai Sag (see Figure  1), with a shallow burial depth interval of 2732.3-3035.0 m for the upper third number of the Shahejie Formation. This causes poor physical properties of the shale oil. As shown in Figure 12, as the depth increases, the viscosity and density of shale oil changes significantly. When the depth is greater than 3500 m, the viscosity and density decrease sharply. Relatively, the interval from 2732.3 to 3035.0 m in the Qing21 well has a high viscosity and density, which is not favorable for shale oil seepage. In addition, shale oil contains a high colloidal constituent and low saturated and aromatic hydrocarbons at the studied interval ( Figure 13). Usually, saturated and aromatic hydrocarbons have strong flowability, while colloidal is composed by heavier hydrocarbons and has poor flowability. Therefore, shale oil extraction will present a great challenge for the Qiang21 well. Structurally, Well Qiang21 is located at the edge of the Yangwuzhai Sag (see Figure  1), with a shallow burial depth interval of 2732.3-3035.0 m for the upper third number of the Shahejie Formation. This causes poor physical properties of the shale oil. As shown in Figure 12, as the depth increases, the viscosity and density of shale oil changes significantly. When the depth is greater than 3500 m, the viscosity and density decrease sharply. Relatively, the interval from 2732.3 to 3035.0 m in the Qing21 well has a high viscosity and density, which is not favorable for shale oil seepage. In addition, shale oil contains a high colloidal constituent and low saturated and aromatic hydrocarbons at the studied interval ( Figure 13). Usually, saturated and aromatic hydrocarbons have strong flowability, while colloidal is composed by heavier hydrocarbons and has poor flowability. Therefore, shale oil extraction will present a great challenge for the Qiang21 well.

Conclusions
This paper presented a new method to quantitatively evaluate the shale oil sweet spot. A comprehensive sweet spot of shale oil is a superposition of geological and engineering sweet spots. The former can be expressed using the four parameters of oil content, movable oil ratio, pressure coefficient of reservoir fluid, and shale permeability; the latter can be reflected by a brittleness index. Based on these five primary parameters, a comprehensive weighting factor was established to quantitatively evaluate favorable shale oil reservoirs. A detailed methodology for calculating the comprehensive weighting factor and various primary parameters was provided.
The proposed method was applied to one key exploratory well (Qiang21) from Raoyang Sag of Bohai Bay basin, China. A sweet spot interval of the shale oil reservoir in the upper third member of the Shahejie Formation was determined at 2926-2955 m, where the interval has high oil content, a high mobility ratio, and a high brittleness index, which are favorable factors for shale oil extraction. However, the shale reservoir also has a normal formation pressure coefficient and relatively poor physical properties, which may increase the difficulty of fracturing the shale reservoir and mining the shale oil. Overall, the parameters of oil content and mobility ratio provide the main contributions to the comprehensive weighting factor.
The interval of 2732.3-3035.0 m in the Qiang21 well has a high viscosity, high density, high colloidal constituent, and low saturated and aromatic hydrocarbons, which are not favorable for shale oil seepage. Thus, shale oil extraction will present a great challenge for the Qiang21 well.  Data Availability Statement: The study did not report any data.

Conclusions
This paper presented a new method to quantitatively evaluate the shale oil sweet spot. A comprehensive sweet spot of shale oil is a superposition of geological and engineering sweet spots. The former can be expressed using the four parameters of oil content, movable oil ratio, pressure coefficient of reservoir fluid, and shale permeability; the latter can be reflected by a brittleness index. Based on these five primary parameters, a comprehensive weighting factor was established to quantitatively evaluate favorable shale oil reservoirs. A detailed methodology for calculating the comprehensive weighting factor and various primary parameters was provided.
The proposed method was applied to one key exploratory well (Qiang21) from Raoyang Sag of Bohai Bay basin, China. A sweet spot interval of the shale oil reservoir in the upper third member of the Shahejie Formation was determined at 2926-2955 m, where the interval has high oil content, a high mobility ratio, and a high brittleness index, which are favorable factors for shale oil extraction. However, the shale reservoir also has a normal formation pressure coefficient and relatively poor physical properties, which may increase the difficulty of fracturing the shale reservoir and mining the shale oil. Overall, the parameters of oil content and mobility ratio provide the main contributions to the comprehensive weighting factor.
The interval of 2732.3-3035.0 m in the Qiang21 well has a high viscosity, high density, high colloidal constituent, and low saturated and aromatic hydrocarbons, which are not favorable for shale oil seepage. Thus, shale oil extraction will present a great challenge for the Qiang21 well.