1 Introduction
The deformation structure of soft sediments refers to the deformation structure formed by external force acting on loose sediments, which is mainly due to the sudden and substantial increase of pore pressure in water-saturated sediments confined below the impermeable layer (
Liu and Yu, 1990;
Qiao et al., 2017;
Su et al., 2022;
Guo et al., 2023), which are widely developed in sedimentary basins (
Van Loon, 2009;
Li et al., 2019;
Chen, 2020), especially in fluvial and lacustrine sedimentary basins (
Moretti, 2000). They are of great significance to the analysis of the basin evolution background and oil and gas exploration (
Du and Han, 2000;
Chen et al., 2003;
Hurst et al., 2003;
Du et al., 2005;
Qiao and Gao, 2007;
Wei et al., 2007;
Xia and Tian, 2007;
Li et al., 2012;
Chen, 2020). The formation mechanism of soft sediment deformation structure (SSDS) mainly includes external factors such as seismic action and internal factors such as gravity and tide (
Owen and Moretti, 2011). Accurate identification of its formation mechanism is of great significance for understanding the type, occurrence time and intensity of tectonic activities in a basin (
Guo et al., 2023), and clarifying the influence of SSDS on oil and gas reservoirs in the basin (
Qiao et al., 2017;
Liu et al., 2021). The sediments formed by deformation and redeposition due to seismic action transformation are called seismite (
Törő and Pratt, 2016;
Qiao et al., 2017). The original sedimentary structure of sediments is destroyed by the seismic action, formed a large number of interconnected structures, such as microfractures, thus improving the physical properties of the reservoir (
Shao et al., 2014;
Leila et al., 2022), but there is a lack of research on the seismite reservoir.
The upper Es4 member in the western sag of the Liaohe depression is one of the key exploration targets for shale oil in China, and high-yield industrial oil flow has been obtained in many wells. SSDS are widely developed in fine-grained sedimentary rocks in this area. In the current study, the SSDS in the Leijia region of Liaohe depression were studied by means of core observation and thin section identification. The formation mechanism of SSDS in this area was judged, and the influence of SSDS on shale oil reservoirs was discussed, which provides new evidence for the exploration and development of shale oil by studying the relationship between the formation of weak fine-grained sedimentary rocks and seismic events.
2 Geological setting
Liaohe depression is a Cenozoic rift faulted basin in the north-east of Bohai Bay Basin. It starts from the Liaodong depression in the east, with Yanshan subsidence zone in the west and Jiyang depression in the south, showing a structural pattern of “two grabens with one base” (
Yu et al., 2007;
Qi et al., 2013;
Ji, 2016) (Fig.1). The western sag is a secondary narrow and long half-graben-shaped sag which is distributed in the Liaohe depression in the north-east direction with an area of 2560 km
2. The western part is a slope landform, and the eastern side is a boundary fault (
Chen et al., 2002;
Shan et al., 2005). Leijia region is located in the northern part of the western sag, close to Taian-Dawa fault (Fig.2(a)) (
Li et al., 2022). Multistage extensional and transtensional tectonic activities controlled the formation and evolution of the Cenozoic Bohai Bay Basin (
Zhang et al., 2017). The Eocene Shahejie (Es) Formation has experienced three stages of tectonic subsidence in the Leijia region (Fig.2). The target layer of this study matches the first stage subsidence (The study layer is shown in Fig.2(c)). During this period, volcanic and seismic activities were strong, and a large number of NNE-NE normal faults were developed and extended for a long time in the Leijia region. (
Shan et al., 2005;
Li et al., 2010;
Qi et al., 2013) (Fig.2(b)).
The strata in the western sag of Liaohe depression are relatively complex. The basement is mainly Archean strata, and the lithology is mainly granite and gneiss. Proterozoic and Paleozoic strata are developed in some areas (
Wang, 2017;
Wang et al., 2019). The Cenozoic sediments in the study area have a large sedimentation rate and many sedimentary cycles. The Eocene lacustrine sediments are mainly composed of Shahejie Formation, including Es1, Es2, Es3 and Es4 (
Yang et al., 2006;
Shan et al., 2014;
Shan et al., 2017). In this study, the lowest member of Es4 is selected as research object. During this period, a large set of organic rich fine-grained sedimentary rocks were developed in the deep lake area of Leijia region, and oolitic carbonate rocks, conglomerates and sandstones were developed in some areas.
3 Methods
In this study, the characteristics of SSDS in core samples collected from wells of L36, L37, L39, L57, L84, T16, and G25-21 are described in detail. The characteristics of SSDS in single well section and connected well profile are compared, which is used to study the vertical and lateral distribution characteristics of SSDS. The cores collected in this study are 10 cm in diameter. For selected cores, approximately 700 high-resolution photos were taken, and for important cores, photos were taken from different angles to cover the entire core. Moreover, some identified SSDS were sketched for better observation.
4 Deformation structures
4.1 Classification of SSDS
Many deformation sedimentary structures are preserved in the sedimentary rocks of the upper Es4 member in the western sag. According to their development characteristics, they are divided into two categories: soft sediment deformation and brittle deformation. The research area includes six types of SSDS: liquefied dikes, liquefied breccia, convoluted laminae, annular bedding, synsedimentary faults, and vein structures.
4.1.1 Liquefied dikes
Liquefied dikes refer to the dikes in which the liquefied flow of sand layer in the sediment intrudes into adjacent sedimentary layers, which are generally connected with the liquefied lithofacies strata (
Qiao et al., 2017). Liquefied dikes in the study area are generally developed in the upper of Es4. The dikes are mostly sandstone-bearing coarse-grained clay to fine-grained sandstone, with colors ranging from light gray to gray. They are typically composed of soft sediment injections that crosscut the laminated strata. The orientation of the dike body is inclined or horizontal to the long axis and the shape includes radial, banded, tadpole, etc (Fig.3(a)−Fig.3(d)). The dick length varies from 0.5 cm to 3 cm and extends vertically more than 10 cm. Its tail end is gradually pinched out, and the surrounding rock is in traction bending shape with liquefaction flow (Fig.3(e) and 3(f)).
4.1.2 Liquefied breccia
Liquefied breccia is a deformation structure formed by the liquefied sand layer penetrating the adjacent argillaceous layer (
Qiao et al., 2017). The lithology of liquefied breccia in the study area is mainly siltstone or argillaceous carbonate rock, which is intruded and cut by mudstone. Gray siltstone and gray white argillaceous carbonate rock are occur in breccia. The shape is mostly ellipsoidal, the angular angle of breccia is not obvious, and the size varies from several millimeters to several centimeters (Fig.4). Some breccias show obvious drag marks or are separated into reticular patterns under microscope (Fig.4(d)). Distinctive liquefied deformation features can be seen below the edge slide, and the micro breccia formed after thin layer liquefaction of argillaceous calcite can be seen (Fig.4(d)−Fig.4(f)).
4.1.3 Convoluted laminae
The convoluted laminae in the study area is developed in laminated shale. The convoluted laminae layer is sandwiched between the undeformation layers, the thickness of the deformation layer is about 1−30 cm (Fig.5(a)−Fig.5(g)), and small related folds are developed at the top. The convoluted laminae are irregular in shape (Fig.5(a)), with some folds and some undulations (Fig.5(c)−Fig.5(g)). The flank of the anticline is gentle, with no obvious sliding structure and injection vein. The orientation of fold axes varies from inclined to nearly horizontal and the fold amplitudes are about 2 cm, and partial dislocation is visible. Recumbent folds are developed in the horizontal direction, with an average fold amplitude of 5−10 cm, gentle wings, and no obvious sliding structures and injection veins (Fig.5(h) and Fig.5(i)).
4.1.4 Annular bedding
Annular bedding is a circular structure connected by the curling deformation of multiple thin rock layers (
Qiao et al., 2017). Annular bedding is mostly distributed in bedded argillaceous siltstone, showing the characteristics of multiple annular layers nested. The diameter of the inner sheath is about 1 cm, and the that of the outer ring sheath is 2–10 cm (Fig.6). Light crystalline calcite veins are developed along the annular laminae and a small number of micro faults are associated with the annular layer (Fig.6(c) and Fig.6(d)).
4.1.5 Synsedimentary faults
Synsedimentary faults are usually developed locally in shale, showing a group or a series of parallel normal faults, and moving to loose strata vertically or at a certain angle. Most synsedimentary faults are normal faults developed between unaffected upper and lower strata. The lithology of deformation sediments is generally composed of fine to medium-grained calcareous clastic rocks and local argillaceous limestone, some of which carbonate gravels of pebble size. The synsedimentary fault has a large dip of about 50°– 80° with an offset between 0.2 cm and 2 cm, and there is a certain gap in the thick layer of the fault (Fig.7). The combination of graben and barrier is also developed in the synsedimentary fault (Fig.7(e) and 7(f)). The dip angle of the fault plane is large and the fault scale is small. There is a strong liquefaction flow along the fault, and the thickness of the upper plate of the contemporaneous sedimentary layer is obviously greater than that of the lower plate (Fig.7(c) and 7(f)).
4.1.6 Vein structures
Vein structure is a group of fine veins arranged in parallel, which is small in size, several microns to several millimeters in width and several microns to several centimeters in height. Each vein is filled with fractures commented by dark fine particles, mainly clay or mica minerals (Fig.8(a)−Fig.8(c)). They appear as a group of S- or X-shaped (double S-crossing pattern) veins, usually without cutting the bedding plane (Fig.8(d)). The veins are regularly arranged, and the interval between them is generally one to several millimeters wide, which is generally in high angle contact or sub-parallel arrangement with the bedding (Fig.8(d)).
4.2 Temporal and spatial characteristics of SSDS
Based on the cores collected from six wells of L37, L36, L39, L84, T16, and G25-21, the well profile is established (the location is shown in Fig.2). The plane and vertical distribution characteristics of deformation structures and sedimentary sequences are analyzed through a cross-section.
4.2.1 Lateral characteristics
The deformation layer is widely distributed along the fault zone and various types of SSDS are developed (Fig.9). The deformation layer is controlled by the Tanlu Fault Zone and extends north-eastward along the fault zone. In areas with strong tectonic activity, the deformation structure types are diverse and the thickness of soft sediments is large. In areas with weak tectonic activity, the deformation structure types of soft sediments are single and the thickness is small. This large-scale distribution of SSDS is likely to be caused by seismic action. Affected by the fault zone, many high-intensity and multi-period episodic earthquakes occur, which can form a large-scale distribution of deformation layer. The SSDS caused by gravity can only be developed locally, and the extension is not far, without large-scale distribution. Therefore, it can be inferred that the SSDS in the study area may be caused by seismic action.
4.2.2 Vertical characteristics
It can be seen from the single well profile that the deformation layer and the undeformation layer appear alternately in the vertical direction, showing the characteristics of multi-stage superposition (Fig.10). The adjacent deformation layers are separated by the undeformation layers in a certain degree of circulation. In the deformation layer, the distribution of various deformation structures has certain regularity. The top of the deformation layer is convoluted laminae, and the bottom is seismic breccia or synsedimentary fault, which also has certain cyclicity. Strata with the same lithology above and below the deformation layer, one part deforms under the influence of deformation structure, and the other part does not deform (Fig.5(d) and Fig.10). When the seismic activity ends, the overlying strata are usually deposited on the deformation layer, resulting in the same lithology of the upper and lower strata, but only the lower strata deform.
5 Discussion
5.1 Formation mechanism of SSDS
Morphological characteristics of deformation structures in soft sediments depend on lithology, physical properties and driving force systems of sediments (
Ezquerro et al., 2015;
Törő et al., 2015). Most of the deformation structures undergo fluidization or liquefaction deformation due to the decrease of sedimentary strength (
Owen et al., 2011), or brittle deformation of the sedimentary layer under the effect of external stress. However, both the decrease of sedimentary strength and the generation of external stress require the function of triggering mechanism. The triggering mechanism of SSDS mainly includes external triggering factors such as seismic action, and internal triggering factors such as storm, tide and gravity (
Du et al., 2001,
2007). However, it is difficult to distinguish between seismic and non-seismic deformation structures only from the external morphology (
Owen and Moretti, 2011;
Feng et al., 2016). Morphological analysis of deformation structures combined with sedimentary and structural settings help to determine the mechanism of sediment deformation (
Törő et al., 2015;
Törő and Pratt, 2015) (Tab.1).
According to previous studies, the characteristics of SSDS caused by seismic action generally have the following characteristics: 1) it is a small deformation structure between the undeformation layers, and the deformation layers and the undeformation layers have similar lithology and sedimentary facies; 2) its spatial distribution has the characteristics of lateral continuity and vertical periodicity, and its distribution area is often large; 3) its deformation intensity and deformation frequency are obviously related to the fault. In general, the closer to the fault, the greater the deformation strength. (
Hilbert-Wolf et al., 2009;
Van Loon, 2009;
Owen et al., 2011;
Qiao et al., 2017).
From the tectonic setting, the optimal triggering mechanism of these deformation structures can be interpreted as synsedimentary to post-depositional seismic action. The following factors indicate that these deformation structures may be triggered by seismic action. 1) The liquefied breccia in the study area is splicable, indicating that the breccia has no obvious displacement and belongs to
in situ or near-
in situ deposition. It is obviously different from the off-site sedimentary liquefied breccia caused by gravity flow or storms, which is consistent with the characteristics of liquefied breccia formed by seismic action (
Kahle, 2002;
el Taki and Pratt, 2012;
Wallace and Eyles, 2015) (Fig.4); In addition, there is no clear direction of the fold axial plane of the convoluted laminae layer in the study area, which does not conform to characteristic of the convoluted laminae (has the same axial plane) that formed due to gravity. In conclusion, the liquefied deformation structure in the study area is formed by seismic action. 2) There is a systematic and constant spacing between the continuous vein structures in the study area (Fig.8). This phenomenon is caused by the shock wave generated by resonance (
Mazumder et al., 2016). The shock wave can be generated by seismic action and gravity action, but no evidence of gravity action is found in the study interval. Therefore, the vein structure in the study area is most likely to be formed by seismic action. 3) The SSDS is widely distributed in wells of L36, L37, L39, L84, L86, G25-21, and T16, which is consistent with the large-scale distribution of deformation structures caused by seismic action. This is obviously different from the deformation structure caused by gravity flow which only develops locally (
Rossetti et al., 2011;
Li et al., 2012). 4) In the vertical direction, multiple sets of superimposed SSDS are separated by undeformation strata, and the deformation structures show a certain degree of cyclicity (Fig.10). Its top is generally convoluted laminae, and the bottom is generally fractured breccia or liquefied breccia. Meanwhile, the lithology of the deformation layer is the same as that of the adjacent undeformation layer (Fig.5(d) and Fig.5(h)). This phenomenon is most likely caused by seismic action which leads to the destruction of the first deposited rock layer, forming a deformation layer. When the seismic activity ends, the normal sedimentary load will continue to be above the deformation layer. 5) During the sedimentary period of Es4 in the western sag, strong unbalanced fault block movements and a large number of active faults were developed, indicating that these faults may control the topographic relief in the basin and have a great influence on the distribution and deformation of sediments. In this period, the western sag was the main source supply area, forming a lake basin, which had the prerequisite for SSDS caused by seismic action. 6) Paleogene seismites have also been found in the eastern sag of Liaohe depression, including drainage structure, load structure, flame structure, spherical pillow structure and synsedimentary faults. The structure of the western sag is similar to that of the eastern sag (
Liu et al., 2015). Therefore, the simultaneous development of SSDS in the western sag is likely to be the cause of seismic action; 7) In the study area, the tidal effect was weak during the sedimentary period of Es4, and large-scale deformation structure layers could not be formed, so the tidal effect could be excluded to form deformation structures. And during this period, the basin is in a deep lake environment, and the deformation structure caused by groundwater intrusion is generally formed in a river environment, so it does not exist in the lake environment of this study.
Similar to modern earthquakes, paleoearthquakes may be episodic, and there are usually multiple continuous events during seismic activity, and sedimentary deformation structures are one of the good records (
Qiao et al., 2006;
Qiao and Guo, 2013). Several soft sediment deformation layers are separated by undeformation layers, which can be regarded as an earthquake episode. When the time interval between two earthquake episodes is longer than the time interval between the respective earthquake sedimentary records, many earthquake episodes constitute an active period. The thicker undeformation sediment layer between earthquake episodes indicate a longer time interval. In core observation, if multiple soft sediment deformation layers and brittle deformation layers occur continuously, the deformation is most likely caused by paleoearthquakes rather than other factors (
Qiao et al., 2017).
In summary, a small part of the SSDS in the study area may be causeed by internal triggering factors, but extensive deformation structures are mainly caused by seismic action.
5.2 Effect of SSDS on reservoir quality
5.2.1 Pore characteristics of rocks developed in SSDS
By comparing the permeability and porosity of the undeformation and deformation layers in the study area, it is found that the porosity and permeability of the deformation layer are significantly higher than those of the undeformation layer. The porosity of the undeformation layer in the study area is 0.64%−17.55% (mean 10.32%), and the permeability is 0.0008−0.73 × 10−3 μm2 (mean 0.125 × 10−3 μm2). However, the porosity of the deformation layer is 3.73%−29.25% (mean 18.94%), and the permeability of the deformation zone is 0.0016−0.96 × 10−3 μm2 (mean 0.175 × 10−3 μm2) (Tab.2).
A large number of identifiable pores and fractures can be seen in the deformation layers. The width of fractures ranges from tens of microns to hundreds of millimeters, including both brittle fracture and plastic deformation (Fig.11). A large number of pores formed after brecciation of fine sediments can be seen under thin sections, with pore diameters of several hundred microns. This shows that the seismic action has a direct transformation effect on the reservoir, and the formed fractures have a significant effect on connecting pores and improving reservoir permeability. At the same time, it is also conducive to the activities of pore water and groundwater, thus promoting the development of dissolution pores and formed a unified pore, hole and fracture system.
In the core, fractures and liquefaction channels can connect with pores, thereby improving the permeability of reservoirs. Cores fractures are developed in seismite strata of wells L36 and L57, and there are different degrees of oil and gas reservoirs on the surface of cores, and even asphalt occurs along the fractures on the surface of dark shale. The reservoir shows heterogeneity after transformation, and the liquefied dikes can invade and pierce the surrounding rock. For example, after the carbonate veins invade the shale strata with poor permeability, the effective porosity and permeability of the reservoir increase, making the originally dense fine-grained sediments a good reservoir after transformation (Fig.12). These phenomena suggest that seismic action improves the reservoir properties of fine-grained reservoirs in the study area.
The micro-fractures formed by the seismic action provide channels for the later fluid seepage, which in turn to form dissolution pores and improve the porosity of the reservoir. In thin section observation, a series of dissolution pores are developed along the fractures in shale, which effectively increases the porosity of reservoir and provides a good place for oil and gas accumulation.
The upper member of Es4 in the study area belongs to deep lake facies, and a set of fine-grained sediments containing carbonate shale is deposited. In theory, the reservoir conditions should be poor. However, seismic action can make tight reservoirs form fractures and secondary pores, thus improving the physical properties of fine-grained reservoirs.
5.2.2 Significance of shale reservoir exploration
The shale oil resources of fine-grained reservoirs are abundant (
Jiang et al., 2022). However, its ultra-low porosity and ultra-low permeability are the bottleneck that restrict the exploration and development of shale oil. This study found that the porosity and permeability of the deformation layer in the study area are better than those of the undeformation layer, indicating that the seismic SSDS has an important influence on the physical properties of shale oil reservoirs, which is of great significance for the exploration and development of shale oil (reservoir physical contrast see Section 5.2.1). In China, seismites have been found and reported in the fine-grained reservoir strata of several major petroliferous basins, such as Yanchang Formation in Ordos Basin, Xujiahe Formation in Sichuan Basin, Qingshankou Formation in Songliao Basin, and Shahejie Formation in Bohai Bay Basin (
Wei et al., 2007;
Zhang et al., 2007;
Yang et al., 2008a,
2008b;
Shao et al., 2014;
Zhao et al., 2019). The above basins are also the key targets of shale oil exploration in China. The relationship between seismites and shale reservoirs is worthy of attention in future research. Identifying this relationship directly determines exploration deployment and dessert prediction. For shale oil reservoirs controlled by seismic action, favorable exploration targets may be fine-grained reservoirs developed along seismic zones.
6 Conclusions
1) A variety of soft sedimentary deformation structures such as liquefied dikes, liquefied breccia, convoluted laminae, annular bedding, synsedimentary faults and vein structures were identified in the study area. These structures are widely distributed in the horizontal direction, and have the characteristics of multi-stage superposition in the vertical direction. Most of the deformation structures are in situ deformation structures, and their structural characteristics conform to the characteristics of SSDS caused by seismic action, indicating that seismic action is the main driving force for the formation of these structures.
2) The micro-fractures and faults generated by seismic action can become good oil and gas migration channels, and provide channels for later fluid seepage, formed dissolution pores in the formation, and improving the permeability and porosity of tight fine-grained sedimentary rock reservoirs. Therefore, the seismic SSDS is of great significance for the exploration of shale oil reservoirs.