Carbon number determination and migration analysis method, device, equipment and medium for shale oil

Abstract

The present disclosure relates to a method and device for determining the carbon number and migration analysis of shale oil, an electronic device and a storage medium, and relates to the field of oilfield technology. The method for determining the carbon number of shale oil comprises: performing a full hydrocarbon gas chromatography analysis on the shale oil to obtain the contents corresponding to multiple different carbon numbers; performing a laser spectral analysis on the shale oil to obtain the fluorescence intensity of the shale oil in different bands; determining the first proportion of the first component and / or the second proportion of the second component corresponding to the set receiving band according to the fluorescence intensity in the different bands; determining the critical carbon number of the first component and the second component corresponding to the shale oil based on the first proportion / second proportion and the contents corresponding to the multiple different carbon numbers. The embodiments of the present disclosure can realize the determination of the carbon number and migration analysis of shale oil.

Description

Technical Field

[0001] The present disclosure relates to the field of oil field technology, and in particular to a method and device for determining the carbon number and migration analysis of shale oil, an electronic device and a storage medium, and specifically to a quantitative analysis method and device, an electronic device and a storage medium for the differential distribution of different crude oil components in micro-regions in unconventional reservoirs represented by mud shale during oil field exploration and development. Background Art

[0002] At present, due to the shortage of energy and the growing demand for energy in my country, major oil fields have carried out a lot of research work in the field of shale oil and gas exploration, and shale oil has also received more and more attention. Oil content evaluation, as an important attribute of shale oil reservoir evaluation, needs to be studied in depth. The traditional analysis methods for oil content evaluation mainly include rock pyrolysis, oil saturation, nuclear magnetic resonance analysis, etc., but the above methods can only give the total amount of shale oil, and it is urgent to develop quantitative analysis technology for the differential distribution of different crude oil components in micro-areas. The existing laser confocal technology abroad can realize full oil observation, but it cannot give the results of the spatial distribution differences of various crude oil components, which restricts the research on the law of crude oil enrichment. Summary of the invention

[0003] The present disclosure proposes a method and device for determining the carbon number and migration analysis of shale oil, an electronic device, and a storage medium technical solution.

[0004] According to one aspect of the present disclosure, a method for determining the carbon number of shale oil is provided, comprising:

[0005] The shale oil is subjected to total hydrocarbon gas chromatography analysis to obtain the contents corresponding to multiple different carbon numbers;

[0006] Performing laser spectroscopy analysis on the shale oil to obtain the fluorescence intensity of the shale oil in different wavebands;

[0007] Determine, according to the fluorescence intensity in the different wavelength bands, a first proportion of the first component and / or a second proportion of the second component corresponding to the set receiving wavelength band;

[0008] Based on the first proportion / second proportion and the contents corresponding to the multiple different carbon numbers, the critical carbon numbers of the first component and the second component corresponding to the shale oil are determined.

[0009] Preferably, the method for determining the critical carbon number of the light component and the heavy component corresponding to the shale oil based on the first proportion / the second proportion and the contents corresponding to the multiple different carbon numbers includes:

[0010] According to the contents corresponding to the multiple different carbon numbers, respectively determine the first content proportion corresponding to each carbon number;

[0011] Taking multiple set carbon numbers as the critical values, and determining multiple second content proportions corresponding to the multiple set carbon numbers and the carbon numbers below the multiple set carbon numbers according to the first content proportion at each carbon number; and / or

[0012] Based on the first proportion / second proportion and multiple second content proportions, the critical carbon number of the light component and the heavy component corresponding to the shale oil is determined; and / or, the method for determining the first proportion of the first component and / or the second proportion of the second component corresponding to the set receiving band according to the fluorescence intensity under the different bands includes: constructing the fluorescence intensity curve under the different bands; calculating the area ratio corresponding to the set receiving band and the fluorescence intensity curve under all different bands, and configuring the area ratio as the first proportion of the first component and / or the second proportion of the second component corresponding to the set receiving band; and / or,

[0013] The method for calculating the area ratio corresponding to the fluorescence intensity curves in the set receiving band and all different bands includes: integrating the fluorescence intensity curves in all different bands to obtain the total area of ​​the fluorescence intensity curves; integrating the fluorescence intensity curves in the set receiving band to obtain the area of ​​the fluorescence intensity curve region; calculating the ratio of the area of ​​the fluorescence intensity curve region to the total area of ​​the fluorescence intensity curve to obtain the corresponding area ratio; and / or,

[0014] The method for determining the critical carbon number of the light component and the heavy component corresponding to the shale oil based on the first proportion / the second proportion and the plurality of second content proportions includes:

[0015] Calculating the differences between the plurality of second content proportions and the first proportion / second proportion respectively to obtain a plurality of proportion differences;

[0016] The set carbon number corresponding to the minimum proportion difference among the multiple proportion differences is determined as the critical carbon number of the light component and the heavy component corresponding to the shale oil.

[0017] According to one aspect of the present disclosure, a method for analyzing the migration of shale oil is provided, comprising: the carbon number determination method as described above; and

[0018] The first component and the second component are respectively configured as the light component and the heavy component of the shale oil, and a first set receiving band and a second set receiving band corresponding to the critical carbon numbers of the light component and the heavy component are obtained;

[0019] Performing laser scanning on the rock of the shale oil, collecting reflected light signals and transmitted light signals corresponding to the set wavelength, and collecting corresponding first fluorescence signals and second fluorescence signals based on the first set receiving band and the second set receiving band, respectively;

[0020] Determine the distribution of the light component and / or the heavy component based on the reflected light signal, the transmitted light signal, the first fluorescence signal, and the second fluorescence signal;

[0021] Based on the distribution of the light components and / or heavy components, the migration rate of the light components and / or heavy components of the shale oil is determined.

[0022] Preferably, the method for determining the distribution of the light component and / or the heavy component based on the reflected light signal, the transmitted light signal, the first fluorescent signal and the second fluorescent signal, respectively, comprises:

[0023] Determine a first distribution characteristic map corresponding to matrix pores of the rock and a second distribution characteristic map corresponding to cracks based on the reflected light signal and the transmitted light signal respectively;

[0024] The first fluorescent signal and the second fluorescent signal are respectively arranged in the first distribution characteristic diagram and the second distribution characteristic diagram to respectively determine the distribution of the light component and the heavy component; and / or,

[0025] The method for determining the migration rate of the light component and / or the heavy component of the shale oil based on the distribution of the light component and / or the heavy component comprises:

[0026] Based on the distribution of the light component and / or the heavy component, the first volume V of the light component and / or the heavy component in the seam is determined respectively. 轻1 and the second volume V 重1 , and respectively determine the third volume V of the light component and / or heavy component in the matrix pores 轻2 and the fourth volume V 重2 ;

[0027] Based on the first volume V 轻1 and the third volume V 轻2 , determining the migration rate of the light component of the shale oil;

[0028] Based on the second volume V 重1 and the fourth volume V 重2 , determining the migration rate of the heavy components of the shale oil; and / or,

[0029] Before the laser scanning of the shale oil rock, one or more of the energy intensity of the laser and / or the laser wavelength and / or the pinhole value and / or the scanning mode and / or the line accumulation and / or the surface average and / or the scanning speed and / or the top and bottom of the scanning and / or the scanning thickness and / or the number of scanning layers and / or the interlayer spacing are configured; and / or,

[0030] Before laser scanning the shale oil rock, the shale oil rock is prepared, including:

[0031] Obtaining an oil-bearing shale rock sample, and cutting the shale rock sample into rock blocks according to a set first size or into rock columns according to a set second size;

[0032] Gluing the rock block or rock column onto a glass slide, and cutting the rock block or rock column into straight rock slices according to a first set thickness;

[0033] The straight rock slice is ground to a second set thickness to obtain the shale oil rock for laser scanning.

[0034] Preferably, before configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic diagram and the second distribution characteristic diagram, respectively, and determining the distribution of the light component and the heavy component, respectively, the first distribution characteristic diagram and the second distribution characteristic diagram are corrected to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram, respectively, and the first fluorescence signal and the second fluorescence signal are respectively configured in the first corrected distribution characteristic diagram and the second corrected distribution characteristic diagram, respectively, to determine the distribution of the light component and the heavy component; wherein the method of correcting the first distribution characteristic diagram and the second distribution characteristic diagram to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram, respectively, includes:

[0035] Superimposing the first distribution characteristic graph and the second distribution characteristic graph to obtain a third distribution characteristic graph;

[0036] Correcting the first distribution characteristic map and the second distribution characteristic map based on the overlapping area of ​​the third distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map; and / or,

[0037] based on the first volume V 轻1 and the third volume V 轻2 The method for determining the migration rate of the light component of the shale oil comprises:

[0038] The first volume V 轻1 Divide by the first volume V 轻1 and the third volume V 轻2 The sum of the above is used to obtain the migration rate of the light component of the shale oil; and / or,

[0039] The second volume V 重1 and the fourth volume V 重2 The method for determining the migration rate of the heavy component of the shale oil comprises:

[0040] The first volume V 重1Divide by the first volume V 重1 and the third volume V 重2 The migration rate of the heavy components of the shale oil is obtained by summing the above-mentioned values; and / or,

[0041] The method of correcting the first distribution characteristic map and the second distribution characteristic map based on the overlapping area of ​​the third distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map includes:

[0042] Performing edge detection on the overlapping area to obtain edge features;

[0043] If the edge feature is a line, the edge feature that is a line is deleted from the first distribution feature map to obtain a first corrected distribution feature map;

[0044] If the edge feature is a point, the edge feature that is a point is deleted from the second distribution feature map to obtain a second corrected distribution feature map.

[0045] According to one aspect of the present disclosure, a shale oil analysis method is provided, comprising: the carbon number determination method as described above; and / or, the migration analysis method as described above; and

[0046] The first component and the second component are respectively configured as the light component and the heavy component of the shale oil, and a first set receiving band and a second set receiving band corresponding to the critical carbon numbers of the light component and the heavy component are obtained;

[0047] Performing laser scanning on the rock of the shale oil, collecting reflected light signals and transmitted light signals corresponding to the set wavelength, and collecting corresponding first fluorescence signals and second fluorescence signals based on the first set receiving band and the second set receiving band, respectively;

[0048] Determine the distribution of the light component and / or the heavy component based on the reflected light signal, the transmitted light signal, the first fluorescence signal, and the second fluorescence signal;

[0049] Based on the distribution of the light components and / or heavy components, the light-to-weight ratio of crude oil in the seams of the shale oil and / or the light-to-weight ratio of crude oil in the matrix pores and / or the light-to-weight ratio of the shale oil and / or the volume percentage content in the shale oil seams and / or the volume percentage content in the shale oil matrix pores and / or the oil saturation are determined.

[0050] Preferably, the method for determining the distribution of the light component and / or the heavy component based on the reflected light signal, the transmitted light signal, the first fluorescent signal and the second fluorescent signal, respectively, comprises:

[0051] Determine a first distribution characteristic map corresponding to matrix pores of the rock and a second distribution characteristic map corresponding to cracks based on the reflected light signal and the transmitted light signal respectively;

[0052] The first fluorescent signal and the second fluorescent signal are respectively arranged in the first distribution characteristic diagram and the second distribution characteristic diagram to respectively determine the distribution of the light component and the heavy component; and / or,

[0053] The method for determining the crude oil light-to-weight ratio in the shale oil seams and / or the crude oil light-to-weight ratio in the matrix pores and / or the shale oil light-to-weight ratio and / or the volume percentage content in the shale oil seams and / or the volume percentage content in the shale oil matrix pores and / or the oil saturation based on the distribution of the light component and / or the heavy component comprises:

[0054] Based on the distribution of the light component and / or the heavy component, the first volume V of the light component and / or the heavy component in the seam is determined respectively. 轻1 and the second volume V 重1 , and respectively determine the third volume V of the light component and / or heavy component in the matrix pores 轻2 and the fourth volume V 重2 ;

[0055] Based on the first volume V 轻1 and the second volume V 重1 , determining the crude oil light-to-heavy ratio in the seams of the shale oil;

[0056] Based on the third volume V 轻2 and the fourth volume V 重2 , determining the crude oil weight ratio in the matrix pores of the shale oil;

[0057] Based on the first volume V 轻1 , the second volume V 重1 , the third volume V 轻2 and the fourth volume V 重2 , determining the light-to-heavy ratio of the shale oil;

[0058] Based on the first volume V 轻1 , the second volume V 重1 and the rock volume of the shale oil, to determine the volume percentage content in the shale oil fracture;

[0059] Based on the third volume V 轻2 、The fourth volume V 重2 and the rock volume of the shale oil, to determine the volume percentage in the pores of the shale oil matrix;

[0060] Determining the oil saturation of the shale oil based on the volume percentage content in the shale oil seams, the volume percentage content in the shale oil seams and the porosity of the rock of the shale oil; and / or,

[0061] Before the laser scanning of the shale oil rock, one or more of the energy intensity of the laser and / or the laser wavelength and / or the pinhole value and / or the scanning mode and / or the line accumulation and / or the surface average and / or the scanning speed and / or the top and bottom of the scanning and / or the scanning thickness and / or the number of scanning layers and / or the interlayer spacing are configured; and / or,

[0062] Before laser scanning the shale oil rock, the shale oil rock is prepared, including:

[0063] Obtaining an oil-bearing shale rock sample, and cutting the shale rock sample into rock blocks according to a set first size or into rock columns according to a set second size;

[0064] Gluing the rock block or rock column onto a glass slide, and cutting the rock block or rock column into straight rock slices according to a first set thickness;

[0065] and grinding the straight rock slice to a second set thickness to obtain the shale oil rock for laser scanning; and / or,

[0066] Before configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic diagram and the second distribution characteristic diagram, respectively, and determining the distribution of the light component and the heavy component, respectively, the first distribution characteristic diagram and the second distribution characteristic diagram are corrected to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram, respectively, and the first fluorescence signal and the second fluorescence signal are respectively configured in the first corrected distribution characteristic diagram and the second corrected distribution characteristic diagram, respectively, to determine the distribution of the light component and the heavy component; wherein the method of correcting the first distribution characteristic diagram and the second distribution characteristic diagram to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram, respectively, comprises:

[0067] Superimposing the first distribution characteristic graph and the second distribution characteristic graph to obtain a third distribution characteristic graph;

[0068] Correcting the first distribution characteristic map and the second distribution characteristic map based on the overlapping area of ​​the third distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map; and / or,

[0069] based on the first volume V 轻1 and the second volume V 重1 The method for determining the crude oil weight ratio in the seam of the shale oil comprises: the first volume V 轻1 Divide by the second volume V重1 , obtaining the crude oil weight ratio in the seam of the shale oil; and / or,

[0070] The third volume V 轻2 and the fourth volume V 重2 , a method for determining the crude oil weight ratio in the matrix pores of the shale oil, comprising: the third volume V 轻2 Divide by the fourth volume V 重2 , obtaining the crude oil weight ratio in the matrix pores of the shale oil; and / or,

[0071] based on the first volume V 轻1 , the second volume V 重1 , the third volume V 轻2 and the fourth volume V 重2 , a method for determining the weight ratio of the shale oil comprises: 轻1 and the third volume V 轻2 Sum, get the first total; for the second volume V 重1 and the fourth volume V 重2 Sum them to obtain a second sum; divide the first sum by the second sum to obtain the light-to-heavy ratio of the shale oil; and / or,

[0072] based on the first volume V 轻1 , the second volume V 重1 and the rock volume of the shale oil, and a method for determining the volume percentage content in the shale oil fracture, comprising: 轻1 and the second volume V 重1 After summing, divide by the rock volume of the shale oil to obtain the volume percentage content in the shale oil fracture; and / or,

[0073] The third volume V 轻2 、The fourth volume V 重2 and the rock volume of the shale oil, and a method for determining the volume percentage content in the pores of the shale oil matrix, comprising: 轻2 and the second volume V 重2 After summing, divide by the rock volume of the shale oil to obtain the volume percentage content in the pores of the shale oil matrix; and / or,

[0074] The method for determining the oil saturation of the shale oil based on the volume percentage content in the shale oil fracture, the volume percentage content in the shale oil fracture and the porosity of the rock of the shale oil comprises: summing the volume percentage content in the shale oil fracture and the volume percentage content in the shale oil fracture, and dividing the sum by the porosity of the rock of the shale oil to obtain the oil saturation of the shale oil; and / or,

[0075] The method of correcting the first distribution characteristic map and the second distribution characteristic map based on the overlapping area of ​​the third distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map includes:

[0076] Performing edge detection on the overlapping area to obtain edge features;

[0077] If the edge feature is a line, the edge feature that is a line is deleted from the first distribution feature map to obtain a first corrected distribution feature map;

[0078] If the edge feature is a point, the edge feature that is a point is deleted from the second distribution feature map to obtain a second corrected distribution feature map.

[0079] According to one aspect of the present disclosure, there is provided a device, comprising: a shale oil carbon number determination device and / or a shale oil migration analysis device and / or a shale oil analysis device;

[0080] The device for determining the carbon number of shale oil comprises:

[0081] A full hydrocarbon gas chromatography analysis unit is used to perform full hydrocarbon gas chromatography analysis on shale oil to obtain contents corresponding to multiple different carbon numbers;

[0082] A laser spectrum analysis unit, used for performing laser spectrum analysis on the shale oil to obtain the fluorescence intensity of the shale oil in different wavebands;

[0083] A first determining unit, configured to determine a first proportion of the first component and / or a second proportion of the second component corresponding to a set receiving band according to the fluorescence intensities in the different bands;

[0084] A second determination unit is configured to determine the critical carbon numbers of the first component and the second component corresponding to the shale oil based on the first proportion / the second proportion and the contents corresponding to the multiple different carbon numbers; and / or,

[0085] The shale oil migration analysis device includes: a shale oil carbon number determination device; and,

[0086] A first acquisition unit, configured to configure the first component and the second component as the light component and the heavy component of the shale oil, respectively, and to acquire a first set receiving band and a second set receiving band corresponding to the critical carbon numbers of the light component and the heavy component;

[0087] a first laser scanning unit, configured to perform laser scanning on the rock of the shale oil, collect reflected light signals and transmitted light signals corresponding to a set wavelength, and collect corresponding first fluorescence signals and second fluorescence signals based on the first set receiving band and the second set receiving band, respectively;

[0088] a first distribution determination unit, configured to determine the distribution of the light component and / or the heavy component based on the reflected light signal, the transmitted light signal, the first fluorescent signal, and the second fluorescent signal;

[0089] A migration rate analysis unit, for determining the migration rate of the light component and / or the heavy component of the shale oil based on the distribution of the light component and / or the heavy component; and / or,

[0090] The shale oil analysis device includes: a shale oil carbon number determination device and / or a shale oil migration analysis device; and,

[0091] A second acquisition unit is used to configure the first component and the second component as the light component and the heavy component of the shale oil respectively, and obtain a first set receiving band and a second set receiving band corresponding to the critical carbon number of the light component and the heavy component;

[0092] a second laser scanning unit, configured to perform laser scanning on the rock of the shale oil, collect reflected light signals and transmitted light signals corresponding to a set wavelength, and collect corresponding first fluorescence signals and second fluorescence signals based on the first set receiving band and the second set receiving band, respectively;

[0093] a second distribution determination unit, configured to determine the distribution of the light component and / or the heavy component based on the reflected light signal, the transmitted light signal, the first fluorescent signal, and the second fluorescent signal;

[0094] An analysis unit, for determining, based on the distribution of the light component and / or the heavy component, the crude oil light-to-weight ratio in the shale oil seams and / or the crude oil light-to-weight ratio in the matrix pores and / or the shale oil light-to-weight ratio and / or the volume percentage content in the shale oil seams and / or the volume percentage content in the shale oil matrix pores and / or the oil saturation. According to one aspect of the present disclosure, an electronic device is provided, comprising:

[0095] processor;

[0096] a memory for storing processor-executable instructions;

[0097] Wherein, the processor is configured to: execute the above-mentioned carbon number determination method; and / or, the above-mentioned migration analysis method; and / or, the above-mentioned analysis method.

[0098] According to one aspect of the present disclosure, a computer-readable storage medium is provided, on which computer program instructions are stored, and when the computer program instructions are executed by a processor, the carbon number determination method described above; and / or, the migration analysis method described above; and / or, the analysis method described above are implemented.

[0099] In the embodiments disclosed herein, the proposed method and device for determining the carbon number and migration analysis of shale oil, electronic equipment and storage medium can realize the determination of the critical carbon number of the first component and the second component corresponding to the shale oil, so as to solve the problem that the current results of the spatial distribution differences of the various components of crude oil cannot be given, which restricts the research on the enrichment law of crude oil.

[0100] It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the disclosure.

[0101] Further features and aspects of the present disclosure will become apparent from the following detailed description of exemplary embodiments with reference to the attached drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0102] The drawings herein are incorporated into the specification and constitute a part of the specification. These drawings illustrate embodiments consistent with the present disclosure and are used to illustrate the technical solutions of the present disclosure together with the specification.

[0103] Figure 1 A flow chart showing a method for determining the carbon number of shale oil according to an embodiment of the present disclosure;

[0104] Figure 2 The full hydrocarbon gas chromatogram and laser spectrum analysis diagram corresponding to crude oil samples of different densities according to the embodiment of the present disclosure are shown;

[0105] Figure 3 A flow chart showing a method for analyzing the migration of shale oil according to an embodiment of the present disclosure;

[0106] Figure 4 A distribution characteristic diagram showing light components and / or heavy components of crude oil in corresponding matrix pores and / or fractures according to an embodiment of the present disclosure;

[0107] Figure 5 is a block diagram of an electronic device 800 according to an exemplary embodiment;

[0108] Figure 6 is a block diagram of an electronic device 1900 according to an exemplary embodiment. DETAILED DESCRIPTION

[0109] Various exemplary embodiments, features and aspects of the present disclosure will be described in detail below with reference to the accompanying drawings. The same reference numerals in the accompanying drawings represent elements with the same or similar functions. Although various aspects of the embodiments are shown in the accompanying drawings, the drawings are not necessarily drawn to scale unless otherwise specified.

[0110] The word “exemplary” is used exclusively herein to mean “serving as an example, example, or illustration.” Any embodiment described herein as “exemplary” is not necessarily to be construed as preferred or advantageous over other embodiments.

[0111] The term "and / or" herein is only a description of the association relationship of the associated objects, indicating that there may be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the term "at least one" herein represents any combination of at least two of any one or more of a plurality of. For example, including at least one of A, B, and C can represent including any one or more elements selected from the set consisting of A, B, and C.

[0112] In addition, in order to better illustrate the present disclosure, numerous specific details are given in the following specific embodiments. It should be understood by those skilled in the art that the present disclosure can also be implemented without certain specific details. In some examples, methods, means, components and circuits well known to those skilled in the art are not described in detail in order to highlight the subject matter of the present disclosure.

[0113] It can be understood that the above-mentioned embodiments of the shale oil carbon number determination and migration analysis methods mentioned in the present disclosure can be combined with each other to form combined embodiments without violating the principle logic. Due to space limitations, the present disclosure will not go into details.

[0114] In addition, the present disclosure also provides a shale oil carbon number determination, migration analysis device, electronic device, computer-readable storage medium, and program, all of which can be used to implement any shale oil carbon number determination, migration analysis method provided by the present disclosure. The corresponding technical solutions and descriptions are referred to the corresponding records in the method part and will not be repeated here.

[0115] Figure 1 A flow chart showing a method for determining the carbon number of shale oil according to an embodiment of the present disclosure is shown as follows: Figure 1 As shown, the method for determining the carbon number of shale oil includes: step 101: performing full hydrocarbon gas chromatography analysis on the shale oil to obtain the contents corresponding to multiple different carbon numbers; step 102: performing laser spectral analysis on the shale oil to obtain the fluorescence intensity of the shale oil in different bands; step 103: determining the first proportion of the first component and / or the second proportion of the second component corresponding to the set receiving band according to the fluorescence intensity in the different bands; step 104: determining the critical carbon number of the first component and the second component corresponding to the shale oil based on the first proportion / second proportion and the contents corresponding to the multiple different carbon numbers. The critical carbon number of the first component and the second component corresponding to the shale oil can be determined to solve the problem that the current results of the spatial distribution differences of the components of crude oil cannot be given, which restricts the research on the enrichment law of crude oil.

[0116] Step 101: Perform a total hydrocarbon gas chromatography analysis on the shale oil to obtain contents corresponding to multiple different carbon numbers.

[0117] In the embodiments of the present disclosure and other possible embodiments, the crude oil full hydrocarbon gas chromatograph can determine the range and content of normal alkanes from the lowest carbon number to the highest carbon number contained in the sample. For example, the crude oil full hydrocarbon gas chromatograph can be used to perform full hydrocarbon gas chromatographic analysis on shale oil to obtain the contents corresponding to multiple different carbon numbers.

[0118] In the embodiments of the present disclosure and other possible embodiments, selected shale oil crude oil is taken, and a set volume (for example, 4 ml to 6 ml) of shale oil crude oil sample is taken with a glass container and sealed; if the shale oil crude oil contains water, sampling is performed according to the oil-water ratio, and the sampling volume should ensure that the sample after dehydration is a set volume (for example, 4, ml to 6 ml), and sealed. Before chromatographic analysis, the set temperature (for example, 40 to 50 ° C) is kept constant for a set time (for example, 1 to 2 hours), and after the oil and water are separated, samples are taken from the oil layer for analysis, and the shale oil is subjected to full hydrocarbon gas chromatography analysis using a crude oil full hydrocarbon gas chromatograph to obtain the contents corresponding to multiple different carbon numbers.

[0119] Step 102: Perform laser spectral analysis on the shale oil to obtain the fluorescence intensity of the shale oil in different bands (different wavelengths).

[0120] In the embodiments of the present disclosure and other possible embodiments, a laser spectrometer in the field of chemical engineering may be used to perform laser spectral analysis on the shale oil to obtain the fluorescence intensity of the shale oil in different bands.

[0121] In the embodiments of the present disclosure and other possible embodiments, a laser spectrum analyzer with a set wavelength can be used to perform laser spectrum analysis on the shale oil to obtain the fluorescence intensity of the shale oil in different bands. The set wavelength can be configured to 488 nm. Similarly, those skilled in the art can configure the set wavelength according to actual needs.

[0122] Step 103: Determine a first proportion of the first component and / or a second proportion of the second component corresponding to a set receiving band according to the fluorescence intensities in the different bands.

[0123] In the embodiments of the present disclosure and other possible embodiments, the first component and the second component can be configured as the light component and the heavy component of the shale oil, respectively. In addition, in the embodiments of the present disclosure and other possible embodiments, the set receiving band can be configured as 500-550nm. Similarly, those skilled in the art can configure the set receiving band according to actual needs. Among them, the set receiving band is configured as the receiving band corresponding to the first component or the second component.

[0124] For example, the set receiving band of the first component configuration can be 500-550nm, and a laser spectrometer in the field of chemical engineering is used to perform laser spectral analysis on the shale oil to obtain the fluorescence intensity of the shale oil in different bands, and then the first proportion of the first component corresponding to the set receiving band is determined based on the fluorescence intensity in different bands.

[0125] In an embodiment of the present disclosure, the method for determining the first proportion of the first component and / or the second proportion of the second component corresponding to the set receiving band based on the fluorescence intensity under the different bands includes: constructing fluorescence intensity curves under the different bands; calculating the area ratio corresponding to the set receiving band and the fluorescence intensity curves under all different bands, and configuring the area ratio as the first proportion of the first component and / or the second proportion of the second component corresponding to the set receiving band.

[0126] In the embodiments of the present disclosure and other possible embodiments, a laser spectrometer in the field of chemical engineering is used to perform laser spectral analysis on the shale oil to obtain the fluorescence intensity of the shale oil in different bands (different wavelengths). That is, the different bands (different wavelengths) can be used as the horizontal coordinate, and the fluorescence intensity of the shale oil corresponding to the different bands (different wavelengths) can be used as the vertical coordinate to construct the fluorescence intensity curves in the different bands (different wavelengths).

[0127] In an embodiment of the present disclosure, based on the fluorescence intensity curve, the method for calculating the area ratio corresponding to the set receiving band and the fluorescence intensity curves in all different bands includes: integrating the fluorescence intensity curves in all different bands to obtain the total area of ​​the fluorescence intensity curve; integrating the fluorescence intensity curve in the set receiving band to obtain the regional area of ​​the fluorescence intensity curve; calculating the ratio of the regional area of ​​the fluorescence intensity curve to the total area of ​​the fluorescence intensity curve to obtain the corresponding area ratio.

[0128] For example, the set receiving band of the first component configuration can be 500-550nm, and the fluorescence intensity curves in all different bands of 500-800nm ​​are integrated to obtain the total area of ​​the fluorescence intensity curve (for example, 100); the fluorescence intensity curve in the set receiving band of 500-550nm is integrated to obtain the area of ​​the fluorescence intensity curve region (for example, 63.55); the ratio of the area of ​​the fluorescence intensity curve region to the total area of ​​the fluorescence intensity curve is calculated to obtain the corresponding area ratio of 63.55%.

[0129] Step 104: Based on the first proportion / second proportion and the contents corresponding to the multiple different carbon numbers, determine the critical carbon numbers of the first component and the second component corresponding to the shale oil, thereby providing a corresponding basis for the division of light components and heavy components of shale oil by laser confocal microscopy.

[0130] In an embodiment of the present disclosure, the method for determining the critical carbon number of the light components and heavy components corresponding to the shale oil based on the first proportion / second proportion and the contents corresponding to the multiple different carbon numbers includes: determining the first content proportion corresponding to each carbon number according to the contents corresponding to the multiple different carbon numbers; taking multiple set carbon numbers as the critical number, and determining the multiple second content proportions corresponding to the multiple set carbon numbers and the carbon numbers below them according to the first content proportion under each carbon number; determining the critical carbon number of the light components and heavy components corresponding to the shale oil based on the first proportion / second proportion and the multiple second content proportions.

[0131] For example, in the embodiments of the present disclosure and other possible embodiments, multiple set carbon numbers can be respectively configured as carbon 0 (C1) to carbon 60 (C60), with multiple set carbon numbers as critical (i.e., carbon 1 (C1) to carbon 60 (C60)), according to the first content proportion under each carbon number, for example: the first content proportion under C0 is 0, the first content proportion under C1 is 0, the first content proportion under C2 is 0.02, the first content proportion under C3 is 0.30, the first content proportion under C4 is 1.12, the first content proportion under C5 is 2.29, the first content proportion under C6 is 3.68,…

[0132] In the embodiments of the present disclosure and other possible embodiments, a method for determining multiple second content proportions corresponding to the multiple set carbon numbers and carbon numbers below, based on the first content proportion at each carbon number, includes: summing the first content proportions at each carbon number of the multiple set carbon numbers and carbon numbers below, to obtain multiple second content proportions corresponding to the multiple set carbon numbers and carbon numbers below. For example, multiple set carbon numbers can be respectively configured as carbon 0 (C1) to carbon 60 (C60), and the multiple second content ratios corresponding to the first set carbon number C1 and the carbon numbers below the multiple set carbon numbers are configured as 0; the multiple second content ratios corresponding to the second set carbon number C2 and the carbon number (C1) below the multiple set carbon numbers are configured as 0.02 (0.02+0); the multiple second content ratios corresponding to the third set carbon number C3 and the carbon numbers (C1 and C2) below the multiple set carbon numbers are configured as 0.32 (0.30+0.02+0); the multiple second content ratios corresponding to the fourth set carbon number C4 and the carbon numbers (C1, C2 and C3) below the multiple set carbon numbers are configured as 1.44 (1.12+0.30+0.02+0);….

[0133] In an embodiment of the present disclosure, the method for determining the critical carbon number of the light components and heavy components corresponding to the shale oil based on the first proportion / second proportion and multiple second content proportions includes: respectively calculating the difference between the multiple second content proportions and the first proportion / second proportion to obtain multiple proportion differences; determining the set carbon number corresponding to the minimum proportion difference among the multiple proportion differences as the critical carbon number of the light components and heavy components corresponding to the shale oil.

[0134] In the embodiments of the present disclosure and other possible embodiments, after respectively calculating the difference between the multiple second content proportions and the first proportion / second proportion, the absolute value of the difference is taken to obtain multiple proportion differences. For example, in the embodiments of the present disclosure and other possible embodiments, the first proportion / second proportion is configured as 1, and the second content proportions corresponding to the carbon numbers below C1, C2, C3, and C4 are respectively configured as 0, 0.02, 0.32, and 1.44, and the difference between the second content proportion and the first proportion / second proportion is respectively calculated, and the absolute value of the difference is taken to obtain the proportion differences of 0, 0.98, 0.68, and 0.44; among them, the proportion difference between the carbon number below C4 and the first proportion / second proportion is the minimum proportion difference, so C4 is determined as the critical carbon number of the light component and the heavy component corresponding to the shale oil.

[0135] In the embodiments of the present disclosure and other possible embodiments, Figure 2 The full hydrocarbon gas chromatogram and laser spectrum analysis diagram corresponding to crude oil samples of different densities according to the embodiment of the present disclosure are shown as follows: Figure 2 As shown, Figure 2 (a)-(e) have a density of 0.7787 g / cm 3 、0.9389g / cm 3 、0.7438g / cm 3 、0.8903g / cm 3 , 0.8588g / cm 3 The total hydrocarbon gas chromatogram and laser spectrum analysis diagram corresponding to the crude oil (shale oil) sample. Figure 2 The left sides of (a)-(e) have a density of 0.7787 g / cm 3 、0.9389g / cm 3 、0.7438g / cm 3 、0.8903g / cm 3 , 0.8588g / cm 3 The total hydrocarbon gas chromatogram corresponding to the crude oil (shale oil) sample (the ordinate represents the peak intensity pA, and the abscissa represents the C number.); while, Figure 2The right sides of (a)-(e) have densities of 0.7787 g / cm 3 、0.9389g / cm 3 、0.7438g / cm 3 、0.8903g / cm 3 , 0.8588g / cm 3 Laser spectrum analysis diagram corresponding to the crude oil (shale oil) sample.

[0136] like Figure 2 (a) shows that the density is 0.8588 g / cm 3 The results of the gas chromatography analysis of all hydrocarbons in crude oil samples were obtained. Spectral scanning was performed. Using laser excitation with a set wavelength (e.g., 488nm) and scanning with a set step length (e.g., 5nm), the results showed that the integrated area of ​​the set receiving band (e.g., 500nm to 550nm) accounted for 63.55%, and the integrated area of ​​the part above the boundary of the set receiving band of 550nm accounted for 36.45%, which is equivalent to the gas chromatography data of all hydrocarbons. As shown in Table 1 and Figure 2 (a) shows that the density is 0.8588 g / cm 3 The results of gas chromatography analysis of crude oil total hydrocarbons showed that: 15 The components with C numbers below account for 61.59% of the total hydrocarbon components. 15 And above C arrays account for 38.41%; the C number corresponding to the boundary value of the receiving band of 550nm is set to C15.

[0137] like Figure 2 (b) shows that the density is 0.9389 g / cm 3 The results of the gas chromatography analysis of all hydrocarbons in crude oil samples were obtained. Spectral scanning was performed. Using laser excitation with a set wavelength (e.g., 488nm) and scanning with a set step length (e.g., 5nm), the results showed that the integrated area of ​​the set receiving band (e.g., 500nm to 550nm) accounted for 28.15%, and the integrated area of ​​the part above the boundary of the set receiving band of 550nm accounted for 71.85%, which is equivalent to the gas chromatography data of all hydrocarbons. As shown in Table 2 and Figure 2 (b) shows that the density is 0.9389 g / cm 3 The results of gas chromatography analysis of crude oil total hydrocarbons showed that: 15 The components with C numbers below account for 30.09% of the total hydrocarbon components. 15 And above C arrays account for 69.91%; the C number corresponding to the boundary value of the receiving band of 550nm is set to C15.

[0138] like Figure 2 (c) shows that the density is 0.7438 g / cm 3The results of the gas chromatography analysis of all hydrocarbons in crude oil samples were obtained. Spectral scanning was performed. Laser excitation with a set wavelength (e.g., 488nm) was used and scanning was performed with a set step length (e.g., 5nm). The results showed that the integrated area of ​​the set receiving band (e.g., 500nm to 550nm) accounted for 100%, and the integrated area of ​​the part above the boundary of the set receiving band of 550nm accounted for 0%, which was equivalent to the gas chromatography data of all hydrocarbons. As shown in Table 3 and Figure 2 (c) shows that the density is 0.7438 g / cm 3 The results of gas chromatography analysis of crude oil total hydrocarbons showed that: 15 The components with C numbers below account for 94.89% of the total hydrocarbon components. 15 And above C arrays account for 5.11%; the C number corresponding to the boundary value of the receiving band of 550nm is set to C15.

[0139] like Figure 2 As shown in (d), the density is 0.8903 g / cm 3 The results of the gas chromatography analysis of all hydrocarbons in crude oil samples were obtained. Spectral scanning was performed. Using laser excitation with a set wavelength (e.g., 488nm) and scanning with a set step length (e.g., 5nm), the results showed that the integrated area of ​​the set receiving band (e.g., 500nm to 550nm) accounted for 36.80%, and the integrated area of ​​the part above the boundary of the set receiving band of 550nm accounted for 63.20%, which is equivalent to the gas chromatography data of all hydrocarbons. As shown in Table 4 and Figure 2 As shown in (d), the density is 0.8903 g / cm 3 The results of gas chromatography analysis of crude oil total hydrocarbons showed that: 15 The components with C numbers below account for 24.74% of the total hydrocarbon components. 15 And above C arrays account for 75.26%; the C number corresponding to the boundary value of the receiving band of 550nm is set to C15.

[0140] like Figure 2 As shown in (e), the density is 0.8588 g / cm 3 The results of the gas chromatography analysis of all hydrocarbons in crude oil samples were obtained. Spectral scanning was performed. Using laser excitation with a set wavelength (e.g., 488nm) and scanning with a set step length (e.g., 5nm), the results showed that the integrated area of ​​the set receiving band (e.g., 500nm to 550nm) accounted for 42.77%, and the integrated area of ​​the part above the boundary of the set receiving band of 550nm accounted for 57.23%, which is equivalent to the gas chromatography data of all hydrocarbons. As shown in Table 5 and Figure 2 As shown in (e), the density is 0.8588 g / cm 3 The results of gas chromatography analysis of crude oil total hydrocarbons showed that: 15 The components with C numbers below account for 27.43% of the total hydrocarbon components.15 And above C arrays account for 72.57%; the C number corresponding to the boundary value of the receiving band of 550nm is set to C15.

[0141] Table 1 Density is 0.8588 g / cm 3 The results of gas chromatography analysis of total hydrocarbons in crude oil samples

[0142]

[0143] Table 2 Density is 0.9389 g / cm 3 The results of gas chromatography analysis of total hydrocarbons in crude oil samples

[0144]

[0145] Table 3 density is 0.7438g / cm 3 The results of gas chromatography analysis of total hydrocarbons in crude oil samples

[0146]

[0147] Table 4 density is 0.8903g / cm 3 The results of gas chromatography analysis of total hydrocarbons in crude oil samples

[0148]

[0149] Table 5 density is 0.8588g / cm 3 The results of total hydrocarbon gas chromatography analysis of crude oil samples.

[0150]

[0151] Figure 3 A flow chart of a shale oil migration analysis method according to an embodiment of the present disclosure is shown. Figure 3 As shown, the migration analysis method of shale oil includes: the carbon number determination method as mentioned above; and step S1101: respectively configuring the first component and the second component as the light component and the heavy component of the shale oil, and obtaining the first set receiving band and the second set receiving band corresponding to the critical carbon number of the light component and the heavy component; step S1102: performing laser scanning on the rock of the shale oil, collecting the reflected light signal and the transmitted light signal corresponding to the set wavelength, and collecting the corresponding first fluorescence signal and the second fluorescence signal based on the first set receiving band and the second set receiving band respectively; step S1103: respectively determining the distribution of the light component and / or the heavy component based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal; step S1104: determining the migration rate of the light component and / or the heavy component of the shale oil based on the distribution of the light component and / or the heavy component.

[0152] Step S1101: respectively configuring the first component and the second component as the light component and the heavy component of the shale oil, and obtaining the first set receiving band and the second set receiving band corresponding to the critical carbon number of the light component and the heavy component.

[0153] For example, in the embodiments of the present disclosure and other possible embodiments, the critical carbon number of the light component and the heavy component can be configured as C15; wherein the first setting receiving wave can be configured as 500nm~550nm, and the second setting receiving wave band can be configured as 550nm~800nm.

[0154] Step S1102: perform laser scanning on the shale oil rock, collect reflected light signals and transmitted light signals corresponding to the set wavelength, and collect corresponding first fluorescence signals and second fluorescence signals based on the first set receiving band and the second set receiving band, respectively.

[0155] In an embodiment of the present disclosure, before laser scanning is performed on the shale oil rock, the shale oil rock is prepared, including: obtaining an oil-containing shale rock sample, cutting the shale rock sample into rock blocks according to a set first size or cutting the rock columns according to a set second size; sticking the rock blocks or rock columns to a glass slide, and cutting the rock blocks or rock columns into straight rock slices according to a first set thickness; and grinding the straight rock slices to a second set thickness to obtain the shale oil rock for laser scanning.

[0156] In the embodiments of the present disclosure and other possible embodiments, an oil-bearing shale rock sample is selected, and rock blocks or rock columns are cut. The size of the rock block should be 25mm×25mm×5mm (set the first size), and the diameter of the rock column should be 25mm×5mm (set the second size). The sample cannot be soaked in organic solvents before preparation. Use non-fluorescent α-cyanoacrylate ethyl instant strong adhesive to stick the polished rock block or rock column to a glass slide, and cut it into a straight rock slice with a thickness of 1.5mm (the first set thickness). Then, according to the rock type and particle size, the rock slice is finely ground to a thickness of 50μm to 100μm (the second set thickness).

[0157] In the embodiments of the present disclosure, before the laser scanning of the shale oil rock is performed, one or more of the energy intensity of the laser and / or the laser wavelength and / or the pinhole value and / or the scanning mode and / or the line accumulation and / or the surface average and / or the scanning speed and / or the top and bottom of the scan and / or the scanning thickness and / or the number of scanning layers and / or the interlayer spacing are configured.

[0158] In the embodiments of the present disclosure and other possible embodiments, the prepared oil-bearing rock slice is placed on the stage of a laser confocal microscope at a set room temperature (e.g., 20° C.), the laser energy intensity is set to 30%, a 488 nm wavelength laser is selected as a light source to excite the oil-bearing rock slice sample (a rock slice of shale oil or a rock of shale oil), a filter corresponding to the 488 nm wavelength laser is selected, a default pinhole value is selected, an XYZ scanning mode is selected, line accumulation is set to 4, surface average is set to 2, the scanning speed is set to 100 Hz, and a unidirectional scan is performed. According to the lithology and particle size of the rock slice, at least 10 μm above and below the maximum focal plane is used as the top and bottom of the scan, the scanning thickness is preferably greater than 20 μm, the number of scanning layers is preferably greater than 20 layers, and the interlayer spacing is preferably less than 1 μm. The rock slice is scanned in the order of point by point, line by line, surface by surface, and layer by layer.

[0159] In the embodiments of the present disclosure and other possible embodiments, the reflected light signal and the transmitted light signal corresponding to the set wavelength (for example, 488nm) are collected, and the corresponding first fluorescence signal and second fluorescence signal are collected based on the first set receiving band and the second set receiving band, respectively. Specifically, the reflected light signal at a wavelength of 488nm is collected to obtain the matrix pores (for example, the gaps between the skeleton particles) corresponding to the rock slice; the transmitted light signal at a wavelength of 488nm is collected to obtain the distribution characteristics of the cracks (for example, reticular cracks, foliation cracks, etc.) corresponding to the rock slice. At the same time, the first fluorescence signal corresponding to the light crude oil receiving band range (the first set receiving band) determined above to represent the light components of crude oil is collected, and the second fluorescence signal of the heavy crude oil receiving band range (the second set receiving band) determined above to represent the heavy components of crude oil is collected.

[0160] In the embodiments of the present disclosure and other possible embodiments, when collecting data (reflected light signal, transmitted light signal, first fluorescence signal and second fluorescence signal), adjust the signal intensity gain knob of each channel so that the signal intensity of each channel is at the maximum value below overexposure, and the gain intensity of the receiving channel of the light and heavy components of crude oil should be consistent. After the scanning is completed, save the three-dimensional data volume (the three-dimensional data volume corresponding to the matrix pores and cracks), and record the scanned field of view length μm, width μm and height μm. More specifically, collect the reflected light signal at a wavelength of 488nm to obtain the matrix pores of the rock, and adjust the receiving gain knob appropriately to ensure that the matrix pores are fully displayed and collected completely; at the same time, collect the transmitted light signal to obtain the distribution characteristics of microcracks (reticular cracks, lamellae cracks), and adjust the receiving gain knob appropriately during collection to ensure that the cracks (reticular cracks, lamellae cracks, etc.) are all displayed and collected completely.

[0161] Step S1103: determining the distribution of the light component and / or heavy component based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal respectively.

[0162] In the embodiments of the present disclosure and other possible embodiments, light and heavy components of crude oil in fractures (reticular fractures, foliation fractures, etc.) and matrix pores are quantitatively analyzed (separation of matrix pores and fractures, separation of light and heavy components of crude oil).

[0163] In an embodiment of the present disclosure, the method for determining the distribution of the light component and / or heavy component based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal, respectively, includes: determining a first distribution characteristic map corresponding to the matrix pores of the rock and a second distribution characteristic map corresponding to the cracks based on the reflected light signal and the transmitted light signal, respectively; configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic map and the second distribution characteristic map, respectively, to determine the distribution of the light component and the heavy component, respectively.

[0164] In the embodiments of the present disclosure and other possible embodiments, the saved three-dimensional data volume (the three-dimensional data volume corresponding to the matrix pores and fractures) is imported into the multi-dimensional microscopic image analysis software, and the three-dimensional surface modeling is performed on the reflected light signal to obtain a first distribution characteristic map of the rock matrix pores; the surface modeling is performed on the transmitted light signal to obtain a second distribution characteristic map of the fracture (reticular fractures, foliation fractures, etc.) micro-area.

[0165] In the embodiments of the present disclosure and other possible embodiments, the method of respectively configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic map and the second distribution characteristic map to respectively determine the distribution of the light component and the heavy component includes: superimposing the first fluorescence signal with the first distribution characteristic map to obtain a first superimposed fluorescence distribution characteristic map corresponding to the slit; superimposing the second fluorescence signal with the second distribution characteristic map to obtain a second superimposed fluorescence distribution characteristic map corresponding to the matrix pores; and modeling the first superimposed fluorescence distribution characteristic map and the second superimposed fluorescence distribution characteristic map to obtain the distribution of the light component and the heavy component.

[0166] In the embodiments of the present disclosure and other possible embodiments, the first fluorescence signal corresponding to the light component of crude oil is superimposed with the first distribution characteristic map after surface modeling corresponding to the lamellae fractures and reticular fractures, and the first superimposed fluorescence distribution characteristic map corresponding to the fractures is obtained by using the co-localization mode to characterize the distribution of the light component. Similarly, the second fluorescence signal corresponding to the heavy component of crude oil is superimposed with the second distribution characteristic map after surface modeling corresponding to the matrix pores, and the second superimposed fluorescence distribution characteristic map corresponding to the matrix pores is obtained by using the co-localization mode to characterize the distribution of the heavy component.

[0167] More specifically, in the embodiments of the present disclosure and other possible embodiments, the first fluorescence signal corresponding to the light component of crude oil is superimposed with the first distribution characteristic map after surface modeling corresponding to the lamina fractures and reticular fractures, and the co-localization mode is used to obtain the first superimposed fluorescence distribution characteristic map corresponding to the fracture, and the fluorescence signal of the light component of crude oil in the lamina fractures, reticular fractures and other fractures in the first superimposed fluorescence distribution characteristic map is marked, and all the fluorescence signals of the light component of crude oil outside the lamina fractures, reticular fractures and other fractures are removed, so that the fluorescence signal of the light component of crude oil in the lamina fractures, reticular fractures and other fractures can be extracted, and the surface modeling is performed on the fluorescence signal to obtain the first superimposed fluorescence distribution characteristic map. Similarly, the second fluorescence signal corresponding to the heavy component of crude oil is superimposed with the second distribution characteristic map after surface modeling corresponding to the matrix pores, and the co-localization mode is used to obtain the second superimposed fluorescence distribution characteristic map corresponding to the matrix pores, and the fluorescence signal of the light component of crude oil in the matrix pores in the second superimposed fluorescence distribution characteristic map is marked, and all the fluorescence signals of the light component of crude oil outside the matrix pores are removed, so that the fluorescence signal of the light component of crude oil in the matrix pores can be extracted, and the surface modeling is performed on the fluorescence signal to obtain the first superimposed fluorescence distribution characteristic map.

[0168] In an embodiment of the present disclosure, before configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic diagram and the second distribution characteristic diagram, respectively, and determining the distribution of the light component and the heavy component, respectively, the first distribution characteristic diagram and the second distribution characteristic diagram are corrected to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram, respectively; configuring the first fluorescence signal and the second fluorescence signal in the first corrected distribution characteristic diagram and the second corrected distribution characteristic diagram, respectively, and determining the distribution of the light component and the heavy component, respectively; wherein, the method of correcting the first distribution characteristic diagram and the second distribution characteristic diagram to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram, respectively, includes: superimposing the first distribution characteristic diagram and the second distribution characteristic diagram to obtain a third distribution characteristic diagram; correcting the first distribution characteristic diagram and the second distribution characteristic diagram based on the overlapping area of ​​the third distribution characteristic diagram to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram.

[0169] In an embodiment of the present disclosure, the method of correcting the first distribution characteristic map and the second distribution characteristic map based on the overlapping area of ​​the third distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map includes: performing edge detection on the overlapping area to obtain edge features; if the edge features are lines, deleting the edge features that are lines in the first distribution characteristic map to obtain a first corrected distribution characteristic map; if the edge features are points, deleting the edge features that are points in the second distribution characteristic map to obtain a second corrected distribution characteristic map.

[0170] In the embodiments of the present disclosure and other possible embodiments, the algorithm for edge detection of the overlapping area can be configured as one or more of the Roberts Cross operator, Prewitt operator, Sobel operator, Kirsch operator, compass operator, Marr-Hildreth, the zero crossing of the second-order derivative in the gradient direction, Canny operator, and Laplacian operator.

[0171] In the embodiments of the present disclosure and other possible embodiments, the first distribution characteristic map and the second distribution characteristic map after surface modeling of the reflected light signal and the transmitted light signal are superimposed to obtain a third distribution characteristic map; the overlapping area of ​​the first distribution characteristic map and the second distribution characteristic map in the third distribution characteristic map is determined; the first distribution characteristic map and the second distribution characteristic map are corrected based on the overlapping area to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map; the first distribution characteristic map and the second distribution characteristic map are further analyzed and split and reorganized to respectively split out the first corrected distribution characteristic map and the second corrected distribution characteristic map corresponding to the foliation cracks, reticular cracks and matrix pores, and then respectively obtain the micro-area distribution of the foliation cracks, reticular cracks and matrix pores.

[0172] In the embodiments of the present disclosure and other possible embodiments, the method of respectively configuring the first fluorescence signal and the second fluorescence signal in the first corrected distribution characteristic diagram and the second corrected distribution characteristic diagram to respectively determine the distribution of the light component and the heavy component includes: superimposing the first corrected fluorescence signal with the first distribution characteristic diagram to obtain a first superimposed fluorescence distribution characteristic diagram corresponding to the slit; superimposing the second fluorescence signal with the second corrected distribution characteristic diagram to obtain a second superimposed fluorescence distribution characteristic diagram corresponding to the matrix pores; and modeling the first superimposed fluorescence distribution characteristic diagram and the second superimposed fluorescence distribution characteristic diagram to obtain the distribution of the light component and the heavy component.

[0173] In the embodiments of the present disclosure and other possible embodiments, the first fluorescence signal corresponding to the light component of crude oil is superimposed with the first corrected distribution characteristic map after surface modeling corresponding to the lamellae fractures and reticular fractures, and the first superimposed fluorescence distribution characteristic map corresponding to the fractures is obtained by using the co-localization mode to characterize the distribution of the light component. Similarly, the second fluorescence signal corresponding to the heavy component of crude oil is superimposed with the second corrected distribution characteristic map after surface modeling corresponding to the matrix pores, and the second superimposed fluorescence distribution characteristic map corresponding to the matrix pores is obtained by using the co-localization mode to characterize the distribution of the heavy component.

[0174] More specifically, in the embodiments of the present disclosure and other possible embodiments, the first fluorescence signal corresponding to the light component of crude oil is superimposed with the first corrected distribution characteristic map after surface modeling corresponding to the lamina fractures and reticular fractures, and the first superimposed fluorescence distribution characteristic map corresponding to the fractures is obtained by using the co-localization mode, and the fluorescence signal of the light component of crude oil in the lamina fractures, reticular fractures and other fractures in the first superimposed fluorescence distribution characteristic map is marked, and all the fluorescence signals of the light component of crude oil outside the lamina fractures, reticular fractures and other fractures are removed, so that the fluorescence signal of the light component of crude oil in the lamina fractures, reticular fractures and other fractures can be extracted, and the surface modeling is performed on the fluorescence signal to obtain the first superimposed fluorescence distribution characteristic map. Similarly, the second fluorescence signal corresponding to the heavy component of crude oil is superimposed with the second corrected distribution characteristic map after surface modeling corresponding to the matrix pores, and the co-localization mode is used to obtain the second superimposed fluorescence distribution characteristic map corresponding to the matrix pores, and the fluorescence signal of the light component of crude oil in the matrix pores in the second superimposed fluorescence distribution characteristic map is marked, and all the fluorescence signals of the light component of crude oil outside the matrix pores are removed, so that the fluorescence signal of the light component of crude oil in the matrix pores can be extracted, and the surface modeling is performed on the fluorescence signal to obtain the first superimposed fluorescence distribution characteristic map.

[0175] Step S1104: Based on the distribution of the light components and / or heavy components, determine the migration rate of the light components and / or heavy components of the shale oil.

[0176] In the embodiment of the present disclosure and other possible embodiments, the fluorescence signals of the light components of crude oil in the lamina fractures, reticular fractures and other fractures in the first superimposed fluorescence distribution characteristic diagram are marked, and all the fluorescence signals of the light components of crude oil outside the lamina fractures, reticular fractures and other fractures are removed, so that the fluorescence signals of the light components of crude oil in the lamina fractures, reticular fractures and other fractures can be extracted, and the surface modeling of the fluorescence signals is performed to obtain the first superimposed fluorescence distribution characteristic diagram, and the volume data (first volume V ) of the light components of crude oil in the lamina fractures, reticular fractures and other fractures are obtained and exported by using the data statistics function of the Imaris multi-dimensional microscopic image analysis software. 轻1 Similarly, the volume data of the light components of crude oil in the matrix pores (the third volume V 轻2 ). Obtain and export the volume data of the heavy components of crude oil in the lamina fractures, reticular fractures, etc. (the second volume V 重1 ). Obtain and derive the volume data of the heavy components of crude oil in the matrix pores (the fourth volume V 重2 ).

[0177] In an embodiment of the present disclosure, the method for determining the migration rate of the light component and / or the heavy component of the shale oil based on the distribution of the light component and / or the heavy component comprises: determining the first volume V of the light component and / or the heavy component in the fracture based on the distribution of the light component and / or the heavy component, respectively. 轻1 and the second volume V重1 , and respectively determine the third volume V of the light component and / or heavy component in the matrix pores 轻2 and the fourth volume V 重2 Based on the first volume V 轻1 and the third volume V 轻2 , determine the migration rate of the light component of the shale oil; based on the second volume V 重1 and the fourth volume V 重2 , determine the migration rate of the heavy components of the shale oil.

[0178] In an embodiment of the present disclosure, the first volume V 轻1 and the third volume V 轻2 The method for determining the migration rate of the light component of the shale oil comprises: the first volume V 轻1 Divide by the first volume V 轻1 and the third volume V 轻2 The migration rate of the light component of the shale oil is obtained by summing the second volume V 重1 and the fourth volume V 重2 The method for determining the migration rate of the heavy component of the shale oil comprises: the first volume V 重1 Divide by the first volume V 重1 and the third volume V 重2 The sum of the above mentioned components can be used to obtain the migration rate of the heavy components of the shale oil.

[0179] In the embodiments of the present disclosure and other possible embodiments, the migration rates of the light and heavy components of the crude oil can be obtained by dividing the volume of the light (heavy) components of the crude oil in the lamina fractures, network fractures and other fractures by the sum of the volumes of the light (heavy) components of the crude oil in the lamina fractures, network fractures and other fractures and in the matrix pores.

[0180] The migration rate of light components in crude oil is calculated according to formula (1):

[0181]

[0182] Where:

[0183] η 1 ——Migration rate of light components in crude oil, in % .

[0184] The migration rate of heavy components in crude oil is calculated according to formula (2):

[0185]

[0186] Where:

[0187] η 2 ——Migration rate of heavy components of crude oil, in % .

[0188] The present disclosure also proposes a shale oil analysis method, including: the carbon number determination method as described above; and / or, the migration analysis method as described above; and respectively configuring the first component and the second component as the light component and the heavy component of the shale oil, and obtaining a first set receiving band and a second set receiving band corresponding to the critical carbon number of the light component and the heavy component; performing laser scanning on the rock of the shale oil, collecting the reflected light signal and the transmitted light signal corresponding to the set wavelength, and respectively based on the first set receiving band and the second set receiving band The receiving band collects the corresponding first fluorescence signal and second fluorescence signal; based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal, the distribution of the light component and / or the heavy component is respectively determined; based on the distribution of the light component and / or the heavy component, the light-to-weight ratio of the crude oil in the fractures of the shale oil and / or the light-to-weight ratio of the crude oil in the matrix pores and / or the light-to-weight ratio of the shale oil and / or the volume percentage content in the shale oil fractures and / or the volume percentage content in the matrix pores of the shale oil and / or the oil saturation are determined.

[0189] In an embodiment of the present disclosure, the method for determining the distribution of the light component and / or heavy component based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal, respectively, includes: determining a first distribution characteristic map corresponding to the matrix pores of the rock and a second distribution characteristic map corresponding to the cracks based on the reflected light signal and the transmitted light signal, respectively; configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic map and the second distribution characteristic map, respectively, to determine the distribution of the light component and the heavy component, respectively.

[0190] In the embodiments of the present disclosure, before the laser scanning of the shale oil rock is performed, one or more of the energy intensity of the laser and / or the laser wavelength and / or the pinhole value and / or the scanning mode and / or the line accumulation and / or the surface average and / or the scanning speed and / or the top and bottom of the scan and / or the scanning thickness and / or the number of scanning layers and / or the interlayer spacing are configured.

[0191] In the embodiments of the present disclosure and other possible embodiments, the prepared oil-bearing rock slice is placed on the stage of a laser confocal microscope at a set room temperature (e.g., 20° C.), the laser energy intensity is set to 30%, a 488 nm wavelength laser is selected as a light source to excite the oil-bearing rock slice sample (a rock slice of shale oil or a rock of shale oil), a filter corresponding to the 488 nm wavelength laser is selected, a default pinhole value is selected, an XYZ scanning mode is selected, line accumulation is set to 4, surface average is set to 2, the scanning speed is set to 100 Hz, and a unidirectional scan is performed. According to the lithology and particle size of the rock slice, at least 10 μm above and below the maximum focal plane is used as the top and bottom of the scan, the scanning thickness is preferably greater than 20 μm, the number of scanning layers is preferably greater than 20 layers, and the interlayer spacing is preferably less than 1 μm. The rock slice is scanned in the order of point by point, line by line, surface by surface, and layer by layer.

[0192] In an embodiment of the present disclosure, before laser scanning is performed on the shale oil rock, the shale oil rock is prepared, including: obtaining an oil-containing shale rock sample, cutting the shale rock sample into rock blocks according to a set first size or cutting the rock columns according to a set second size; sticking the rock blocks or rock columns to a glass slide, and cutting the rock blocks or rock columns into straight rock slices according to a first set thickness; and grinding the straight rock slices to a second set thickness to obtain the shale oil rock for laser scanning.

[0193] In the embodiments of the present disclosure and other possible embodiments, an oil-bearing shale rock sample is selected, and rock blocks or rock columns are cut. The size of the rock block should be 25mm×25mm×5mm (set the first size), and the diameter of the rock column should be 25mm×5mm (set the second size). The sample cannot be soaked in organic solvents before preparation. Use non-fluorescent α-cyanoacrylate ethyl instant strong adhesive to stick the polished rock block or rock column to a glass slide, and cut it into a straight rock slice with a thickness of 1.5mm (the first set thickness). Then, according to the rock type and particle size, the rock slice is finely ground to a thickness of 50μm to 100μm (the second set thickness).

[0194] Before configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic diagram and the second distribution characteristic diagram, respectively, and determining the distribution of the light component and the heavy component, respectively, the first distribution characteristic diagram and the second distribution characteristic diagram are corrected to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram, respectively, and the first fluorescence signal and the second fluorescence signal are configured in the first corrected distribution characteristic diagram and the second corrected distribution characteristic diagram, respectively, to determine the distribution of the light component and the heavy component, respectively; wherein, the method of correcting the first distribution characteristic diagram and the second distribution characteristic diagram to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram, respectively, includes: superimposing the first distribution characteristic diagram and the second distribution characteristic diagram to obtain a third distribution characteristic diagram; correcting the first distribution characteristic diagram and the second distribution characteristic diagram based on the overlapping area of ​​the third distribution characteristic diagram to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram.

[0195] The method of correcting the first distribution feature map and the second distribution feature map based on the overlapping area of ​​the third distribution feature map to obtain a first corrected distribution feature map and a second corrected distribution feature map includes: performing edge detection on the overlapping area to obtain edge features; if the edge features are lines, deleting the edge features that are lines in the first distribution feature map to obtain the first corrected distribution feature map; if the edge features are points, deleting the edge features that are points in the second distribution feature map to obtain the second corrected distribution feature map.

[0196] In the embodiments of the present disclosure and other possible embodiments, the algorithm for edge detection of the overlapping area can be configured as one or more of the Roberts Cross operator, Prewitt operator, Sobel operator, Kirsch operator, compass operator, Marr-Hildreth, the zero crossing of the second-order derivative in the gradient direction, Canny operator, and Laplacian operator.

[0197] In an embodiment of the present disclosure, the method for determining the light-to-weight ratio of crude oil in the fracture of the shale oil and / or the light-to-weight ratio of crude oil in the matrix pores and / or the light-to-weight ratio of the shale oil and / or the volume percentage content in the shale oil fracture and / or the volume percentage content in the shale oil matrix pores and / or the oil saturation based on the distribution of the light component and / or the heavy component comprises: determining the first volume V of the light component and / or the heavy component in the fracture based on the distribution of the light component and / or the heavy component, respectively. 轻1 and the second volume V 重1 , and respectively determine the third volume V of the light component and / or heavy component in the matrix pores 轻2 and the fourth volume V 重2Based on the first volume V 轻1 and the second volume V 重1 , determine the crude oil weight ratio in the seam of the shale oil; based on the third volume V 轻2 and the fourth volume V 重2 , determine the crude oil weight ratio in the matrix pores of the shale oil; based on the first volume V 轻1 , the second volume V 重1 , the third volume V 轻2 and the fourth volume V 重2 , determine the weight ratio of the shale oil; based on the first volume V 轻1 , the second volume V 重1 and the rock volume of the shale oil, determine the volume percentage content in the shale oil fracture; based on the third volume V 轻2 、The fourth volume V 重2 and the rock volume of the shale oil, determine the volume percentage in the pores of the shale oil matrix; based on the volume percentage in the shale oil fractures, the volume percentage in the shale oil fractures and the porosity of the shale oil rock, determine the oil saturation of the shale oil.

[0198] In an embodiment of the present disclosure, the first volume V 轻1 and the second volume V 重1 The method for determining the crude oil weight ratio in the seam of the shale oil comprises: the first volume V 轻1 Divide by the second volume V 重1 , and obtain the weight of the crude oil in the seams of the shale oil.

[0199] In the embodiments of the present disclosure and other possible embodiments, in the embodiments of the present disclosure and other possible embodiments, the volume of the light components of the crude oil in the lamina and reticular fractures is divided by the volume of the heavy components of the crude oil in the lamina and reticular fractures to obtain the light-to-heavy ratio of the crude oil in the lamina and reticular fractures.

[0200] The weight ratio of crude oil in lamination fractures and reticular fractures is calculated according to formula (3):

[0201]

[0202] Where:

[0203] ω 1 ——The weight ratio of crude oil in the lamina fractures and network fractures (the weight ratio of crude oil in the fractures of the shale oil);

[0204] V 轻1 ——The volume of the light components of crude oil in the interlaminar fracture network (the first volume V 轻1 ), the unit is cubic micrometer (μm 3 );

[0205] V 重1 ——The volume of the heavy components of crude oil in the interlaminar fracture network (the second volume V 重1 ), the unit is cubic micrometer (μm 3 ).

[0206] In an embodiment of the present disclosure, the third volume V 轻2 and the fourth volume V 重2 , a method for determining the crude oil weight ratio in the matrix pores of the shale oil, comprising: the third volume V 轻2 Divide by the fourth volume V 重2 , and obtain the crude oil weight ratio in the matrix pores of the shale oil.

[0207] In the embodiment of the present disclosure and other possible embodiments, the light-to-heavy ratio of the crude oil in the matrix pores may be obtained by dividing the volume of the light components of the crude oil in the matrix pores by the volume of the heavy components of the crude oil in the matrix pores.

[0208] The weight ratio of crude oil in the matrix pores is calculated according to formula (4):

[0209]

[0210] Where:

[0211] ω 2 ——The weight ratio of crude oil in the matrix pores (the weight ratio of crude oil in the matrix pores of the shale oil);

[0212] V 轻2 ——The volume of light components of crude oil in the matrix pores (three volumes V 轻2 ), the unit is cubic micrometer (μm 3 );

[0213] V 重2 ——The volume of the heavy components of crude oil in the matrix pores (V 重2 ), the unit is cubic micrometer (μm 3 ).

[0214] In an embodiment of the present disclosure, the first volume V 轻1 , the second volume V 重1 , the third volume V 轻2 and the fourth volume V 重2 , a method for determining the weight ratio of the shale oil comprises: 轻1 and the third volume V 轻2 Sum, get the first total; for the second volume V 重1 and the fourth volume V 重2 The sum is calculated to obtain a second sum; the first sum is divided by the second sum to obtain the light-to-heavy ratio of the shale oil.

[0215] In the embodiment of the present disclosure and other possible embodiments, the light-to-heavy ratio of the crude oil of the sample can be obtained by dividing the volume of the light components of the crude oil in the lamina, reticular fractures, and matrix pores by the volume of the heavy components of the crude oil.

[0216] The crude oil weight ratio of the sample is calculated according to formula (5):

[0217]

[0218] Where:

[0219] ω 3 ——The light-to-heavy ratio of the crude oil (shale oil) of the sample;

[0220] For example, the crude oil weight ratio (ω 1 ) is greater than the light-to-heavy ratio of crude oil in the matrix pores (ω 2 ), it indicates that the sample kerogen is in the mature hydrocarbon expulsion stage, and the light components of the generated crude oil preferentially migrate and gradually accumulate from the matrix pores to the adjacent lamellae and reticular fractures.

[0221] In an embodiment of the present disclosure, the first volume V 轻1 , the second volume V 重1 and the rock volume of the shale oil, and a method for determining the volume percentage content in the shale oil fracture, comprising: 轻1 and the second volume V 重1 After summing up, the values ​​are divided by the rock volume of the shale oil to obtain the volume percentage content in the shale oil fractures.

[0222] In the embodiments of the present disclosure and other possible embodiments, the volume percentage of the crude oil in the lamina fractures, network fractures, etc. can be obtained by adding the volumes of the light components and the heavy components of the crude oil in the lamina fractures, network fractures, etc. and dividing the sum by the rock volume corresponding to the field of view.

[0223] The volume percentage of crude oil in lamination fractures and reticular fractures is calculated according to formula (6):

[0224]

[0225] Where:

[0226] δ 1 ——The volume percentage of crude oil in fractures such as lamellae and reticular fractures (the volume percentage of crude oil in the shale oil fractures);

[0227] L——The length of rock corresponding to the field of view, in micrometers (μm);

[0228] W——The width of the rock corresponding to the field of view, in micrometers (μm);

[0229] H is the height of the rock corresponding to the field of view, in micrometers (μm).

[0230] In an embodiment of the present disclosure, the third volume V 轻2 、The fourth volume V 重2 and the rock volume of the shale oil, and a method for determining the volume percentage content in the pores of the shale oil matrix, comprising: 轻2 and the second volume V 重2 After the sum is calculated, it is divided by the rock volume of the shale oil to obtain the volume percentage content in the pores of the shale oil matrix.

[0231] In the embodiment of the present disclosure and other possible embodiments, the volume percentage of the crude oil in the matrix pores can be obtained by adding the volumes of the light components and the heavy components of the crude oil in the matrix pores and dividing the sum by the rock volume corresponding to the field of view.

[0232] The volume percentage of crude oil in the matrix pores is calculated according to formula (7):

[0233]

[0234] Where:

[0235] δ 2 ——The volume percentage of crude oil in the matrix pores (the volume percentage of the shale oil in the matrix pores).

[0236] In an embodiment of the present disclosure, the method for determining the oil saturation of the shale oil based on the volume percentage content in the shale oil fracture, the volume percentage content in the shale oil fracture and the porosity of the rock of the shale oil comprises: summing the volume percentage content in the shale oil fracture and the volume percentage content in the shale oil fracture, and dividing the sum by the porosity of the rock of the shale oil to obtain the oil saturation of the shale oil.

[0237] In the embodiment of the present disclosure and other possible embodiments, the oil saturation of the sample can be obtained by adding the volume percentages of light and heavy components of crude oil in lamellae, reticular fractures and matrix pores and dividing the sum by the porosity of the sample.

[0238] The oil saturation of the sample is calculated according to formula (8):

[0239]

[0240] Where:

[0241] τ——oil saturation of the sample;

[0242] φ——the porosity of the sample.

[0243] The execution subject of the method for determining the carbon number and migration analysis of shale oil may be a shale oil carbon number determination and migration analysis device, for example, the method for determining the carbon number and migration analysis of shale oil may be executed by a terminal device or a server or other processing device, wherein the terminal device may be a user equipment (UE), a mobile device, a user terminal, a terminal, a cellular phone, a cordless phone, a personal digital assistant (PDA), a handheld device, a computing device, a vehicle-mounted device, a wearable device, etc. In some possible implementations, the method for determining the carbon number and migration analysis of shale oil may be implemented by a processor calling computer-readable instructions stored in a memory.

[0244] Those skilled in the art will understand that in the above-mentioned method for determining the carbon number and analyzing the migration of shale oil in a specific implementation manner, the writing order of the steps does not mean a strict execution order and does not constitute any limitation on the implementation process. The specific execution order of each step should be determined by its function and possible internal logic.

[0245] In addition, the present disclosure also proposes a device, including: a shale oil carbon number determination device and / or a shale oil migration analysis device and / or a shale oil analysis device; wherein the shale oil carbon number determination device includes: a full hydrocarbon gas chromatography analysis unit, used to perform full hydrocarbon gas chromatography analysis on the shale oil to obtain contents corresponding to multiple different carbon numbers; a laser spectral analysis unit, used to perform laser spectral analysis on the shale oil to obtain the fluorescence intensity of the shale oil in different bands; a first determination unit, used to determine the first proportion of the first component and / or the second proportion of the second component corresponding to the set receiving band according to the fluorescence intensity in the different bands; a second determination unit, used to determine the first proportion / second proportion and the contents corresponding to the multiple different carbon numbers based on the first proportion / second proportion, Determine the critical carbon number of the first component and the second component corresponding to the shale oil; and / or, wherein the migration analysis device of shale oil includes: a shale oil carbon number determination device; and, a first acquisition unit, used to configure the first component and the second component as the light component and the heavy component of the shale oil, respectively, and obtain the first set receiving band and the second set receiving band corresponding to the critical carbon number of the light component and the heavy component; a first laser scanning unit, used to perform laser scanning on the rock of the shale oil, collect the reflected light signal and the transmitted light signal corresponding to the set wavelength, and collect the corresponding first fluorescence signal and the second fluorescence signal based on the first set receiving band and the second set receiving band respectively; a first distribution determination unit, used to determine the first fluorescence signal and the second fluorescence signal based on the reflected light signal and the second set receiving band a first light signal, a second light signal, a second fluorescence signal, and a second fluorescence signal for respectively determining the distribution of the light component and / or the heavy component; a migration rate analysis unit for determining the migration rate of the light component and / or the heavy component of the shale oil based on the distribution of the light component and / or the heavy component; and / or, wherein the shale oil analysis device comprises: a shale oil carbon number determination device and / or a shale oil migration analysis device; and a second acquisition unit for respectively configuring the first component and the second component as the light component and the heavy component of the shale oil, and acquiring the first set receiving band and the second set receiving band corresponding to the critical carbon number of the light component and the heavy component; a second laser scanning unit for performing a rock scan on the shale oil; Laser scanning is performed to collect reflected light signals and transmitted light signals corresponding to the set wavelength, and the corresponding first fluorescence signals and second fluorescence signals are collected based on the first set receiving band and the second set receiving band respectively; a second distribution determination unit is used to determine the distribution of the light component and / or the heavy component based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal respectively; an analysis unit is used to determine the light-to-weight ratio of crude oil in the fractures of the shale oil and / or the light-to-weight ratio of crude oil in the matrix pores and / or the light-to-weight ratio of the shale oil and / or the volume percentage content in the fractures of the shale oil and / or the volume percentage content in the matrix pores of the shale oil and / or the oil saturation based on the distribution of the light component and / or the heavy component.The critical carbon numbers of the first and second components corresponding to shale oil can be determined to solve the problem that the current results of the spatial distribution differences of crude oil components cannot be given, which restricts the research on crude oil enrichment laws.

[0246] In some embodiments, the functions or modules included in the device provided in the embodiments of the present disclosure can be used to execute the method for determining the carbon number and migration analysis of shale oil described in the above method embodiments. The specific implementation thereof can refer to the description of the above embodiment of the method for determining the carbon number and migration analysis of shale oil, and for the sake of brevity, it will not be repeated here.

[0247] The disclosed embodiment also proposes a computer-readable storage medium on which computer program instructions are stored, and when the computer program instructions are executed by a processor, the above-mentioned shale oil carbon number determination and migration analysis method is implemented. The computer-readable storage medium may be a non-volatile computer-readable storage medium. The critical carbon number of the first component and the second component corresponding to the shale oil can be determined to solve the problem that the current results of the spatial distribution differences of the components of crude oil cannot be given, which restricts the research on the enrichment law of crude oil.

[0248] The disclosed embodiment also proposes an electronic device, comprising: a processor; a memory for storing instructions executable by the processor; wherein the processor is configured to perform the above-mentioned carbon number determination and migration analysis method for shale oil. The electronic device can be provided as a terminal, a server or other form of equipment. The critical carbon number of the first component and the second component corresponding to the shale oil can be determined to solve the problem that the current results of the spatial distribution differences of the components of crude oil cannot be given, which restricts the research on the enrichment law of crude oil.

[0249] The present disclosure has at least the following beneficial effects: it solves the shortcoming that the fluorescent signal channels of light and heavy components of crude oil cannot be separated into different positions of lamina, reticular fractures and matrix pores, and can realize intuitive visualization observation of the micro-region distribution of different crude oil components, and can also realize quantitative analysis of the differential distribution of different crude oil components in micro-regions, determine the migration rate of light and heavy components of crude oil, and then finely characterize the oil content of mud shale samples, and quantitatively evaluate the fluidity of crude oil. Therefore, this technical method and device, electronic equipment and storage medium can realize the fine evaluation of the oil content and fluidity of shale oil, which is of great significance to the calculation of shale oil resources and recoverable reserves and the effective use of shale oil, and provide data support for the evaluation of mud shale reservoirs.

[0250] Figure 4 The distribution characteristics of light components and / or heavy components of crude oil in the corresponding matrix pores and / or fractures according to the embodiment of the present disclosure are shown. Figure 4 As shown, Figure 4 (a) is the distribution characteristic diagram of light components corresponding to matrix pores and cracks; Figure 4(b) is the distribution characteristic diagram of heavy components corresponding to matrix pores and cracks; Figure 4 (c) is the (corrected) distribution characteristic diagram of the light components of crude oil in the fracture; Figure 4 (d) is the (corrected) distribution characteristic diagram corresponding to the light components of heavy oil in the fracture; Figure 4 (e) is the (corrected) distribution characteristic diagram of the light components of crude oil in the matrix pores; Figure 4 (f) is the (corrected) distribution characteristic diagram corresponding to the light components of heavy oil in the matrix pores.

[0251] Comprehensive analysis shows that the light-to-heavy ratio (ω1) of crude oil in the sample's lamina fractures, reticular fractures, etc. is 1:0.67, the light-to-heavy ratio (ω2) of crude oil in the matrix pores is 1:2.17, the light-to-heavy ratio (ω3) of crude oil in this sample is 1:0.96, the migration rate of the light component of crude oil (η1) is 80.51%, the migration rate of the heavy component of crude oil (η2) is 56.11%, the volume percentage of crude oil in lamina fractures, reticular fractures, etc. (δ1) is 1.93%, the volume percentage of crude oil in the matrix pores (δ2) is 0.89%, and the oil saturation of this sample is 74.21%. From the above data, it can be seen that the shale sample has good overall oil content, the light-to-heavy ratio of crude oil in fractures such as lamina fractures and network fractures is greater than the first set value (for example, 1), the migration rate of the light component of crude oil exceeds or is greater than the second set value (for example, 80%), and the migration rate of the heavy component of crude oil exceeds or is greater than the third set value (for example, 56%), so the crude oil also has good fluidity.

[0252] In summary, the present disclosure provides a method, device, equipment and medium for determining the carbon number and migration analysis of shale oil in order to solve the problem of visualizing the micro-area distribution of different components of shale oil and fine quantitative analysis of crude oil components. This effective analysis method and device are used to observe and analyze shale samples of a certain size. The present invention mainly uses a core cutter and a polisher to quickly cut and polish rock samples to prepare oil-bearing rock slices, uses a laser confocal microscope to collect reflected light signals, and obtains three-dimensional data of rock morphology; collects transmitted light signals to obtain three-dimensional data of the distribution characteristics of reticular fractures, lamella fractures and other fractures and matrix pores, and collects different fluorescence signals to obtain three-dimensional data of light and heavy components of crude oil respectively. LASAF software and Imaris multidimensional microscopic image analysis software are used to quantitatively analyze the light and heavy components of crude oil in lamina, reticular fractures, and matrix pores, and then the volume percentage of light and heavy components of crude oil in lamina, reticular fractures, and matrix pores, the migration rate of light and heavy components of crude oil, oil saturation, etc. are obtained, and finally a quantitative analysis method for the differential distribution of different crude oil components in micro-regions is formed.

[0253] Figure 58 is a block diagram of an electronic device 800 according to an exemplary embodiment. For example, the electronic device 800 may be a mobile phone, a computer, a digital broadcast terminal, a messaging device, a game console, a tablet device, a medical device, a fitness device, a personal digital assistant, or the like.

[0254] Reference Figure 5 , the electronic device 800 may include one or more of the following components: a processing component 802 , a memory 804 , a power component 806 , a multimedia component 808 , an audio component 810 , an input / output (I / O) interface 812 , a sensor component 814 , and a communication component 816 .

[0255] The processing component 802 generally controls the overall operation of the electronic device 800, such as operations associated with display, phone calls, data communications, camera operations, and recording operations. The processing component 802 may include one or more processors 820 to execute instructions to complete all or part of the steps of the above-mentioned method. In addition, the processing component 802 may include one or more modules to facilitate the interaction between the processing component 802 and other components. For example, the processing component 802 may include a multimedia module to facilitate the interaction between the multimedia component 808 and the processing component 802.

[0256] The memory 804 is configured to store various types of data to support operations on the electronic device 800. Examples of such data include instructions for any application or method operating on the electronic device 800, contact data, phone book data, messages, pictures, videos, etc. The memory 804 can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic disk or optical disk.

[0257] The power supply component 806 provides power to the various components of the electronic device 800. The power supply component 806 may include a power management system, one or more power supplies, and other components associated with generating, managing, and distributing power for the electronic device 800.

[0258] The multimedia component 808 includes a screen that provides an output interface between the electronic device 800 and the user. In some embodiments, the screen may include a liquid crystal display (LCD) and a touch panel (TP). If the screen includes a touch panel, the screen may be implemented as a touch screen to receive input signals from the user. The touch panel includes one or more touch sensors to sense touch, slide, and gestures on the touch panel. The touch sensor may not only sense the boundaries of the touch or slide action, but also detect the duration and pressure associated with the touch or slide operation. In some embodiments, the multimedia component 808 includes a front camera and / or a rear camera. When the electronic device 800 is in an operating mode, such as a shooting mode or a video mode, the front camera and / or the rear camera may receive external multimedia data. Each front camera and rear camera may be a fixed optical lens system or have a focal length and optical zoom capability.

[0259] The audio component 810 is configured to output and / or input audio signals. For example, the audio component 810 includes a microphone (MIC), and when the electronic device 800 is in an operating mode, such as a call mode, a recording mode, and a voice recognition mode, the microphone is configured to receive an external audio signal. The received audio signal can be further stored in the memory 804 or sent via the communication component 816. In some embodiments, the audio component 810 also includes a speaker for outputting audio signals.

[0260] I / O interface 812 provides an interface between processing component 802 and peripheral interface modules, such as keyboards, click wheels, buttons, etc. These buttons may include but are not limited to: a home button, a volume button, a start button, and a lock button.

[0261] The sensor assembly 814 includes one or more sensors for providing various aspects of status assessment for the electronic device 800. For example, the sensor assembly 814 can detect the open / closed state of the electronic device 800, the relative positioning of the components, such as the display and keypad of the electronic device 800, and the sensor assembly 814 can also detect the position change of the electronic device 800 or a component of the electronic device 800, the presence or absence of contact between the user and the electronic device 800, the orientation or acceleration / deceleration of the electronic device 800, and the temperature change of the electronic device 800. The sensor assembly 814 may include a proximity sensor configured to detect the presence of a nearby object without any physical contact. The sensor assembly 814 may also include an optical sensor, such as a CMOS or CCD image sensor, for use in imaging applications. In some embodiments, the sensor assembly 814 may also include an acceleration sensor, a gyroscope sensor, a magnetic sensor, a pressure sensor, or a temperature sensor.

[0262] The communication component 816 is configured to facilitate wired or wireless communication between the electronic device 800 and other devices. The electronic device 800 can access a wireless network based on a communication standard, such as WiFi, 2G or 3G, or a combination thereof. In an exemplary embodiment, the communication component 816 receives a broadcast signal or broadcast-related information from an external broadcast management system via a broadcast channel. In an exemplary embodiment, the communication component 816 also includes a near field communication (NFC) module to facilitate short-range communication. For example, the NFC module can be implemented based on radio frequency identification (RFID) technology, infrared data association (IrDA) technology, ultra-wideband (UWB) technology, Bluetooth (BT) technology and other technologies.

[0263] In an exemplary embodiment, the electronic device 800 may be implemented by one or more application-specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field programmable gate arrays (FPGAs), controllers, microcontrollers, microprocessors, or other electronic components to perform the above methods.

[0264] In an exemplary embodiment, a non-volatile computer-readable storage medium is also provided, such as a memory 804 including computer program instructions, which can be executed by a processor 820 of an electronic device 800 to perform the above method.

[0265] Figure 6 1 is a block diagram of an electronic device 1900 according to an exemplary embodiment. For example, the electronic device 1900 may be provided as a server. Figure 6 , the electronic device 1900 includes a processing component 1922, which further includes one or more processors, and a memory resource represented by a memory 1932 for storing instructions executable by the processing component 1922, such as an application. The application stored in the memory 1932 may include one or more modules, each corresponding to a set of instructions. In addition, the processing component 1922 is configured to execute instructions to perform the above method.

[0266] The electronic device 1900 may also include a power supply component 1926 configured to perform power management of the electronic device 1900, a wired or wireless network interface 1950 configured to connect the electronic device 1900 to a network, and an input / output (I / O) interface 1958. The electronic device 1900 may operate based on an operating system stored in the memory 1932, such as Windows Server™, Mac OS X™, Unix™, Linux™, FreeBSD™ or the like.

[0267] In an exemplary embodiment, a non-volatile computer-readable storage medium is also provided, such as a memory 1932 including computer program instructions, which can be executed by the processing component 1922 of the electronic device 1900 to perform the above method.

[0268] The present disclosure may be a system, a method and / or a computer program product. The computer program product may include a computer-readable storage medium carrying computer-readable program instructions for causing a processor to implement various aspects of the present disclosure.

[0269] A computer-readable storage medium may be a tangible device that can hold and store instructions used by an instruction execution device. A computer-readable storage medium may be, for example, but not limited to, an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. More specific examples of computer-readable storage media (a non-exhaustive list) include: a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a static random access memory (SRAM), a portable compact disk read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanical encoding device, such as a punch card or a raised structure in a groove on which instructions are stored, and any suitable combination of the foregoing. As used herein, a computer-readable storage medium is not to be interpreted as a transient signal per se, such as a radio wave or other freely propagating electromagnetic wave, an electromagnetic wave propagating through a waveguide or other transmission medium (e.g., a light pulse through a fiber optic cable), or an electrical signal transmitted through a wire.

[0270] The computer-readable program instructions described herein can be downloaded from a computer-readable storage medium to each computing / processing device, or downloaded to an external computer or external storage device via a network, such as the Internet, a local area network, a wide area network, and / or a wireless network. The network can include copper transmission cables, optical fiber transmissions, wireless transmissions, routers, firewalls, switches, gateway computers, and / or edge servers. The network adapter card or network interface in each computing / processing device receives the computer-readable program instructions from the network and forwards the computer-readable program instructions for storage in the computer-readable storage medium in each computing / processing device.

[0271] The computer program instructions for performing the operation of the present disclosure may be assembly instructions, instruction set architecture (ISA) instructions, machine instructions, machine-related instructions, microcode, firmware instructions, state setting data, or source code or object code written in any combination of one or more programming languages, including object-oriented programming languages, such as Smalltalk, C++, etc., and conventional procedural programming languages, such as "C" language or similar programming languages. Computer-readable program instructions may be executed completely on a user's computer, partially on a user's computer, as an independent software package, partially on a user's computer, partially on a remote computer, or completely on a remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer via any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., using an Internet service provider to connect via the Internet). In some embodiments, an electronic circuit, such as a programmable logic circuit, a field programmable gate array (FPGA), or a programmable logic array (PLA), may be customized by utilizing the state information of the computer-readable program instructions, and the electronic circuit may execute the computer-readable program instructions, thereby realizing various aspects of the present disclosure.

[0272] Various aspects of the present disclosure are described herein with reference to the flowcharts and / or block diagrams of the methods, devices (systems) and computer program products according to the embodiments of the present disclosure. It should be understood that each box in the flowchart and / or block diagram and the combination of each box in the flowchart and / or block diagram can be implemented by computer-readable program instructions.

[0273] These computer-readable program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, thereby producing a machine, so that when these instructions are executed by the processor of the computer or other programmable data processing device, a device that implements the functions / actions specified in one or more boxes in the flowchart and / or block diagram is generated. These computer-readable program instructions can also be stored in a computer-readable storage medium, and these instructions cause the computer, programmable data processing device, and / or other equipment to work in a specific manner, so that the computer-readable medium storing the instructions includes a manufactured product, which includes instructions for implementing various aspects of the functions / actions specified in one or more boxes in the flowchart and / or block diagram.

[0274] Computer-readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device so that a series of operating steps are performed on the computer, other programmable data processing apparatus, or other device to produce a computer-implemented process, thereby causing the instructions executed on the computer, other programmable data processing apparatus, or other device to implement the functions / actions specified in one or more boxes in the flowchart and / or block diagram.

[0275] The flow chart and block diagram in the accompanying drawings show the possible architecture, function and operation of the system, method and computer program product according to multiple embodiments of the present disclosure. In this regard, each square box in the flow chart or block diagram can represent a part of a module, program segment or instruction, and a part of the module, program segment or instruction includes one or more executable instructions for realizing the specified logical function. In some alternative implementations, the function marked in the square box can also occur in a sequence different from that marked in the accompanying drawings. For example, two continuous square boxes can actually be executed substantially in parallel, and they can sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each square box in the block diagram and / or flow chart, and the combination of the square boxes in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs the specified function or action, or can be implemented with a combination of special hardware and computer instructions.

[0276] The embodiments of the present disclosure have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The selection of terms used herein is intended to best explain the principles of the embodiments, practical applications, or technical improvements to the technology in the market, or to enable other persons of ordinary skill in the art to understand the embodiments disclosed herein.

Claims

1. A method for determining the carbon number of shale oil, It is characterized in that include: The shale oil is subjected to total hydrocarbon gas chromatography analysis to obtain the contents corresponding to multiple different carbon numbers; Performing laser spectroscopy analysis on the shale oil to obtain the fluorescence intensity of the shale oil in different wavebands; According to the fluorescence intensity under the different bands, determine the first proportion of light components and the second proportion of heavy components corresponding to the set receiving band; wherein, according to the fluorescence intensity under the different bands, determine the first proportion of light components and the second proportion of heavy components corresponding to the set receiving band, including: constructing fluorescence intensity curves under the different bands; calculating the area ratio corresponding to the set receiving band and the fluorescence intensity curves under all different bands, and configuring the area ratio as the first proportion of light components and the second proportion of heavy components corresponding to the set receiving band; Based on the ratio of the first proportion to the second proportion and the contents corresponding to the multiple different carbon numbers, the critical carbon number of the light component and the heavy component corresponding to the shale oil is determined; wherein, the method of determining the critical carbon number of the light component and the heavy component corresponding to the shale oil based on the ratio of the first proportion to the second proportion and the contents corresponding to the multiple different carbon numbers includes: determining the first content proportion corresponding to each carbon number according to the contents corresponding to the multiple different carbon numbers; taking multiple set carbon numbers as the critical value, summing up the first content proportion under each carbon number of the multiple set carbon numbers and the carbon numbers below, respectively, to obtain multiple second content proportions corresponding to the multiple set carbon numbers and the carbon numbers below; respectively calculating the difference between the multiple second content proportions and the ratio of the first proportion to the second proportion, to obtain multiple proportion differences; determining the set carbon number corresponding to the minimum proportion difference among the multiple proportion differences as the critical carbon number of the light component and the heavy component corresponding to the shale oil.

2. The determination method according to claim 1, It is characterized in that The calculating the area ratio of the set receiving band to the fluorescence intensity curves of all different bands includes: Integrate the fluorescence intensity curves in all different bands to obtain the total area of ​​the fluorescence intensity curve; Integrating the fluorescence intensity curve under the set receiving band to obtain the area of ​​the fluorescence intensity curve; The ratio of the area of ​​the fluorescence intensity curve region to the total area of ​​the fluorescence intensity curve is calculated to obtain a corresponding area ratio.

3. A shale oil migration analysis method, include: The method for determining carbon number according to any one of claims 1 to 2, characterized in that: Obtaining a first set receiving band and a second set receiving band corresponding to the critical carbon numbers of the light component and the heavy component; Performing laser scanning on the rock of the shale oil, collecting reflected light signals and transmitted light signals corresponding to the set wavelength, and collecting corresponding first fluorescence signals and second fluorescence signals based on the first set receiving band and the second set receiving band, respectively; The distribution of the light component and the heavy component are determined based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal respectively; wherein, the distribution of the light component and the heavy component are determined based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal respectively, including: determining a first distribution characteristic map corresponding to the matrix pores of the rock and a second distribution characteristic map corresponding to the cracks based on the reflected light signal and the transmitted light signal respectively; configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic map and the second distribution characteristic map respectively, and determining the distribution of the light component and the heavy component respectively; Based on the distribution of the light components and the heavy components, the migration rates of the light components and the heavy components of the shale oil are determined; wherein, the method for determining the migration rates of the light components and the heavy components of the shale oil based on the distribution of the light components and the heavy components comprises: based on the distribution of the light components and the heavy components, respectively determining the first volume of the light components in the fractures and the second volume of the heavy components in the fractures, and respectively determining the third volume of the light components in the matrix pores and the fourth volume of the heavy components in the matrix pores; based on the first volume and the third volume, determining the migration rate of the light components of the shale oil; based on the second volume and the fourth volume, determining the migration rate of the heavy components of the shale oil.

4. The migration analysis method according to claim 3, It is characterized in that Before laser scanning the shale oil rock, one or more of the laser energy intensity, laser wavelength, pinhole value, scanning mode, line accumulation, surface average, scanning speed, scanning top and bottom, scanning thickness, number of scanning layers, and interlayer spacing are configured.

5. The migration analysis method according to any one of claims 3 or 4, It is characterized in that Before laser scanning the shale oil rock, the shale oil rock is prepared, including: Obtaining an oil-bearing shale rock sample, and cutting the shale rock sample into rock blocks according to a set first size or into rock columns according to a set second size; Gluing the rock block or rock column onto a glass slide, and cutting the rock block or rock column into straight rock slices according to a first set thickness; The straight rock slice is ground to a second set thickness to obtain the shale oil rock for laser scanning.

6. The migration analysis method according to any one of claims 3 or 4, It is characterized in that Before respectively configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic diagram and the second distribution characteristic diagram to respectively determine the distribution of the light component and the heavy component, the first distribution characteristic diagram and the second distribution characteristic diagram are corrected to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram, respectively, and the first fluorescence signal and the second fluorescence signal are respectively configured in the first corrected distribution characteristic diagram and the second corrected distribution characteristic diagram to respectively determine the distribution of the light component and the heavy component; Among them, the correcting the first distribution characteristic map and the second distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map, respectively, includes: superimposing the first distribution characteristic map and the second distribution characteristic map to obtain a third distribution characteristic map; correcting the first distribution characteristic map and the second distribution characteristic map based on the overlapping area of ​​the third distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map.

7. The migration analysis method according to claim 5, It is characterized in that Before respectively configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic diagram and the second distribution characteristic diagram to respectively determine the distribution of the light component and the heavy component, the first distribution characteristic diagram and the second distribution characteristic diagram are corrected to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram, respectively, and the first fluorescence signal and the second fluorescence signal are respectively configured in the first corrected distribution characteristic diagram and the second corrected distribution characteristic diagram to respectively determine the distribution of the light component and the heavy component; Among them, the correcting the first distribution characteristic map and the second distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map, respectively, includes: superimposing the first distribution characteristic map and the second distribution characteristic map to obtain a third distribution characteristic map; correcting the first distribution characteristic map and the second distribution characteristic map based on the overlapping area of ​​the third distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map.

8. The migration analysis method according to any one of claims 3, 4 or 7, It is characterized in that The determining the migration rate of the light component of the shale oil based on the first volume and the third volume includes: The first volume is divided by the sum of the first volume and the third volume to obtain the migration rate of the light component of the shale oil.

9. The migration analysis method according to claim 5, It is characterized in that The determining the migration rate of the light component of the shale oil based on the first volume and the third volume includes: The first volume is divided by the sum of the first volume and the third volume to obtain the migration rate of the light component of the shale oil.

10. The migration analysis method according to claim 6, It is characterized in that The determining the migration rate of the light component of the shale oil based on the first volume and the third volume includes: The first volume is divided by the sum of the first volume and the third volume to obtain the migration rate of the light component of the shale oil.

11. The migration analysis method according to any one of claims 3, 4, 7, 9 or 10, It is characterized in that The determining the migration rate of the heavy component of the shale oil based on the second volume and the fourth volume includes: The second volume is divided by the sum of the second volume and the fourth volume to obtain the migration rate of the heavy component of the shale oil.

12. The migration analysis method according to claim 5, It is characterized in that The determining the migration rate of the heavy component of the shale oil based on the second volume and the fourth volume includes: The second volume is divided by the sum of the second volume and the fourth volume to obtain the migration rate of the heavy component of the shale oil.

13. The migration analysis method according to claim 6, It is characterized in that The determining the migration rate of the heavy component of the shale oil based on the second volume and the fourth volume includes: The second volume is divided by the sum of the second volume and the fourth volume to obtain the migration rate of the heavy component of the shale oil.

14. The migration analysis method according to claim 8, It is characterized in that The determining the migration rate of the heavy component of the shale oil based on the second volume and the fourth volume includes: The second volume is divided by the sum of the second volume and the fourth volume to obtain the migration rate of the heavy component of the shale oil.

15. The migration analysis method according to claim 6, It is characterized in that The method of correcting the first distribution characteristic map and the second distribution characteristic map based on the overlapping area of ​​the third distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map includes: Performing edge detection on the overlapping area to obtain edge features; If the edge feature is a line, the edge feature that is a line is deleted from the first distribution feature map to obtain a first corrected distribution feature map; If the edge feature is a point, the edge feature that is a point is deleted from the second distribution feature map to obtain a second corrected distribution feature map.

16. A method for analyzing shale oil. include: The method for determining carbon number according to any one of claims 1 to 2, characterized in that: Obtaining a first set receiving band and a second set receiving band corresponding to the critical carbon numbers of the light component and the heavy component; Performing laser scanning on the rock of the shale oil, collecting reflected light signals and transmitted light signals corresponding to the set wavelength, and collecting corresponding first fluorescence signals and second fluorescence signals based on the first set receiving band and the second set receiving band, respectively; The distribution of the light component and the heavy component are determined based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal respectively; wherein, the distribution of the light component and the heavy component are determined based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal respectively, including: determining a first distribution characteristic map corresponding to the matrix pores of the rock and a second distribution characteristic map corresponding to the cracks based on the reflected light signal and the transmitted light signal respectively; configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic map and the second distribution characteristic map respectively, and determining the distribution of the light component and the heavy component respectively; Based on the distribution of the light components and the heavy components, determine one or more of the crude oil light-to-weight ratio in the fractures of the shale oil, the crude oil light-to-weight ratio in the matrix pores, the light-to-weight ratio of the shale oil, the volume percentage content in the fractures of the shale oil, the volume percentage content in the matrix pores of the shale oil, and the oil saturation; Wherein, based on the distribution of the light component and the heavy component, determining one or more of the crude oil light-to-weight ratio in the fracture of the shale oil, the crude oil light-to-weight ratio in the matrix pores, the light-to-weight ratio of the shale oil, the volume percentage content in the fracture of the shale oil, the volume percentage content in the matrix pores of the shale oil, and the oil saturation includes: based on the distribution of the light component and the heavy component, respectively determining the first volume and the second volume of the light component and the heavy component in the fracture, and respectively determining the third volume and the fourth volume of the light component and the heavy component in the matrix pores; based on the first volume and the second volume, determining the crude oil light-to-weight ratio in the fracture of the shale oil; Based on the third volume and the fourth volume, the weight ratio of crude oil in the matrix pores of the shale oil is determined; based on the first volume, the second volume, the third volume and the fourth volume, the weight ratio of the shale oil is determined; based on the first volume, the second volume and the rock volume of the shale oil, the volume percentage content in the cracks of the shale oil is determined; based on the third volume, the fourth volume and the rock volume of the shale oil, the volume percentage content in the matrix pores of the shale oil is determined; based on the volume percentage content in the cracks of the shale oil, the volume percentage content in the matrix pores of the shale oil and the porosity of the shale oil rock, the oil saturation of the shale oil is determined.

17. The analysis method according to claim 16, It is characterized in that Before laser scanning the shale oil rock, one or more of the laser energy intensity, laser wavelength, pinhole value, scanning mode, line accumulation, surface averaging, scanning speed, scanning top and bottom, scanning thickness, number of scanning layers, or interlayer spacing are configured.

18. The analytical method according to any one of claims 16 or 17, It is characterized in that Before laser scanning the shale oil rock, the shale oil rock is prepared, including: Obtaining an oil-bearing shale rock sample, and cutting the shale rock sample into rock blocks according to a set first size or into rock columns according to a set second size; Gluing the rock block or rock column onto a glass slide, and cutting the rock block or rock column into straight rock slices according to a first set thickness; The straight rock slice is ground to a second set thickness to obtain the shale oil rock for laser scanning.

19. The analysis method according to claim 16, It is characterized in that Before respectively configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic diagram and the second distribution characteristic diagram to respectively determine the distribution of the light component and the heavy component, the first distribution characteristic diagram and the second distribution characteristic diagram are corrected to obtain a first corrected distribution characteristic diagram and a second corrected distribution characteristic diagram, respectively, and the first fluorescence signal and the second fluorescence signal are respectively configured in the first corrected distribution characteristic diagram and the second corrected distribution characteristic diagram to respectively determine the distribution of the light component and the heavy component; Among them, the correcting the first distribution characteristic map and the second distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map, respectively, includes: superimposing the first distribution characteristic map and the second distribution characteristic map to obtain a third distribution characteristic map; correcting the first distribution characteristic map and the second distribution characteristic map based on the overlapping area of ​​the third distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map.

20. The analytical method according to claim 16, It is characterized in that The step of determining the crude oil weight ratio in the seam of the shale oil based on the first volume and the second volume includes: The first volume is divided by the second volume to obtain the crude oil weight ratio in the seam of the shale oil.

21. The analytical method according to any one of claims 16 or 20, It is characterized in that The step of determining the crude oil weight ratio in the matrix pores of the shale oil based on the third volume and the fourth volume includes: The third volume is divided by the fourth volume to obtain the crude oil weight ratio in the matrix pores of the shale oil.

22. The analytical method according to claim 16, It is characterized in that The determining the light-to-heavy ratio of the shale oil based on the first volume, the second volume, the third volume, and the fourth volume includes: The first volume and the third volume are summed to obtain a first total; the second volume and the fourth volume are summed to obtain a second total; The first sum is divided by the second sum to obtain the light-to-heavy ratio of the shale oil.

23. The analytical method according to claim 21, It is characterized in that The determining the light-to-heavy ratio of the shale oil based on the first volume, the second volume, the third volume, and the fourth volume includes: The first volume and the third volume are summed to obtain a first total; the second volume and the fourth volume are summed to obtain a second total; The first sum is divided by the second sum to obtain the light-to-heavy ratio of the shale oil.

24. The analytical method according to any one of claims 16, 20, 22 or 23, It is characterized in that The determining the volume percentage content in the fracture of the shale oil based on the first volume, the second volume and the rock volume of the shale oil includes: The first volume and the second volume are summed and then divided by the rock volume of the shale oil to obtain the volume percentage of the shale oil in the fracture.

25. The analytical method according to claim 21, It is characterized in that The determining the volume percentage content in the fracture of the shale oil based on the first volume, the second volume and the rock volume of the shale oil includes: The first volume and the second volume are summed and then divided by the rock volume of the shale oil to obtain the volume percentage of the shale oil in the fracture.

26. The analytical method according to any one of claims 16 or 20 or 22 or 23 or 25, It is characterized in that The determining the volume percentage content in the matrix pores of the shale oil based on the third volume, the fourth volume and the rock volume of the shale oil comprises: The third volume and the fourth volume are summed and then divided by the rock volume of the shale oil to obtain the volume percentage of the shale oil in the matrix pores.

27. The analytical method according to claim 21, It is characterized in that The determining the volume percentage content in the matrix pores of the shale oil based on the third volume, the fourth volume and the rock volume of the shale oil comprises: The third volume and the fourth volume are summed and then divided by the rock volume of the shale oil to obtain the volume percentage of the shale oil in the matrix pores.

28. The analytical method according to claim 24, It is characterized in that The determining the volume percentage content in the matrix pores of the shale oil based on the third volume, the fourth volume and the rock volume of the shale oil comprises: The third volume and the fourth volume are summed and then divided by the rock volume of the shale oil to obtain the volume percentage of the shale oil in the matrix pores.

29. The analytical method according to any one of claims 16 or 20 or 22 or 23 or 25 or 27 or 28, It is characterized in that The determining of the oil saturation of the shale oil based on the volume percentage content in the cracks of the shale oil, the volume percentage content in the matrix pores of the shale oil and the porosity of the rock of the shale oil comprises: The oil saturation of the shale oil is obtained by summing the volume percentage content in the cracks of the shale oil and the volume percentage content in the matrix pores of the shale oil and dividing it by the porosity of the rock of the shale oil.

30. The analytical method according to claim 21, It is characterized in that The determining of the oil saturation of the shale oil based on the volume percentage content in the cracks of the shale oil, the volume percentage content in the matrix pores of the shale oil and the porosity of the rock of the shale oil comprises: The oil saturation of the shale oil is obtained by summing the volume percentage content in the cracks of the shale oil and the volume percentage content in the matrix pores of the shale oil and dividing it by the porosity of the rock of the shale oil.

31. The analytical method according to claim 24, It is characterized in that The determining of the oil saturation of the shale oil based on the volume percentage content in the cracks of the shale oil, the volume percentage content in the matrix pores of the shale oil and the porosity of the rock of the shale oil comprises: The oil saturation of the shale oil is obtained by summing the volume percentage content in the cracks of the shale oil and the volume percentage content in the matrix pores of the shale oil and dividing it by the porosity of the rock of the shale oil.

32. The analytical method according to claim 26, It is characterized in that The determining of the oil saturation of the shale oil based on the volume percentage content in the cracks of the shale oil, the volume percentage content in the matrix pores of the shale oil and the porosity of the rock of the shale oil comprises: The oil saturation of the shale oil is obtained by summing the volume percentage content in the cracks of the shale oil and the volume percentage content in the matrix pores of the shale oil and dividing it by the porosity of the rock of the shale oil.

33. The analytical method according to claim 19, It is characterized in that The method of correcting the first distribution characteristic map and the second distribution characteristic map based on the overlapping area of ​​the third distribution characteristic map to obtain a first corrected distribution characteristic map and a second corrected distribution characteristic map includes: Performing edge detection on the overlapping area to obtain edge features; If the edge feature is a line, the edge feature that is a line is deleted from the first distribution feature map to obtain a first corrected distribution feature map; If the edge feature is a point, the edge feature that is a point is deleted from the second distribution feature map to obtain a second corrected distribution feature map.

34. A method for analyzing shale oil. include: The migration analysis method according to any one of claims 3 to 15, characterized in that: Obtaining a first set receiving band and a second set receiving band corresponding to the critical carbon numbers of the light component and the heavy component; Performing laser scanning on the rock of the shale oil, collecting reflected light signals and transmitted light signals corresponding to the set wavelength, and collecting corresponding first fluorescence signals and second fluorescence signals based on the first set receiving band and the second set receiving band, respectively; The distribution of the light component and the heavy component are determined based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal respectively; wherein, the distribution of the light component and the heavy component are determined based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal respectively, including: determining a first distribution characteristic map corresponding to the matrix pores of the rock and a second distribution characteristic map corresponding to the cracks based on the reflected light signal and the transmitted light signal respectively; configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic map and the second distribution characteristic map respectively, and determining the distribution of the light component and the heavy component respectively; Based on the distribution of the light components and the heavy components, determine one or more of the crude oil light-to-weight ratio in the fractures of the shale oil, the crude oil light-to-weight ratio in the matrix pores, the light-to-weight ratio of the shale oil, the volume percentage content in the fractures of the shale oil, the volume percentage content in the matrix pores of the shale oil, and the oil saturation; Wherein, based on the distribution of the light component and the heavy component, determining one or more of the crude oil light-to-weight ratio in the fracture of the shale oil, the crude oil light-to-weight ratio in the matrix pores, the light-to-weight ratio of the shale oil, the volume percentage content in the fracture of the shale oil, the volume percentage content in the matrix pores of the shale oil, and the oil saturation includes: based on the distribution of the light component and the heavy component, respectively determining the first volume and the second volume of the light component and the heavy component in the fracture, and respectively determining the third volume and the fourth volume of the light component and the heavy component in the matrix pores; based on the first volume and the second volume, determining the crude oil light-to-weight ratio in the fracture of the shale oil; Based on the third volume and the fourth volume, the weight ratio of crude oil in the matrix pores of the shale oil is determined; based on the first volume, the second volume, the third volume and the fourth volume, the weight ratio of the shale oil is determined; based on the first volume, the second volume and the rock volume of the shale oil, the volume percentage content in the cracks of the shale oil is determined; based on the third volume, the fourth volume and the rock volume of the shale oil, the volume percentage content in the matrix pores of the shale oil is determined; based on the volume percentage content in the cracks of the shale oil, the volume percentage content in the matrix pores of the shale oil and the porosity of the shale oil rock, the oil saturation of the shale oil is determined.

35. A device, It is characterized in that include: A device for determining the carbon number of shale oil; wherein the device for determining the carbon number of shale oil comprises: A full hydrocarbon gas chromatography analysis unit is used to perform full hydrocarbon gas chromatography analysis on shale oil to obtain contents corresponding to multiple different carbon numbers; A laser spectrum analysis unit, used for performing laser spectrum analysis on the shale oil to obtain the fluorescence intensity of the shale oil in different wavebands; A first determination unit is used to determine a first proportion of light components and a second proportion of heavy components corresponding to a set receiving band according to the fluorescence intensities under the different bands; wherein, the determination of the first proportion of light components and the second proportion of heavy components corresponding to the set receiving band according to the fluorescence intensities under the different bands includes: constructing fluorescence intensity curves under the different bands; calculating an area ratio corresponding to the set receiving band and the fluorescence intensity curves under all different bands, and configuring the area ratio as the first proportion of light components and the second proportion of heavy components corresponding to the set receiving band; The second determination unit is used to determine the critical carbon number of the light components and heavy components corresponding to the shale oil based on the ratio of the first proportion to the second proportion and the contents corresponding to the multiple different carbon numbers; wherein, the determination of the critical carbon number of the light components and heavy components corresponding to the shale oil based on the ratio of the first proportion to the second proportion and the contents corresponding to the multiple different carbon numbers includes: determining the first content proportion corresponding to each carbon number according to the contents corresponding to the multiple different carbon numbers; taking multiple set carbon numbers as the critical, summing up the first content proportion under each carbon number of the multiple set carbon numbers and below, respectively, to obtain multiple second content proportions corresponding to the multiple set carbon numbers and below; respectively calculating the difference between the multiple second content proportions and the ratio of the first proportion to the second proportion, to obtain multiple proportion differences; determining the set carbon number corresponding to the minimum proportion difference among the multiple proportion differences as the critical carbon number of the light components and heavy components corresponding to the shale oil.

36. The device according to claim 35, It is characterized in that Also includes: A shale oil migration analysis device; wherein the shale oil migration analysis device comprises: A first acquisition unit is used to acquire a first set receiving band and a second set receiving band corresponding to the critical carbon numbers of the light component and the heavy component; wherein the critical carbon numbers of the light component and the heavy component are determined by the carbon number determination device of the shale oil; a first laser scanning unit, configured to perform laser scanning on the rock of the shale oil, collect reflected light signals and transmitted light signals corresponding to a set wavelength, and collect corresponding first fluorescence signals and second fluorescence signals based on the first set receiving band and the second set receiving band, respectively; A first distribution determination unit is used to determine the distribution of the light component and the heavy component based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal, respectively; wherein, the distribution of the light component and the heavy component is determined based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal, respectively, including: determining a first distribution characteristic map corresponding to the matrix pores of the rock and a second distribution characteristic map corresponding to the cracks based on the reflected light signal and the transmitted light signal, respectively; configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic map and the second distribution characteristic map, respectively, to determine the distribution of the light component and the heavy component, respectively; A migration rate analysis unit is used to determine the migration rates of the light components and heavy components of the shale oil based on the distribution of the light components and heavy components; wherein the method for determining the migration rates of the light components and heavy components of the shale oil based on the distribution of the light components and heavy components comprises: determining the first volume of the light components in the fractures and the second volume of the heavy components in the fractures, and determining the third volume of the light components in the matrix pores and the fourth volume of the heavy components in the matrix pores, respectively, based on the distribution of the light components and the heavy components; determining the migration rate of the light components of the shale oil based on the first volume and the third volume; and determining the migration rate of the heavy components of the shale oil based on the second volume and the fourth volume.

37. The device according to any one of claims 35 or 36, It is characterized in that Also includes: An analysis device for shale oil; wherein the analysis device for shale oil comprises: A second acquisition unit is used to acquire a first set receiving band and a second set receiving band corresponding to the critical carbon numbers of the light component and the heavy component; wherein the critical carbon numbers of the light component and the heavy component are determined by the carbon number determination device of the shale oil; a second laser scanning unit, configured to perform laser scanning on the rock of the shale oil, collect reflected light signals and transmitted light signals corresponding to a set wavelength, and collect corresponding first fluorescence signals and second fluorescence signals based on the first set receiving band and the second set receiving band, respectively; A second distribution determination unit is used to determine the distribution of the light component and the heavy component based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal, respectively; wherein, the distribution of the light component and the heavy component is determined based on the reflected light signal, the transmitted light signal, the first fluorescence signal and the second fluorescence signal, respectively, including: determining a first distribution characteristic map corresponding to the matrix pores of the rock and a second distribution characteristic map corresponding to the cracks based on the reflected light signal and the transmitted light signal, respectively; configuring the first fluorescence signal and the second fluorescence signal in the first distribution characteristic map and the second distribution characteristic map, respectively, to determine the distribution of the light component and the heavy component, respectively; An analysis unit is used to determine one or more of the light-to-weight ratio of crude oil in the fractures of the shale oil, the light-to-weight ratio of crude oil in the matrix pores, the light-to-weight ratio of the shale oil, the volume percentage content in the fractures of the shale oil, the volume percentage content in the matrix pores of the shale oil, and oil saturation based on the distribution of the light components and the heavy components; wherein, the determination of one or more of the light-to-weight ratio of crude oil in the fractures of the shale oil, the light-to-weight ratio of crude oil in the matrix pores, the light-to-weight ratio of the shale oil, the volume percentage content in the fractures of the shale oil, the volume percentage content in the matrix pores of the shale oil, and oil saturation based on the distribution of the light components and the heavy components includes: determining the first volume and the second volume of the light component and the heavy component in the fracture, respectively, and determining the volume of the light component and the heavy component in the fracture, respectively. The invention relates to a method for determining the light-to-weight ratio of crude oil in the fractures of the shale oil based on the first volume and the second volume; the light-to-weight ratio of crude oil in the matrix pores of the shale oil based on the third volume and the fourth volume; the light-to-weight ratio of shale oil based on the first volume, the second volume, the third volume and the fourth volume; the volume percentage of shale oil in fractures based on the first volume, the second volume and the rock volume of the shale oil; the volume percentage of shale oil in matrix pores based on the third volume, the fourth volume and the rock volume of the shale oil; the oil saturation of the shale oil based on the volume percentage of shale oil in fractures, the volume percentage of shale oil in matrix pores and the porosity of the shale oil rock.

38. An electronic device, It is characterized in that include: processor; a memory for storing processor-executable instructions; The processor is configured to call the instructions stored in the memory to execute the carbon number determination method according to any one of claims 1 to 2.

39. An electronic device, It is characterized in that include: processor; a memory for storing processor-executable instructions; The processor is configured to call the instructions stored in the memory to execute the migration analysis method as described in any one of claims 3-15.

40. An electronic device, It is characterized in that include: processor; a memory for storing processor-executable instructions; Wherein, the processor is configured to call the instructions stored in the memory to execute the analysis method as described in any one of claims 16-33.

41. A computer readable storage medium having computer program instructions stored thereon, It is characterized in that When the computer program instructions are executed by a processor, the carbon number determination method according to any one of claims 1 to 2 is implemented.

42. A computer readable storage medium having computer program instructions stored thereon, It is characterized in that When the computer program instructions are executed by a processor, the migration analysis method according to any one of claims 3 to 15 is implemented.

43. A computer readable storage medium having computer program instructions stored thereon, It is characterized in that When the computer program instructions are executed by a processor, the analysis method according to any one of claims 16 to 33 is implemented.

Images