Effect of intensity of sedimentary cover deformation on hydrocarbon accumulation in Dongying Sag, Bohai Bay Basin, China

The developmental phase of the fault deformation zone denotes the zone of weak deformation (with strong concealment) that evolves within the sedimentary cover of the basin. Recent studies have unveiled the objectively existing tectonic phenomenon of weakly deformed tectonic belts within the sedimentary basin cover, closely intertwined with oil and gas accumulation. To elucidate the deformation intensity and hydrocarbon accumulation scale within the cap cover deformation zone, a pivotal concern in oil and gas geology, this study focuses on the Dongying Sag. The structural physical simulation experiment method, incorporating variable caprock thickness and variable shear strength, is employed to scrutinize the impact of basement fault strike-slip activity on the development of faults in the sedimentary caprock of the basin and dyed oil is charged. In conjunction with sag examples, Early R shear single-channel migration-isolated aggregation, Early and mid-term R shear main channel migration-geese and beaded aggregation, P shear main channel migration-intermittent zonal aggregation, Full channel migration-continuous belt aggregation accumulation models of basement faults are established. It is emphasized that the R shear pressurized deformation section and the R and P shear intersection section in the deformation zone are favorable target areas for oil and gas exploration.


Effect of intensity of sedimentary cover deformation on hydrocarbon accumulation in Dongying Sag, Bohai Bay Basin, China
Weiwei Zhou 1,2* , Changqi Zhao 1 & He Chang 1* The developmental phase of the fault deformation zone denotes the zone of weak deformation (with strong concealment) that evolves within the sedimentary cover of the basin.Recent studies have unveiled the objectively existing tectonic phenomenon of weakly deformed tectonic belts within the sedimentary basin cover, closely intertwined with oil and gas accumulation.To elucidate the deformation intensity and hydrocarbon accumulation scale within the cap cover deformation zone, a pivotal concern in oil and gas geology, this study focuses on the Dongying Sag.The structural physical simulation experiment method, incorporating variable caprock thickness and variable shear strength, is employed to scrutinize the impact of basement fault strike-slip activity on the development of faults in the sedimentary caprock of the basin and dyed oil is charged.In conjunction with sag examples, Early R shear single-channel migration-isolated aggregation, Early and mid-term R shear main channel migration-geese and beaded aggregation, P shear main channel migration-intermittent zonal aggregation, Full channel migration-continuous belt aggregation accumulation models of basement faults are established.It is emphasized that the R shear pressurized deformation section and the R and P shear intersection section in the deformation zone are favorable target areas for oil and gas exploration.
Non-renewable energy sources, such as methane gas, petroleum, and coal, are crucial for sustaining modern society 1 .With consistent growth in global hydrocarbon energy demand, reaching 116 million tons of standard coal equivalent from 2010 to 2019, accounting for 79% of global energy consumption 2 , meeting this demand requires refined exploration processes and structural analysis of the Dongying Sag fault trend.The study aims to establish a comprehensive model for the Dongying Sag fault trend, providing new theoretical foundations for oil and gas exploration 1, 3,4 .The geological term "fault" refers to the rupture of crustal rock layers induced by stress, resulting in significant relative movement along the fault plane 5 .Concealed faults, found in both the basement and subsurface, exhibit distinctive characteristics, including evident fault surfaces (zones) and identifiable features.The developmental phase of deformation zones within faults refers to the presence of weak deformation structures within the sedimentary cover, signifying early and intermediate stages of fault zone formation.Recognition is achieved through the systematic arrangement of various geological units (secondary faults, oil reservoirs, traps, facies belts, depressions, rock masses, buried hills, etc.).To elucidate the deformation intensity and hydrocarbon accumulation scale within the cap cover deformation zone.In recent years, scholars globally have extensively investigated the early and intermediate activities of basement faults and the development of deformation zones in sedimentary basin covers.Referred to by various names such as "immature strike-slip faults," "deformation bands," and "penetrative structures," these phenomena have been associated with thermal evolution, source rock distribution, reservoir properties, trap development, and hydrocarbon accumulation [6][7][8] .While previous research has conducted quantitative studies on structural deformation using structural sandbox models [9][10][11][12] and experimental investigations on hydrocarbon charging processes 8,[13][14][15] , such as Atilla Agdin 16 , who identified the significant impact of the heterogeneous internal structure of fault zones on the fluid flow of oil and gas through field observations and simulation experiments.Scholars like Fu Xiaofei, Luo Qun, and others [6][7][8][17][18][19] further explored the influence of the internal structure of fault zones on fault sealing and opening, as well as their control over the migration and accumulation of oil and gas. This udy represents a pioneering effort in

Deformation zones and hydrocarbon distribution
The Dongying Depression is a Middle and Cenozoic fault basin developed on top of the North China Craton and Paleozoic basement.Two sets of effective hydrocarbon source rocks and multiple reservoirs are developed in the depression, and there are two hydrocarbon-bearing systems in the upper Sha Si-Sha 4 and Sha 2 sections and in the Sha 3-Sha 2 and Sha 3 sections 34 .The development of faults has made the longitudinal transport of hydrocarbons very active, which is matched with the raw storage and cover condition, resulting in multi-phase reservoir formation that has made it as the most oil-rich formation in the Jiyang Depression.The Dongying depression has multiple hydrocarbon generation and discharge centers with multiple phases of hydrocarbon generation and reservoir formation.The hydrocarbons generated in each depression are transported and assembled nearby, forming a relatively independent hydrocarbon system 35 .The basement faults in the Dongying depression were distributed in north-east and north-west directions during the basin-forming period (Fig. 2).
The Bamianhe and Wangjiagang fault zones are trap-level deformation of sedimentary cover zones (Table 1), which are basal fault-slip type in genesis 31,32 .They are manifested as a series of sedimentary cover deformation zones formed by a linear arrangement of en-echelon arrayed blocks, arc-shaped blocks or composite blocks (Fig. 3a and b), with a fault width less than 1 km and a length less than 15 km, generated by the right-lateral slip activity of the north-east-trending basement fault paleo syncline.If the conditions of oil source, storage and transmission are available, it can form oil reservoirs distributed in a string of beads, en-echelon columns and strips, thus identifying the existence of the deformation of sedimentary cover zone.The fault system of Wang jiagang tectonic zone is mainly composed of a series of near east-west and north-east spreading positive faults.The north-east faults are mainly distributed in the flanks of the tectonic zone and are south-east dipping.The near east-west faults are distributed in the main tectonic zone, mainly north-dipping, with the characteristic of arc curvature, and the overall en-echelon column spreading, and the broken nose traps are independently distributed in a string of beads on the deformation of sedimentary cover zone (Fig. 3a).This feature indicates that the main fault surface of the fault zone is not penetrated, and is in the early stage of the deformation zone development, and that the tectonic deformation is weak, and the seismic section shows negative flower-like and semi-flower-like structures (Fig. 3c,d,e).The Bamianhe tectonic zone is a near east-west directional, nearly linear, small fault Yan-row type spreading, and the reverse roof block traps are in a belt-like gathering (Fig. 3b), indicating that the main fault surface has been gradually formed, and the seismic section shows negative flowerlike structure (Fig. 3f).Its evolution degree is greater than that of the Wangjiagang deformation zone."High": In the context of basin structures, "high" refers to a positive or uplifted area, typically compared to the surrounding lower areas."Unconformity-bounded stratigraphic unit": An unconformity is a gap or interruption in the geologic record, indicating a period of erosion or non-deposition.A stratigraphic unit refers to a distinct layer of rock or sediment.An unconformity-bounded stratigraphic unit refers to a specific unit of rock or sediment that is bounded by unconformities above and below it.It is a term used to describe a type of stratigraphic unit that has formed between periods of erosion or non-deposition."Controlling reserves": This term refers to the factors or elements that have a significant influence on the distribution, accumulation, and availability of natural resources, particularly oil and gas reserves.It relates to understanding and identifying the geological and structural features that control the location and extent of reserves within a basin.

Quantitative study on deformation strength of deformation zones
Previous studies have revealed an uncertain functional relationship between basement uplift and basal fault slip 36 .While the angle between associated faults and the principal basement fault (PDZ) demonstrates a functional relationship, its connection to the activity of the basal fault remains unclear [37][38][39][40] .The length of Richard shear exhibits synchronous increase with basal fault slip until the fault zone becomes fully connected 41 .To investigate the activity of the overlying layer deformation zone, the Riedel shear length is extracted from numerous structural physical analogue experiments and normalized with the thickness of the overlying layer (Z).A multiple linear regression model is established using the SPSS software, with the Riedel shear length/overlying layer thickness (y) Table 1.Summary table of identification of deformation of sedimentary cover zones in the Bohai Bay Basin.Note on the Levels in Geological Structural Units: The levels in geological structural units form a hierarchical framework that delineates the relative size and hierarchy of various geological features.At the highest level, the first-order structural unit encompasses expansive geological forms such as basins or mountain ranges.Within the realm of hydrocarbons, these first-order units are commonly referred to as depressions or uplifts, signifying the primary areas of crustal subsidence or uplift.Transitioning to an intermediate level, we encounter the subfirst-order structural unit.This level acts as a bridge between the first-order and second-order structural units.Specifically in the context of hydrocarbons, sub-first-order units denote secondary regions within first-order structural units, exemplified by sag basins.Situated below the first-order level, the second-order structural unit represents comparatively smaller geological structures, including domes or anticline Zones.These secondorder units play an integral role as constituent elements within first-order structural units, contributing to their overall composition.Descending further, we encounter the third-order structural unit, capturing smallerscale structural entities such as sub-sags or uplifts.These third-order units are nested within second-order structural units, accentuating localized geological deformations and structural features.In summary, the levels in geological structural units establish a systematic framework for understanding the scale and arrangement of geological features, enabling comprehensive analyses within the fields of marine and petroleum geology.as the dependent variable and the basal strike slip/overlying layer thickness ( x 1 ) and basal lateral slip/overlying layer thickness ( x 2 ) as the independent variables.Corrected basal slip values for various periods in the Palaeogene of the Bamenhe and Wangjiagang deformation zones are then calculated.

Approach for analogue experiments
A series of physical analogue experiments is conducted to explore the development and evolution of secondary faults during the basal slip process.The experiments vary the thickness of the overlying layer and shear strength to collect relevant parameters 38,42,43

Analysis and processing of experimental data
The study collects extensive data from structural physical analogue experiments, including basal fault slip displacement, lateral displacement, and the length of en-echelon faults (Riedel shear) at various angles of basement fault oblique shearing, overlying layer thicknesses, and shear strengths.To eliminate the influence of overlying layer thickness, these variables are divided by the overlying layer thickness.The length of en-echelon faults/ overlying layer thickness (L NBD = L/Z), basal fault slip displacement/overlying layer thickness (D NBD = D/Z), and lateral displacement/overlying layer thickness (H NBD = H/Z) are chosen as the dependent variable (y) and independent variables ( x 1 and x 2 ), respectively.Bivariate linear regression analysis is performed using these variables.
The model fit analysis yields the R-squared value, which indicates the proportion of the dependent variable's variation explained by the independent variables and serves as a measure of the equation's goodness of fit.A higher R-squared value is desirable, with a value above 0.8 indicating a good fit to the sample points (Table 2).Tables 2, 3, 4 present the results of the multiple linear regression analysis conducted with the SPSS software.Table 1 demonstrates an adjusted R-squared value of 0.809, indicating a good fit to the sample points.The Sig value, representing the significance level, should be less than 0.05.In this study, a Sig value less than 0.05 is considered significant.Table 2 shows that the variables x 1 (basal fault slip displacement/overlying layer thickness) and x 2 (lateral displacement/overlying layer thickness) have a significant impact on the dependent variable y (length  www.nature.com/scientificreports/ of en-echelon faults/overlying layer thickness), with a Sig value of 0. The F-value, representing the goodness of fit, is relatively large with a value of 187.268 (Table 3).
The quantitative relationship between strike-slip displacement of basement faults, basement stretching and twisting angles, and the length of bead-like faults in the cover can be determined by employing the non-standardized coefficients (B values) from Table 4 as the coefficients for the corresponding independent variables: Among them, y = L NBD = L/Z (The ratio of en-echelon cracks length to sedimentary cover thickness), x 1 = D NBD = D/Z (strike-slip displacement of basement faults divided by sedimentary cover thickness), and x 2 = H NBD = H/Z (lateral displacement of basement faults divided by sedimentary cover thickness).Let the angle of the basal fault tension-torsion motion be α, then × 2/ × 1 = tanα, × 2 = × 1*tanα, and Eq.(1) can be transformed into: where y is the ratio of en-echelon crack length to sedimentary cover thickness, x 1 is the ratio of basal fault dis- placement to sedimentary cover thickness, and α is the angle of basal fault rotation.

Results of deformation of sedimentary cover zone activity (D NBD )
Based on the regression analysis of experimental data using a multivariate quadratic function, the thickness of geological strata in the Bamianhe and Wangjiagang areas (after erosion restoration) during different periods of the Upper Cenozoic, Riedel shear length measured from the structural diagrams for each period, and estimated torsional angles, the basement fault slip values were calculated for the Bamianhe and Wangjiagang deformation zones.Using this data, a quantitative analysis of the tectonic activity intensity of the overlying deformation zones was conducted, and the critical value DCV (Critical Value) at which the overlying deformation zone transitions to a strike-slip fault was determined.
The evolution stage of the overlying deformation zone depends on the ratio of the basement fault slip displacement to the thickness of the overlying strata, which we define as the Normal Basement Displacement (D NBD ).A higher D NBD value indicates a higher evolution stage of the overlying deformation zone.The critical value at which the overlying deformation zone transitions to a strike-slip fault (explicit fault) is defined as DCV.It represents the point at which the main fault plane in the overlying strata becomes apparent.As the thickness of the overlying strata increases, a greater basement fault slip displacement is required for the main fault plane to penetrate.When comparing the results, it can be observed that the D NBD (Normal basement displacement) values for each period in the Bamianhe area are greater than those in the Wangjiagang area (Tables 5, 6).This indicates that the evolution stages of the Bamianhe deformation zone during the Upper Cenozoic are higher than those of the Wangjiagang (1) y = 1.71x 1 + 1.840x 2 + 0.912 When DCV < 2.53, the fault zone is in the stage of deformation zone development.When D NBD > 2.89, significant displacement occurs along the fault, and the fault nature becomes apparent, exhibiting explicit strike-slip fault characteristics.At this stage, the fault zone is connected along the strike and exhibits good transmissibility.This explains why the hydrocarbons in the Bamianhe deformation zone exhibit a zone-shaped block aggregation pattern, while in the Wangjiagang deformation zone, they exhibit a linear block aggregation pattern.

Overburden deformation intensity and hydrocarbon charging analogue model setup
The objective of this study is to design an analogue experiment device for investigating the hydrocarbons charging effect of fault trend zones in the Bohai Bay Basin.The device comprises four main components: the experimental model, the fluid charging system (Fig. 4), the formation condition simulation system (Fig. 5), and the data image acquisition system (SONY Alpha A7M2).The experimental model is specifically designed for conducting physical simulations of hydrocarbons charging.It consists of two parallel wooden boards filled with granular media.By introducing simulated hydrocarbons using a mixture of red ink and paste at a volume ratio of 4:1 at different stages, we examined the impact of fault deformation characteristics on oil transport and accumulation.Furthermore, we conducted a detailed analysis of the evolution and reservoir control of fault trend zones in the cap rock of the basin, using the case study of the Dongying Sag as a reference.

Analogue model setup
The experimental model setup consisted of two solid wooden boards measuring 80 cm × 30 cm.One board was securely fixed to the edge of the device frame, while the other board could be moved horizontally at a constant speed of 0.4 mm per minute.This discontinuous motion of the wooden boards simulated dextral strike-slip motion in the basement, enabling us to observe the distribution of cracks in the overburden.The clay model (the physical property parameters of materials are shown in Table 7) was prepared using a mixture of quartz sand (7 cm), a combination of quartz sand and kaolin (10 cm), and kaolin (8 cm) to simulate a formation water www.nature.com/scientificreports/content of 13.5%.In the wooden model, the sand-soil mixture was filled to a specific height, the surface was leveled, and a 0.5 cm thick PDMS layer was applied to prevent clay leakage and ensure the uniform transmission of shear deformation to the upper clay block.During the experiment, we initially conducted the slip test, followed by applying glue to the bottom of the plank to prevent substrate leakage.We simulated hydrocarbons using a mixture of red ink and paste at a volume ratio of 4:1.The fluid pump switch was turned on, and the mixture was charged at a constant speed and a predetermined distance from the bottom.We recorded the pressure gauge and flow count values.

Analogue result analysis
In the high oil charging experiment, during the initial stage of structural deformation, the basal fault exhibited right-lateral displacement with a slip distance (d) of 1.56 cm (D NBD = 0.39) (Fig. 6b.I, b.AB).The overburden deformation zone was in the early phase, showing a weak staggered pattern of fault structures.With a charging pressure of 90 kPa and a charging volume of 12 ml, only a small amount of hydrocarbons accumulated in the fault traps near the source, with a trap area fill ratio of 20%.The basal fault continued to slip, and the cumulative corrected slip distance (D NBD ) reached 3.76 cm (D NBD = 0.94) (Fig. 6b.II).The overburden deformation zone evolved into the early to middle stages, with small cracks expanding and forming larger cracks in a staggered distribution.Small cracks also appeared on the edges of Riedel shear, indicating a charge pressure of 90 kPa and a charge volume of 20 ml.Under sufficient oil supply, hydrocarbons accumulated in a staggered block pattern along the entire overburden deformation zone (Fig. 6b.CD).The section near the oil source was more enriched in hydrocarbons, with a trap area fill ratio of 55%.When the basal fault slip distance reached a cumulative value of d = 7.32 cm (D NBD = 1.83) (Fig. 6b.III), secondary co-directional and reverse shear faults (P) began to appear in the overburden deformation zone.The charge pressure was 70 kPa, and the charge volume was 24 ml.Localized permeability along the strike of the deformation zone facilitated the migration of hydrocarbons (Fig. 6b.III, Fig. 6b.EF), and the trap area fill ratio reached 75%.When the basal fault slip distance reached d = 11.6 cm (D NBD = 2.9) (Fig. 6b.IV), the overburden deformation zone was penetrated by Y-shear faults (Fig. 6.GH), and the fault nature became apparent, exhibiting strong permeability along the strike.The charging volume was 24 ml, and the charging pressure started to decrease to 45 kPa.Hydrocarbons loss decreased, and the entire fault zone became enriched with hydrocarbons, with a trap area fill ratio as high as 90%.Arc-shaped, intersecting, and composite oil-bearing fault traps were formed in the overburden deformation zone (Fig. 6b.IV).By comparing the physical simulation of the structural features and hydrocarbons distribution patterns with those of the Baimianhe and Wangjiagang overburden deformation zones, it was revealed that Wangjiagang corresponds to the "b" evolution stage with a D NBD of 1.83, and Baimianhe corresponds to the "d" evolution stage with a D NBD of 2.93.This corresponds to the slip distances of the latent and manifest stages of the overburden deformation zones calculated based on Formula ①.Therefore, the multivariate quadratic function fitted by the SPSS software is a reasonable representation of the actual geological processes.It can be used to calculate the slip distance of the basal fault, quantitatively evaluate the activity intensity of the overburden deformation zone, and analyze and predict the distribution patterns of hydrocarbons.This method is worth exploring.Under the condition of low oil charging, with the same charging pressure as in the high oil charging experiment, in the early stage of the experiment (Fig. 6a.I,II), when D NBD < 0.95, the trap area fill ratio was less than 20%, and hydrocarbons accumulated more in the traps on the side of Riedel shear.The distribution pattern of hydrocarbons was a staggered block accumulation mode (Fig. 6a.AB) in the early to middle stages of the  7. Physical properties of granular materials: grain size (D), cohesion (C), angle of internal friction (ϕ), bulk density ( ρ b ), porosity (θ)and permeability (K).overburden deformation zone evolution when D NBD reached 1.9.In the sections where P-shear locally penetrated, the hydrocarbons appeared as a string-like distribution, while in the sections where P-shear did not penetrate, the distribution was staggered or point-like (Fig. 6a.III,CD), and the trap area fill ratio reached 45%.As the slip distance of the basal fault continued to increase to D NBD = 2.98, the Y-shear penetrated and formed the main fault plane (PDZ), and hydrocarbons appeared as a belt-like block accumulation pattern (Fig. 6a.IV,EF), with a trap area fill ratio of 70%.Compared to the high oil charging experiment, the degree of hydrocarbons enrichment was lower under the same D NBD conditions in the case of low oil charging.The tendency of hydrocarbons enrichment under the condition of low oil charging was more pronounced in the early to middle stages of the Riedel shear structure in the overburden deformation zone.Although we did not have access to CT equipment to observe the distribution of cracks and fluid within the clay model, Riedel has already completed the analysis of the Riedel shear structure internationally.Riedel's analysis shows that CT scan slices reveal Riedel shear as a helical shape rooted in the basal fault, and the Riedel shear surface on the active side forms a convex shape consistent with rotation.In other words, during the early to middle stages of the development of the overburden deformation zone, the trend surface of Riedel shear is always connected to the basal fault.When the D NBD value reaches a certain threshold, it provides the conditions for vertical migration of hydrocarbons.This is the main reason for the staggered block accumulation pattern of hydrocarbons in the www.nature.com/scientificreports/overburden deformation zone during the early to middle stages of structural evolution.During the structural deformation process, the convex surface of the Riedel shear on the active side acts as a preferred migration channel according to the theory of section accumulation of hydrocarbons.When the hydrocarbons charging volume is limited, hydrocarbons migrate along this preferred channel and preferentially fill the traps on the active side of the main displacement zone (PDZ), resulting in a distribution pattern of hydrocarbons enriched on one side.
In addition to the charging volume, the charging location also significantly affects the distribution pattern of hydrocarbons.In the experiment, the charging location was set in the center of the model, and it was found that during the early to middle stages of the overburden deformation zone when D NBD < 0.95, hydrocarbons preferentially filled the larger faults at a moderate distance from the oil source, rather than the smaller faults nearby or the larger faults farther away.Small faults closer to the oil source can accumulate hydrocarbons but have a lower trap area fill ratio.Under constant pressure experimental conditions, regardless of the charging volume, the overburden deformation zones located farther from the charging location did not accumulate hydrocarbons.
Overall, the experimental results provide valuable insights into the structural deformation and distribution patterns of hydrocarbons in the overburden deformation zone.The sealing techniques employed in the wooden model effectively prevented hydrocarbons channeling along the contact surface, allowing for a more accurate simulation of the natural conditions.The observations and data obtained from the experiments support the use of the multivariate quadratic function fitted by the SPSS software as a tool for quantitative analysis and prediction of hydrocarbons distribution.Further research can explore variations in charging parameters and geological conditions to enhance our understanding of the complex processes involved in the migration and accumulation of hydrocarbons.
Based on physical analogue experiments and research examples on fault-controlled reservoirs, four types of reservoir formation modes in different evolutionary stages of the Dongying Depression's cap layer deformation zone have been established.
1. Riedel shear single-channel migration-isolated accumulation model (Fig. 7a): Under weak deformation conditions, fewer faults appear as discontinuous linear or isolated and dispersed formations.They form a few small-scale fault-block traps with non-connected orientations and poor vertical channel connectivity.These traps require higher charging pressure to form fault-block hydrocarbons reservoirs with low oil columns and small reserves.2. Riedel shear main channel migration-en-line beaded accumulation model (Fig. 7b): Under moderate to strong deformation conditions, there are relatively fewer faults in the cap layer deformation zone, arranged in a serial bead-like pattern.A series of small-scale fault noses and fault-block closures constitute a closure zone with non-connected orientations but good vertical channel connectivity.The resulting fault-block hydrocarbons reservoirs have low oil columns but larger reserves.This mode corresponds to the reservoir formation mode in the Wangjiagang cap layer deformation zone.3. P shear main channel migration-discontinuous zone accumulation model (Fig. 7c): Under moderate to strong deformation conditions, small faults or fissures develop locally in an intermittent belt-like distribution.This forms a large number of larger fault blocks that constitute closure zones.The orientations of the channel connections are intermittent in the vertical direction but have good vertical channel connectivity.
The resulting hydrocarbons reservoirs have medium oil column heights and reserves.4. Full-channel migration-continuous belt-like accumulation model (Fig. 7d): Under high deformation intensity, small faults are densely distributed in a belt-like pattern, or the main fault surfaces are essentially interconnected.Fault block groups are formed in the cap layer deformation zone, and a three-dimensional fully connected migration channel is formed within the deformation zone.The oil charging pressure for these reservoirs is low, and they exhibit a rich variety of reservoir types, large oil column heights, and high reserves.This mode corresponds to the reservoir formation mode in the Baminhe cap layer deformation zone.
Therefore, the blank sections near bead-like reservoirs in the cap layer deformation zone and the intersecting positions of grid-like patterns in multiple orientations are favorable target areas.This understanding provides a new theoretical basis for further hydrocarbons exploration, indicating new directions for exploration.The experiments have preliminarily revealed the quantitative relationship between the deformation intensity of the cap layer deformation zone and the hydrocarbons charging volume.As the deformation intensity increases, the channel conditions improve, resulting in larger oil column heights and greater closure area saturation, leading to higher hydrocarbons enrichment levels (Tables 8 and 9).This understanding lays the foundation for subsequent research to reveal the quantitative relationship between basement strike-slip displacement (deformation intensity of the cap layer deformation zone) and the hydrocarbons charging volume (degree of hydrocarbons enrichment).

Matching relationship between fault activity and oil and gas migration and accumulation periods
Research on the pool-controlling effect of fault trend zones indicates that, beyond adhering to source control theory, the accumulation of oil and gas occurs in annular belts and is influenced by various factors.This comprehension is governed by the basement faults of the formation, demonstrating grid-like aggregation characteristics.The Dongying Depression exhibits the characteristics of late-stage reservoir formation and has two hydrocarbon generation and expulsion periods: the Dongying Late period and the Guantao-Minghuazhen period (Fig. 8), with the latter being the main period.It can be observed that before the large-scale migration of hydrocarbons, the fault activity in the Baminhe and Wangjiagang regions is at its peak, generating a significant number of structural traps and structurally related composite traps within the cap layer deformation zone.During this period, the main fault surface of the Wangjiagang fault zone is not fully connected and has non-connected orientations.According to experiments, it is known that the Riedel shear faults are always connected to the basement, providing relatively good vertical channel connectivity.However, due to incomplete channel connectivity, the volume of hydrocarbons charging is limited during migration.They accumulate and form reservoirs only at the positions of convex-shaped ridge sections of the Riedel shear faults, exhibiting a bead-like block pattern of hydrocarbons accumulation.In contrast, the main fault surface of the Baminhe fault zone is fully connected.www.nature.com/scientificreports/Apart from vertical migration, hydrocarbons can also migrate along the active plate's side, forming dominant migration channels.Hydrocarbons migrate along these channels and preferentially charge the closures on the active side of the principal displacement zone (PDZ), resulting in intermittent-continuous patterns of hydrocarbons migration accumulation.Scholars such as Luo Qun (2010), Fu Xiaofei (2012), Jiang Youlu (2019), and many others 7,17,18,44,45 have recognized that fault zones actually have a binary structure as a three-dimensional geological body.Efficient migration occurs in the section where the fault is fully connected, acting as a "seismic pump."In the development areas of inactive or weakly active deformation zones, induced faults and pores can still undergo continuous migration (prior to fault filling) or slow migration (after fault filling).This understanding aligns with the concept that slow hydrocarbons accumulation can still occur during the early and middle stages of fault evolution when the deformation intensity is low.
In the future research landscape, there are critical avenues to explore in the context of sandbox modeling.Attention should be directed towards refining internal structure and fluid distribution imaging, addressing challenges like sandbox sealing, material compatibility with imaging devices, and the representation of metallic materials.Additionally, future models should consider density differences in internal structural deformation and fluid charging, incorporating parameters like reservoir temperature and pressure.Investigating the impact of

Figure 1 .
Figure 1. Outline of the distribution of basement faults in the Pre-Paleogene in the Bohai Bay Basin.
. The influence of overlying layer thickness and shear strength on the development and evolution of secondary faults in the overlying layer is analyzed.Sensitive parameters are identified based on the trends observed, which serve as indicators of the activity intensity of the overlying layer deformation zone.A tectonic evolution model of the overlying layer deformation zone is established, supported by structural examples.

Figure 2 .
Figure 2. The distribution of oil fields and deformation of sedimentary cover zones in the Dongying Sag.

Figure 3 .
Figure 3. Structural characteristics of the Bamianhe and Wangjiagang deformation zone.(a) Planar structural characteristics of Wangjiagang deformation (WD); (b) Planar structural characteristics of Bamianhe deformation (BD); (c) Cross-sectional structural characteristics of the southern segment of Wangjiagang deformation (WD); (d) Cross-sectional structural characteristics of the central segment of Wangjiagang deformation (WD); (e) Cross-sectional structural characteristics of the northern segment of Wangjiagang deformation (WD); (f) Cross-sectional structural characteristics of Bamianhe deformation (BD).

Figure 5 .
Figure 5. Schematic diagram of the initial experimental setup for the layered sandbox model.

Figure 6 .
Figure 6.Physical analogue experiment of structural deformation and hydrocarbon charging.(a) The analogue experiment results of the total oil charging volume of 43 ml and the hydrocarbons accumulation mode; (b) The analogue experiment results of the oil charging volume of 80 ml and the oil gas accumulation mode.

Figure 7 .
Figure 7.The structure, evolution, and reservoir control model of the deformation of sedimentary cover zone (DSCZ).(a) Riedel shear single-channel migration-isolated accumulation model; (b) Riedel shear main channel migration-en-line beaded accumulation model; (c) P shear main channel migration-discontinuous zone accumulation model; (d) Full-channel migration-continuous belt-like accumulation model.

Figure 8 .
Figure 8. Matching between faulting stages and the main stages of hydrocarbon generation and expulsion (wangjiagang and bamianhe deformation zone).

Table 2 .
Model fit analysis results.

Table 4 .
Model coefficient analysis results.Note that the go-slip displacement in the table = D NBD (base fault go-slip displacement/cover thickness), and transverse displacement = H NBD (base fault transverse displacement/ cover thickness).deformation zone.The slip displacement estimation results show that the critical value for the development of Upper Cenozoic faults in the southern slope of the Dongying Sag falls within the range of 2.53 < DCV < 2.89.

Table 5 .
Strike-slip calculation results of the Bamianhe deformation of sedimentary cover zone.

Table 6 .
Calculation results of strike-slip in Wangjiagang deformation of sedimentary cover zone.

Table 8 .
Quantitative parameter table of deformation strength and hydrocarbon enrichment degree of fault trend zone (cumulative oil charging volume 43 ml).I, II, III, and IV represent the four stages in the development of the low oil charging deformation of sedimentary cover zone.• The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Table 9 .
Quantitative parameter table of deformation strength and hydrocarbon enrichment degree of fault trend zone (cumulative oil charging volume 80 ml).I, II, III, and IV represent the four stages in the development of high oil charging covering rock deformation zones.