• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Unconventional Petroleum Sedimentology: A Key to Understanding Unconventional Hydrocarbon Accumulation

    2022-02-13 09:18:08CnnZouZnQuJqnZnZynHnyLuJnuZoTnYnZuHuToFnyunZnYumnWnQnZnWnLuHnnLuZqnFnDnLuJnnGoRonLuYn
    Engineering 2022年11期

    Cnn Zou, Zn Qu*, Jqn Zn, Zyn L, Hny W, B Lu, Jnu Zo,Tn Yn,S Zu,Hu To,Fnyun Zn,Yumn WnQn ZnWn Lu,Hnn Lu,Zqn Fn, Dn Lu, Jnn Go, Ron Lu, Yn L

    a Research Institute of Petroleum Exploration and Development, China National Petroleum Corporation, Beijing 100083, China

    b National Energy Shale Gas Research and Development (Experiment) Center, Beijing 100083, China

    c State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China

    d Department of Geology, Colorado College Colorado Springs, Colorado, CO 80903, USA

    e School of Geoscience and Technology, Southwest Petroleum University, Chengdu 610500, China

    f Key Laboratory of Tectonics and Petroleum Resources of Ministry of Education, China University of Geosciences, Wuhan 430074, China

    g College of Geosciences, China University of Petroleum, Qingdao 266580, China

    h State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Institute of Sedimentary Geology Chengdu University of Technology, Chengdu 610059, China

    i China University of Petroleum, Beijing 102249, China

    j Oil and Gas Research Center, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China

    k College of Earth Sciences, Jilin University, Changchun 130061, China

    l School of Energy Resources, China University of Geosciences (Beijing), Beijing 100083, China

    Keywords:Sedimentology Black shales Fine-grained sediments Organic matter accumulation Extraordinarily high organic matter Unconventional petroleum

    A B S T R A C T The commercial exploitation of unconventional petroleum resources(e.g.,shale oil/gas and tight oil/gas)has drastically changed the global energy structure within the past two decades. Sweet-spot intervals(areas),the most prolific unconventional hydrocarbon resources,generally consist of extraordinarily high organic matter (EHOM) deposits or closely associated sandstones/carbonate rocks. The formation of sweet-spot intervals(areas)is fundamentally controlled by their depositional and subsequent diagenetic settings,which result from the coupled sedimentation of global or regional geological events,such as tectonic activity, sea level (lake level) fluctuations, climate change, bottom water anoxia, volcanic activity,biotic mass extinction or radiation,and gravity flows during a certain geological period.Black shales with EHOM content and their associated high-quality reservoir rocks deposited by the coupling of major geological events provide not only a prerequisite for massive hydrocarbon generation but also abundant hydrocarbon storage space. The Ordovician-Silurian Wufeng-Longmaxi shale of the Sichuan Basin,Devonian Marcellus shale of the Appalachian Basin, Devonian-Carboniferous Bakken Formation of the Williston Basin, and Triassic Yanchang Formation of the Ordos Basin are four typical unconventional hydrocarbon systems selected as case studies herein. In each case, the formation of sweet-spot intervals for unconventional hydrocarbon resources was controlled by the coupled sedimentation of different global or regional geological events, collectively resulting in a favorable environment for the production,preservation,and accumulation of organic matter,as well as for the generation,migration,accumulation,and exploitation of hydrocarbons. Unconventional petroleum sedimentology, which focuses on coupled sedimentation during dramatic environmental changes driven by major geological events, is key to improve the understanding of the formation and distribution of sweet-spot intervals(areas)in unconventional petroleum systems.

    1. Introduction

    Unconventional petroleum resources are abundant and widely distributed(Fig.S1(a)in Appendix A),represented by shale oil/gas,tight oil/gas, coalbed methane/oil, heavy oil, and oil sand [1]. The emergence of new production technologies(e.g.,horizontal drilling and multi-stage hydraulic fracturing) has allowed the commercial exploitation and production of oil and gas from tight (lowpermeability) unconventional reservoirs that were unprofitable decades ago. In the United States, between 2009 and 2019, the annual production of shale gas and tight oil (including shale oil)increased from 1.4 × 1011to 7.2 × 1011m3and from 3.2 × 107to 3.9 × 108t, respectively [2]. Outside North America, the exploitation of unconventional petroleum resources is the greatest in China. In 2020, China’s annual shale gas production exceeded 2.0 × 1010m3, and tight gas and tight oil production exceeded 4.5 × 1010m3and 3.0 × 106t, respectively [3]. Although several technical challenges remain, the discovery of unconventional petroleum resources has fundamentally altered the global energy structure.

    According to the US Energy Information Administration, the technically recoverable resources of shale gas and shale oil discovered in the continental United States are 2.1 × 1013m3and 3.3 × 109t, respectively, which are mainly distributed over 22 shale plays[4](Fig.S1(b)in Appendix A).Regarding geological conditions, these shale plays are distributed over the Canadian Shield and its periphery, and accumulated in foreland or cratonic basins(Table S1 in Appendix A). China has 35 organic-rich shale strata(intervals) spanning from the Mesoproterozoic to the Cenozoic[5]. These shale strata were deposited under drastically different environmental settings (marine, marine-continental, and lacustrine) in different types of sedimentary basins (Fig. S1(c) and Table S1 in Appendix A). Shale gas is one of the first unconventional resources exploited industrially in China,with an estimated technically recoverable resource volume of 2.2×1013m3,whereof the largest proportion (1.3 × 1013m3) originates from marine shales[6].Conversely,shale oil resources in China mainly comprise lacustrine organic-rich shale and mudstone. The estimated total mass of the technically recoverable resources of shale oil in China is between 7.4 × 109and 3.7 × 1010t.

    Unconventional petroleum exploration and development mainly target certain sweet-spot intervals and areas [1,6]. The term ‘‘sweet-spot interval” refers to a stratigraphic interval enriched in unconventional oil and gas that can be economically exploited under current technical conditions.The sweet-spot interval typically accounts for a small portion of the total thickness of the shale strata [1,6], and its accumulation area in geography is called the ‘‘sweet-spot area” [6]. Accurately determining sweet-spot intervals (areas) for shale oil and gas has been historically difficult because their formation requires a suite of geological factors, including a high total organic carbon (TOC)content, porosity, gas or oil content, brittle mineral content, high strata pressure coefficient, and abundant lamination fissures and microfractures [6,7]. Therefore, improving the understanding of the accumulation mechanisms of unconventional oil and gas is critical to improve the accuracy of the prediction of the spatial distribution of sweet-spot intervals (areas) in any unconventional reservoir.

    A prerequisite for the accumulation of unconventional petroleum in sweet-spot intervals(areas)is a source rock with substantial amounts of retained hydrocarbons[6,8].The TOC content is the primary parameter when evaluating source rocks, and a threshold concentration of 0.5 wt% TOC characterizes mudstone as a hydrocarbon source rock in conventional petroleum geology [9]. However, shale oil sweet-spot intervals in immature to low-mature stages of organic matter (OM) are usually characterized by TOC≥6.0 wt%, and shale gas sweet-spot intervals in high-mature to over-mature stages of OM are characterized by TOC ≥3.0 wt%[6,10]. OM-rich intervals with a TOC content ≥3.0 wt% were first proposed as extraordinarily high OM (EHOM) deposits [11],wherein OM enrichment is a fundamental factor controlling the formation and distribution of sweet-spot intervals (areas) of shale oil/gas. Similar to tight oil, oil-producing or sweet-spot intervals correspond to sandstones or carbonates with better reservoir properties than those of shales,and are within or immediately adjacent to organic-rich shales. The key to understanding unconventional hydrocarbon accumulation is resolving the mechanism of enrichment of the EHOM and the formation of coexisting high-quality reservoir rocks.

    OM in sediments is mainly derived from primary producers,such as algae, in overlying surface waters. Under normal marine conditions, approximately 10% of the photosynthesized OM migrates from the photic zone and settles in the underlying sediment [12]. Upon the settling of OM at the sediment-water interface, active degradation by aerobic and anaerobic remineralization processes causes less than 1%of the photosynthesized OM to be retained and buried in sediments[13]. The enrichment of OM in sediments and the formation of EHOM deposits are thus likely controlled collectively by high primary productivity,anoxic bottom water conditions, and appropriate sedimentation rates [14], signifying that EHOM deposits only occur in specific depositional environments.For instance,large-scale EHOM deposits are usually closely related to several global or regional geological events, such as sea (or lake) level fluctuations (e.g., Ref. [15]),climate change (e.g., Ref. [16]), oceanic anoxic events (e.g., Ref.[17]), volcanic activity (e.g., Ref. [18]), and biotic events (e.g., Ref.[19]). Similarly, the formation of tight oil sweet-spot intervals(areas) is usually caused by a dramatic change in the sedimentary system associated with tectonic events(e.g.,Ref.[20]),climate,sea level changes (e.g., Ref. [21]), gravity flow events (e.g., Ref. [22]),and other geological processes.Therefore,the formation of unconventional hydrocarbon sweet-spot intervals (areas) may occur from the coupling of sedimentation controlled by major geological events [6,23].

    The quality of petroleum reservoir rocks is controlled by the original sediment composition and diagenetic overprinting.Diagenetic processes in sandstones and carbonates, including compaction, cementation, and recrystallization, also occur in fine-grained mudstones/shales. However, the diagenesis of organic-rich black shales involves a complex mass transfer, such as OM thermal maturation (e.g., Ref. [24]) and clay mineral transformation (e.g., Ref. [25]). The diagenetic evolutionary pathway of shales is controlled by shale lithofacies, which are fundamentally governed by depositional processes, and the environment because sedimentation provides the materials that diagenesis acts upon during burial. Integrating the deposition and diagenesis of unconventional reservoir rocks can help determine potential sweet-spot intervals(areas)in unconventional petroleum systems.

    In this study, first, we review recent progresses in the understanding of depositional and diagenetic processes in the fine-grained sediments of unconventional petroleum systems.Subsequently, we discuss the role of geological events in controlling the formation of unconventional petroleum sweet-spot intervals(areas) by presenting four case studies: the Late Ordovician-early Silurian Wufeng-Longmaxi shale of the Sichuan Basin, the Middle Devonian Marcellus shale of the Appalachian Basin,the Late Devonian-Early Carboniferous Bakken Formation of the Williston Basin,and the Middle-Late Triassic Yanchang Formation of the Ordos Basin. This study aims to review relevant sedimentological data obtained from these shale strata to improve the understanding of unconventional hydrocarbon accumulations from a new perspective—unconventional petroleum sedimentology—and hopefully define new perspectives for the future exploration and development of unconventional petroleum resources.

    2.Depositional processes in unconventional petroleum systems

    2.1. OM-rich shales or mudstones

    2.1.1. Deposition of shales or mudstones

    Shales and mudstones are fine-grained sedimentary rocks composed of ≥50%particles with a grain size<62.5 μm.Although they are both widely used terms for fine-grained terrigenous clastic rocks(Fig.1),the word‘‘shale”is commonly used in the petroleum industry. Nevertheless, the term ‘‘shale” should be used with caution because it implies fissility,usually a byproduct of weathering,which is not always observed in freshly exposed rocks [26]. Typically, mudstones contain components of varying grain sizes and compositions derived from weathering, primary production, and diagenesis [27]. Because of the large cryptic variation in grain size and sedimentary features, mudstones have historically been considered to be mostly deposited through suspension fallout under low-energy hydrodynamic conditions. The well-developed fine laminations(undisrupted by macrobenthos)[28,29](Fig.1),which are commonly present in many organic-matter-rich mudstones,have been used as key evidence to suggest that these sediments were deposited under persistently stagnant (and anoxic) bottom water conditions through the continuous suspension settling of mineral grains and OM(e.g.,Ref.[30]).However,this simple model of mud deposition is being increasingly challenged.

    Over the past few decades, a paradigm shift has occurred for clarifying the processes that govern the transport, deposition,and erosion of fine-grained sediments, owing to studies on mud transport and deposition from multiple aspects. The dispersal and accumulation of mud by several hydrodynamic processes(e.g.,river floods and oceanic currents)in relatively energetic environments have been directly documented in many modern muddy shelves [31]. Mud in rivers and seas—regardless of salinity—is likely to form flocculated aggregates, termed ‘‘flocs” [32]. The transport and deposition of mud as flocculated aggregates have also been commonly observed in many modern mud-dominated depositional systems [33].

    Advances in the understanding of the depositional processes of fine-grained sediments have also been complemented by findings from flume experiments. These findings indicate that flocculated muds can be transported in bedload and form floccule ripples at sufficiently strong flow velocities(~25 cm·s-1at 5 cm flow depth),thereby facilitating the bedload transport of sand [32,34]. The increasing number of careful facies analyses of ancient mudstone successions has also revealed many sedimentary and biogenic features in mudstones. The common presence of low-angle ripplescale cross lamination and localized erosional features in mudstones indicates the influence of flows that carry sediment across the seafloor or lake bottom [35]. The wide range of sedimentary structures (e.g., current and wave ripple cross lamination, parallel lamination,and grading and loading structures)identified in mudstones indicate that the deposition of fine-grained sediments is subject to various types of gravity flows,such as turbidity currents,muddy hyperpycnal flows, storm-induced currents, tide-related currents, or a combination of two or more flow types [28,29,36-39](Fig.1).Features such as erosional surfaces,lag layers(of detrital, biogenic, or diagenetic origin), and enriched early diagenetic products (e.g., carbonate concretions, pyrite nodules, and phosphatic particles)indicate erosion,hiatuses,or condensation within fine-grained sedimentary successions [22].

    The accumulation of mud, as well as the facies variation in mudstones, is governed by the rate and characteristics of the mud supply and energy levels in the depositional environment.Recent multifaceted studies have proven that the model of suspension settling of mud under completely stable and anoxic conditions is an exception rather than the norm. In general, mud accumulation appears to be much more episodic and instantaneous than conventional considerations. Many mudstone successions, even those that are OM-rich, have been re-interpreted to have been deposited under the influence of bottom currents in oxygen-rich conditions after sedimentary facies analyses and petrographic studies (e.g., Ref. [40]).

    2.1.2. OM accumulation

    OM accumulation in sediments is mainly controlled by primary productivity, bottom water redox conditions, and sedimentation rate[14].Primary productivity represents the rate of the photosynthetic conversion of solar energy to chemical energy stored in the plant biomass [41]. Bottom water redox conditions are typically determined between 0.5 and 1.0 m above the sediment-water interface [42]. Redox conditions can be categorized into oxic(> 2.0 mL·L-1O2), dysoxic (2.0-0.2 mL·L-1O2), suboxic (0.2-0 mL·L-1O2), anoxic (0 mL·L-1O2), and euxinic (anoxic and free H2S)conditions[12,43].However,suboxic conditions are generally considered to be anoxic.Thus,we subdivided the redox conditions into oxic, dysoxic, anoxic, and euxinic conditions. In principle,redox conditions in a water column are determined by physical and biogeochemical processes. Physical processes involve oxygen solubility and oxygen resupply from vertical mixing and horizontal ocean circulation(including bottom currents),and biogeochemical processes involve oxygen demand during OM decomposition.Although more than 90% of the OM derived from primary productivity can be decomposed under both oxic and anoxic conditions,anoxic environments preserve more residual OM than oxic environments because anoxic conditions can decrease the degradation rate of some metabolizable materials [42]. The gradual buildup of anoxic conditions can lead to a TOC content of 1-3 wt%,and a TOC content higher than 3-6 wt% mostly occurs because of primary productivity [11,42]. As elevated levels of primary productivity lead to a high in situ oxygen demand, marine anoxia is often accompanied by high primary productivity.

    The sedimentation rates often control the dilution effect and decomposition time of OM in different bacterial zones. High sedimentation rates (e.g., >100 cm per thousand years (kyr)) in carbonate rocks or sandstones dilute OM concentrations in sediments. In modern sediments, a strong correlation is often observed between high productivity and high sediment accumulation rate[42].OM associated with high sedimentation rates is generally reactive.The decomposition of OM in sediments can rapidly consume oxygen in the porewater.Thus,rapidly accumulating sediments tend to have anoxic pore waters regardless of the redox conditions in the overlying bottom water [42]. Water column redox conditions play an important role in OM accumulation when sedimentation rates are<100 cm·kyr-1,especially at<10 cm·kyr-1.Comparing the effects of primary productivity and redox conditions on OM enrichment under low sedimentation rates is thus necessary. The most important factors for OM enrichment are enhanced primary productivity and suppressed OM decomposition.

    2.2. Deposition of siltstones and fine sandstones

    Fig.1. Some typical sedimentary structures in shale/mudstone.(a)Parallel lamination,small ripple lamination,Wufeng Formation,Ordovician,China;(b)parallel lamination,Marcellus shale,Devonian,southwestern Pennsylvania,USA[28];(c)parallel lamination and ripple lamination,Bakken Formation,Mississippian,Montana,USA;(d)current ripple lamination,Qiongzhusi Formation,Cambrian,Sichuan Basin, China;(e)current ripple lamination,Yanchang Formation,Triassic, Ordos Basin,China;(f)current ripple cross lamination, Lucaogou Formation, Permian, Junggar Basin, China; (g) wave ripple lamination, Eagle Ford shale, Cretaceous, Texas, USA [29]; (h) lag layer consisting of fossil fragments and reworked pyrite(black-colored particles),Wufeng Formation,Ordovician,Sichuan Basin,China;and(i)lag layers of fossil fragments,Bakken Formation,Mississippian, Montana, USA.

    Another important type of unconventional hydrocarbon resource comprises siltstone and fine sandstone reservoirs; generally, tight sandstone reservoirs are defined as low-permeability sandstone reservoirs with an in situ matrix permeability to gas of less than 0.1 millidarcys (MD, 1 mD = 1 × 10-3μm2), exclusive of the natural fracture permeability [44]. From a depositional perspective, tight sandstones can be formed by the co-deposition of sand and mud:The excessive muddy matrix fills the pore spaces,significantly decreasing the permeability of sandstones.Thus,tight sandstones may occur as immature sandstones deposited close to the source area or as muddy debris flows, hybrid flows, or hyperpycnal flows [45]. However, tight sandstone reservoirs have been found in a wide range of depositional environments (e.g., fluvial,shallow marine,and deep-water environments)[46].This indicates that rather than specific depositional processes, sediment composition and depositional conditions exert a first-order control on subsequent diagenetic alterations, which further alters porosity and permeability characteristics and leads to the formation of tight sandstones. For example, tight sandstone reservoirs in the United States are generally clean sandstones deposited in high-energy depositional environments, wherein intergranular pores are occluded by authigenic cement rather than immature muddy sandstones [47]. The formation of tight sandstone reservoirs also requires relatively quiescent tectonics and the deposition of regionally continuous sandstones with low permeability.Later tectonic alterations, such as faulting and fracturing, as well as overpressure caused by hydrocarbon generation, can then enhance the quality of tight sandstone reservoirs [48].

    Deep-water gravity flows represent one of the most important sediment transport mechanisms on a global scale [49] and can yield large, thick sand bodies of coarse-grained and associated fine-grained clastic sediments on submarine slopes,abyssal plains,and deep lakes [50]. Thick-bedded coarse-grained clastic rocks,such as tight sandstone reservoirs, are usually close to highquality source rocks. The hydrocarbons generated from the source rocks readily migrate to these sandstones and form high-yield lithologic reservoirs [50]. Fine-grained deposits produced by deep-water gravity flows are usually rich in organic material and are potential sweet-spot intervals (areas) for unconventional petroleum systems [50].

    3. Diagenetic processes in unconventional petroleum systems

    3.1. Evolution of OM and organic pores

    OM in black shales is a source of oil and gas in conventional and unconventional petroleum systems[24,51].Dispersed OM in black shales is composed of macerals of different origins[14,24,52].During thermal maturation, oil-prone macerals, such as amorphous OM and alginite,first transform to pre-oil bitumen through a bituminization process with little or no migration involved and subsequently transform into oil and post-oil bitumen in the oil window and condensate-wet gas window [24,51,53]. The post-oil bitumen continues to transform into gas and pyrobitumen in the dry gas window. The generated oil can also secondarily crack into gas and pyrobitumen in the dry gas window [51,53]. However, gasprone macerals (e.g., vitrinite), inert macerals (e.g., inertinite),and zooclasts (e.g., graptolites) do not exhibit significant changes in morphology and occurrence during thermal maturation because of their low hydrocarbon generation potential[24,53]. These macerals can still be observed in high-maturity shales.

    Secondary maceral solid bitumen becomes the dominant OM in black shales when the thermal maturity reaches the peak oil window maturity (vitrinite reflectance (Ro) 0.8%-1.0%) [24,51,53](Fig.S2 in Appendix A).Solid bitumen occurs as a speckled,wispy,and interconnected OM. It occupies the interparticle spaces between mineral grains (Fig. 2) and intraparticle spaces within mineral grains and fossil shells, such as in foraminifera tests[54,55]. As solid bitumen was once liquid, it may form a threedimensional interconnected OM pore network with organic pores at high maturity, preserving methane and improving porosity[24,51].

    In unconventional petroleum systems, OM-hosted pores are an important constituent of the pore network (Fig. 3) and can be the dominant pore type in some gas shales [56]. The gas content and methane adsorption capacity of gas shales are positively correlated with the TOC content,indicating that organic pores play an important role in shale gas preservation [6,57]. The development of organic pores is controlled by both thermal maturity and OM type,and their preservation is controlled by thermal maturity,OM content, shale composition, and pore pressure [53,56,58].

    Organic pores in black shales may have primary or secondary origins. Primary organic pores occur within structured OM, such as cellular pores derived from the cell lumen of higher terrestrial plants [53,58]. Secondary organic pores are generally thought to form when gaseous hydrocarbons are generated and expelled from residual OM (solid bitumen or pyrobitumen) [51,53,55,56,58,59].Secondary organic pores are more abundant in the gas window than in the oil window because ①gas generation and expulsion mainly occur in the gas window, and ②the migration of bitumen and oil in the oil window can infill newly formed organic pores and mask their presence [24,53].

    Fig.2. Photomicrographs(reflected white light and oil immersion)of pyrobitumen(red arrows) in black shales. (a) Wufeng-Longmaxi shale (Ro 3.07%). Sample from Changning County, Sichuan, China. (b) Marcellus shale (Ro 2.41%). Sample from Canastota, New York, USA.

    3.2. Diagenesis of minerals and mineral-associated pores in shales

    The brittle mineral content controls the fracturing effectiveness of shale and is essential for the commercial exploitation of shale oil and gas. As the main brittle component in organic-matter-rich shales, quartz not only has a detrital origin but also an authigenic origin [54,60] (Fig. 4). The dissolution and reprecipitation of siliceous biological skeletons have been proposed to explain the formation of microcrystalline authigenic quartz infill found in intergranular pores/biocavities in shales, such as the Chattanooga,Barnett, and Mowry shales in the United States [61]. The pressure dissolution of clastic quartz particles and the transformation of clay minerals can also lead to the extensive development of quartz overgrowth, as documented for Haynesville shale in the United States [62]. Therefore, authigenic quartz from the dissolution and reprecipitation of siliceous biological skeletons and transformation of clay minerals are the main mechanisms for the formation of authigenic quartz in shale oil and gas reservoirs.

    Fig. 3. Nano-scale organic pores in the Wufeng-Longmaxi shale, Sichuan Basin.Sample from Changning County, Sichuan, China. Organic pores in panel (a) are protected by a rigid framework composed of diagenetic quartz (Qtz). Organic matter in panel (b) has a spongy texture.

    Calcite cementation, replacement, and recrystallization are common in shale and mudstone [63] (Fig. S3 in Appendix A).Calcareous biological components (e.g., foraminifera); seasonal carbonate deposition formed by periodic algal blooms;and biological, chemical, and mechanical compaction typically provide the materials for calcite recrystallization during diagenesis (Fig. S3)[64]. Sparry calcite usually occurs in organic-matter-rich mudstones.It is mainly distributed at the edge of the OM laminae,indicating a close relationship between sparry calcite and hydrocarbon generation. The release of CO2and organic acids during hydrocarbon generation causes the dissolution of carbonate minerals. The resulting pressure or concentration gradient can then induce transportation in the pores and fractures over a short distance via infiltration or expansion, thus precipitating granular carbonate minerals when the binding Mg2+is removed from fluids [65].

    Clay minerals are closely related to oil and gas generation and exploration. Acidic clay minerals have been demonstrated to have catalytic effects on hydrocarbon generation reactions of OM.Therefore,the diagenetic evolution of clay minerals and their interactions with OM have attracted extensive attention[66].The transformation of smectite into illite through the illite/smectite (I/S)mixed layer is the most common type of clay mineral diagenesis;however, the underlying mechanism remains controversial [67].Nevertheless, the temperature of each stage of smectite transformation to illite is generally consistent among studies.The transformation of smectite to the I/S mixed layer occurs at 70-95 °C [68],and approximately 95%of smectite can be converted to illite before the temperature reaches 200°C[69];this process is also controlled by the chemical composition of the interlayer solution, formation pressure, and time.

    The primary sedimentary components control the diagenetic pathway of shale/mudstone and may affect the pore evolution.Although authigenic minerals formed in the early diagenetic stage can reduce porosity, they effectively inhibit mechanical compaction in later stages[68,70,71](Fig.5).Shales rich in quartz,particularly diagenetic silica, contain more organic pores than clayrich shales because the former can form a rigid silica framework that protects organic pores from mechanical compaction[72].Dissolution pores and OM-hosted pores typically develop in the middle and late diagenetic stages. Additionally, under the influence of pressure and tectonic stress, fractures of different scales formed during diagenesis can store oil and gas and act as the main conduits for hydrocarbon migration,greatly improving the permeability of shale and mudstone reservoirs.

    4. Sedimentological implications for unconventional hydrocarbon accumulation

    The formation of unconventional hydrocarbon systems is the result of the depositional and diagenetic processes of finegrained sediments coupled with hydrocarbon accumulation.Depositional environments control OM enrichment (forming source rocks and shale oil/gas reservoirs) or the concurrent deposition of coarse-grained clastics and carbonates(forming tight oil/gas reservoirs). Diagenesis involves mineral alteration, reservoir porethroat evolution,and the thermal maturation of OM. Hydrocarbon accumulation mainly refers to the formation and distribution of sweet-spot intervals(areas)in unconventional petroleum systems,and its prerequisites include large-scale hydrocarbon expulsion and migration, effective trapping of hydrocarbons, and sealing by roof and floor strata. Fundamentally, the sedimentological history controls the development of source rocks, reservoirs, sealed roof and floor strata, diagenetic evolution, and the hydrocarbon accumulation process, which collectively control the formation of unconventional petroleum systems [6,23].

    Unconventional oil and gas exploration and development target sweet-spot intervals (areas) within shale strata, that is, the EHOM intervals and coexisting tight sandstone or carbonate intervals[6,11]. The formation of the EHOM and sealed strata is closely related to their depositional settings and is the result of sedimentation coupled with several global or regional geological events,such as tectonic activity, volcanic activity, sea level (lake level)fluctuations, climate change, water column anoxia, biotic mass extinction, and gravity flow events. The coupled spatiotemporal deposition of these geological events was recorded in detail during several major transition intervals in the Phanerozoic (Fig. S4 in Appendix A). Four typical examples of unconventional hydrocarbon reservoirs are discussed in the context of how coupled sedimentation and geological histories influence the formation of their sweet-spot intervals, as well as hydrocarbon accumulation.

    4.1. Typical unconventional petroleum reservoirs

    4.1.1. Late Ordovician-Early Silurian Wufeng-Longmaxi shale gas

    Fig. 4. Typical authigenic quartz in shales/mudstones. (a) Radiolarian fossils in Longmaxi shale, Weiyuan,Sichuan Basin. (b) Sponge spicules (s)were replaced by quartz in Barnett shale, Texas, Fort Worth Basin. For classification purposes, dolomite replacement (d) is added to the bio-siliceous grain volume [54]. (c) Silica nanospheres (Si) and plate-like authigenic quartz in Woodford shale, Texas, Permian Basin [60]. (d) Micron-scale authigenic quartz (small red triangles) aggregates in Longmaxi shale, Weiyuan,Sichuan Basin.

    From the Late Ordovician to the early Silurian, Wufeng-Longmaxi shale was widely deposited throughout the Yangtze Shelf Sea, spanning across 13 graptolite zones, during five stages[23] (Fig. 6). Although its total thickness generally exceeds 300 m, only the graptolitic black shales are excellent source rocks and the major target for shale gas exploration and development,with a sweet-spot area of approximately 2 × 104km2in south China [6] (Table S1). Stratigraphically, these graptolite shales (i.e.,sweet-spot intervals)are mainly located at the bottom of the shale strata and comprise the most organic-rich intervals(TOC >3.0 wt%)in the Wufeng and Long-1 member shales [23] (Fig. 6). These sweet-spot intervals are generally 10-40 m thick and have a high porosity(>4.0%),gas content(>3.0 m3·t-1),brittle mineral content(>70%),and strata pressure coefficient(>1.2),with relatively abundant lamination fissures and microfractures [6]. In 2021, the annual production of Wufeng-Longmaxi shale gas was 2.3 × 1010m3, and the cumulative proven reserves of shale gas were over 2.7 × 1012m3, making China the largest shale gas producer outside North America.

    4.1.2. Middle Devonian Marcellus shale hale gas

    The largest unconventional gas-producing formation worldwide, the Marcellus shale,was deposited during the Middle Devonian in an asymmetrical northeast-trending Appalachian foreland basin in southeastern Laurentia (Fig. 7 [73-81]). It contains three members:Union Springs shale,Cherry Valley limestone,and Oatka Creek shale[74].According to the TOC content,the Marcellus shale can be divided into two intervals: an organic-rich interval and an organic-poor interval. The organic-rich interval comprises the entire Union Springs shale and the lower part of the Oatka Creek shale, characterized by high TOC contents (ranging from 5 to 13 wt%, with an average of 7 wt%) [74,78-80]. The organic-poor interval comprises the upper part of the Oatka Creek shale, with a comparatively lower TOC content (ranging from 2 to 4 wt%, with an average of 3 wt%) [74,78,80]. The organic-rich interval, which primarily consists of black shales with a low clay mineral content, is the most important gas-producing part and a sweet-spot interval in the Marcellus shale. The reservoir space in these shales consists of interparticle pores and open fractures,with a porosity of 3%-10%[82],which is significantly affected by the OM content and thermal maturity[83]. In shale gas industrial production areas, the thickness of the sweet-spot intervals generally exceeds 15 m[84].The Marcellus shale is the largest single natural gas-producing formation in the United States,with a daily production of over 6.7 × 108m3in 2020, accounting for more than onethird of total daily US production [2].

    4.1.3. Late Devonian-Early Carboniferous Bakken tight oil/shale oil

    The Bakken Formation was deposited in the Williston Basin,an elliptical intracratonic basin in western Laurentia,during the Devonian-Carboniferous transitional interval (Fig. 8 [19,21,48,85-96].Stratigraphically, the Bakken Formation can be subdivided into four stratigraphic members in ascending order: the Pronghorn,Lower Bakken, Middle Bakken, and Upper Bakken members [86](Fig.8).Of these four members,the Lower and Upper Bakken members have average thicknesses of 3 and 2 m,respectively[97].However, the Lower and Upper Bakken members are world-class petroleum source rocks. They are characterized by organic-rich,massive to parallel-laminated black shales,with average TOC contents between 8 and 10 wt% [21,92,93,97].

    Fig.5. Mineral,OM,and pore evolution during diagenesis[68,70,71].Vol.:volatile;bit.:bituminous;anthrac.:anthracite;Opal-A:amorphous opal;Opal-CT:crystalline opal.

    Fig.6. Major geological events during the Ordovician and Silurian transition and characteristics of shale gas sweet-spot intervals for Wufeng-Longmaxi shale in the Sichuan Basin, south China (modified from Ref. [23]). Ma: million year; KB: Kuanyinchiao bed; LOMEI: Late Ordovician mass extinction interval.

    Fig. 7. Middle Devonian geological records and characteristics of shale gas sweet-spot intervals of Marcellus shale in Appalachian Basin, North America. Generalized stratigraphic column of the Marcellus Formation is from Ref. [74], with the star indicating the age of the volcanic ash [73]. Tectonic activity intensity is modified from Ref.[74].Bottom water redox conditions are based on Ref.[75].The global sea level change curve is from Ref.[76].The relative sea level change curve is from Ref.[77].The clastic input curve is from Ref. [74]. The data on TOC contents are collected from Refs. [74,78-81].

    The Bakken Formation is a representative unconventional hydrocarbon play in the United States that produces more than 1.6 × 105tonnes of tight oil and 5.5 × 107m3of shale gas per day in 2020 [2]. Most of the hydrocarbon resources in the Bakken Formation were extracted from the Middle Bakken and Pronghorn members.The Middle Bakken member has an average thickness of 13 m[98].Tight oil sweet-spot intervals within the Middle Bakken comprise clean fine-grained sandstones,sandy skeletal-lime grainstones, and microcrystalline dolomites, with an average matrix porosity of approximately 9%, and maximum up to 14%[48,87,94-96]. The estimated prospective area of the Bakken tight oil/shale oil sweet-spot area is over 3.4 × 104km2in the United States and more than 2.3 × 104km2in Canada [99] (Table S1).

    4.1.4. Middle-Late Triassic Yanchang tight oil/shale oil

    During the Triassic period,the Yanchang Formation was deposited in the Ordos Basin, a rectangular basin on the southwestern margin of the north China Block. The Yanchang Formation can be subdivided into ten members based on lithofacies associations,logging characteristics, and tuff marker layers: Chang 1-Chang 10 members with increasing stratigraphic depths. The Chang 7 member mainly comprises black shales and mudstones intercalated with siltstones and sandstones [23] (Fig. 9). Black shale intervals at the bottom of the Chang 7 member (Chang 73)are the sweet-spot intervals for shale oil and have high TOC contents (up to 30 wt%), with a thickness of approximately 30-60 m. Sandstones in the middle and upper parts of the Chang 7 member (Chang 72and Chang 71), which were deposited by gravity flows, are sweet-spot intervals for tight oil (Fig. 9).Sweet-spot areas for tight oil and shale oil over approximately 1.06 × 104and 2.60 × 104km2, respectively [100] (Table S1). In 2020, 9.31 × 105t of tight oil was produced from the Chang 7 member [101]. In 2021, the cumulative proven reserves of tight oil in the Chang 7 member were over 1 × 109t, and prospective shale oil resources were approximately 6 × 109t [100].

    4.2. Geological events for EHOM accumulation

    4.2.1. High nutrient influx promotes primary productivity

    Fig. 8. Late Devonian-Early Carboniferous geological records and characteristics of the tight oil/shale oil sweet-spot intervals in the Bakken Formation in Williston Basin,North America.Generalized stratigraphic column is modified from Ref.[86].Tectonic event columns are based on Refs.[19,88,89].Global sea level change curve is from Ref.[76].Climate change information is from Refs.[19,93].Bioturbation is redrawn from Ref.[91];the width of bars represents the relative strength of bioturbation.Bottom water redox conditions are based on Ref.[89].Data on TOC contents are collected from Refs. [21,92,93].Data on porosity is collected from Refs.[48,94-96].Distribution of sweetspot intervals is from Ref.[87].Biodiversity affected by the Hangenberg Crisis is from Ref.[19];crosses represent extinction,and the width of the bars represents the relative abundance of taxa. PD: present day; Iceh.: icehouse; DCB: Devonian-Carboniferous boundary.

    Fig. 9. Major geological events during the Middle-Late Triassic and characteristics of tight oil/shale oil sweet-spot intervals of Yanchang Formation in Ordos Basin, north China (modified from Ref. [23]).

    OM in sediments is derived from primary producers, either in the overlying surface waters or on land. In open marine environments, primary producers in the ocean surface are dominated by various types of phytoplankton, and the total biomass determines the OM input. A high nutrient supply fuels primary productivity and ultimately determines the degree of OM production in the water column [102]. The production of EHOM generally requires additional nutrient sources derived from enhanced weathering on land, upwelling of nutrient-replete deep-water masses, or hydrothermal fluids. These additional nutrient sources are closely related to important geological events, including tectonic activity,climate change, sea level (lake level) fluctuations, and volcanic activity [15,18]. The production of EHOM sediments associated with climate change controls the flux and pathways of terrestrial nutrients delivered to sedimentary basins.Regional climate change from arid to more humid conditions typically enhances continental weathering and runoff,resulting in an increased input of terrestrial nutrients that induce suitable conditions for phytoplankton blooms in surface waters. For example, EHOM deposits (TOC content up to 35 wt%)of the Lower and Upper Bakken were deposited during humid climate conditions, and those in the Middle Bakken with poor OM corresponded to semiarid climatic conditions with little rainfall [19,90] (Fig. 8).

    Ancient OM-rich deposits associated with upwelling are widely distributed globally[21,103].Passive upwelling driven by sea level rise supplied large amounts of nutrients and greatly increased primary productivity, leading to EHOM enrichment in the Lower and Upper Bakken shales. Indeed, large-scale sea level changes have been invoked as a driving mechanism for enhanced nutrient input via upwelling during the Devonian-Carboniferous transitional interval (Fig. 8). A study [21] demonstrated that a relatively rapid increase in sea level may have caused an estuary-like pattern of marine circulation between the intracratonic Williston Basin and open-ocean conditions, wherein nutrient-rich upwelling seawater was injected into the semi-enclosed basin at intermediate depths(100-200 m), leading to the proliferation of phytoplankton and high primary productivity.

    Volcanic and hydrothermal activities, which are closely associated with tectonic events,are also important for EHOM accumulation(Figs.6,7,9,and 10).Volcanic and hydrothermal activities can enhance primary productivity via two pathways: Volcanic ash deposition and hydrothermal fluids release large amounts of nutrients directly into seawater,and volcanic eruptions and hydrothermal activities release large amounts of CO2,resulting in high partial pressure of CO2(pCO2), eventually inducing climate warming and intense terrestrial chemical weathering, accelerating the influx of nutrients into oceans and lakes [104-108]. The marine deposition in the EHOM interval in the Wufeng Formation on the Yangtze Shelf Sea is related to volcanic eruptions (Fig. 6). Volcanic ash was extensively deposited in the Yangtze Shelf Sea of south China during the Ordovician-Silurian transition[18,23,104,109,110].Fertilization associated with volcanism during this period of active tectonics resulted in a significant increase in the content of bioavailable nitrogen and elemental nutrients (P, K, Si, Fe, and Zn) [104] in the Yangtze Shelf Sea. These processes enhanced the primary productivity and eventually formed EHOM deposits in the upper part of the Wufeng Formation (Fig. 6). The deposition of the EHOM interval at the bottom of the Marcellus shale in the Appalachian Basin further strengthened this relationship (Fig. 7).

    The EHOM interval(TOC up to 30 wt%)of the Yanchang Formation in the Ordos Basin was deposited in an ancient lacustrine environment(Fig.9[23]),closely related to both volcanic eruptions and hydrothermal activity,and probably driven by the Qinling orogenic event [111]. The EHOM deposits in the Chang 7 member contain 150 ash layers, with a single layer being up to 1 m thick [112].Hydrothermal activity during the Chang 7 depositional period was established by several petrological and geochemical characteristics, such as deposits of bedded chert and laminar ankerite,abnormally high abundances of Fe and Mn, and the heavy sulfur isotope composition of framboidal pyrite [111]. Hydrothermal input contributes to primary productivity via significant nutrient input(N,P,Fe,Si,and Zn).This can lead to the enhanced preservation of OM through the degassing of H2S and SO2, resulting in water column stratification with anoxia at depth [113]. Moreover,the peak volcanic and hydrothermal activity occurred during the Chang 73sub-member depositional period, which corresponds well to the deposition of the EHOM [105] (Fig. 9).

    Additionally, biotic extinction/radiation events may also be linked to increased primary productivity via the significant restructuring of marine/lacustrine ecosystems and their associated food webs [19]. Such an extinction event, the ‘‘Hangenberg Crisis,”was recorded in the Devonian-Carboniferous Bakken Formation.This extinction event affected both marine and terrestrial ecosystems throughout 100-300 kyr, resulting in the extinction of about 16% of all marine families and about 21% of all marine genera, as well as the global disappearance of Archaeopteris forests on land[19](Fig.8).Two significant second/third-order extinctions(Annulata Event and Dasberg Event) occurred in the lower and middle parts of Lower Bakken [114]. These extinction events were associated with surface seawater warming and increased nutrient supply caused by sea level rise,leading to the proliferation of primary producers and thus EHOM accumulation in the Lower and Upper Bakken formations in the Williston Basin[83,114].The Late Ordovician mass extinction(ca.445 Ma),the first of the‘‘Big Five”P(pán)hanerozoic extinction events [115], resulted in the extinction of marine benthic animals and nekton species, and survivors, including some graptolites, flourished because of reduced competition and increased niche space [116]. Gradual global warming and the removal of higher predators during the extinction event were conducive to the widespread proliferation of plankton and algae,probably resulting in the widespread deposition of EHOM and graptolite-rich intervals in the lower part of the Longmaxi Formation [23] (Fig. 6).

    4.2.2. Extensive anoxia facilitates OM preservation

    OM enrichment is related to the high primary productivity of OM and the effective preservation of OM due to anoxic or euxinic bottom water conditions [117]. The Late Ordovician to Early Silurian interval witnessed two global anoxic/euxinic events (Fig. 6)that contributed to the global deposition of EHOM(North America,Europe, and south China). The first of these events occurred from the Late Katian to Early Hirnantian, mainly influenced by deepwater shelf settings, and the second event was related to a rapid sea level rise and the upwelling of deep euxinic waters from the mid-late Hirnantian to the Early Rhudanian, causing distribution throughout the shelf (including the shallow-water inner shelf)[17]. The deposition of the EHOM intervals of Wufeng-Longmaxi shale in the Yangtze Shelf Sea corresponded well to these two global euxinic events (Fig. 6). The EHOM interval in the Wufeng Formation during the first euxinic event has a relatively limited spatial distribution,whereas the EHOM interval related to the second euxinic event is spatially more extensive [17,23]. These two euxinic events also promoted nutrient element recycling (e.g.,phosphorus) [118], thereby further enhancing primary productivity.

    A globally widespread oceanic anoxic event also occurred during the Devonian-Carboniferous transition [19] and was recorded in the Lower Bakken black shales of the Williston Basin[89,99].These shales are characterized by high OM content,enrichment of pyrite,and an absence of benthic fauna or related ichnofossils(Fig.6).Similar black shales were also widely deposited in other contemporaneous basins in Europe, North America, southeastern Asia, North Africa, Russia, Thailand, and south China (e.g., Ref. [119]). Black shales covered about 21% of the global depositional area during the Late Famennian[120].Although the cause of this oceanic anoxic event remains controversial, the widespread deposition of black shale corresponds well to a high sea level and may be associated with an expansion of anoxic water onto shelves owing to sea level rise under redox-stratified global ocean conditions[117].

    Bottom water anoxia during the deposition of Marcellus EHOM was likely caused by multiple factors (Fig. 7). In the early stage of the deposition of the Marcellus Formation, bottom water anoxia was caused by the weakening or stagnation of vertical ocean circulation,which may be related to factors such as the semi-closed nature of the Appalachian epicontinental sea, rapid tectonic subsidence, sea level rise, seasonal temperature stratification of the water column, and the rain shadow effect [74,77,121]. Subsequently, primary productivity increased, owing to a massive input of terrestrial nutrients, increasing the sinking flux of OM. OM decomposition resulted in water column oxygen depletion,leading to reducing conditions in the water column and underlying sediment pore waters [52,77]. This phenomenon is also known as the ‘‘productivity-anoxia feedback” mechanism [15,122].

    Fig. 10. Schematic of a typical sedimentary route for unconventional petroleum systems. Nutrients supplied by three sources—rivers, wind, and upwelling—are the key determinants of primary productivity in surface waters. Volcanism and hydrothermal activities can enhance the influx of nutrients into oceans and lakes. Anoxic or even euxinic bottom waters are favorable for the preservation of EHOM.Rapid or sudden changes in the sedimentary system associated with tectonic events,climate and sea level changes, and subaqueous gravity flow activity yield juxtaposed reservoirs. Many organic pores generated during the thermal maturation of the EHOM deposit provide abundant reservoir space for hydrocarbons.Abundant intergranular pores developed in the silty laminae associated with short-term geological events significantly improve the reservoir quality of mudrocks. Scanning electron microscopy images are from Yanchang shale (left) and Wufeng-Longmaxi shale (right) in China.

    4.2.3. Low sedimentation rate reduces dilution of sediments

    Enhanced clastic inputs and associated high sedimentation rates (>100 cm·kyr-1) generally dilute the OM content in sediments. In the Marcellus shale, the terrigenous input gradually increased(based on increased Al and Ti contents),and the OM content decreased from the Union Springs member to the Oatka Greek member (Fig. 7). Most ancient OM enrichment occurs in environments with low sedimentation rates. The average sedimentation rates for marine EHOM sediments are much lower than 100 cm·kyr-1: ranging from 0.2 to 0.4 cm·kyr-1(average TOC =2-8 wt%) in the Wufeng-Longmaxi Formation [5], 0.1 to 0.3 cm·kyr-1(average TOC = 8-11 wt%) in the Bakken Formation[21],and 0.17 to 1.3 cm·kyr-1(average TOC=8-10 w%)in the basal Marcellus Formation[73,123].The sedimentation rate of lacustrine EHOM sediments was close to 5 cm·kyr-1based on the sedimentation rate of the Chang 73black shale, yielding an average TOC concentration of 15 wt%[100,112].Changes in sedimentation rates resulting from terrigenous clastic sediment inputs have a major impact on the dilution or condensation of OM,and their processes are mainly controlled by the coupling of major geological events,including tectonic activity, climate change, and sea level fluctuation.

    4.3. Geological events that induced the formation of high-quality unconventional reservoirs

    4.3.1. Dramatic changes in depositional environments promote the formation of tight reservoirs

    The Bakken Formation and the Chang 7 member of the Yanchang Formation are representative strata of marine and lacustrine tight oil/shale oil, respectively. The features they have in common include closely coexisting source and reservoir rocks and oil accumulation without long-distance migration. Nearby sandstone and carbonate reservoirs have significantly higher matrix porosities than shale reservoirs and are the main oil-producing intervals[20,124](Figs.8 and 9).These reservoir rocks were deposited during dramatic changes in the depositional environment because of the coupling of sedimentology and geological events, such as climate cooling,rapid sea level fall,and subaqueous sediment gravity flow[97,125],providing abundant reservoir space for hydrocarbon accumulation [126].

    The Bakken Formation comprises two world-class petroleum source rock members (i.e., the Lower and Upper Bakken) and two tight reservoir members(the Middle Bakken and Pronghorn members). Hydrocarbons generated from the Upper and Lower Bakken organic-rich shales primarily charged the Middle Bakken and—locally—the Pronghorn member [93] (Fig. 8). Fine sandstone,bioclastic grainstone, and microcrystalline dolostone intervals within the Middle Bakken member have high matrix porosity(mean value of about 9%,maximum of 14%)and are promising targets for tight oil exploration and development in the Bakken play[127].Sedimentological studies indicate that the Upper and Lower Bakken were deposited under relatively deep marine anoxic conditions (>200 m depth), whereas the Middle Bakken was deposited under shallow-water high-energy conditions (< 10 m depth)[124].Dramatic changes in the depositional environment between adjacent members in the Bakken Formation are closely related to the climate-driven sea level fluctuation, reflected by lithofacies variations from deep-water fine-grained sediments(Lower Bakken member) with a low sedimentation rate (0.1 to 0.3 cm·kyr-1) to shallow-water coarse clastic deposits (Middle Bakken member)with a high sedimentation rate [97,124]. Subsequently, sea level rise resulted in the slow deposition of deep-water fine-grained lithofacies (Upper Bakken member) (Fig. 8). During the Devonian-Carboniferous transition, abrupt climate cooling (and sea level fall)may be linked to high organic carbon burial in the ocean(represented by widespread black shale deposition)and the proliferation of land plants, which enhanced chemical weathering and accelerated CO2drawdown from the atmosphere[19].Subsequent global warming(and attendant sea level rise)may be related to an increase in atmospheric CO2, linked to volcanic activity, increased carbonate weathering, and OM oxidation during sea level fall[128].

    The formation of lacustrine tight reservoirs of the Chang 7 member in the Ordos Basin is related to gravity flow activities(Fig. 9). Some geological events, such as volcanic eruptions, earthquakes, and storm waves, can trigger gravitational destabilization and the sliding of earlier sediments to form subaqueous sediment gravity flows,which can transport large amounts of sediment into semi-deep and deep lakes to form sandstone reservoirs [6,125].These deep-water sand bodies are laterally continuous and in close vertical proximity to EHOM shales. This pattern of source-reservoir contact can shorten the migration distance of oil and improve the efficiency of oil accumulation [6,129]. Although these gravity flow sandstone reservoirs are subject to strong diagenetic alteration (e.g., compaction and cementation) and represent typical tight reservoirs, they are relatively high-quality reservoirs with an average porosity of 4%-8% and permeability of 0.02-0.10 mD,which is significantly better than those of mudstones/shales with an average porosity less than 2.5% and permeability less than 0.01 mD [130]. Therefore, tight oil sweet-spot intervals within the Chang 7 member are dominated by these sandstones (Figs. 9 and 10), the formation and distribution of which are controlled by river-fed hyperpycnal flows and intrabasinal sediment gravity flows.

    4.3.2. Evolution of EHOM improves shale reservoir quality

    Shales, long regarded as hydrocarbon source rocks and seals,have become economic hydrocarbon reservoirs largely because of the abundance of pores in the OM.OM-hosted pores have irregular,bubble-like, elliptical shapes, and generally have sizes ranging from a few to several hundred nanometers (Fig. 3) [6,56]. These pores are widely developed in types II and I kerogen when the thermal maturity reaches an Roof about 0.6%or higher,accounting for up to 40%-50%of the total porosity and significantly improving shale reservoir quality [56]. The porosity of the sweet-spot intervals with EHOM in shales is comparatively high because the volume of OM pores is negligible in other intervals with TOC contents less than 3.0 wt%. For example, the porosity (≥4.0%) of the sweet-spot intervals in Wufeng-Longmaxi shale is generally higher than that of other intervals [6,56]. The enrichment of OM is conducive to increasing OM pores when the TOC content is less than 6 wt% [72,131]. However, when the TOC content is too high,the development of OM pores is inhibited because the lower structural rigidity of sediments with high TOC contents leads to the collapse and closure of OM pores during compaction [72,131].

    The OM in black shales plays an important role in controlling methane adsorption capacity and shale gas storage. Although the relationship between the total amount of methane adsorbed onto shale and its thermal maturity remains controversial,the adsorbed gas content in shales worldwide is directly proportional to their TOC content [132]. A stronger gas adsorption capacity of shale reservoirs implies a higher gas content, for example, the total gas content of the Wufeng-Longmaxi shale is positively correlated with its TOC content.Therefore,the evolution of the thermal maturity of the EHOM during diagenesis could improve the quality of shale reservoirs via two pathways: increasing total porosity and enhancing methane adsorption capacity.

    4.3.3. Well-developed silica-rich laminae allow efficient development of shale reservoir

    Shales usually exhibit heterogeneity and anisotropy, owing to their complex mineral composition and abundant lamina structure(Fig.1[32,33]).According to their mineral composition,laminae in shales can be classified into four types: organic-rich, clay-rich,carbonate-rich,and silica-rich[133].Silica-rich laminae are mainly composed of silt quartz and feldspar grains, also known as silty laminae(Fig.10)[134].The development of silty laminae improves the porosity and permeability of shale reservoirs and thus has a positive impact on the transport and exploitation of hydrocarbons[135]. Moreover, the high silica content makes the shale reservoir more brittle,which is conducive to the formation of complex fracture networks during reservoir fracturing reconstruction and improves the recovery efficiency of shale oil and gas [136,137].

    Silty laminae in shales originate from biogenic and abiogenic sources.Biogenic silty laminae are closely associated with the proliferation of siliceous organisms, such as radiolarians and sponge spicules in surface waters,which is caused by an increase in nutrient input as a likely consequence of an increase in sea level or volcanic and hydrothermal activities [14]. A high abundance of siliceous organisms has led to the deposition of numerous organic-rich black shales worldwide (e.g. Refs. [24,136]). Conversely,abiogenic silty laminae represent short-term and relatively high siliciclastic mineral deposition rates associated with turbidity currents, hyperpycnal flows, or episodic inputs of volcanic ash[138].Changes in the depositional setting and sedimentation rates accompanying these short-term geological events are conducive to the development of laminae in shales/mudstones. Thus, the response of geological events in sediments may improve the shale reservoir quality through various mechanisms,including increased organic pore development,reservoir adsorption capacity,reservoir percolation capacity, and brittleness. All these processes have important implications for the accumulation and industrial recovery of unconventional oil and gas.

    5. Conclusions

    This study provides a holistic review of the deposition and diagenesis of fine-grained sediments in unconventional petroleum systems and discusses the role of coupled sedimentation of major geological events in controlling the formation of unconventional petroleum sweet-spot intervals(areas)by considering four unconventional plays (with Ordovician to Triassic ages) from China and North America.

    The following conclusions were drawn:

    (1) The sweet-spot intervals (areas) of unconventional petroleum exploration and development are generally composed of EHOM deposits(TOC ≥3 wt%)or adjacent tight sandstone and carbonate rocks. EHOM deposits in unconventional petroleum systems provide not only a prerequisite for the massive generation and accumulation of hydrocarbons by OM enrichment during deposition but also abundant storage space for hydrocarbons by the evolution of OM-hosted pores and dissolution pores of minerals during diagenesis. Although tight sandstone reservoirs, including siltstones and fine sandstones, can form in a wide range of depositional environments,they are closely associated with EHOM deposits in shale strata because of dramatic changes in the depositional environment, such as deep-water gravity flows, providing high-quality reservoirs for unconventional petroleum resources.

    (2) The sweet-spot intervals of four typical examples of unconventional hydrocarbon reservoirs, namely the Wufeng-Longmaxi shale of the Sichuan Basin, Marcellus shale of the Appalachian Basin, Bakken Formation of the Williston Basin, and Yanchang Formation of the Ordos Basin, were identified and discussed to examine the roles of major geological events in their formation.The results established that their formation is the result of the spatial and temporal coupled sedimentation of global or regional geological events, such as tectonic activity, sea level (lake level)fluctuation,climate change,bottom water anoxia,volcanic activity,biotic mass extinction or radiation, and gravity flows during key geological transitional periods.

    (3) The coupled geological events collectively provide the fundamental geological settings for EHOM enrichment, including a high nutrient influx, enhanced bottom water anoxia/euxinia, and suitable sedimentation rate. Increased nutrient influx is mainly derived from enhanced weathering on land, the upwelling of nutrient-replete deep-water masses and hydrothermal fluids closely related to tectonic activity, climate change, sea level (lake level)fluctuations,or volcanic activity,which promote the proliferation of primary producers and further increase primary productivity. Extensive bottom water anoxia/euxinia facilitates OM preservation by decreasing OM decomposition. Suitable sedimentation rates controlled by the coupling of major geological events,including tectonic activity, climate change, and sea level fluctuations, are favorable for the enrichment of OM as they reduce the dilution of sediments.

    (4) Coupled geological events provide favorable depositional environments for the formation of high-quality unconventional reservoirs. Dramatic environmental changes driven by major geological events, such as deep-water gravity flows, are conducive to the deposition of high-quality tight oil/gas reservoirs closely associated with EHOM deposits in shale strata.The deposition of EHOM intervals related to geological events also promotes shale oil/gas reservoir quality through the evolution of OM-hosted pores and dissolution pores of minerals during diagenesis. In addition, welldeveloped silica-rich laminae in shales related to the proliferation of siliceous organisms or bottom currents can cause shale reservoirs to become brittle,allowing the efficient development of shale oil/gas.

    (5)Unconventional petroleum sedimentology,which focuses on coupled sedimentation during dramatic environmental changes driven by major geological events, is key to improve the understanding of the formation and distribution of unconventional hydrocarbon sweet-spot intervals (areas). Further research should aim to define new perspectives for the future exploration and development of unconventional petroleum resources.

    Acknowledgments

    This work was jointly supported by the Scientific Research and Technological Development Programs of CNPC (2021yjcq02 and 2021DJ2001). We thank Germán Otharán and anonymous reviewer(s) for helpful reviews and Editors Shasha Zhao and Lin Zhu for constructive comments.

    Compliance with ethical guidelines

    Caineng Zou,Zhen Qiu,Jiaqiang Zhang,Zhiyang Li,Hengye Wei,Bei Liu, Jianhua Zhao, Tian Yang, Shifa Zhu, Huifei Tao, Fengyuan Zhang,Yuman Wang,Qin Zhang,Wen Liu,Hanlin Liu,Ziqing Feng,Dan Liu,Jinliang Gao,Rong Liu,and Yifan Li declare that they have no conflicts of interest or financial conflicts to disclose.

    Appendix A. Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.eng.2022.06.016.

    日韩欧美一区二区三区在线观看| 神马国产精品三级电影在线观看 | 国产麻豆69| 久久人妻熟女aⅴ| 亚洲久久久国产精品| 天天躁狠狠躁夜夜躁狠狠躁| 精品福利观看| 曰老女人黄片| 国语自产精品视频在线第100页| 天天一区二区日本电影三级 | 美女高潮到喷水免费观看| 亚洲第一青青草原| 色综合亚洲欧美另类图片| 巨乳人妻的诱惑在线观看| 97人妻精品一区二区三区麻豆 | 日本五十路高清| 91精品三级在线观看| 亚洲一码二码三码区别大吗| 久久久久久亚洲精品国产蜜桃av| 色综合婷婷激情| 亚洲成人久久性| 天堂动漫精品| 黄色片一级片一级黄色片| 脱女人内裤的视频| 精品电影一区二区在线| 90打野战视频偷拍视频| 女人被狂操c到高潮| 9191精品国产免费久久| 午夜福利18| 亚洲av电影在线进入| 国产欧美日韩综合在线一区二区| 日日摸夜夜添夜夜添小说| 欧美亚洲日本最大视频资源| 亚洲精品国产精品久久久不卡| 成年人黄色毛片网站| av天堂久久9| 久热爱精品视频在线9| 午夜a级毛片| 亚洲七黄色美女视频| 免费不卡黄色视频| netflix在线观看网站| 久久欧美精品欧美久久欧美| or卡值多少钱| 欧美一级毛片孕妇| 国产视频一区二区在线看| 日韩欧美一区视频在线观看| 色精品久久人妻99蜜桃| 成人18禁高潮啪啪吃奶动态图| av免费在线观看网站| 在线观看一区二区三区| 一区二区三区精品91| 成人18禁高潮啪啪吃奶动态图| 亚洲成人免费电影在线观看| 亚洲精品久久成人aⅴ小说| 男女下面插进去视频免费观看| 久久久久国产一级毛片高清牌| 老司机深夜福利视频在线观看| 国产亚洲精品第一综合不卡| 欧美成人性av电影在线观看| 在线播放国产精品三级| 精品少妇一区二区三区视频日本电影| 亚洲成av片中文字幕在线观看| 亚洲成av人片免费观看| 亚洲情色 制服丝袜| 午夜两性在线视频| 欧美日本中文国产一区发布| 国产精品秋霞免费鲁丝片| 成人精品一区二区免费| 亚洲成人精品中文字幕电影| 夜夜爽天天搞| 欧美一级a爱片免费观看看 | 亚洲熟女毛片儿| 久久天堂一区二区三区四区| 黄片大片在线免费观看| 99国产精品99久久久久| 人人妻人人澡欧美一区二区 | 久久香蕉精品热| 1024香蕉在线观看| 18禁裸乳无遮挡免费网站照片 | 91九色精品人成在线观看| 久久久久亚洲av毛片大全| 大陆偷拍与自拍| 性少妇av在线| 亚洲国产精品久久男人天堂| 在线观看www视频免费| 亚洲熟妇中文字幕五十中出| 精品国产一区二区久久| 香蕉久久夜色| 成人三级黄色视频| 日本 欧美在线| 国产片内射在线| 亚洲熟女毛片儿| 国产一区二区三区综合在线观看| 日本vs欧美在线观看视频| 精品久久久精品久久久| 午夜福利18| 国产一区二区三区在线臀色熟女| 老鸭窝网址在线观看| 给我免费播放毛片高清在线观看| 久久久久久亚洲精品国产蜜桃av| 国产精品永久免费网站| 色尼玛亚洲综合影院| 国产精品一区二区免费欧美| 91精品国产国语对白视频| 国产激情欧美一区二区| 欧美日韩瑟瑟在线播放| 国产精品亚洲一级av第二区| 大码成人一级视频| 正在播放国产对白刺激| 18禁黄网站禁片午夜丰满| 首页视频小说图片口味搜索| www国产在线视频色| 欧美最黄视频在线播放免费| 久久国产精品影院| 亚洲视频免费观看视频| 亚洲av第一区精品v没综合| 欧美丝袜亚洲另类 | 免费看十八禁软件| 看黄色毛片网站| 亚洲少妇的诱惑av| 在线观看午夜福利视频| 黄色毛片三级朝国网站| 男女下面进入的视频免费午夜 | 极品人妻少妇av视频| 一本综合久久免费| 久久中文看片网| 亚洲色图 男人天堂 中文字幕| 国产精品电影一区二区三区| 亚洲国产精品999在线| 久久久久久久精品吃奶| 成人欧美大片| 精品一品国产午夜福利视频| 中文字幕人妻熟女乱码| 咕卡用的链子| 国产高清激情床上av| 黄片大片在线免费观看| 美国免费a级毛片| 国产一区在线观看成人免费| 成年人黄色毛片网站| 热re99久久国产66热| 日本欧美视频一区| 男人操女人黄网站| 香蕉久久夜色| 午夜免费观看网址| 国内精品久久久久久久电影| 一级a爱视频在线免费观看| 男女床上黄色一级片免费看| 老司机午夜十八禁免费视频| 亚洲专区国产一区二区| 黄色视频,在线免费观看| 精品久久久久久成人av| 99久久久亚洲精品蜜臀av| 国产一区二区在线av高清观看| 免费不卡黄色视频| 亚洲欧美精品综合一区二区三区| 亚洲无线在线观看| www.熟女人妻精品国产| 日本五十路高清| 久久久久久久久久久久大奶| 真人做人爱边吃奶动态| 他把我摸到了高潮在线观看| 很黄的视频免费| 亚洲一区中文字幕在线| 国产一区二区在线av高清观看| 亚洲一区高清亚洲精品| 激情视频va一区二区三区| 成人亚洲精品一区在线观看| 中文字幕最新亚洲高清| 亚洲欧美日韩无卡精品| 后天国语完整版免费观看| 亚洲一区二区三区不卡视频| 国产成年人精品一区二区| 亚洲人成伊人成综合网2020| 国产视频一区二区在线看| 欧美激情 高清一区二区三区| 亚洲一码二码三码区别大吗| 亚洲精品国产一区二区精华液| 成人精品一区二区免费| 天天躁夜夜躁狠狠躁躁| 女性生殖器流出的白浆| 久9热在线精品视频| 亚洲男人的天堂狠狠| 免费一级毛片在线播放高清视频 | 一级片免费观看大全| 97超级碰碰碰精品色视频在线观看| 成人18禁高潮啪啪吃奶动态图| 熟妇人妻久久中文字幕3abv| 欧美日本亚洲视频在线播放| 极品教师在线免费播放| aaaaa片日本免费| 成人18禁高潮啪啪吃奶动态图| 脱女人内裤的视频| 欧美另类亚洲清纯唯美| 999久久久国产精品视频| 黄色 视频免费看| 麻豆久久精品国产亚洲av| 精品高清国产在线一区| 成人精品一区二区免费| 久99久视频精品免费| 久久久水蜜桃国产精品网| 午夜两性在线视频| 色哟哟哟哟哟哟| 女人高潮潮喷娇喘18禁视频| 亚洲成人国产一区在线观看| 精品国产美女av久久久久小说| 欧美激情 高清一区二区三区| 手机成人av网站| 中文亚洲av片在线观看爽| 中文字幕最新亚洲高清| www.熟女人妻精品国产| 最近最新中文字幕大全电影3 | 国产精品1区2区在线观看.| av天堂在线播放| 丰满的人妻完整版| 久久精品国产亚洲av高清一级| 国产国语露脸激情在线看| 91在线观看av| 色在线成人网| 大型av网站在线播放| 久久伊人香网站| 免费无遮挡裸体视频| 两个人免费观看高清视频| 制服丝袜大香蕉在线| 很黄的视频免费| 波多野结衣巨乳人妻| 少妇粗大呻吟视频| 亚洲激情在线av| 亚洲黑人精品在线| 国产区一区二久久| 咕卡用的链子| 亚洲国产精品999在线| 搡老妇女老女人老熟妇| 日本vs欧美在线观看视频| 韩国av一区二区三区四区| 成人三级黄色视频| 久久久久久久精品吃奶| 多毛熟女@视频| 国产亚洲av高清不卡| 精品国产美女av久久久久小说| 久久天躁狠狠躁夜夜2o2o| 日韩欧美国产在线观看| 搡老岳熟女国产| 成年人黄色毛片网站| 国产精品乱码一区二三区的特点 | 亚洲国产日韩欧美精品在线观看 | 亚洲精品久久国产高清桃花| 欧美性长视频在线观看| 久久久国产精品麻豆| 国产精品 欧美亚洲| 18禁观看日本| 中国美女看黄片| 亚洲午夜理论影院| 天堂√8在线中文| 成熟少妇高潮喷水视频| 婷婷精品国产亚洲av在线| 日韩大码丰满熟妇| 亚洲最大成人中文| 久久久精品欧美日韩精品| www.熟女人妻精品国产| 亚洲精品在线观看二区| 国产熟女午夜一区二区三区| 777久久人妻少妇嫩草av网站| 欧美日韩中文字幕国产精品一区二区三区 | 国产极品粉嫩免费观看在线| 一区二区三区激情视频| 大型av网站在线播放| 看片在线看免费视频| 国产免费av片在线观看野外av| 国产成人一区二区三区免费视频网站| 黄网站色视频无遮挡免费观看| 久久中文字幕人妻熟女| 亚洲情色 制服丝袜| 日本欧美视频一区| av片东京热男人的天堂| 男人舔女人的私密视频| 十八禁人妻一区二区| 欧美 亚洲 国产 日韩一| 国产精品一区二区三区四区久久 | 精品一区二区三区av网在线观看| 国产欧美日韩一区二区三| 在线av久久热| 日韩欧美国产一区二区入口| 可以在线观看的亚洲视频| 国产精品亚洲一级av第二区| 一级,二级,三级黄色视频| 中文字幕精品免费在线观看视频| 久久精品人人爽人人爽视色| 50天的宝宝边吃奶边哭怎么回事| 成人三级黄色视频| 首页视频小说图片口味搜索| 最好的美女福利视频网| 在线观看日韩欧美| 黄网站色视频无遮挡免费观看| 69av精品久久久久久| 久久精品国产综合久久久| 88av欧美| 丰满人妻熟妇乱又伦精品不卡| 欧美av亚洲av综合av国产av| 又黄又粗又硬又大视频| 少妇粗大呻吟视频| 麻豆国产av国片精品| 黄色成人免费大全| 男女下面插进去视频免费观看| 久久青草综合色| 他把我摸到了高潮在线观看| 在线播放国产精品三级| 午夜福利一区二区在线看| 99久久国产精品久久久| 好男人电影高清在线观看| 给我免费播放毛片高清在线观看| 在线观看免费视频日本深夜| 1024视频免费在线观看| 在线观看66精品国产| 日日干狠狠操夜夜爽| 色播在线永久视频| 成人免费观看视频高清| 一区二区三区国产精品乱码| 国产亚洲欧美在线一区二区| 久久国产亚洲av麻豆专区| 国产欧美日韩一区二区三| 人妻久久中文字幕网| 男女午夜视频在线观看| 亚洲免费av在线视频| 久久人人97超碰香蕉20202| 成人手机av| 国产精品一区二区精品视频观看| av超薄肉色丝袜交足视频| 搞女人的毛片| 男女之事视频高清在线观看| 无限看片的www在线观看| 90打野战视频偷拍视频| 黄片播放在线免费| 午夜影院日韩av| av视频在线观看入口| 熟女少妇亚洲综合色aaa.| 黑人巨大精品欧美一区二区mp4| 国产欧美日韩一区二区三| 中文字幕另类日韩欧美亚洲嫩草| 一区在线观看完整版| 91精品国产国语对白视频| 国产精品99久久99久久久不卡| 国产精品亚洲av一区麻豆| 日韩一卡2卡3卡4卡2021年| 色综合婷婷激情| 女同久久另类99精品国产91| 午夜福利,免费看| 欧美人与性动交α欧美精品济南到| 亚洲五月色婷婷综合| 1024香蕉在线观看| 亚洲色图 男人天堂 中文字幕| 窝窝影院91人妻| 午夜福利18| 欧美日韩一级在线毛片| 欧美不卡视频在线免费观看 | 免费高清在线观看日韩| 色综合欧美亚洲国产小说| 精品国产亚洲在线| 欧美日韩中文字幕国产精品一区二区三区 | 久久久久久大精品| 国产精品秋霞免费鲁丝片| 男男h啪啪无遮挡| 此物有八面人人有两片| а√天堂www在线а√下载| 一卡2卡三卡四卡精品乱码亚洲| 91精品国产国语对白视频| 在线观看免费视频日本深夜| 久久人人精品亚洲av| 99国产精品免费福利视频| 色播亚洲综合网| 两人在一起打扑克的视频| 久久国产精品人妻蜜桃| 超碰成人久久| www.自偷自拍.com| 琪琪午夜伦伦电影理论片6080| 亚洲成av人片免费观看| 久久人人爽av亚洲精品天堂| 久久精品91蜜桃| 黄色视频,在线免费观看| 久久午夜综合久久蜜桃| 久久狼人影院| 国产97色在线日韩免费| 午夜精品国产一区二区电影| 嫁个100分男人电影在线观看| 久久精品亚洲精品国产色婷小说| 国产午夜精品久久久久久| 日韩欧美在线二视频| 久久欧美精品欧美久久欧美| 欧美色视频一区免费| 色婷婷久久久亚洲欧美| 色综合站精品国产| 日本在线视频免费播放| 国产在线观看jvid| 午夜福利成人在线免费观看| 国产欧美日韩一区二区三| 久久人妻av系列| 欧美黑人欧美精品刺激| 国产精品电影一区二区三区| 成在线人永久免费视频| 欧美乱码精品一区二区三区| 性色av乱码一区二区三区2| 少妇熟女aⅴ在线视频| 免费搜索国产男女视频| 亚洲人成网站在线播放欧美日韩| 欧美成狂野欧美在线观看| 麻豆久久精品国产亚洲av| 久99久视频精品免费| 国产精品久久久久久亚洲av鲁大| 中文字幕色久视频| 少妇被粗大的猛进出69影院| 高潮久久久久久久久久久不卡| 亚洲熟妇熟女久久| 女人爽到高潮嗷嗷叫在线视频| 男女之事视频高清在线观看| 一二三四在线观看免费中文在| 高清毛片免费观看视频网站| 国产高清视频在线播放一区| 夜夜夜夜夜久久久久| 无遮挡黄片免费观看| 午夜福利成人在线免费观看| 亚洲成人国产一区在线观看| 日日摸夜夜添夜夜添小说| 久久久精品欧美日韩精品| 亚洲自偷自拍图片 自拍| 人成视频在线观看免费观看| tocl精华| 老司机在亚洲福利影院| 两性夫妻黄色片| 女人被躁到高潮嗷嗷叫费观| 国产区一区二久久| 日韩欧美一区视频在线观看| 50天的宝宝边吃奶边哭怎么回事| or卡值多少钱| 欧美精品啪啪一区二区三区| 99riav亚洲国产免费| 亚洲精品中文字幕在线视频| 日本vs欧美在线观看视频| 亚洲第一av免费看| 最近最新中文字幕大全电影3 | 黑丝袜美女国产一区| 国产精品乱码一区二三区的特点 | 91字幕亚洲| 97超级碰碰碰精品色视频在线观看| 日本欧美视频一区| 最新美女视频免费是黄的| 日韩三级视频一区二区三区| 国产精品亚洲美女久久久| 免费人成视频x8x8入口观看| 999久久久精品免费观看国产| 99久久国产精品久久久| 欧美在线黄色| 久久人人爽av亚洲精品天堂| 亚洲狠狠婷婷综合久久图片| 曰老女人黄片| 好男人在线观看高清免费视频 | 久久中文字幕一级| 亚洲欧美精品综合一区二区三区| 电影成人av| 欧美 亚洲 国产 日韩一| 亚洲七黄色美女视频| 麻豆成人av在线观看| 国产精华一区二区三区| 亚洲伊人色综图| 首页视频小说图片口味搜索| 国产在线观看jvid| 成年女人毛片免费观看观看9| 欧美日本亚洲视频在线播放| 色尼玛亚洲综合影院| 婷婷六月久久综合丁香| 亚洲国产精品久久男人天堂| 在线视频色国产色| 国产区一区二久久| 国产亚洲精品第一综合不卡| x7x7x7水蜜桃| 国产精品影院久久| 精品无人区乱码1区二区| 欧美一区二区精品小视频在线| 日韩高清综合在线| a在线观看视频网站| 国产亚洲欧美精品永久| 午夜免费成人在线视频| 99久久国产精品久久久| www.自偷自拍.com| 欧美日韩福利视频一区二区| 国产一级毛片七仙女欲春2 | av视频在线观看入口| 曰老女人黄片| 中亚洲国语对白在线视频| 久热爱精品视频在线9| 嫁个100分男人电影在线观看| 精品国产亚洲在线| 精品熟女少妇八av免费久了| ponron亚洲| 欧美色欧美亚洲另类二区 | 久久天堂一区二区三区四区| 国产精品1区2区在线观看.| 嫩草影院精品99| 一级毛片精品| 久久久久久久午夜电影| 精品国产国语对白av| 久久中文字幕一级| 亚洲一区中文字幕在线| 日本 欧美在线| 少妇 在线观看| 高清毛片免费观看视频网站| 亚洲中文字幕一区二区三区有码在线看 | 黄片大片在线免费观看| 夜夜看夜夜爽夜夜摸| 国产av又大| 日韩免费av在线播放| 一进一出抽搐动态| 神马国产精品三级电影在线观看 | 国产一区二区在线av高清观看| 一级毛片女人18水好多| 99在线人妻在线中文字幕| 9色porny在线观看| 国产精品一区二区三区四区久久 | 欧美日韩亚洲国产一区二区在线观看| 久久婷婷人人爽人人干人人爱 | 深夜精品福利| 中文字幕av电影在线播放| 免费高清在线观看日韩| 757午夜福利合集在线观看| 久久天躁狠狠躁夜夜2o2o| 一区二区三区激情视频| 午夜精品在线福利| 久久香蕉精品热| 亚洲人成电影观看| 麻豆久久精品国产亚洲av| 日本一区二区免费在线视频| 人人妻,人人澡人人爽秒播| aaaaa片日本免费| 久久久国产欧美日韩av| 国产97色在线日韩免费| 99re在线观看精品视频| 亚洲精品国产一区二区精华液| 一a级毛片在线观看| 久久精品国产99精品国产亚洲性色 | 日本 av在线| 久久久久久久久免费视频了| 国产精品日韩av在线免费观看 | 精品国产乱码久久久久久男人| 欧美精品啪啪一区二区三区| 日韩大码丰满熟妇| 18禁观看日本| 中出人妻视频一区二区| 国产又爽黄色视频| 国产精品久久电影中文字幕| 亚洲片人在线观看| 午夜久久久在线观看| 黄色女人牲交| 9色porny在线观看| 19禁男女啪啪无遮挡网站| 国产精品亚洲美女久久久| 免费人成视频x8x8入口观看| 亚洲成人免费电影在线观看| 黄色片一级片一级黄色片| 亚洲国产精品999在线| 国产激情久久老熟女| 99热只有精品国产| 亚洲色图综合在线观看| 嫩草影院精品99| 免费无遮挡裸体视频| 久久午夜亚洲精品久久| 乱人伦中国视频| 欧美另类亚洲清纯唯美| 国产精品,欧美在线| 看免费av毛片| 午夜亚洲福利在线播放| 亚洲欧美日韩另类电影网站| 免费在线观看黄色视频的| 成人三级黄色视频| 国产aⅴ精品一区二区三区波| 国产亚洲欧美在线一区二区| 亚洲国产欧美网| 免费人成视频x8x8入口观看| 日本黄色视频三级网站网址| av福利片在线| 国产野战对白在线观看| 色尼玛亚洲综合影院| 变态另类丝袜制服| 久久精品亚洲熟妇少妇任你| 日韩高清综合在线| 97碰自拍视频| 亚洲欧美精品综合久久99| 99精品欧美一区二区三区四区| 一级毛片精品| 国产伦人伦偷精品视频| 亚洲专区字幕在线| 母亲3免费完整高清在线观看| 一级a爱视频在线免费观看| 欧美成狂野欧美在线观看| 伊人久久大香线蕉亚洲五| 午夜老司机福利片| 久久国产乱子伦精品免费另类| 老司机午夜福利在线观看视频| 一进一出好大好爽视频| 一级毛片精品| а√天堂www在线а√下载| 久热爱精品视频在线9| 99在线视频只有这里精品首页| ponron亚洲| 国产区一区二久久| 色综合亚洲欧美另类图片| 99国产精品99久久久久| 日本a在线网址| 国产视频一区二区在线看| 麻豆国产av国片精品| 村上凉子中文字幕在线| 日韩欧美在线二视频| 99精品欧美一区二区三区四区| 亚洲中文日韩欧美视频| av在线天堂中文字幕| 老司机在亚洲福利影院| 一区二区三区精品91| 99国产极品粉嫩在线观看|