# Quantitative Method for Evaluating Shale Oil Resources Based on Movable Oil Content

Authors: Zhao Xianzheng; Pu Xiugang; Jin Fengming; Han Wenzhong; Shi Zhannan; Zhang Wei; Dong Xiongying; Wang Hu; Li Haodong
Journal: Geofluids
 (2021)

---

## Abstract

The light hydrocarbon content (S1) of shale oil is lost in aboveground experimental measurements, which do not accurately reflect actual underground light hydrocarbon content and cannot meet the demands of resource quantity calculation. Based on field and laboratory experimental data from the second member of the Kongdian Formation in the Cangdong Sag, Bohai Bay Basin, the retained oil and movable oil contents in shale were calculated using a mathematical formula, and the total resources were quantified. The correction coefficient of S1 from pyrolysis and the adsorption of oil by the total organic carbon (TOC) were determined to be 1.25 and 100 mg/g, respectively. The calculated parameter movable oil content (Smovable) and the corresponding calculation formula of S1 are proposed. The lower limit of S1 is 100 mg HC/g rock, the TOC content is 4 wt % and 6 wt %, and the corrected movable oil content of 3 mg/g and 6 mg/g, respectively, will be of great significance in shale oil geology and engineering. The optimal geological and engineering settings are divided into three categories and eight subcategories on the basis of these findings. Methods for evaluating total, movable, and recoverable oil resources from shale are discussed, which elucidate a new method for quantitative evaluation and ranking of shale oil resources. This approach is suitable for application in other shale oil exploration and development areas globally.


---

## Body

## 1. Introduction

Compared with North American shale oil, China’s lacustrine oil shales are quite different in terms of reservoir and fluid properties and reservoir transformation characteristics [1, 2]. Lacustrine shales in China are mainly grayish and argillaceous shales [3]. The total organic carbon (TOC) content varies greatly from 0.7% to 28%; the vitrinite reflectance is generally lower than 1.0%; the light oil content is also relatively small; the porosity and permeability are relatively poor, and the pore throat diameter is less than 300 nm; the fluid viscosity is up to 350.0 MPa·s; the density is between 0.67 and 0.94 g/cm3 [4–6]. It is of great significance to understand shale oil resources, especially movable resources, using current engineering technology [7]. There are three principal evaluation methods for shale oil resource evaluation: analogy, volumetric, and genetic methods [8]. The analogy method is used to estimate oil and gas resources by comparing it with oil and gas fields with similar geological characteristics. The volumetric method calculates the oil volume in a unit volume of rock and has higher requirements for data in the study area. The genetic method estimates hydrocarbon generation resources by simulating hydrocarbon generation methods, constrained by experimental objects and methods. S1 is the key parameter for shale oil bearing evaluation, which is consistent with the concept of oil saturation in sandstone reservoirs. S1 characterizes shale-free oil [9, 10]. Pyrolysis S1 combined with the oil saturation index OSI (S1/TOC) is the key parameter for shale oil bearing and mobility evaluation [11, 12]. Generally, the parameters are obtained from the wellhead cuttings, but the cuttings may be lost in the process of formation breaking from the bottom hole to the wellhead and placement. The saturation index should be the index of oil and gas mobility rather than the actual momentum. Methods for accurately correcting S1 parameters to reflect the quantity of free oil in the underground shale reservoir and quantitatively characterizing the free oil parameters of shale oil are required. Based on the data of closed cores and cuttings, combined with field and indoor pyrolysis experiments and swelling experiments, this study discussed the calculation of the S1 recovery coefficient, shale oil potential, and shale oil resources with the aim of establishing a suitable method for evaluating China’s lacustrine shale oil resources.


## 2. Geological Setting

The Cangdong Sag is one of several oil-rich sag areas in the Bohai Bay Basin. It is located in the southwest of the Huanghua Depression. It is sandwiched between the Cangxian uplift in the west, Xuhei bulge in the east, Kongdian bugle in the north, and Dongguan uplift in the south (Figure1). It is a continental faulted lake basin formed by regional extension during the Cenozoic. During the sedimentary period of the second member of the Kongdong Formation, the Cangdong Sag contained an inland freshwater to brackish water elliptical closed lake and the input of clastic materials from the four major provenances (Kongdian, Cangxian, Dongguan, and Xuhei uplifts) around the basin, the sedimentary facies showed regular changes from the edge of the lake basin to the middle of the lake basin, and the outer ring comprises extensive braided river delta deposits (yellow-orange colors in Figure 1) [13, 14]. The middle ring comprises a dense zone composed of silty mudstone, siltstone, and dolomite, and the inner ring consists of fine-grained sediments composed of thick dark shale with thin layers of siltstone and dolomite. These sediments are mainly fine-grained sedimentary rocks (>50%) with particle sizes less than 0.0625 mm and are rich in preserved organic matter, which is the source of oil in the shale (Figure 1).Figure 1
Structural subdivision and sedimentary facies map of the Kong 2 member in Cangdong Sag, Bohai Bay Basin.


## 3. Samples and Experiments

All samples were obtained from the Kong 2 member in the Cangdong Sag. To determine the amount of oil adsorbed by kerogens with different maturities, 12 samples were taken from the core of G19 well for the swelling adsorption experiment, which is described below. To determine the rate of loss of light hydrocarbons, 46 pieces of cuttings from the G25 well closed coring section and corresponding depth wellhead cuttings were subjected to thermal extraction. Experiments were conducted to analyze the distribution of light hydrocarbon content and calculate the movable light hydrocarbon content. A total of 808 cutting samples were taken, including 192 at 3939 m depth in horizontal well GY1-3-1, 274 at 3714 m depth in horizontal well GY2-1-1, 207 at 3702 m depth in horizontal well GY2-1-2, and 135 at 3777 m depth in horizontal well GY2-1-4 (Table1; wells are marked in Figure 1).Table 1
Experimental sample details.
WellSampleExperimentTypeNumberNameResultG19Core12Swelling adsorptionAdsorbed hydrocarbonG25Core23TOC, thermal extraction and pyrolysisTOC,S1G25Cuttings23TOC, thermal extraction and pyrolysisTOC,S1GY1-3-1Cuttings192TOC, thermal extraction and pyrolysisTOC,S1GY2-1-1Cuttings274TOC, thermal extraction and pyrolysisTOC,S1GY2-1-2Cuttings207TOC, thermal extraction and pyrolysisTOC,S1GY2-1-4Cuttings135TOC, thermal extraction and pyrolysisTOC,S1


### 3.1. Swelling Adsorption Experiment

The inner diameter of the screen tube was 60 mm, and the height of the screen tube was 60 mm. In the experimental process, the sample was crushed and ground to 120 mesh and dried in an oven; 200 mg of sample was placed into the test tube and centrifuged for 30 min at 8000 rpm, and the height (H1) of the sample was measured. The glass tube with a sieve hole at the lower end was placed into a centrifuge tube containing an organic solvent, and the sample was centrifuged for 30 min. The kerogen was soaked at 30°C for 24 h and then centrifuged under the same conditions. The height after swelling (H2) was measured, and the swelling ratio (Qv=H2/H1) was calculated. The total hydrocarbon retention and swelling ratio of kerogen can be converted according to the following formula [15]:
(1)M=Qv−1.0ρ1ρ2,where M is the amount of hydrocarbon retained in unit kerogen (g/g); Qv is the swelling ratio; ρ1 is the weighted average density of n-tetradecane, o-xylene, acetic acid, isopropanol, and ethanol retained in kerogen (0.9 g/cm3; and ρ2 is the density of kerogen (1.28–1.35 g/cm3), corresponding to a vitrinite reflectance (Ro) 0.5–1.3%.


### 3.2. TOC, Thermal Extraction, and Pyrolysis

The TOC values were measured using a CS-600 organic carbon analyzer, following the Chinese National Standards GB/T 19145–2003 [3] after the removal of carbonate minerals with acid.A Vinci Technologies Rock-Eval 6 instrument was used for thermal extraction (S1) and pyrolysis (S2) following a procedure similar to that of Lin et al. [16] and Hood et al. [17]. For S1, the oven was kept isothermally at 300°C, and for S2, the oven temperature was increased to 550°C at a rate of 25°C/min. The experiment was completed at the Wuxi Institute of Petroleum Geology, Research Institute of Petroleum Exploration and Development, Sinopec. The steps are as follows: Soxhlet extraction to remove soluble organic matter from source rocks, acid dissolution to remove inorganic minerals, ultrasound for 30 min to fully saturate adsorption, and high-speed centrifugation for 30 min to remove the nonadsorbed part and yield the saturated adsorption capacity.


## 3.1. Swelling Adsorption Experiment

The inner diameter of the screen tube was 60 mm, and the height of the screen tube was 60 mm. In the experimental process, the sample was crushed and ground to 120 mesh and dried in an oven; 200 mg of sample was placed into the test tube and centrifuged for 30 min at 8000 rpm, and the height (H1) of the sample was measured. The glass tube with a sieve hole at the lower end was placed into a centrifuge tube containing an organic solvent, and the sample was centrifuged for 30 min. The kerogen was soaked at 30°C for 24 h and then centrifuged under the same conditions. The height after swelling (H2) was measured, and the swelling ratio (Qv=H2/H1) was calculated. The total hydrocarbon retention and swelling ratio of kerogen can be converted according to the following formula [15]:
(1)M=Qv−1.0ρ1ρ2,where M is the amount of hydrocarbon retained in unit kerogen (g/g); Qv is the swelling ratio; ρ1 is the weighted average density of n-tetradecane, o-xylene, acetic acid, isopropanol, and ethanol retained in kerogen (0.9 g/cm3; and ρ2 is the density of kerogen (1.28–1.35 g/cm3), corresponding to a vitrinite reflectance (Ro) 0.5–1.3%.


## 3.2. TOC, Thermal Extraction, and Pyrolysis

The TOC values were measured using a CS-600 organic carbon analyzer, following the Chinese National Standards GB/T 19145–2003 [3] after the removal of carbonate minerals with acid.A Vinci Technologies Rock-Eval 6 instrument was used for thermal extraction (S1) and pyrolysis (S2) following a procedure similar to that of Lin et al. [16] and Hood et al. [17]. For S1, the oven was kept isothermally at 300°C, and for S2, the oven temperature was increased to 550°C at a rate of 25°C/min. The experiment was completed at the Wuxi Institute of Petroleum Geology, Research Institute of Petroleum Exploration and Development, Sinopec. The steps are as follows: Soxhlet extraction to remove soluble organic matter from source rocks, acid dissolution to remove inorganic minerals, ultrasound for 30 min to fully saturate adsorption, and high-speed centrifugation for 30 min to remove the nonadsorbed part and yield the saturated adsorption capacity.


## 4. Results

### 4.1. Thermal Extraction and Pyrolysis

There were 46 samples in total; the TOC values were evenly distributed, the minimum value was 0.67% wt, the maximum value was 7.29% wt, and the average value was 3.49% wt. The minimum pyrolysisS1 value measured by closed coring well samples was 0.34 mg HC/g rock, and the maximum pyrolysis S1 value was 6.84 mg HC/g rock, with an average of 1.19 mg HC/g rock. Corresponding to the depth of cutting samples, the minimum pyrolysis S1 value is 0.30 mg HC/g rock, and the maximum pyrolysis S1 value is 4.89 mg HC/g rock, with an average of 0.91 mg HC/g rock. The S1 loss value of the sample was 0.04 mg HC/g rock, the maximum was 2.14 mg HC/g rock, and the average was 0.39 mg HC/g rock. The minimum loss rate was 11.3%, the maximum was 52.7%, and the average loss rate was 29.4% (Table 2).Table 2
Original hydrocarbon content, loss amount, and loss rate with conventional pyrolysis of closed core and wellhead cuttings.
No.123456789101112WellG25G25G25G25G25G25G25G25G25G25G25G25Depth (m)3365.93353.73359.33355.73351.03342.53376.43339.53347.03346.03344.53352.7TOC (%)1.020.671.712.281.842.982.405.663.252.062.463.84Core (mg/g)0.340.600.840.911.141.211.221.241.241.631.671.79Cutting (mg/g)0.300.380.570.680.950.681.080.920.870.921.271.03Loss (mg/g)0.040.220.270.230.190.530.140.320.370.710.400.76Loss ratio (%)11.835.932.625.216.044.111.326.029.843.624.142.3N.O1314151617181920212223AverageWellG25G25G25G25G25G25G25G25G25G25G25Depth (m)3374.43367.93369.93351.63377.43375.43363.33366.93361.33356.73345.0TOC (%)4.691.853.297.292.963.463.534.666.565.975.72Core (mg/g)1.871.922.732.983.123.383.563.754.484.526.842.30Cutting (mg/g)1.020.912.332.431.772.312.673.163.972.384.891.63Loss (mg/g)0.851.011.401.551.351.071.890.590.512.141.950.67Loss ratio (%)45.452.714.618.443.331.724.915.911.347.328.529.4


### 4.2. TOC andS1∗ Features

In general, the TOC value was positively correlated withS1. When TOC is between 1 and 3 wt %, the high-value part of S1 increases rapidly with the increase in TOC. When TOC is greater than 3 wt %, the high-value part of S1 changes or disappears in a horizontal shape, and there is an obvious upper envelope. As TOC increases, S1 increases. When TOC is less than 1% wt, S1 does not exceed 4 mg/g. When TOC is greater than 3% wt, S1 approaches the maximum value, generally 12–14 mg/g, and GY2-1-2 reaches 20 mg/g (Figure 2).Figure 2
Distribution chart of TOC andS1∗ characteristics in the horizontal section of key wells: (a) well GY2-1-1, equivalent to a vertical depth of 3714 m; (b) well GY1-3-1, equivalent to a vertical depth of 3939 m; (c) well GY2-1-2, equivalent to a vertical depth of 3702 m; (d) well GY2-1-4, equivalent to a vertical depth of 3777 m. Samples were obtained from horizontal well cuttings of the second member, Kongdian Formation.
(a)(b)(c)(d)


### 4.3. Adsorbed Hydrocarbon Content

The content of light hydrocarbons adsorbed by low-maturity shale was high. With an increase in maturity, the amount of hydrocarbon adsorbed gradually decreases. At the low maturity stage (Ro less than 0.6%), the adsorption capacity is more than 150 mg HC/g TOC; at the mature stage (Ro>0.6%), the amount of adsorbed hydrocarbon is 82–155/109 mg HC/g TOC. When Ro is greater than 0.8% (corresponding to a depth of 3300 m), the average amount of adsorbed hydrocarbon was 95 mg HC/g TOC (Figure 3).Figure 3
Quantitative relationship between shale maturity and hydrocarbon adsorption determined by swelling test.


## 4.1. Thermal Extraction and Pyrolysis

There were 46 samples in total; the TOC values were evenly distributed, the minimum value was 0.67% wt, the maximum value was 7.29% wt, and the average value was 3.49% wt. The minimum pyrolysisS1 value measured by closed coring well samples was 0.34 mg HC/g rock, and the maximum pyrolysis S1 value was 6.84 mg HC/g rock, with an average of 1.19 mg HC/g rock. Corresponding to the depth of cutting samples, the minimum pyrolysis S1 value is 0.30 mg HC/g rock, and the maximum pyrolysis S1 value is 4.89 mg HC/g rock, with an average of 0.91 mg HC/g rock. The S1 loss value of the sample was 0.04 mg HC/g rock, the maximum was 2.14 mg HC/g rock, and the average was 0.39 mg HC/g rock. The minimum loss rate was 11.3%, the maximum was 52.7%, and the average loss rate was 29.4% (Table 2).Table 2
Original hydrocarbon content, loss amount, and loss rate with conventional pyrolysis of closed core and wellhead cuttings.
No.123456789101112WellG25G25G25G25G25G25G25G25G25G25G25G25Depth (m)3365.93353.73359.33355.73351.03342.53376.43339.53347.03346.03344.53352.7TOC (%)1.020.671.712.281.842.982.405.663.252.062.463.84Core (mg/g)0.340.600.840.911.141.211.221.241.241.631.671.79Cutting (mg/g)0.300.380.570.680.950.681.080.920.870.921.271.03Loss (mg/g)0.040.220.270.230.190.530.140.320.370.710.400.76Loss ratio (%)11.835.932.625.216.044.111.326.029.843.624.142.3N.O1314151617181920212223AverageWellG25G25G25G25G25G25G25G25G25G25G25Depth (m)3374.43367.93369.93351.63377.43375.43363.33366.93361.33356.73345.0TOC (%)4.691.853.297.292.963.463.534.666.565.975.72Core (mg/g)1.871.922.732.983.123.383.563.754.484.526.842.30Cutting (mg/g)1.020.912.332.431.772.312.673.163.972.384.891.63Loss (mg/g)0.851.011.401.551.351.071.890.590.512.141.950.67Loss ratio (%)45.452.714.618.443.331.724.915.911.347.328.529.4


## 4.2. TOC andS1∗ Features

In general, the TOC value was positively correlated withS1. When TOC is between 1 and 3 wt %, the high-value part of S1 increases rapidly with the increase in TOC. When TOC is greater than 3 wt %, the high-value part of S1 changes or disappears in a horizontal shape, and there is an obvious upper envelope. As TOC increases, S1 increases. When TOC is less than 1% wt, S1 does not exceed 4 mg/g. When TOC is greater than 3% wt, S1 approaches the maximum value, generally 12–14 mg/g, and GY2-1-2 reaches 20 mg/g (Figure 2).Figure 2
Distribution chart of TOC andS1∗ characteristics in the horizontal section of key wells: (a) well GY2-1-1, equivalent to a vertical depth of 3714 m; (b) well GY1-3-1, equivalent to a vertical depth of 3939 m; (c) well GY2-1-2, equivalent to a vertical depth of 3702 m; (d) well GY2-1-4, equivalent to a vertical depth of 3777 m. Samples were obtained from horizontal well cuttings of the second member, Kongdian Formation.
(a)(b)(c)(d)


## 4.3. Adsorbed Hydrocarbon Content

The content of light hydrocarbons adsorbed by low-maturity shale was high. With an increase in maturity, the amount of hydrocarbon adsorbed gradually decreases. At the low maturity stage (Ro less than 0.6%), the adsorption capacity is more than 150 mg HC/g TOC; at the mature stage (Ro>0.6%), the amount of adsorbed hydrocarbon is 82–155/109 mg HC/g TOC. When Ro is greater than 0.8% (corresponding to a depth of 3300 m), the average amount of adsorbed hydrocarbon was 95 mg HC/g TOC (Figure 3).Figure 3
Quantitative relationship between shale maturity and hydrocarbon adsorption determined by swelling test.


## 5. Discussion

### 5.1. Correction Coefficient of Light Hydrocarbon Content

The pyrolysis parameterS1 is the key parameter for characterizing the oil content of shale oil. Shale oil resources may change by 2–3 times before and after the correction of light hydrocarbon loss [18, 19]. Currently, most correction coefficients are obtained by indoor pyrolysis and gold tube experiments (Li et al. 2012). However, owing to the influence of core storage conditions, experimental test and analysis technology, and swelling effects, there is a large error between the number of S1 recovery and the actual underground value [12]. According to the analysis of chromatographic results of crude oils with different properties under different storage conditions, Cooler et al. [20] and Hunt et al. [21] suggested that light hydrocarbons account for nearly 1/3 of the total oil, and almost all are lost. Zhang et al. [8] compared the results of direct pyrolysis with those of pyrolysis after different storage times; it was concluded that the longer the storage time, the more light hydrocarbon was lost. The light hydrocarbon recovery was based on 50% loss. These experiments provide a lot of good references for us to understand light hydrocarbon recovery, but there are always differences between field application and experimental analysis to provide the applicability of the method. In this study, the recovery coefficient of light hydrocarbon content was determined by comparing the pyrolysis S1 value of wellhead cuttings with that of closed coring at the same depth. The correction coefficient was 1.42, based on the loss of pyrolysis S1 value of 29.4%. It is worth noting that the coefficient is closely related to the maturity of organic matter, which is suitable for this study area, but corresponding experiments should be conducted in other maturity sections or other study areas.


### 5.2. Calculation Method of Movable Hydrocarbon

At present, the commonly used characterization index of the movable hydrocarbon content in shale oil is the oil saturation index (S1/TOC). When the saturation index is greater than 100 mg/g, the movable hydrocarbon content is greater, and the larger the value, the greater the movable hydrocarbon content [12, 22]. According to S1/TOC, the index reflects the movable efficiency rather than the actual momentum. To quantitatively analyze the real amount of mobile hydrocarbons, a method for subtracting the adsorption amount from the total amount of light hydrocarbon is proposed. When the maturity is 0.5–1.0%, the amount of hydrocarbon adsorbed by shale in the study area is 95 mg/g (Figure 4(a)), the total hydrocarbon amount is the experimental data, and the absolute amount of real movable hydrocarbon is the distance between the total hydrocarbon amount and the adsorption baseline (Figure 4(a)).Figure 4
Characterization and formula derivation of movable hydrocarbon in shale oil: (a) schematic diagram of TOC andS1 outer envelope in horizontal section of key well; (b) parameter diagram of shale movable hydrocarbon calculation formula.
(a)(b)The derivation process is as follows.M is the value of TOC and S1∗ at any point, the straight line S1∗=a∗TOC is the lower limit of the movable oil, and a is the coefficient of organic matter adsorption movable oil content. The movable oil content can be expressed as the distance from any point m to the straight line S1∗=a∗TOC (Figure 4(b)).
(2)D=S1∗−A∗TOC12+A2=Smovable12+A2,(3)cosθ=112+A2=DL.Substitute Equation (3) into (2):
(4)Smoveable=L.


### 5.3. Shale Oil Resource Evaluation

#### 5.3.1. Resource Classification

The normal development of shale oil resources requires fracturing in addition to good geological conditions [23]. High-quality geological conditions mainly include good physical properties and oil-bearing properties [24]. Currently, the average porosity of shale is 3–8%, and the permeability is lower than 0.1 mDa. The difference in oil-bearing properties is key to the success of shale oil wells. Zhao et al. (2020) considered that the probability of high yield is high when S1 is greater than 3 mg/g, and S1∗ increases linearly with TOC when TOC is less than 3% wt. When the TOC is approximately 3% wt, the linear distance from the corner of the outer loop to the lower adsorption limit was the largest, indicating that the mobile hydrocarbon content was large. Brittle mineral content is an important reference index for fracturing. TOC is a plastic mineral, and its content has a significant influence on brittleness (Figure 5). When the TOC content was 4 wt %, the content of brittle minerals was only 70%. When the TOC value reached 6 wt %, the brittle mineral content was close to 60%. Furthermore, the brittleness of shale with low organic matter (TOC<4 wt %) is the best for S1∗, followed by that of shale with high organic matter (TOC>4 wt %). Concurrently, when TOC is greater than 4 wt %, S1∗ decreases with an increase in TOC.Figure 5
Effect of TOC on the content of brittle minerals.To summarize, the geological and engineering evaluation template of shale oil was established based on these findings. The lower limit of movable oil is 100 mg HC/g rock; TOC content of 4% wt and 6% wt and corrected movable oil content of 3 mg/g and 6 mg/g divide the optimal geological and engineering settings of shale oil into three categories and eight subcategories (Figure6).Figure 6
Shale oil resource evaluation template.


#### 5.3.2. Calculation Method

The specific steps of quantitative calculation of shale oil resources are as follows: (1) obtain the correction coefficient of movable oil content (S1) using an experimental method, obtain the movable oil content (S1∗) close to the formation state, calculate the total amount of retained oil in shale, and evaluate the total shale oil resources; (2) determine the recoverable amount of shale oil by combining the correction of the pyrolysis movable oil content (S1∗) and residual TOC content. The key parameter can be used to calculate the geological movable resources of shale oil; and (3) establish an evaluation template for shale oil geology and engineering to calculate the recoverable resources of shale oil.(1) Total Light Hydrocarbon Resources of Shale Oil. The correction coefficient of S1 was obtained by physical experimentation.The free oil loss rate was determined using the free oil content of fresh samples from sealed coring and cutting samples from the corresponding depth of the wellhead, and the free oil correction coefficient (Kc) was determined.
(5)Qtotal=0.1×Sshale×hshale×ρ×S1×Kc.Qtotal is shale-free oil resource (104 t); Sshale is total shale area (km2); hshale is total shale effective thickness (m); ρ is shale density (2.5 g/cm3); S1 is light hydrocarbon content of residual oil (mg HC/g rock); and Kc is the dimensionless correction coefficient of shale free oil loss.(2) Movable Shale Oil Resources.
(1)
The movable oil content (S1∗) can be obtained by using the correction coefficient to correct the residual oil(2)
Based on the experimental data of kerogen swelling adsorption, the content of residual TOC adsorbed oil and amount of adsorbed oil per unit organic matter (Kabs) were determined(3)
Calculate the movable oil content (Smovable) with Equation (5)(6)Smovable=S1∗−Kabs×TOC,S1∗=S1×Kc,where Smovable is the amount of residual oil (mg HC/g rock); S1∗ is the moveable oil content after correction, unit: mg HC/g rock; S1 is the content of conventional thermally extractable oil, unit: 1 mg HC/g rock; and Kabs is the adsorption capacity of unit organic matter, mg HC/g TOC.
(7)Qmovable=0.1×Sshale×hshale×ρ×S1∗,where Qmoveable is the movable resources of shale-free oil (104 t); Sshale is total shale area (km2); hshale is total shale effective thickness (m); ρ is shale density (2.5 g/cm3); and S1∗ is the content of conventional pyrolysis movable oil (mg HC/g rock).(3) Recoverable Resources.
(1)
The TOC content andS1∗ chart of corrected movable oil content were established to quickly evaluate optimal geological settings(2)
Combined with the brittleness index template, the optimal engineering template settings can be evaluated(3)
Determine the lower limit of recoverable resources and determine the range of recoverable resourcesThe template was used to classify the study area horizontally and vertically to calculate the recoverable shale oil resources. The formula for calculating the shale-oil movable resources is as follows:(8)Qrecoverable=0.1×Ssweet×hsweet×ρ×Smovable,where Qrecoverable is shale-free oil resource (104 t); Ssweet is shale sweet spot area (km2); hsweet is shale sweet spot thickness (m); ρ is shale density (2.5 g/cm3); and Smovable is the amount of residual movable oil (mg HC/g rock).Vertically, single wells are classified by TOC content, residual moveable oil, and the brittleness index. The distribution characteristics of the plane TOC content, residual oil moveable, and brittleness index in the study area were determined by seismic inversion. Superposition classification was conducted, and shale oil moveable resources were calculated using Equation (8) (Table 3).Table 3
Shale oil movable resources of XX layer in XX block.
TypeThickness (m)Favorable area (km2)Smovable (kg oil/t rock)Resource (104 t)II15023.610.23009.0I24624.28.42337.7IIII13534.54.81449.0II22515.64.2409.5IIIIII14552.42.51473.8III25039.82.51243.8Total9922.7Average thickness, favorable area, and measurable parameters are calculated by the weighted average method.


## 5.1. Correction Coefficient of Light Hydrocarbon Content

The pyrolysis parameterS1 is the key parameter for characterizing the oil content of shale oil. Shale oil resources may change by 2–3 times before and after the correction of light hydrocarbon loss [18, 19]. Currently, most correction coefficients are obtained by indoor pyrolysis and gold tube experiments (Li et al. 2012). However, owing to the influence of core storage conditions, experimental test and analysis technology, and swelling effects, there is a large error between the number of S1 recovery and the actual underground value [12]. According to the analysis of chromatographic results of crude oils with different properties under different storage conditions, Cooler et al. [20] and Hunt et al. [21] suggested that light hydrocarbons account for nearly 1/3 of the total oil, and almost all are lost. Zhang et al. [8] compared the results of direct pyrolysis with those of pyrolysis after different storage times; it was concluded that the longer the storage time, the more light hydrocarbon was lost. The light hydrocarbon recovery was based on 50% loss. These experiments provide a lot of good references for us to understand light hydrocarbon recovery, but there are always differences between field application and experimental analysis to provide the applicability of the method. In this study, the recovery coefficient of light hydrocarbon content was determined by comparing the pyrolysis S1 value of wellhead cuttings with that of closed coring at the same depth. The correction coefficient was 1.42, based on the loss of pyrolysis S1 value of 29.4%. It is worth noting that the coefficient is closely related to the maturity of organic matter, which is suitable for this study area, but corresponding experiments should be conducted in other maturity sections or other study areas.


## 5.2. Calculation Method of Movable Hydrocarbon

At present, the commonly used characterization index of the movable hydrocarbon content in shale oil is the oil saturation index (S1/TOC). When the saturation index is greater than 100 mg/g, the movable hydrocarbon content is greater, and the larger the value, the greater the movable hydrocarbon content [12, 22]. According to S1/TOC, the index reflects the movable efficiency rather than the actual momentum. To quantitatively analyze the real amount of mobile hydrocarbons, a method for subtracting the adsorption amount from the total amount of light hydrocarbon is proposed. When the maturity is 0.5–1.0%, the amount of hydrocarbon adsorbed by shale in the study area is 95 mg/g (Figure 4(a)), the total hydrocarbon amount is the experimental data, and the absolute amount of real movable hydrocarbon is the distance between the total hydrocarbon amount and the adsorption baseline (Figure 4(a)).Figure 4
Characterization and formula derivation of movable hydrocarbon in shale oil: (a) schematic diagram of TOC andS1 outer envelope in horizontal section of key well; (b) parameter diagram of shale movable hydrocarbon calculation formula.
(a)(b)The derivation process is as follows.M is the value of TOC and S1∗ at any point, the straight line S1∗=a∗TOC is the lower limit of the movable oil, and a is the coefficient of organic matter adsorption movable oil content. The movable oil content can be expressed as the distance from any point m to the straight line S1∗=a∗TOC (Figure 4(b)).
(2)D=S1∗−A∗TOC12+A2=Smovable12+A2,(3)cosθ=112+A2=DL.Substitute Equation (3) into (2):
(4)Smoveable=L.


## 5.3. Shale Oil Resource Evaluation

### 5.3.1. Resource Classification

The normal development of shale oil resources requires fracturing in addition to good geological conditions [23]. High-quality geological conditions mainly include good physical properties and oil-bearing properties [24]. Currently, the average porosity of shale is 3–8%, and the permeability is lower than 0.1 mDa. The difference in oil-bearing properties is key to the success of shale oil wells. Zhao et al. (2020) considered that the probability of high yield is high when S1 is greater than 3 mg/g, and S1∗ increases linearly with TOC when TOC is less than 3% wt. When the TOC is approximately 3% wt, the linear distance from the corner of the outer loop to the lower adsorption limit was the largest, indicating that the mobile hydrocarbon content was large. Brittle mineral content is an important reference index for fracturing. TOC is a plastic mineral, and its content has a significant influence on brittleness (Figure 5). When the TOC content was 4 wt %, the content of brittle minerals was only 70%. When the TOC value reached 6 wt %, the brittle mineral content was close to 60%. Furthermore, the brittleness of shale with low organic matter (TOC<4 wt %) is the best for S1∗, followed by that of shale with high organic matter (TOC>4 wt %). Concurrently, when TOC is greater than 4 wt %, S1∗ decreases with an increase in TOC.Figure 5
Effect of TOC on the content of brittle minerals.To summarize, the geological and engineering evaluation template of shale oil was established based on these findings. The lower limit of movable oil is 100 mg HC/g rock; TOC content of 4% wt and 6% wt and corrected movable oil content of 3 mg/g and 6 mg/g divide the optimal geological and engineering settings of shale oil into three categories and eight subcategories (Figure6).Figure 6
Shale oil resource evaluation template.


### 5.3.2. Calculation Method

The specific steps of quantitative calculation of shale oil resources are as follows: (1) obtain the correction coefficient of movable oil content (S1) using an experimental method, obtain the movable oil content (S1∗) close to the formation state, calculate the total amount of retained oil in shale, and evaluate the total shale oil resources; (2) determine the recoverable amount of shale oil by combining the correction of the pyrolysis movable oil content (S1∗) and residual TOC content. The key parameter can be used to calculate the geological movable resources of shale oil; and (3) establish an evaluation template for shale oil geology and engineering to calculate the recoverable resources of shale oil.(1) Total Light Hydrocarbon Resources of Shale Oil. The correction coefficient of S1 was obtained by physical experimentation.The free oil loss rate was determined using the free oil content of fresh samples from sealed coring and cutting samples from the corresponding depth of the wellhead, and the free oil correction coefficient (Kc) was determined.
(5)Qtotal=0.1×Sshale×hshale×ρ×S1×Kc.Qtotal is shale-free oil resource (104 t); Sshale is total shale area (km2); hshale is total shale effective thickness (m); ρ is shale density (2.5 g/cm3); S1 is light hydrocarbon content of residual oil (mg HC/g rock); and Kc is the dimensionless correction coefficient of shale free oil loss.(2) Movable Shale Oil Resources.
(1)
The movable oil content (S1∗) can be obtained by using the correction coefficient to correct the residual oil(2)
Based on the experimental data of kerogen swelling adsorption, the content of residual TOC adsorbed oil and amount of adsorbed oil per unit organic matter (Kabs) were determined(3)
Calculate the movable oil content (Smovable) with Equation (5)(6)Smovable=S1∗−Kabs×TOC,S1∗=S1×Kc,where Smovable is the amount of residual oil (mg HC/g rock); S1∗ is the moveable oil content after correction, unit: mg HC/g rock; S1 is the content of conventional thermally extractable oil, unit: 1 mg HC/g rock; and Kabs is the adsorption capacity of unit organic matter, mg HC/g TOC.
(7)Qmovable=0.1×Sshale×hshale×ρ×S1∗,where Qmoveable is the movable resources of shale-free oil (104 t); Sshale is total shale area (km2); hshale is total shale effective thickness (m); ρ is shale density (2.5 g/cm3); and S1∗ is the content of conventional pyrolysis movable oil (mg HC/g rock).(3) Recoverable Resources.
(1)
The TOC content andS1∗ chart of corrected movable oil content were established to quickly evaluate optimal geological settings(2)
Combined with the brittleness index template, the optimal engineering template settings can be evaluated(3)
Determine the lower limit of recoverable resources and determine the range of recoverable resourcesThe template was used to classify the study area horizontally and vertically to calculate the recoverable shale oil resources. The formula for calculating the shale-oil movable resources is as follows:(8)Qrecoverable=0.1×Ssweet×hsweet×ρ×Smovable,where Qrecoverable is shale-free oil resource (104 t); Ssweet is shale sweet spot area (km2); hsweet is shale sweet spot thickness (m); ρ is shale density (2.5 g/cm3); and Smovable is the amount of residual movable oil (mg HC/g rock).Vertically, single wells are classified by TOC content, residual moveable oil, and the brittleness index. The distribution characteristics of the plane TOC content, residual oil moveable, and brittleness index in the study area were determined by seismic inversion. Superposition classification was conducted, and shale oil moveable resources were calculated using Equation (8) (Table 3).Table 3
Shale oil movable resources of XX layer in XX block.
TypeThickness (m)Favorable area (km2)Smovable (kg oil/t rock)Resource (104 t)II15023.610.23009.0I24624.28.42337.7IIII13534.54.81449.0II22515.64.2409.5IIIIII14552.42.51473.8III25039.82.51243.8Total9922.7Average thickness, favorable area, and measurable parameters are calculated by the weighted average method.


## 5.3.1. Resource Classification

The normal development of shale oil resources requires fracturing in addition to good geological conditions [23]. High-quality geological conditions mainly include good physical properties and oil-bearing properties [24]. Currently, the average porosity of shale is 3–8%, and the permeability is lower than 0.1 mDa. The difference in oil-bearing properties is key to the success of shale oil wells. Zhao et al. (2020) considered that the probability of high yield is high when S1 is greater than 3 mg/g, and S1∗ increases linearly with TOC when TOC is less than 3% wt. When the TOC is approximately 3% wt, the linear distance from the corner of the outer loop to the lower adsorption limit was the largest, indicating that the mobile hydrocarbon content was large. Brittle mineral content is an important reference index for fracturing. TOC is a plastic mineral, and its content has a significant influence on brittleness (Figure 5). When the TOC content was 4 wt %, the content of brittle minerals was only 70%. When the TOC value reached 6 wt %, the brittle mineral content was close to 60%. Furthermore, the brittleness of shale with low organic matter (TOC<4 wt %) is the best for S1∗, followed by that of shale with high organic matter (TOC>4 wt %). Concurrently, when TOC is greater than 4 wt %, S1∗ decreases with an increase in TOC.Figure 5
Effect of TOC on the content of brittle minerals.To summarize, the geological and engineering evaluation template of shale oil was established based on these findings. The lower limit of movable oil is 100 mg HC/g rock; TOC content of 4% wt and 6% wt and corrected movable oil content of 3 mg/g and 6 mg/g divide the optimal geological and engineering settings of shale oil into three categories and eight subcategories (Figure6).Figure 6
Shale oil resource evaluation template.


## 5.3.2. Calculation Method

The specific steps of quantitative calculation of shale oil resources are as follows: (1) obtain the correction coefficient of movable oil content (S1) using an experimental method, obtain the movable oil content (S1∗) close to the formation state, calculate the total amount of retained oil in shale, and evaluate the total shale oil resources; (2) determine the recoverable amount of shale oil by combining the correction of the pyrolysis movable oil content (S1∗) and residual TOC content. The key parameter can be used to calculate the geological movable resources of shale oil; and (3) establish an evaluation template for shale oil geology and engineering to calculate the recoverable resources of shale oil.(1) Total Light Hydrocarbon Resources of Shale Oil. The correction coefficient of S1 was obtained by physical experimentation.The free oil loss rate was determined using the free oil content of fresh samples from sealed coring and cutting samples from the corresponding depth of the wellhead, and the free oil correction coefficient (Kc) was determined.
(5)Qtotal=0.1×Sshale×hshale×ρ×S1×Kc.Qtotal is shale-free oil resource (104 t); Sshale is total shale area (km2); hshale is total shale effective thickness (m); ρ is shale density (2.5 g/cm3); S1 is light hydrocarbon content of residual oil (mg HC/g rock); and Kc is the dimensionless correction coefficient of shale free oil loss.(2) Movable Shale Oil Resources.
(1)
The movable oil content (S1∗) can be obtained by using the correction coefficient to correct the residual oil(2)
Based on the experimental data of kerogen swelling adsorption, the content of residual TOC adsorbed oil and amount of adsorbed oil per unit organic matter (Kabs) were determined(3)
Calculate the movable oil content (Smovable) with Equation (5)(6)Smovable=S1∗−Kabs×TOC,S1∗=S1×Kc,where Smovable is the amount of residual oil (mg HC/g rock); S1∗ is the moveable oil content after correction, unit: mg HC/g rock; S1 is the content of conventional thermally extractable oil, unit: 1 mg HC/g rock; and Kabs is the adsorption capacity of unit organic matter, mg HC/g TOC.
(7)Qmovable=0.1×Sshale×hshale×ρ×S1∗,where Qmoveable is the movable resources of shale-free oil (104 t); Sshale is total shale area (km2); hshale is total shale effective thickness (m); ρ is shale density (2.5 g/cm3); and S1∗ is the content of conventional pyrolysis movable oil (mg HC/g rock).(3) Recoverable Resources.
(1)
The TOC content andS1∗ chart of corrected movable oil content were established to quickly evaluate optimal geological settings(2)
Combined with the brittleness index template, the optimal engineering template settings can be evaluated(3)
Determine the lower limit of recoverable resources and determine the range of recoverable resourcesThe template was used to classify the study area horizontally and vertically to calculate the recoverable shale oil resources. The formula for calculating the shale-oil movable resources is as follows:(8)Qrecoverable=0.1×Ssweet×hsweet×ρ×Smovable,where Qrecoverable is shale-free oil resource (104 t); Ssweet is shale sweet spot area (km2); hsweet is shale sweet spot thickness (m); ρ is shale density (2.5 g/cm3); and Smovable is the amount of residual movable oil (mg HC/g rock).Vertically, single wells are classified by TOC content, residual moveable oil, and the brittleness index. The distribution characteristics of the plane TOC content, residual oil moveable, and brittleness index in the study area were determined by seismic inversion. Superposition classification was conducted, and shale oil moveable resources were calculated using Equation (8) (Table 3).Table 3
Shale oil movable resources of XX layer in XX block.
TypeThickness (m)Favorable area (km2)Smovable (kg oil/t rock)Resource (104 t)II15023.610.23009.0I24624.28.42337.7IIII13534.54.81449.0II22515.64.2409.5IIIIII14552.42.51473.8III25039.82.51243.8Total9922.7Average thickness, favorable area, and measurable parameters are calculated by the weighted average method.


## 6. Conclusion

The light hydrocarbon recovery coefficient of shale obtained by the underground core and wellhead cutting method is more practical than other conventional methods. Affected by the TOC content and thermal maturity of organic matter, the recovery coefficient of light hydrocarbons varies in different areas and depths.The shale oil moveable index is the index of shale oil moveable efficiency, and the moveable oil content index can be used to quantitatively characterize shale oil movability.The corrected light hydrocarbon content was used to calculate total shale oil resources, while the movable oil content index was used to estimate movable shale oil resources. There are abundant shale oil resources in the continental basins of China. The Cangdong Sag, Bohai Bay Basin, has been distributed in 100 million tons of shale oil reserves, which requires further research and methodological development in the field of geology and engineering.



---

*Source: 5554880-2021-05-18.xml*