Method for determining carrying path of black shale material source

Abstract

The invention discloses a method for determining a black shale material source carrying path, and relates to the technical field of shale oil and shale gas exploration and development engineering. According to the method, the contents of clay minerals and feldspar are specially selected to be combined with the granularity of shale particles to establish a series of plane contour maps, and through superposition of the plane contour maps, a carrying path with siliceous clastic rock as a main material source, a carrying path with carbonate rock as a main material source and a carrying path of gravity flow are judged; and the carrying path of the black shale of the target layer series in the research area is determined through superposition of the carrying paths. According to the method for determining the carrying path of the black shale material source, the ancient landform of the sedimentary basin can be effectively recovered, an important means is provided for optimization of a favorable shale oil and gas exploration and development area, and then the success rate of well drilling can be effectively increased.

Description

Technical Field

[0001] The present invention belongs to the technical field of shale oil and shale gas exploration and development engineering, and in particular relates to a method for determining the provenance transportation path of black shale. Background Art

[0002] Provenance analysis includes source area structural analysis, source rock component analysis, and source transport path analysis. Among them, transport path analysis is of great significance. It can not only indicate the basin transport and sedimentation process, but also has important significance for finding oil and gas reservoirs. For coarse-grained clastic rocks, people have always determined the material transport path by compiling conglomerate isopach maps, sandstone isopach maps, and sand-to-ground ratio contour maps (such as the literature: Feng Zengzhao. New Exploration of the Lithofacies Paleogeography of the Lower Ordovician in North China. Journal of East China Petroleum Institute, 1977, 19(3):57-79; Feng Zengzhao. New Exploration of the Lithofacies Paleogeography of the Early Ordovician in North China. Geological Science, 1979(4):302-313). This research method has been widely recognized by the academic community and has been widely used in major oil and gas provinces. However, for black shale, since the composition of the material is mainly fine-grained sediments, conventional particle mapping of various particle sizes is not feasible.

[0003] In recent years, some scholars have proposed a sedimentary microfacies mapping method based on single-factor analysis and multi-factor synthesis for black shale (e.g., Shi Zhensheng, Zhou Tianqi, Guo Wei, Liang Pingping, Cheng Feng. Quantitative paleogeographic mapping of marine shales and division of deep-water shelf sedimentary microfacies: A case study of the Longyi 11-4 sublayer of the Wufeng Formation-Longmaxi Formation in the Luzhou area of ​​southern Sichuan. Acta Sedimentologica Sinica, 2022, 40(6):1728-1744.). The research idea of ​​this method is to determine the planar distribution of each sedimentary microfacies by compiling a planar distribution map of black shale stratum thickness, carbonate mineral content, quartz content, and clay mineral content. The application of this method provides an important way to divide the sedimentary microfacies of black shale and compile a quantitative paleogeographic map of black shale. However, this method has the following problems: (1) The thickness of black shale formations is controlled by multiple factors such as provenance, accommodation space, and the supply of sedimentary materials. Its distribution may be inconsistent with the deposition center, thus bringing about multiple interpretations of sedimentary microfacies; (2) Quartz and carbonate minerals in black shale come from multiple sources, and the planar distribution patterns of minerals from different sources vary greatly, so there may be deviations in the identification and compilation of sedimentary microfacies.

[0004] Therefore, there is still no mature method for determining the provenance and transport path of black shale. To more efficiently identify favorable areas for shale oil and gas exploration and development and improve the success rate of exploration and development, this paper provides a method for determining the provenance and transport path of black shale by using shale mineral particle size and mineral component content. Summary of the Invention

[0005] The main purpose of the present invention is to provide a method for determining the provenance and transportation path of black shale, which can effectively restore the ancient landform of the sedimentary basin, so as to more efficiently find favorable areas for shale oil and gas exploration and development and improve the success rate of exploration and development.

[0006] To achieve the above objectives, the present invention provides a method for determining the provenance transport path of black shale, comprising the following steps:

[0007] Prepare the plane contour map of feldspar content, clay mineral content, average grain size of terrigenous quartz, average grain size of feldspar, average grain size of terrigenous carbonate minerals, and average grain size of clay minerals for the black shale of the target strata in the study area;

[0008] Based on the plane contour map of the average grain size of terrigenous quartz, the plane contour map of the average grain size of feldspar, and the plane contour map of feldspar content, it is determined that the transportation path is mainly sourced from siliciclastic rock; based on the plane contour map of the average grain size of terrigenous carbonate minerals, it is determined that the transportation path is mainly sourced from the carbonate rock; based on the plane contour map of the clay mineral content and the plane contour map of the average grain size of clay minerals, it is determined that the transportation path is a gravity flow;

[0009] By superimposing the transport path of the siliceous clastic rock as the main source, the transport path of the carbonate rock as the main source and the transport path of the gravity flow, the completely superimposed transport path is determined as the black shale source transport path.

[0010] The clay minerals and feldspar in black shale are mainly terrestrial in origin, and the changes in their composition and content are mainly controlled by the transport path. Clay minerals are mainly derived from unstable igneous rocks in terrestrial debris. As the transport distance increases, their content decreases, and the mineral composition also changes. Feldspar is also derived from terrestrial origin, and its content also decreases as the transport distance increases. The particle size of the debris particles in black shale is also mainly controlled by the transport distance. The larger the transport distance, the finer its particle size. Therefore, the present invention particularly selects the parameters of clay mineral and feldspar content in conjunction with the particle size of shale particles (terrestrial quartz average particle size, feldspar average particle size, terrestrial carbonate mineral average particle size and clay mineral average particle size) to establish a series of plane contour maps. By superposition of these plane contour maps, it is judged that the transport path with siliceous clastic rock as the main source, the transport path with terrestrial carbonate minerals as the main source, and the transport path of gravity flow are determined. Then, the transport path of the black shale of the target stratum in the study area is determined by the superposition of these transport paths. The method provided by the present invention for determining the provenance and transport path of black shale can effectively restore the ancient landform of the sedimentary basin, provide an important means for optimizing favorable areas for shale oil and gas exploration and development, and thus effectively improve the success rate of drilling.

[0011] Furthermore, the specific steps of determining the transport path of the siliciclastic rock as the main source are as follows: based on the plane contour map of the average grain size of the terrigenous quartz and the plane contour map of the average grain size of the feldspar, preliminarily determining the location where the average grain sizes of both are high as the source location of the siliciclastic rock;

[0012] Combined with the location of the high-value area on the feldspar content plane contour map, the source area dominated by siliciclastic rock was finally determined;

[0013] Based on the plane contour map of the average grain size of the terrigenous quartz and the feldspar, and the plane contour map of the feldspar content, the path of content and average grain size from high-value areas to low-value areas was determined to be the transportation path with the siliceous clastic rock as the main source.

[0014] Furthermore, the specific steps for determining the transportation path of the carbonate rock as the main source are: based on the plane contour map of the average grain size of the terrestrial carbonate minerals, the area with high average grain size is determined as the source area location dominated by carbonate rock, and the path from the high-value area to the low-value area of ​​the average grain size is determined as the transportation path of the carbonate rock as the main source.

[0015] Furthermore, the specific steps for determining the transport path of the gravity flow are: based on the plane contour map of clay mineral content, determining the high-value area of ​​clay mineral content as the source area of ​​the gravity flow; comprehensively considering the plane contour map of clay mineral content and average particle size, determining the path from the high-value area to the low-value area of ​​content and average particle size as the transport path of the gravity flow.

[0016] Furthermore, the specific steps of compiling the plane contour map of feldspar content and the plane contour map of clay mineral content are as follows: using mineral composition logging data to calculate the arithmetic mean of feldspar content and clay mineral content of each well in the target stratum;

[0017] Based on the location information of the research work area, a work area location map is compiled; the coordinates of the wells drilled into the target layer are imported into the work area location map to form a well location distribution map of the work area;

[0018] Importing the arithmetic mean data of the feldspar content of each well into the well location distribution map of the work area, and using the Kriging interpolation method to compile the feldspar content plane contour map;

[0019] The arithmetic mean value data of the clay mineral content of each well is imported into the well location distribution map of the work area, and the plane contour map of the clay mineral content is compiled by using the Kriging interpolation method.

[0020] Furthermore, the mineral composition logging data are the coordinates of all drilling wells that encounter the target layer in the study area and the corresponding mineral composition logging data that have been processed by data cleaning and outlier elimination.

[0021] Furthermore, the specific steps for compiling the plane contour map of the average grain size of terrigenous quartz are as follows: the average grain size data of terrigenous quartz of each well is imported into the well location distribution map of the work area, and the plane contour map of the average grain size of terrigenous quartz is compiled using the Kriging interpolation method.

[0022] Furthermore, the specific steps of compiling the plane contour map of the average particle size of feldspar are as follows: the average particle size data of feldspar of each well is imported into the well location distribution map of the work area, and the plane contour map of the average particle size of feldspar is compiled by using the Kriging interpolation method.

[0023] Furthermore, the specific steps for compiling the plane contour map of the average particle size of the terrestrial carbonate minerals are as follows: the average particle size data of the terrestrial carbonate minerals of each well is imported into the drilling well location distribution map of the work area, and the plane contour map of the average particle size of the terrestrial carbonate minerals is compiled using the Kriging interpolation method.

[0024] Furthermore, the specific steps of compiling the clay mineral average particle size plane contour map are as follows: importing the average particle size data of clay minerals of each well on the drilling well location distribution map of the work area, and using the Kriging interpolation method to compile the clay mineral average particle size plane contour map.

[0025] Furthermore, the average particle size of the terrigenous quartz, the feldspar, the clay minerals, and the terrigenous carbonate minerals is obtained by using the following method: using mineral composition logging data to determine the drilling depth of each well at which the clay mineral content is at its maximum, and for each well in the study area, collecting a core sample at the drilling depth at which the clay mineral content is at its maximum;

[0026] The core samples were made into argon ion polished slices for scanning electron microscopy analysis, and the particle size of the terrigenous quartz, feldspar, clay minerals and terrigenous carbonate minerals under the scanning electron microscopy images of the argon ion polished slices was analyzed to calculate the average particle size data of the terrigenous quartz, feldspar, clay minerals and terrigenous carbonate minerals for each well.

[0027] Furthermore, 3 to 5 core samples are evenly collected within the depth of the maximum distribution of clay mineral content.

[0028] Furthermore, the sample size of the core sample is not less than 1 cm×1 cm×1 cm.

[0029] Compared with the existing technology, the present invention has the following beneficial effects: the present invention specifically selects the clay mineral and feldspar content in combination with the parameters of the shale particle size (average particle size of terrigenous quartz, average particle size of feldspar, average particle size of terrigenous carbonate minerals, and average particle size of clay minerals) to establish a series of plane contour maps. By superimposing these plane contour maps, the transportation path with siliceous clastic rock as the main source, the transportation path with carbonate rock as the main source, and the transportation path of gravity flow are determined. Then, by superimposing these transportation paths, the transportation path of the black shale of the target stratum in the study area is determined. The method for determining the black shale provenance transportation path provided by the present invention can effectively restore the paleogeomorphology of the sedimentary basin, provide an important means for selecting favorable areas for shale oil and gas exploration and development, and thus effectively improve the success rate of drilling. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] Figure 1 A schematic flow chart of the method for determining the black shale provenance transport path of the present invention is shown;

[0031] Figure 2 The diagram shows the distribution of drilling well locations in the Q well area in the southern Sichuan Basin in an embodiment of the present invention;

[0032] Figure 3 The figure shows the distribution of feldspar content in the well area Q in the southern Sichuan Basin in the embodiment of the present invention;

[0033] Figure 4 The feldspar content plane contour map of the Q well area in the southern Sichuan Basin in an embodiment of the present invention is shown;

[0034] Figure 5 The figure shows the distribution of clay minerals in the well area Q in the southern Sichuan Basin according to the embodiment of the present invention;

[0035] Figure 6 The following is a planar contour map of clay mineral content in the Q well area in the southern Sichuan Basin according to an embodiment of the present invention;

[0036] Figure 7 A cross-sectional diagram of mineral composition logging data of a single well in the Q well area in the southern Sichuan Basin in an embodiment of the present invention is shown.

[0037] Figure 8 The figure shows the distribution of wells containing average grain size of terrigenous quartz in the Q well area in the southern Sichuan Basin in the embodiment of the present invention;

[0038] Figure 9 The plane contour map of the average grain size of terrigenous quartz in the Q well area in the southern Sichuan Basin in the embodiment of the present invention is shown;

[0039] Figure 10The figure shows the distribution of well locations containing average feldspar grain size in the Q well area in the southern Sichuan Basin in an embodiment of the present invention;

[0040] Figure 11 The figure shows the contour map of the average particle size of feldspar in the Q well area in the southern Sichuan Basin in the embodiment of the present invention;

[0041] Figure 12 The figure shows the distribution of the average particle size of clay minerals in the Q well area in the southern Sichuan Basin according to the embodiment of the present invention;

[0042] Figure 13 The plane contour map of average particle size of clay minerals in the Q well area in the southern Sichuan Basin according to the embodiment of the present invention is shown;

[0043] Figure 14 The figure shows the distribution of well locations containing average particle size of terrigenous carbonate minerals in the Q well area in the southern Sichuan Basin in an embodiment of the present invention;

[0044] Figure 15 The plane contour map of the average grain size of terrigenous carbonate minerals in the Q well area in the southern Sichuan Basin in an embodiment of the present invention is shown;

[0045] Figure 16 A schematic diagram showing the provenance transport path of the Q well area in the southern Sichuan Basin in an embodiment of the present invention is shown. DETAILED DESCRIPTION

[0046] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present invention will be described in detail below with reference to the embodiments.

[0047] Example

[0048] Taking the Q well area in the southern Sichuan Basin as an example, the logging data of the Long-1 layer of 20 wells were collected, and the samples were collected and analyzed according to the requirements of the present invention. The specific analysis process is as follows: Figure 1 As shown, the specific implementation steps are as follows:

[0049] 1. Compilation of plane contour maps of feldspar and clay mineral content

[0050] Step 1: Collection and processing of logging and drilling well location data.

[0051] (1) Collect and organize the well location coordinates and corresponding mineral composition logging data of 20 wells drilled into the Long-1 layer, and prepare corresponding statistical data tables.

[0052] (2) Preprocess the collected mineral composition logging data, including data cleaning and outlier removal, to ensure the data is accurate and reliable.

[0053] (3) Use Shiwen software to calculate the arithmetic mean of feldspar and clay minerals in the Long-1 layer, and prepare the corresponding statistical data table (as shown in Table 1).

[0054] Table 1. Arithmetic mean values ​​of well location coordinates, feldspar, and clay mineral logging data for 20 wells in the Q well area

[0055]

[0056]

[0057] Step 2: Preparation of plane contour map

[0058] (1) The basin boundary, location boundary, important fault distribution, and important geographical locations (including provincial capitals and county-level cities) of the research area should be clearly defined, and a location map of the work area should be generated using professional geological software (such as Surfer).

[0059] (2) Import the well location coordinates of the Yulong 1 sub-layer drilling into the work area location map and form a well location distribution map of the work area (such as Figure 2 shown).

[0060] (3) Import the arithmetic mean data of feldspar content (%) of each well into the well location distribution map of the work area to form a well location distribution map containing feldspar content (such as Figure 3 shown).

[0061] (4) Using the Kriging interpolation method, a plane contour map of the feldspar content of the Longyi 1 layer in the study area was compiled (e.g. Figure 4 shown).

[0062] (5) Import the arithmetic mean value of clay mineral content (%) of each well into the well location distribution map of the work area to form a well location distribution map containing clay mineral content (e.g. Figure 5 shown).

[0063] (6) Using the Kriging interpolation method, a plane contour map of clay mineral content in the Longyi 1 sub-stratum system in the study area was compiled (e.g. Figure 6 shown).

[0064] 2. Compilation of shale grain size contour map

[0065] Step 1: Sample selection and preparation

[0066] (1) Use mineral composition logging data to determine the drilling depth corresponding to the maximum clay mineral content of each well, such as Figure 7 The figure shows a profile of the mineral composition logging data of a single well. The dotted line in the figure indicates the drilling depth corresponding to the maximum clay mineral content, which is 3762 meters.

[0067] (2) For each well in the study area, four core samples were evenly selected within the drilling depth corresponding to the maximum clay mineral content, with a sample size of 1 cm × 1 cm × 1 cm.

[0068] (3) In order to obtain high-precision and large-field mineral and organic matter images, the research steps and methods of argon ion polishing slice production, image acquisition and splicing, and mineral and organic matter analysis were adopted. The specific process is as follows: the diameter of the large thin slice sample should be no less than 8 cm. The size of the argon ion polishing slice is 10 mm × 10 mm × 5 mm. The image acquisition is carried out using a Hitachi field emission scanning electron microscope with cold discharge, equipped with a low-to-high secondary electron probe and an X-ray energy dispersive spectrometer (EDS). The magnification of the scanning electron microscope is 30,000 times (the maximum resolution of a single photo is 9 nm). The image acquisition area is perpendicular to the grain surface, and the cumulative acquisition area is 60 μm × 40 μm. After the image acquisition is completed, Microsoft HDView software is used for mineral composition and organic matter analysis.

[0069] Step 2: Statistical analysis of shale particle size

[0070] (1) Observe the argon ion polished wafer using a field emission scanning electron microscope (SEM).

[0071] (2) Nano Measurer 1.2 software was used to calculate the average particle size of terrigenous quartz, feldspar, clay minerals, and terrigenous carbonate minerals in field emission scanning electron microscope (SEM) images, and a statistical table of different mineral particle sizes was generated (as shown in Table 2).

[0072] Table 2. Average particle size (μm) of terrigenous quartz, feldspar, clay minerals, and terrigenous carbonate minerals from 20 wells in the Q well area.

[0073] Table 2

[0074]

[0075] Step 3: Preparation of plane contour maps of average particle size of different mineral particles

[0076] (1) Import the average grain size data of terrigenous quartz of each well into the well location distribution map of the work area to form a well location distribution map containing the average grain size of terrigenous quartz (e.g. Figure 8 shown).

[0077] (2) Using the Kriging interpolation method, the average grain size contour map of terrigenous quartz in the Longyi 1 layer of the study area was compiled (e.g. Figure 9 shown).

[0078] (3) Import the average feldspar particle size data of each well into the well location distribution map of the work area to form a well location distribution map containing the average feldspar particle size (e.g. Figure 10 shown).

[0079] (4) Using the Kriging interpolation method, the average particle size contour map of the feldspar in the Longyi 1 layer of the study area was compiled (e.g. Figure 11 shown).

[0080] (5) Import the average particle size data of clay minerals of each well into the well location distribution map of the work area to form a well location distribution map containing the average particle size of clay minerals (e.g. Figure 12 shown).

[0081] (6) Using the Kriging interpolation method, the average particle size contour map of clay minerals in the Longyi 1 layer of the study area was compiled (e.g. Figure 13 shown).

[0082] (7) Import the average particle size data of terrigenous carbonate minerals of each well into the well location distribution map of the work area (e.g. Figure 14 shown).

[0083] (8) Using the Kriging interpolation method, a plane contour map of the average grain size of terrigenous carbonate minerals in the Longyi 1 layer of the study area was compiled (e.g. Figure 15 shown).

[0084] 3. Comprehensive determination of material source transportation routes

[0085] Step 1: Transportation path of siliceous clastic rock as the main source

[0086] (1) Analyze the average particle size contour map of terrigenous quartz and feldspar, identify the location of the high-value areas of the average particle size of the two, and preliminarily determine that the main source area is siliceous clastic rock. Figure 9 The high-value area of ​​average grain size of terrigenous quartz is the area where wells Q17 and Q18 are located, with average grain sizes of terrigenous quartz being 50.2 and 49.9, respectively.

[0087] (2) Based on the above step, combined with the location of the high-value area on the plane contour map of feldspar content, the main provenance area is finally determined to be siliceous clastic rock. Figure 4 The feldspar contents in wells Q17 and Q18 are 15.4 and 14.2, respectively, which are also in the high value area. Therefore, the areas where wells Q17 and Q18 are located are source areas dominated by siliceous clastic rocks.

[0088] (3) Comprehensive analysis of the average grain size and feldspar content contour maps of terrigenous quartz and feldspar to clarify the path of average grain size from high to low value areas, and preliminarily determine the transport path of siliciclastic rocks as the main source. These are the two provenance transport paths in the contour map: 1) the path from the area where well Q17 is located in the southwest to the direction of well Q15 in the northeast; and 2) the path from the area where well Q18 is located in the southeast to the direction of well Q19 in the northwest.

[0089] Step 2: Transportation path of carbonate rocks as the main source

[0090] (1) Analyze the plane contour map of the average particle size of carbonate minerals and determine the location of the high-value area of ​​the average particle size as the main source area of ​​carbonate rock. Figure 15 The high-value area of ​​average grain size of terrigenous carbonate rock minerals is the area where wells Q17 and Q18 are located, with average grain size values ​​of 53.2 and 52.3 respectively.

[0091] (2) Analyze the plane contour map of the average grain size of carbonate minerals and determine the path from the high-value area to the low-value area of ​​the average grain size as the transport path of carbonate rocks as the main source. The two provenance transport paths in the contour map are: 1. The path from the area where well Q17 is located in the southwest to the direction of well Q15 in the northeast; 2. The path from the area where well Q18 is located in the southeast to the direction of well Q19 in the northwest.

[0092] Step 3: Gravity flow transport path

[0093] (1) Analyze the clay mineral content contour map and determine the location of the high-value area of ​​clay mineral content as the source area of ​​gravity flow. Figure 6 The areas with high clay mineral content are the areas where wells Q15 and Q19 are located, with values ​​of 45.7 and 38.7, respectively.

[0094] (2) Comprehensively analyzing the clay mineral content and average grain size contour maps, the paths from high-value areas to low-value areas of content and average grain size were identified as gravity flow transport paths. The contour maps show two gravity flow transport paths: 1) the path from the Q15 well area to the northeast of the Q4 well; 2) the path from the Q19 well area to the northwest of the Q12 well.

[0095] Step 4: Comprehensive determination of the provenance and transportation path. By superimposing the source areas and transportation paths of the siliciclastic rocks as the main source, the source areas and transportation paths of the terrigenous carbonate minerals as the main source, and the source areas and transportation paths of the gravity flow, the overlapping source areas (i.e. Figure 16 The area where Q17 and Q18 are located) is determined as the provenance area of ​​the Long-1 layer in the study area, and the above-mentioned overlapping transportation paths (i.e., the two transportation paths: 1) the path from the area where well Q17 is located to the northeast direction of well Q15; 2) the path from the area where well Q18 is located to the northwest direction of well Q19) are determined as the provenance transportation paths of the Long-1 layer in the study area. (The provenance transportation paths are shown in Figure 16 ).

[0096] The embodiments described above are merely illustrative of embodiments of the present invention. Although the descriptions thereof are relatively specific and detailed, they should not be construed as limiting the scope of the present invention. The present invention may also be implemented in other specific ways or in other specific forms without departing from the gist or essential features of the present invention. Therefore, the embodiments described should be considered in all respects as illustrative and not restrictive. The scope of the present invention should be set forth in the appended claims, and any variations inconsistent with the intent and scope of the claims should also be included within the scope of the present invention.

Claims

1. A method for determining the provenance and transport path of black shale, characterized in that: The following steps are involved: Prepare the plane contour map of feldspar content, clay mineral content, average grain size of terrigenous quartz, average grain size of feldspar, average grain size of terrigenous carbonate minerals, and average grain size of clay minerals for the black shale of the target strata in the study area; Based on the plane contour map of the average grain size of terrigenous quartz, the plane contour map of the average grain size of feldspar, and the plane contour map of feldspar content, it is determined that the transportation path is mainly sourced from siliciclastic rock; based on the plane contour map of the average grain size of terrigenous carbonate minerals, it is determined that the transportation path is mainly sourced from the carbonate rock; based on the plane contour map of the clay mineral content and the plane contour map of the average grain size of clay minerals, it is determined that the transportation path is a gravity flow; By superimposing the transport path of the siliceous clastic rock as the main source, the transport path of the carbonate rock as the main source and the transport path of the gravity flow, the completely superimposed transport path is determined as the black shale source transport path.

2. The method for determining the black shale provenance transport path according to claim 1, characterized in that: The specific steps for determining the transport path of the siliciclastic rock as the main source are: Based on the plane contour map of average grain size of terrigenous quartz and the plane contour map of average grain size of feldspar, the location where the average grain size of both is high is preliminarily determined to be the source area dominated by siliciclastic rock; Combined with the location of the high-value area on the feldspar content plane contour map, the source area dominated by siliciclastic rock was finally determined; Based on the plane contour map of the average grain size of the terrigenous quartz and the feldspar, and the plane contour map of the feldspar content, the path of content and average grain size from high-value areas to low-value areas was determined to be the transportation path with the siliceous clastic rock as the main source.

3. The method for determining the black shale provenance transport path according to claim 1, characterized in that: The specific steps for determining the transport path of the carbonate rock as the main source are: Based on the plane contour map of the average grain size of terrestrial carbonate minerals, the area with high average grain size is determined as the source area mainly composed of carbonate rocks, and the path from the high average grain size area to the low average grain size area is determined as the transportation path with the carbonate rocks as the main source.

4. The method for determining the black shale provenance transport path according to claim 1, characterized in that: The specific steps for determining the transport path of the gravity flow are: Based on the plane contour map of clay mineral content, the area with high clay mineral content is determined as the source area of ​​the gravity flow; based on the plane contour map of clay mineral content and average particle size, the path from the high value area to the low value area of ​​content and average particle size is determined as the transportation path of the gravity flow.

5. The method for determining the black shale provenance transport path according to any one of claims 1 to 4, characterized in that: The specific steps for compiling the plane contour map of feldspar content and the plane contour map of clay mineral content are as follows: The arithmetic mean of feldspar content and clay mineral content of each well in the target stratum is calculated using the mineral composition logging data. Based on the location information of the research work area, a work area location map is compiled; the coordinates of the wells drilled into the target layer are imported into the work area location map to form a well location distribution map of the work area; Importing the arithmetic mean data of the feldspar content of each well into the well location distribution map of the work area, and using the Kriging interpolation method to compile the feldspar content plane contour map; The arithmetic mean value data of the clay mineral content of each well is imported into the well location distribution map of the work area, and the plane contour map of the clay mineral content is compiled by using the Kriging interpolation method.

6. The method for determining the black shale provenance transport path according to claim 5, characterized in that: The specific steps for compiling the plane contour map of the average grain size of terrigenous quartz are as follows: The average grain size data of the terrigenous quartz of each well is imported into the well location distribution map of the work area, and the plane contour map of the average grain size of the terrigenous quartz is compiled by using the Kriging interpolation method.

7. The method for determining the black shale provenance transport path according to claim 5, characterized in that: The specific steps of compiling the plane contour map of the average particle size of feldspar are as follows: importing the average particle size data of feldspar of each well on the drilling well location distribution map of the work area, and compiling the plane contour map of the average particle size of feldspar using the Kriging interpolation method.

8. The method for determining the black shale provenance transport path according to claim 5, characterized in that: The specific steps of compiling the plane contour map of the average particle size of the terrigenous carbonate minerals are as follows: the average particle size data of the terrigenous carbonate minerals of each well are imported into the well location distribution map of the work area, and the plane contour map of the average particle size of the terrigenous carbonate minerals is compiled by using the Kriging interpolation method.

9. The method for determining the black shale provenance transport path according to claim 5, characterized in that: The specific steps of compiling the clay mineral average particle size plane contour map are as follows: importing the average particle size data of clay minerals of each well on the drilling well location distribution map of the work area, and using the Kriging interpolation method to compile the clay mineral average particle size plane contour map.

10. The method for determining the black shale provenance transport path according to any one of claims 6 to 9, characterized in that: The average particle size of the terrigenous quartz, the feldspar, the clay mineral and the terrigenous carbonate mineral is obtained in the following manner: Determine the drilling depth of each well at which the clay mineral content is at its maximum value using mineral composition logging data, and collect core samples from each well within the study area at the drilling depth at which the clay mineral content is at its maximum value; The core samples were made into argon ion polished slices for scanning electron microscopy analysis, and the particle size of the terrigenous quartz, feldspar, clay minerals and terrigenous carbonate minerals under the scanning electron microscopy images of the argon ion polished slices was analyzed to calculate the average particle size data of the terrigenous quartz, feldspar, clay minerals and terrigenous carbonate minerals for each well.

Images