Moliqing basin is located at the north east of China, between Changchun and Jilin. It is part of three geologic subdivisions that form the Yitong Graben basin, the other two being Chaluhe and Luxiang basins  .
The basin has been subjected to active exploration by the Chinese Jilin Oil Company since 1985 and more than 50 wells have been drilled with discovered certifying reserves of 3317 × 104 tons oil and 13.29 × 108 m3 gas.
The studies  and  identified four cenozoic tectonic evolution divisions and the further three uplift events which resulted in major regional unconformities in the basin. Additionally, the boundary faults experience a process of right lateral extrusion and uplift  .
Thermal history and petroleum generation of the basin was conducted  and the source rocks have also been studied and discussed, with the Shuangyang formation as the principal source of the basin  among others.
Nevertheless, measured vitrinite reflectance considered the major thermal indicator was not used for the calibration because of the questionable measured values. Instead, the basin temperature and the pyrolysis Tmax were used as thermal maturity indicators.
However, in this study, the Tmax based calculated Vitrinite reflectance was used in addition to the BHT (Bottom Hole Temperature) and the pyrolysis Tmax so as to provide an additional verification method to best describe the basin thermal conditions.
2. Geological Setting
The Yitong graben is a petroliferous late-Mesozoic and Cenozoic slip-extensional basin. This basin is part of the Tanlu fault zone considered done of the largest continental strike-slip faults in the world  and extends from north east to south west of the northern Chinese province of Jiling covering approximately 2200 km2 with a length of 140 km trending NW and a width that ranges between 12 and 20 km trending EW  . Furthermore, it borders at the north and south with the Dahei and Nadanhada Mountains respectively.
Based on its geological characteristics and basin style, is divided into three subunits; from NE to SW, the Chaluhe fault depression, the Luxiang fault depression and the Moliqing fault depression   .
The Moliqing basin has been divided in four areas according to its geological features: the Kaoshan sag, the Ma’anshan fault-step belt, the Jianshan uplift, and the Xiaogushan slope belt  (Cao, Ye, et al. 2011) (Figure 1).
The Moliqing basin contains seven different sedimentary rocks formations from the Eocene to the Quaternary, overlying a pre-tertiary basement Figure 2. The distribution of the sediment thicknesses inside the Moliqing depression is deeper at the northern area of the basin, diminishing in depth towards the southern zone. The total thickness of the sediment ranges from 2000 m to 5000 m.
3. Data and Methods
The input data required to run the BasinMod 1-D of our study are derived from three wells in the Moliqing basin: well Y8, Y34, and Y46 (Figure 3).
The eroded thickness  (Table 1), present thicknesses of each stratigraphic unit, top and base wellbore depths, and the lithology mixes data are derived from composite well logs. Biostratigraphy is derived from the ditch cuttings which were washed, separated and then analyzed for micro fauna (Foraminifera) and flora (spores and pollens). The location of wells (Figure 3), kerogen types,
Figure 1. Tectonic units of the Moliqing fault depression.
Figure 2. Stratigraphic columns in the Moliqing basin.
Figure 3. Location of the modeled Moliqing basin wells.
Table 1. Calculated eroded thicknesses for the wells used in the 1D modeling (Modified after  ).
total organic carbon (TOC), kinetic algorithms and results of pyrolysis analyses, and organic thermal indicators include vitrinite reflectance values; the water-sediment interface temperatures after  and subsurface bottom hole temperatures (BHTs) through time, the short term sea level proposed by  geologic time scale  , the density of asthenosphere (ρm) and water (ρw) are 3.40 g/cm3 and 1.03 g/cm3 respectively, the Beta factor (ß) and thickness of lithosphere (a) can be obtained accepting default values of 2 and 125 km, respectively, all mixed parameters, the initial porosity, reciprocal and exponential compaction factors, the matrix density, matrix thermal conductivity, and matrix heat capacity are adopted from the default values in BasinMod Software package.
The data and materials for BasinMod 2-D are 1-D well data and Seismic Lines 2, 5, (Figure 4) maps coordinates.
In summary, the modeling processes, its methodology and references are in Table 2 below.
Figure 4. Locations of the cross sections and wells modeled.
Table 2. Mathematical methods selected to build basin models.
4. Results Interpretation and Discussion
4.1. One Dimensional Modeling
The calculations related to this modeling approach include geohistory/burial history, heat flow, temperature, thermal maturity, timing of petroleum generation and pressure.
The 1D model calibration of the Moliqing basin was based on the comparison between the registered formation temperatures obtained during the Bottom Hole Temperature (BHT) tests against the modeled temperature curves.
The northwest border of the Moliqing basin in the Kaoshan sag is the location for Well Y8.
Erosion and sedimentation rate
The basin experienced three uplift processes, producing erosions and cooling effects on the formations and these took place during 42 Ma, 30 Ma and 10 Ma registering erosions of 300 m, 800 m and 500 m, respectively. The missing eroded thickness at this Well is the highest among the studied locations.
The burial history (Figure 5) shows an increase sedimentation rate tendency from the beginning of the basin formation (57.8 Ma) until the starting point of the second erosion event (30 Ma), with values beginning from 63 m/Ma for the Shuangyang formation, up to 147 m/Ma for the Yongji formation. This increase was closely followed by an increment in temperature, causing the Shuangyang formation sediments to reach temperatures over 100˚C. After the second erosion event, the sedimentation rate changed its tendency and started decreasing the magnitude, being 38 m/Ma for the Yongji formation and dropping to 9 m/Ma for the Q formation. The highest paleotemperature was reached during the Neocene 10 Ma with a magnitude over 110˚C.
The calibration of the basin geothermal models was done using the BHT temperature in order to determine the heat flow. The modeled maturity and temperature results have shown a good match with the measured properties (Figure 6), indicating the reliability of the modeled well to describe the thermal settings of the basin. Furthermore, the Tmax based Vitrinite reflectance has been plotted on the maturity charts, and a good tendency have been appreciated, certifying the capacity of this geothermal indicator to be used for calibration purposes.
Source rock maturity
A current heat flow of 64.1 mW/m2 was calculated for the study area. The maturity and temperature chart suggest that two maturity stages have been reached; early and mid-mature. The maturation formation process described by the burial history established that the Shuangyang formation entered the oil threshold 0.5%Ro during the Paleocene at 33 Ma. Afterwards, during the second uplift process (30 Ma), the deepest sediments of the Sheling formation started the transition from immature to the early maturation stage reaching the oil generation phase (Figure 7). After the third uplift process, the deeper sediments of
Figure 5. Burial history and maturity windows of the well Y8.
Figure 6. Maturity and temperature modeling of the well Y8 using two different maturity indicators. (a) Measured Vitrinite reflectance B. Tmax based Vitrinite reflectance.
Figure 7. Hydrocarbon generation history of the well Y8’s source rock.
the Shuangyang formation reached the mid maturity state. The Shuangyang and Sheling formations source rocks, yield oil and gas, being 34 Ma the estimated generation time, however, only the Shuangyang formation expelled oil during the Neocene at 12 Ma. Mg/g TOC.
The Shuangyang formation has been defined as the source rock with the highest hydrocarbon generation potential based on the oil and gas generated (Figure 7); generating 72 Mg/g TOC of oil and 13.42 Mg/g TOC of gas while the Sheling formation generated 11.24 Mg/g of oil and 2.12 Mg/g of gas. Concerning the oil and gas expulsion, the Shuangyang expelled up to the 60% of the oil generated (43 Mg/g TOC) and more than 68% of the gas (7.67 Mg/g TOC).
The well Y34 is located in the Ma’anshan fault step belt of the basin, between the wells Y3 and Y38.
Erosion and sedimentation rate
At this location, the three registered erosions had a magnitude of 200 m, 600 m and 200 m, showing a clear decrement tendency when compare with calculated missing thickness from wells Y8.
The burial history (Figure 8) shows a decreasing sedimentation rate tendency over the time, reaching it maximum rate during the formation of the Shuangyang formation (150 m/Ma), and diminishing the amount of deposited sediments to 62 m/Ma for the Sheling formation, followed by a slight increase to 65 for the Yongji and a finally declining for the Nc (51 m/Ma) and Q formation (8.75 m/Ma). The highest paleotemperature was registered during the Neocene on the Shuangyang formation the over 120˚C.
The model calibration performed by matching the BHT with the modeled temperature curved, generated a close acceptable match between the modeled maturity curved with the calculated Vitrinite reflectance, suggesting the reliability of the computer reconstructed well model (Figure 9).
Source rock maturity
The heat flow was estimated at 78.1 mW/m2. The Shuangyang formation entered the principal stage of oil generation (0.5 - 1 %Ro) during the Paleocene (37 Ma), and currently all the formation reached the mature state having different levels of maturation (Figure 9): the maturity levels at the upper layer sediments is early mature, and at the lower formation layer is mid mature. The lower part of the Sheling formation achieved the oil generation stage during the Neocene (10 Ma) and is currently on early mature stage.
The Shuangyang formation generated oil and gas since the Paleocene at 43 Ma, in two stages (Figure 10). The first stage took place from the 43 Ma until the Neogene 20 Ma, generating hydrocarbon at constant rate over 20 Mg/g TOC.
Figure 8. Burial history and maturity windows of the well Y34.
Figure 9. Maturity and temperature modeling of the well Y34 using two different maturity indicators: (a) Measured Vitrinite reflectance (b) Tmax based Vitrinite reflectance.
Figure 10. Hydrocarbon generation and expelling history of the well Y34’s source rock.
After this point, the oil generation curved showed an exponential increase yielding up to 205.3 Mg/g TOC, and the gas generation was also incremented to 36 Mg/g TOC.
The hydrocarbon expulsion started during the Neocene 14 Ma, and reached 150 Mg/g TOC of oil and 27 Mg/g TOC of gas. The Sheling formation did not generate or expelled hydrocarbons.
The well Y46 is located at the southwest of the basin in the Jianshan uplift near the tectonic boundary.
Erosion and sedimentation rate
The sedimentation rate at the well’s Y46 location shows the lowest level among the studied wells, reaching a peak of 98.75 m/Ma during the Yongji source rock formation (Figure 11). As a consequence of the low sedimentation rate, the burial depth at this location was also the lowest among the studied basin wells, and lower paleotemperatures were generated. The maximum temperatures developed were experienced by the Shuangyang source rock, with values over 90˚C.
Figure 11. Burial history and maturity windows of the well Y46.
In addition to the low sedimentation rate, the erosions encountered by the formations are comparable in magnitude to the missing thickness calculated for the other wells. Consequently, the cooling effect produced by the rifting affected in a higher magnitude the geothermal settings at this location and the maturation of the source rocks.
The simulated Temperature and the BHT data show a good match. Additionally, the maturity curved reaches the 0.5 %Ro threshold at 2000 m, as it was expected (Figure 12). These two elements suggest the reliability of the simulation.
Source rock maturity
The heat flow was estimated at 68.6 mW/m2. As it was mention before, the low sedimentation rate and erosions experienced by the basin caused low paleotemperatures and low maturity in the basin. The Shuangyang formation was the only source rock that entered the oil generation stage during the Neogene (Figure 13); however, the generated oil and gas quantities are negligible: 1.25 Mg/g TOC of oil and 0.22 Mg/g TOC of gas. This source rock did not expelled hydrocarbons.
The reconstruction of 1D thermal history modeling of the Moliqing basin was perform for 7 different locations that describe the different geological settings of the basin. The calculated basin heat flow magnitudes summarized in the Table 3 with values between 64 and 78.1 mW/m2 suggest the Moliqing Basin is characterized with relative higher geothermal background. This geothermal characteristic of the Moliqing basin could favor the thermal maturation of kerogen, offering good conditions for hydrocarbon generation in the basin. The heat flow distribution of the wells Y38 and Y42 has a tendency to grow from the northwest to the southeast, and may attribute to the concentrating distribution of thicker mudstone or shale in northwest basin, which could block the deep geothermal transferring.
Figure 12. Maturity and temperature modeling of the well Y46 using two different maturity indicators: (a) Measured Vitrinite reflectance (b) Tmax based Vitrinite reflectance.
Figure 13. Hydrocarbon generation history of the well Y46’s source rock.
Table 3. Summary of the 1D modeling heat flow and generation time calculations.
The Shuangyang formation is the main source rock of the basin, having generated the majority of the oil and gas produced. The highest amount of generated and expelled oil was located at the location of the wells 34, where also the highest heat flow was estimated.
The Shuangyang and the Sheling formation reached the oil generation stage (0.5 - 1.3 %Ro). Furthermore, the deepest sediments of the Moliqing Shuangyang formation reached mid maturity stage.
4.2. Two-Dimensional Thermal Modeling
Thermal Maturity Results
Thermal Distribution and Maturity Effects
The geomorphology of the Moliqing basin, revealed by the use of the seismic cross sections, exposes an increase of the burial depth in two different directions. The first depth increment is presented along the basin main axis from the southwest to northeast. Additionally, a more dramatic increase is developed from the northwest towards the southeast boundary with differences up to 1800 m.
The modeled temperature profiles (Figure 14 & Figure 15) indicate that the basin temperatures increase with the increment of the burial depth. Similarly, following the temperature tendency, the maturity levels of the source rocks also encounters a higher maturity level towards the deeper location. This is due to the impact that the temperature has on the maturity of the kerogen.
Regarding the BHT measurements, inaccuracies related to the recollection of the data could be responsible for the inconsistence. If the measuring time is shorter than standard, the gathered results would be below the real values. Nevertheless, the scale of the temperature differences suggests that this is not the case. The other cause could be related to the estimation of the thickness of the Shuangyang formation. During the drilling process, the Shuangyang formation was perforated, however, this stratum was never drilled completely, and its real thickness was never registered. The thickness of this formation was consequently calculated using the seismic spectrum in combination with the conversion equation X, and, as it was described earlier the calculation is subjected to errors. If the calculated thickness was sub estimated, the reduction of the strata thickness would decrease the thermal gradient inside the basin, and higher temperatures will be reached. This scenario matched with what is presented in the seismic lines 2D models (Figure 14 & Figure 15).
Figure 14. (a) Temperature cross section and (b) maturity cross section at present day for the Moliqing line 2 seismic line.
Figure 15. (a) Temperature cross section and (b) maturity cross section at present day for the Moliqing seismic line 5.
Due to the inconsistencies above explain, the assessment of a detailed maturity level distribution was not conclusive.
Nevertheless, from the modeled 2D seismic lines it is observable how the generative potential of the basin increases from the Yongji formation, to the Shuangyang.
The verification of the time-depth conversion equation is enhanced in order to increase the precision of the modeled geometry and faults distribution of the basin.
After a careful study of the Moliqing Basin in the Yitong Graben, The following conclusions are presented:
• The Tmax based calculated vitrinite reflectance provided acceptable results.
• The source rocks’ TOC content varies accordingly to their relative basin location and to the depth. The source rock TOC rises from the northeast towards the southwest of the basin, and from the shallower strata to the deeper. Consequently, the Shuangyang formation is the source rock with the highest generative potential with level.
• The three source rocks present in the Moliqing basin, based on the evaluation using the Rock Eval Tmax, have been assessed as mature. These results agree with the results of previous scientific studies  .
• The 2D modeling results were not conclusive due to the lack of reliable maturity indicators to constrain the model.
The main limitation of this study was related to the insufficient and unreliable data needed to characterize the source rock, and to reconstruct and constrain the basin models.
Further study is required to better understand the basin dynamics, especially the geological and geochemical. To achieve this, adequate and reliable data should be readily available.
The authors are grateful to the Department of Oil and Gas exploration, Faculty of Earth Resources, China University of Geosciences (Wuhan) for their support in terms of data and finance.