Comprehensive classification, equipment and electronic equipment for deep shale oil reservoirs

Abstract

The present invention relates to a reservoir classification technology field, and is a deep shale oil reservoir comprehensive classification, device and electronic equipment, including judging the drilling fluid type of the deep shale oil reservoir to be classified; in response to the deep shale oil reservoir to be classified being a water-based drilling fluid shale oil reservoir, a gas logging evaluation classification method, a rock cuttings pyrolysis logging evaluation classification method and a rock brittleness evaluation classification method are used to classify the deep shale oil reservoir to be classified respectively; in response to the deep shale oil reservoir to be classified being an oil-based drilling fluid shale oil reservoir, a gas logging evaluation classification method, a rock brittleness evaluation classification method and a shale oil reservoir oil content evaluation classification method are used to classify the deep shale oil reservoir to be classified respectively. The present invention uses different comprehensive classification methods for the different characteristics of water-based drilling fluid shale oil reservoirs and oil-based drilling fluid shale oil reservoirs, so that the reservoir classification has stronger adaptability and more accurate classification results.

Description

Technical Field

[0001] The present invention relates to the technical field of reservoir classification, in particular to a comprehensive classification device, device and electronic equipment for deep shale oil reservoirs. Background Art

[0002] Shale oil refers to the oil resources contained in organic-rich mud shale formations, including oil in mud shale pores and cracks, and oil in dense carbonate rocks or clastic rock interlayers in mud shale formations. Shale oil resources are an important type of unconventional natural oil reservoirs. They refer to oil that is found in organic-rich nano-scale pore throat shale formations and are important unconventional oil and gas resources. The exploration and production of unconventional shale oil and gas in China has become a research hotspot in the field of unconventional oil and gas resources.

[0003] In the existing technology, reservoir evaluation is mostly achieved by using gas logging, rock pyrolysis, TOC, physical property index, mineral brittleness index and other parameters to independently or jointly establish deep shale oil reservoir classification standards. For example:

[0004] The first existing open patent document, with the publication number CN110965999B, discloses a method for finely identifying the dominant lithology of shale oil. According to the X-ray diffraction mineral composition data of shale layer cuttings or core samples, three main mineral components are selected to establish a chart to achieve fine division of shale lithology; according to gas logging parameters, rock pyrolysis analysis, and brittleness index data, a standard for dividing the dominant lithology of shale formations is established; the present invention uses X-ray diffraction mineral composition data to establish a chart to distinguish lithology, and combines gas logging parameters, rock pyrolysis analysis, and brittleness index data to distinguish the dominant lithology of shale formations. Firstly, the problem of unclear lithology naming of shale oil can be solved well, and the lithology can be accurately identified. Secondly, the dominant lithology data in the lithology can be identified for oil testing and layer selection. As a result, the engineering service needs are met, and a reliable basis is provided for the efficient development of oil and gas resources in shale formations.

[0005] The second existing public patent document, with the publication number CN115130309B, discloses a method, equipment, medium and terminal for determining the brittleness of shale reservoirs using a bonding degree, including: using the grey correlation method to analyze the correlation between the mineral composition, particle arrangement structure and brittleness index of shale; determining the influence of the mineral composition and particle arrangement structure on the brittleness of shale based on the correlation between the mineral composition, particle arrangement structure and brittleness index of shale; constructing a brittleness evaluation model based on mineral composition and structural directional entropy based on the influence of the mineral composition and particle arrangement structure on the brittleness of shale combined with the total organic carbon content of shale; and using the brittleness evaluation model based on mineral composition and structural directional entropy to evaluate the brittleness of shale. The method for determining the brittleness of shale reservoirs using a bonding degree of the present invention not only retains the influence of mineral components, but also considers the role of the directional arrangement structure of particles, can be better applied to oriented structural development mudstone, and can effectively predict the fracturing ability of shale.

[0006] The existing classification process of deep shale oil reservoirs has the following deficiencies:

[0007] (1) The existing water-based drilling fluid shale oil reservoirs and oil-based drilling fluid shale oil reservoirs use the same reservoir classification method. However, due to the different characteristics of water-based drilling fluid and oil-based drilling fluid, the two react differently to the conventional rock pyrolysis method. In water-based drilling fluid shale oil reservoirs, the conventional rock pyrolysis method can well reflect the oil content of the reservoir. However, since drilling with oil-based drilling fluid easily makes the pyrolysis parameters invalid, the conventional rock pyrolysis method will cause inaccurate evaluation of the oil content of oil-based drilling fluid shale oil reservoirs.

[0008] (2) The high pressure in deep formations causes the gas logging background value to be artificially high. During high-angle and directional drilling, the fluctuating value affects the gas logging accuracy. The density, temperature and bottom hole pressure difference of the drilling fluid affect the amount of formation oil and gas entering the drilling fluid, resulting in low accuracy of gas logging parameters, which affects the reservoir evaluation effect.

[0009] (3) The accuracy of the judgment of rock brittleness depends to a great extent on whether the selected mineral brittleness index calculation method is reasonable. Relevant personnel often choose the mineral brittleness index calculation method based on their work experience, which sometimes leads to deviations between the final result and the actual situation. The commonly used XRD mineral brittleness evaluation and classification method currently uses X-ray diffraction mineral (XRD) acquisition equipment to collect XRD spectra. The collected XRD spectra need to be decomposed to obtain the mineral content. Different geological conditions have different mineral types. Due to the differences in geological understanding among different technicians, the mineral types and contents obtained by decomposition are not accurate, resulting in the accuracy of mineral brittleness index calculation, which affects the reservoir evaluation effect. Summary of the invention

[0010] The present invention provides a comprehensive classification, device and electronic equipment for deep shale oil reservoirs, which overcome the shortcomings of the above-mentioned prior art and can effectively solve the problem that the existing water-based drilling fluid shale oil reservoirs and oil-based drilling fluid shale oil reservoirs use the same reservoir classification method, resulting in inaccurate oil content evaluation of oil-based drilling fluid shale oil reservoirs and inaccurate reservoir classification.

[0011] One of the technical solutions of the present invention is achieved by the following measures: A comprehensive classification method for deep shale oil reservoirs, comprising:

[0012] Determine the drilling fluid type of the deep shale oil reservoir to be classified;

[0013] In response to the deep shale oil reservoir to be classified being a water-based drilling fluid shale oil reservoir, the deep shale oil reservoir to be classified is classified using a gas logging evaluation classification method, a cuttings pyrolysis logging evaluation classification method, and a rock brittleness evaluation classification method, respectively, and all classification results are comprehensively analyzed according to a comprehensive classification standard to obtain a classification result of the deep shale oil reservoir to be classified, wherein the rock brittleness evaluation classification method is to determine an element brittleness index using XRF element logging data, and obtain a reservoir classification result based on the element brittleness index;

[0014] In response to the deep shale oil reservoir to be classified being an oil-based drilling fluid shale oil reservoir, the deep shale oil reservoir to be classified is classified respectively using a gas logging evaluation classification method, a rock brittleness evaluation classification method, and a shale oil reservoir oil content evaluation classification method, and all classification results are comprehensively analyzed according to a comprehensive classification standard to obtain a classification result of the deep shale oil reservoir to be classified, wherein the rock brittleness evaluation classification method is to determine the element brittleness index using XRF element logging data, and obtain the reservoir classification result based on the element brittleness index, and the shale oil reservoir oil content evaluation classification method is to obtain a skeleton resistivity prediction value using a skeleton resistivity prediction model, and obtain the reservoir classification result based on the difference between the skeleton resistivity prediction value and the skeleton resistivity measured value.

[0015] The following are further optimizations and / or improvements to the above technical solutions:

[0016] The above-mentioned gas logging evaluation classification method, rock brittleness evaluation classification method and shale oil reservoir oil content evaluation classification method are used to classify the deep shale oil reservoirs to be classified respectively, and all classification results are comprehensively analyzed according to the comprehensive classification standard to obtain the classification results of the deep shale oil reservoirs to be classified, including:

[0017] Determine the gas logging total hydrocarbon anomaly multiples by using gas logging well logging data, and obtain the first reservoir classification result by combining the first classification standard, wherein the first classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiples threshold intervals and corresponding reservoir types;

[0018] The element brittleness index is determined by using XRF element logging data, and a second reservoir classification result is obtained by combining the second classification standard, wherein the second classification standard is composed of a plurality of element brittleness index threshold intervals and corresponding reservoir types;

[0019] The predicted skeleton resistivity value is obtained by using XRF element logging data, and the corresponding third reservoir classification result is determined based on the difference between the predicted skeleton resistivity value and the measured skeleton resistivity value, wherein the measured skeleton resistivity value is obtained from the resistivity logging data while drilling;

[0020] According to the comprehensive classification standard, the first reservoir classification results, the second reservoir classification results and the third reservoir classification results are comprehensively analyzed to obtain the classification results of the deep shale oil reservoir to be classified.

[0021] The above-mentioned gas logging evaluation classification method, cuttings pyrolysis logging evaluation classification method and rock brittleness evaluation classification method are used to classify the deep shale oil reservoirs to be classified respectively, and all classification results are comprehensively analyzed according to the comprehensive classification standard to obtain the classification results of the deep shale oil reservoirs to be classified, including:

[0022] Determine the gas logging total hydrocarbon anomaly multiples by using gas logging well logging data, and obtain the first reservoir classification result by combining the first classification standard, wherein the first classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiples threshold intervals and corresponding reservoir types;

[0023] The element brittleness index is determined by using XRF element logging data, and a second reservoir classification result is obtained by combining the second classification standard, wherein the second classification standard is composed of a plurality of element brittleness index threshold intervals and corresponding reservoir types;

[0024] The pyrolysis index is determined by using the rock pyrolysis method, and the fourth reservoir classification result is obtained by combining the fourth classification standard, wherein the fourth classification standard is composed of a plurality of pyrolysis index threshold intervals and corresponding reservoir types;

[0025] According to the comprehensive classification standard, the classification results of the first reservoir, the second reservoir and the fourth reservoir are comprehensively analyzed to obtain the classification results of the deep shale oil reservoir to be classified.

[0026] The above-mentioned comprehensive analysis of the first reservoir classification result, the second reservoir classification result and the third reservoir classification result according to the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified is the same as the process of comprehensively analyzing the first reservoir classification result, the second reservoir classification result and the fourth reservoir classification result according to the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified, wherein the comprehensive classification standard comprehensively analyzes the first reservoir classification result, the second reservoir classification result and the third reservoir classification result to obtain the classification result of the deep shale oil reservoir to be classified, including:

[0027] Scoring the first reservoir classification result, the second reservoir classification result and the third reservoir classification result according to the scoring rule, and summing the three scoring results;

[0028] The summation result is brought into the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified, wherein the comprehensive classification standard is composed of multiple summation threshold intervals and corresponding reservoir types.

[0029] The above-mentioned XRF element logging data is used to determine the element brittleness index, and combined with the second classification standard to obtain the second reservoir classification result, including:

[0030] The XRF element logging data was normalized, a mineral analysis sample library was established, and the mineral brittleness index was calculated based on the mineral analysis sample library. The formula for calculating the mineral brittleness index is as follows:

[0031]

[0032] The Pearson correlation analysis method is used to screen out multiple elements in the mineral analysis sample library whose correlation with the mineral brittleness index is greater than the preset value;

[0033] All the screened elements were subjected to multiple linear regression with the mineral brittleness index to obtain the element brittleness index;

[0034] The element brittleness index is introduced into the second classification standard to obtain the second reservoir classification result, wherein the second classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiple threshold intervals and corresponding reservoir types.

[0035] The above-mentioned skeleton resistivity prediction value is obtained by using XRF element logging data, and the corresponding third reservoir classification result is determined based on the difference between the skeleton resistivity prediction value and the skeleton resistivity measured value, including:

[0036] Select multiple sensitive elements in the XRF element logging data, input them into the skeleton resistivity prediction model, and obtain the skeleton resistivity prediction value, wherein the skeleton resistivity prediction model is obtained by training the XGBoost regression model using a number of samples, each of the several samples includes multiple sensitive elements and the identification of the skeleton resistivity measured value, and the multiple sensitive elements are obtained by obtaining the gain value of each element in the XRF element logging data and screening out the elements whose gain value is greater than the threshold value;

[0037] The difference between the predicted skeleton resistivity value and the measured skeleton resistivity value is calculated, and is brought into the third classification standard to determine the corresponding third reservoir classification result, wherein the third classification standard is composed of multiple difference threshold intervals and corresponding reservoir types.

[0038] The above-mentioned method of using gas logging data to determine the total hydrocarbon anomaly multiple of gas logging and combining the first classification standard to obtain the first reservoir classification result includes:

[0039] The total hydrocarbon value is determined using gas logging data and corrected based on drilling environment data;

[0040] The total hydrocarbon abnormality multiple is obtained by using the overall trend line of total hydrocarbon value in the whole well section of the shale oil reservoir and the corrected total hydrocarbon value;

[0041] TGy=TG / y

[0042] Among them, TGy is the total hydrocarbon anomaly multiple; TG is the total hydrocarbon value after correction; y is the overall change trend line of the total hydrocarbon value of the whole well section of the shale oil reservoir;

[0043] The total hydrocarbon anomaly multiple is brought into the first classification standard to determine the corresponding first reservoir classification result, wherein the first classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiple threshold intervals and corresponding reservoir types.

[0044] The above-mentioned rock pyrolysis method is used to determine the pyrolysis index, and combined with the fourth classification standard, the fourth reservoir classification results are obtained, including:

[0045] The pyrolysis index is determined by rock pyrolysis method, and the total pyrolysis amount and light-to-heavy ratio are calculated based on the pyrolysis parameters of the rock cuttings;

[0046] The pyrolysis index is calculated based on the total pyrolysis amount and the light-to-heavy ratio;

[0047] Rzs=S T ×Q

[0048] Where Rzs is the thermal decomposition index; S T is the total amount of pyrolysis; Qzhb is the light-heavy ratio;

[0049] The pyrolysis index is brought into the fourth classification standard to determine the corresponding fourth reservoir classification result, wherein the fourth classification standard is composed of a plurality of pyrolysis index threshold intervals and corresponding reservoir types.

[0050] The second technical solution of the present invention is achieved by the following measures: a deep shale oil reservoir comprehensive classification device, comprising:

[0051] A judgment unit, for judging the drilling fluid type of the deep shale oil reservoir to be classified;

[0052] The water-based classification unit, in response to the deep shale oil reservoir to be classified being a water-based drilling fluid shale oil reservoir, uses a gas logging evaluation classification method, a cuttings pyrolysis logging evaluation classification method, and a rock brittleness evaluation classification method to classify the deep shale oil reservoir to be classified, respectively, and conducts a comprehensive analysis of all classification results according to the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified, wherein the rock brittleness evaluation classification method uses XRF element logging data to determine the element brittleness index, and obtains the reservoir classification result based on the element brittleness index;

[0053] The oil-based classification unit responds to the deep shale oil reservoir to be classified as an oil-based drilling fluid shale oil reservoir, and the gas logging evaluation classification method, the rock brittleness evaluation classification method and the shale oil reservoir oil content evaluation classification method are used to classify the deep shale oil reservoir to be classified respectively, and all classification results are comprehensively analyzed according to the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified, wherein the rock brittleness evaluation classification method is to determine the element brittleness index by using the XRF element logging data, and obtain the reservoir classification result based on the element brittleness index, and the shale oil reservoir oil content evaluation classification method is to obtain the skeleton resistivity prediction value by using the skeleton resistivity prediction model, and obtain the reservoir classification result based on the difference between the skeleton resistivity prediction value and the skeleton resistivity measured value.

[0054] The third technical solution of the present invention is achieved through the following measures: an electronic device, characterized in that it includes a processor and a memory, wherein a computer program is stored in the memory, and the computer program is loaded and executed by the processor to implement the steps in the comprehensive classification of deep shale oil reservoirs.

[0055] The present invention uses different comprehensive classification methods according to the different characteristics of water-based drilling fluid shale oil reservoirs and oil-based drilling fluid shale oil reservoirs, so that the reservoir classification has stronger adaptability and more accurate classification results.

[0056] Among them, in the way of reflecting the oil content of the reservoir, facing the oil-based drilling fluid shale oil reservoir, since the oil-based drilling fluid drilling is easy to make the pyrolysis parameters invalid, it will cause the inaccurate evaluation of the oil content of the shale oil reservoir. In the present invention, the skeleton resistivity prediction value is obtained by using XRF element logging data, and the difference between the skeleton resistivity prediction value and the skeleton resistivity measured value is used to reflect the reservoir oil content, thereby determining the corresponding reservoir classification result, thereby improving the accuracy and adaptability of the overall oil content evaluation of the shale oil reservoir on the basis of not being affected by the pyrolysis parameters. On the contrary, the water-based drilling fluid shale oil reservoir will not have the above problems, so the conventional rock cuttings pyrolysis logging method is used adaptively.

[0057] Among them, in the rock brittleness evaluation and classification method, the commonly used method is to use the relative content of brittle minerals to characterize the rock brittleness, while the present invention directly uses XRF element logging data to determine the element brittleness index, which is faster than the conventional method in analysis speed. It is not affected by the accuracy of XRD spectrum interpretation and experimental analysis, and does not reduce the accuracy of brittleness index calculation, and also effectively reduces the test and detection costs of conventional methods.

[0058] Among them, the present invention corrects the total hydrocarbon value based on the drilling environment data, greatly reduces the influence of drilling fluid density, temperature and bottom hole pressure difference on gas logging parameters, and improves the accuracy of gas logging evaluation and classification. It solves the problem that the gas logging background value is artificially high due to high pressure in deep formations, the fluctuating value during high-angle and directional drilling affects the gas logging accuracy, the drilling fluid density, temperature and bottom hole pressure difference affect the amount of formation oil and gas entering the drilling fluid, resulting in low accuracy of gas logging parameters, and the existing gas logging evaluation and classification methods that do not perform corrections are prone to inaccurate evaluation. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] Attached Figure 1 A schematic flow chart of a comprehensive classification method provided for one embodiment of the present invention.

[0060] Attached Figure 2 A schematic flow chart of a method for classifying shale oil reservoirs using oil-based drilling fluids provided in accordance with an embodiment of the present invention.

[0061] Attached Figure 3 A schematic flow chart of a gas testing evaluation classification method provided in accordance with an embodiment of the present invention.

[0062] Attached Figure 4 A schematic flow chart of a brittleness evaluation classification method provided in one embodiment of the present invention.

[0063] Attached Figure 5 A schematic flow chart of a method for evaluating and classifying oil content provided in accordance with an embodiment of the present invention.

[0064] Attached Figure 6 A schematic diagram of a flow chart of a method for comprehensive analysis of classification results provided by an embodiment of the present invention.

[0065] Attached Figure 7 A schematic diagram of classification results of deep shale oil reservoirs in oil-based drilling provided by an embodiment of the present invention.

[0066] Attached Figure 8 A schematic flow chart of a method for classifying shale oil reservoirs using water-based drilling fluid provided in accordance with an embodiment of the present invention.

[0067] Attached Fig. 9 A schematic flow chart of a pyrolysis index evaluation and classification method provided in accordance with an embodiment of the present invention.

[0068] Attached Fig.10 A schematic diagram of classification results of deep shale oil reservoirs obtained by water-based drilling provided in one embodiment of the present invention.

[0069] Attached Fig.11 A schematic diagram of the structure of a comprehensive classification device provided for one embodiment of the present invention. DETAILED DESCRIPTION

[0070] The present invention is not limited by the following embodiments, and specific implementation methods can be determined based on the technical solution of the present invention and actual conditions.

[0071] Those skilled in the art will understand that, unless otherwise stated, a "module" or "unit" in the embodiments of the present invention refers to a computer program or a part of a computer program with a predetermined function, and works together with other related parts to achieve a predetermined goal, and can be implemented in whole or in part by using software, hardware (such as a processing circuit or memory) or a combination thereof. Similarly, a processor (or multiple processors or memories) can be used to implement one or more modules or units. In addition, each module or unit can be part of an overall module or unit that includes the function of the module or unit.

[0072] In addition, the term “plurality” in the embodiments of the present invention refers to two or more than two, and “first” and “second” are used to distinguish descriptions and should not be understood as implying relative importance.

[0073] The embodiment of the present invention provides a method, device and electronic device for comprehensive classification of deep shale oil reservoirs. The device for comprehensive classification of deep shale oil reservoirs can be integrated in a computer device, which can be a server or a terminal or other device; it can also be executed by a terminal and a server together, and the above examples should not be understood as limiting the present invention.

[0074] The above-mentioned terminal may include a mobile phone, a wearable smart device, a tablet computer, a notebook computer, a personal computer (PC), and a vehicle-mounted computer, etc., and the present invention does not limit this. The present invention does not limit the number of terminal devices.

[0075] The above-mentioned server can be an independent physical server, or a server cluster or distributed system composed of multiple physical servers. It can also be a cloud server that provides basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, CDN, as well as big data and artificial intelligence platforms. The present invention does not impose any restrictions on this.

[0076] For example, the computer equipment first determines the drilling fluid type of the deep shale oil reservoir to be classified; if the deep shale oil reservoir to be classified is a water-based drilling fluid shale oil reservoir, the gas logging evaluation classification method, the cuttings pyrolysis logging evaluation classification method and the rock brittleness evaluation classification method are used to classify the deep shale oil reservoir to be classified respectively, and all classification results are comprehensively analyzed according to the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified, among which the rock brittleness evaluation classification method is to determine the element brittleness index using XRF element logging data, and obtain the reservoir classification result based on the element brittleness index; if the deep shale oil reservoir to be classified is an oil-based drilling fluid shale oil For reservoirs, the gas logging evaluation classification method, rock brittleness evaluation classification method and shale oil reservoir oil content evaluation classification method are used to classify the deep shale oil reservoirs to be classified respectively, and all classification results are comprehensively analyzed according to the comprehensive classification standard to obtain the classification results of the deep shale oil reservoirs to be classified. Among them, the rock brittleness evaluation classification method is to determine the element brittleness index using XRF element logging data, and obtain the reservoir classification result based on the element brittleness index. The shale oil reservoir oil content evaluation classification method is to obtain the skeleton resistivity prediction value using the skeleton resistivity prediction model, and obtain the reservoir classification result based on the difference between the skeleton resistivity prediction value and the skeleton resistivity measured value.

[0077] It should also be noted that the method provided in the embodiment of the present invention may involve artificial intelligence (AI) technology and may be implemented based on artificial intelligence technology, for example, by using machine learning to obtain a corresponding model using sample training.

[0078] Among them, machine learning (ML) is a multi-disciplinary interdisciplinary subject, involving probability theory, statistics, approximation theory, convex analysis, algorithm complexity theory and other disciplines. It specializes in studying how computers simulate or implement human learning behavior to acquire new knowledge or skills and reorganize existing knowledge structures to continuously improve their performance. Machine learning is the core of artificial intelligence and the fundamental way to make computers intelligent. Its applications are in all areas of artificial intelligence.

[0079] The above-mentioned machine learning usually includes technologies such as neural networks, belief networks, reinforcement learning, transfer learning, inductive learning, and self-learning.

[0080] Based on this, the technical solution of the present invention will be introduced and explained in combination with several examples below.

[0081] Embodiment 1: As attached Figure 1 As shown, the embodiment of the present invention discloses a comprehensive classification method for deep shale oil reservoirs, including:

[0082] Step S110, determining the drilling fluid type of the deep shale oil reservoir to be classified;

[0083] Step S120, in response to the deep shale oil reservoir to be classified being a water-based drilling fluid shale oil reservoir, the deep shale oil reservoir to be classified is classified using a gas logging evaluation classification method, a cuttings pyrolysis logging evaluation classification method, and a rock brittleness evaluation classification method, respectively, and all classification results are comprehensively analyzed according to a comprehensive classification standard to obtain a classification result of the deep shale oil reservoir to be classified, wherein the rock brittleness evaluation classification method is to determine an element brittleness index using XRF element logging data, and obtain a reservoir classification result based on the element brittleness index;

[0084] Step S130, in response to the deep shale oil reservoir to be classified being an oil-based drilling fluid shale oil reservoir, the deep shale oil reservoir to be classified is classified using a gas logging evaluation classification method, a rock brittleness evaluation classification method, and a shale oil reservoir oil content evaluation classification method, respectively, and all classification results are comprehensively analyzed according to a comprehensive classification standard to obtain a classification result of the deep shale oil reservoir to be classified, wherein the rock brittleness evaluation classification method is to determine the element brittleness index using XRF element logging data, and obtain a reservoir classification result based on the element brittleness index, and the shale oil reservoir oil content evaluation classification method is to obtain a skeleton resistivity prediction value using a skeleton resistivity prediction model, and obtain a reservoir classification result based on the difference between the skeleton resistivity prediction value and the skeleton resistivity measured value.

[0085] The present invention uses different comprehensive classification methods according to the different characteristics of water-based drilling fluid shale oil reservoirs and oil-based drilling fluid shale oil reservoirs, so that the adaptability of reservoir classification is stronger and the classification results are more accurate. When classifying water-based drilling fluid shale oil reservoirs, a gas logging evaluation classification method, a rock cuttings pyrolysis logging evaluation classification method and a rock brittleness evaluation classification method are introduced to classify the reservoirs respectively, and all classification results are comprehensively analyzed according to the comprehensive classification standard to obtain the classification results of the deep shale oil reservoirs to be classified. When classifying oil-based drilling fluid shale oil reservoirs, a gas logging evaluation classification method, a rock brittleness evaluation classification method and a shale oil reservoir oil content evaluation classification method are introduced to classify the reservoirs respectively, and all classification results are comprehensively analyzed according to the comprehensive classification standard to obtain the classification results of the deep shale oil reservoirs to be classified.

[0086] It should also be noted that in the method of reflecting the oil content of the reservoir, when facing the oil-based drilling fluid shale oil reservoir, the pyrolysis parameters are easily invalidated during drilling with the oil-based drilling fluid, which will cause inaccurate evaluation of the oil content of the shale oil reservoir. In the present invention, the skeleton resistivity prediction value is obtained by using XRF element logging data, and the difference between the skeleton resistivity prediction value and the skeleton resistivity measured value is used to reflect the reservoir oil content, thereby determining the corresponding reservoir classification result, thereby improving the accuracy and adaptability of the overall oil content evaluation of the shale oil reservoir without being affected by the pyrolysis parameters. On the contrary, the water-based drilling fluid shale oil reservoir will not have the above problems, so the conventional rock cuttings pyrolysis logging method is suitable.

[0087] It should also be noted that in the rock brittleness evaluation and classification method, the common method is to use the relative content of brittle minerals to characterize the rock brittleness, while the present invention directly uses XRF element logging data to determine the element brittleness index, which is faster than the conventional method in analysis speed. It is not affected by the accuracy of XRD spectrum interpretation and experimental analysis, and does not reduce the accuracy of brittleness index calculation, and also effectively reduces the test and detection costs of conventional methods.

[0088] Embodiment 2: As attached Figure 2 As shown, the embodiment of the present invention is a further optimization of the above embodiment, wherein the deep shale oil reservoir to be classified is an oil-based drilling fluid shale oil reservoir, and the deep shale oil reservoir to be classified is classified using the gas logging evaluation classification method, the rock brittleness evaluation classification method and the shale oil reservoir oil content evaluation classification method, respectively, and all classification results are comprehensively analyzed according to the comprehensive classification standard to obtain the classification results of the deep shale oil reservoir to be classified, including:

[0089] Step S210, using gas logging data to determine the gas logging total hydrocarbon anomaly multiple, and combining the first classification standard to obtain a first reservoir classification result, wherein the first classification standard is composed of multiple gas logging total hydrocarbon anomaly multiple threshold intervals and corresponding reservoir types.

[0090] In this embodiment, as shown in the attached Figure 3 As shown, step S210 specifically includes:

[0091] Step S211, using gas logging data to determine the total hydrocarbon value, and correcting the total hydrocarbon value based on drilling environment data.

[0092] In this step, the gas logging data is the gas logging data of the deep shale oil reservoir to be classified, including gas logging components C1, C2, C3, iC4, nC4, iC5, and nC5;

[0093] Based on the gas logging data of the above-mentioned deep shale oil reservoir to be classified, the total hydrocarbon value TG is calculated by the following formula:

[0094] TG=C1+C2+C3+iC4+nC4+iC5+nC5, unit: %

[0095] On the basis of obtaining the total hydrocarbon value, the total hydrocarbon value is corrected in combination with the drilling environment data, as described below:

[0096]

[0097] in, , The drilling fluid density collected in real time during logging, in g / cm 3 , is the gravitational acceleration constant, take 9.81m / s 2 , is the well depth collected in real time by logging, in meters; ZS is the drilling time required to drill into 1 meter of the formation, collected in real time by logging, in minutes / meters; TG is the total hydrocarbon value; is the bottom formation temperature, according to the well depth data of adjacent wells Temperature and regional geothermal gradient calculate, ; The outlet drilling fluid temperature collected in real time during logging, in °C; is the atmospheric pressure on the ground, which is 0.1Mpa.

[0098] Step S212, using the overall change trend line of the total hydrocarbon value of the entire well section of the shale oil reservoir and the corrected total hydrocarbon value, obtain the total hydrocarbon abnormality multiple;

[0099] TGy=TG / y

[0100] Among them, TGy is the total hydrocarbon anomaly multiple; TG is the total hydrocarbon value after correction; y is the overall change trend line of the total hydrocarbon value of the entire well section of the shale oil reservoir.

[0101] In this step, the overall change trend line of the total hydrocarbon value of the whole well section of the shale oil reservoir is y=a×h+b, where a and b are obtained by linear fitting using the total hydrocarbon anomaly multiple TGx of the whole well section of the shale oil reservoir and the well depth h.

[0102] Step S213, the total hydrocarbon anomaly multiple is brought into the first classification standard to determine the corresponding first reservoir classification result, wherein the first classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiple threshold intervals and corresponding reservoir types.

[0103] In this step, the first classification standard is set accordingly based on the analysis and summary of the historical data, and may be, but is not limited to, as follows:

[0104] (1) If the total hydrocarbon anomaly multiple is greater than 0.6, it is a Class I reservoir;

[0105] (2) If ≤0.2 total hydrocarbon anomaly multiple ≤0.6, it is a Class II reservoir;

[0106] (3) If the total hydrocarbon anomaly multiple is less than 0.2, it is a Class III reservoir.

[0107] Due to the high pressure in deep formations, the gas logging background value is artificially high. The fluctuating values ​​during high-angle and directional drilling affect the gas logging accuracy. The drilling fluid density, temperature and bottom hole pressure difference affect the amount of formation oil and gas entering the drilling fluid, resulting in low accuracy of gas logging parameters. Therefore, the existing gas logging evaluation method without correction is prone to inaccurate evaluation, and thus inaccurate reservoir classification. Therefore, the embodiment of the present invention corrects the total hydrocarbon value based on the drilling environment data, greatly reduces the influence of drilling fluid density, temperature and bottom hole pressure difference on gas logging parameters, and improves the accuracy of gas logging evaluation classification.

[0108] Step S220, using XRF element logging data to determine the element brittleness index, and combining it with a second classification standard to obtain a second reservoir classification result, wherein the second classification standard consists of a plurality of element brittleness index threshold intervals and corresponding reservoir types.

[0109] In this embodiment, as shown in the attached Figure 4 As shown, step S220 specifically includes:

[0110] Step S221, normalize the XRF element logging data, establish a mineral analysis sample library, and obtain the mineral brittleness index based on the mineral analysis sample library, wherein the formula for obtaining the mineral brittleness index is as follows:

[0111]

[0112] Step S222, using the Pearson correlation analysis method to screen out multiple elements in the mineral analysis sample library whose correlation with the mineral brittleness index is greater than a preset value;

[0113] Step S223, performing multiple linear regression on all the screened elements and the mineral brittleness index to obtain the element brittleness index;

[0114] In this step, the element brittleness index constructed by multivariate linear regression of all the screened elements and the mineral brittleness index is as follows:

[0115]

[0116] in, is the element brittleness index; is a constant, , is the regression coefficient; , In order to use the Pearson correlation analysis method to screen out multiple elements in the mineral analysis sample library whose correlation with the mineral brittleness index is greater than the preset value.

[0117] Step S224, the element brittleness index is brought into the second classification standard to obtain a second reservoir classification result, wherein the second classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiple threshold intervals and corresponding reservoir types.

[0118] In this step, the second classification standard is set accordingly based on the analysis and summary of the historical data, and may be, but is not limited to, as follows:

[0119] (1) If the element brittleness index is greater than 3.0, it is a Class I reservoir;

[0120] (2) If 2.3≦element brittleness index≦3.0, it is a Class II reservoir;

[0121] (3) If the element brittleness index is less than 2.3, it is a Class III reservoir.

[0122] Compared with the commonly used rock brittleness evaluation and classification methods, the above method no longer uses the relative content of brittle minerals to characterize rock brittleness, but directly uses XRF element logging data to determine the element brittleness index. The analysis speed is fast, and it is not affected by the accuracy of XRD spectrum interpretation and experimental analysis. It does not reduce the accuracy of brittleness index calculation, and effectively reduces the test and detection costs of conventional methods.

[0123] Step S230, using XRF element logging data to obtain a predicted skeleton resistivity value, and determining a corresponding third reservoir classification result based on the difference between the predicted skeleton resistivity value and the measured skeleton resistivity value, wherein the measured skeleton resistivity value is obtained from the resistivity logging data while drilling.

[0124] In this embodiment, as shown in the attached Figure 5 As shown, step S230 specifically includes:

[0125] Step S231, select multiple sensitive elements in the XRF element logging data, input them into the skeleton resistivity prediction model, and obtain the skeleton resistivity prediction value, wherein the skeleton resistivity prediction model is obtained by training the XGBoost regression model using a number of samples, each of the several samples includes multiple sensitive elements and identifications of the actual measured values ​​of the skeleton resistivity, and the multiple sensitive elements are obtained by obtaining the gain value of each element in the XRF element logging data, and screening out the elements whose gain value is greater than the threshold.

[0126] In this step, the construction process of the skeleton resistivity prediction model includes:

[0127] Obtain historical XRF element logging data of deep shale oil reservoirs from historical oil-based drilling and preprocess them, including processing missing values, outliers, normalization, etc.

[0128] Use decision tree and other algorithms to determine the gain values ​​of each element in the historical XRF element logging data, and select elements with gain values ​​greater than the threshold (i.e., the greatest importance) as sensitive elements;

[0129] Extract all sensitive elements from historical XRF element logging data to form a sample set, and divide the sample set into a training set and a test set in proportion;

[0130] The XGBoost regression model is trained using the training set. The loss function is introduced during the training. When the value of the loss function is stable, the training is terminated to obtain the skeleton resistivity prediction model.

[0131] The trained skeleton resistivity prediction model is tested and evaluated using the test set. Mean square error (MSE), root mean square error (RMSE), mean absolute error (MAE), and mean absolute percentage error (MAPE) are introduced as evaluation indicators to optimize the model parameters of the skeleton resistivity prediction model and output a skeleton resistivity prediction model that meets the test and evaluation requirements.

[0132] Furthermore, the code corresponding to the construction process of the skeleton resistivity prediction model may be, but is not limited to:

[0133] # Import basic tool library

[0134] import numpy as np

[0135] import pandas as pd

[0136] from sklearn.datasets import make_regression

[0137] from sklearn.model_selection import train_test_split

[0138] from xgboost.sklearn import XGBRegressor

[0139] from sklearn.metrics import mean_squared_error

[0140] from sklearn.metrics import mean_absolute_error

[0141] from sklearn.metrics import mean_absolute_percentage_error

[0142] # Import element data for data analysis

[0143] data = pd.read_table('.data.csv',delimiter=',')

[0144] x = data.iloc[:,:-1]

[0145] y = data.iloc[:,-1]

[0146] X, y = make_regression(n_samples=470, n_features=35)

[0147] # Split into training set and test set

[0148] train_X, test_X, train_y, test_y = train_test_split(X, y, test_size=0.3, shuffle=False)

[0149] #Call the XGBoost model and use the training set data for training (fitting)

[0150] my_model = XGBRegressor(

[0151] max_depth=8,

[0152] learning_rate=0.1,

[0153] n_estimators=100,

[0154] silent=True,

[0155] objective='reg: squarederror',

[0156] booster='gbtree',

[0157] n_jobs=15,

[0158] nthread=None,

[0159] gamma=0,

[0160] min_child_weight=1,

[0161] max_delta_step=0,

[0162] subsample=1,

[0163] colsample_bytree=1,

[0164] colsample_bylevel=1,

[0165] reg_alpha=0,

[0166] reg_lambda=1,

[0167] scale_pos_weight=1,

[0168] base_score=0.5,

[0169] random_state=0,

[0170] seed=None,

[0171] missing=None,

[0172] importance_type = 'gain')

[0173] my_model.fit(train_X, train_y)

[0174] # Use the model to predict the test set data

[0175] predictions = my_model.predict(test_X)

[0176] # Evaluate the training and testing results

[0177] print("mean_squared_error: " + str(mean_squared_error(predictions,test_y)))

[0178] print("rmse: "+ str(np.sqrt(mean_squared_error(predictions, test_y)))

[0179] print("mean_absolute_error: "+ str(mean_absolute_error(predictions,test_y)))

[0180] print("mean_absolute_percentage_error: "+str(mean_absolute_percentage_error (predictions, test_y)))

[0181] Step S232, calculate the difference between the predicted skeleton resistivity value and the measured skeleton resistivity value, and bring it into the third classification standard to determine the corresponding third reservoir classification result, wherein the third classification standard is composed of multiple difference threshold intervals and corresponding reservoir types.

[0182] In this step, the measured value of the skeleton resistivity is the measured shale oil resistivity, which is the sum of the series resistivity of the shale rock mineral skeleton resistivity and the oil and gas resistivity in the pores. The shale rock mineral skeleton resistivity and the oil and gas resistivity in the pores can both be obtained from the XRF element logging data, and the difference between the predicted value of the skeleton resistivity and the measured value of the skeleton resistivity in this step can truly reflect the oil content of the reservoir.

[0183] In this step, the third classification standard is set accordingly based on the analysis and summary of the historical data, and may be, but is not limited to, as follows:

[0184] (1) If the difference is greater than 200, it is a Class I reservoir;

[0185] (2) If 100 ≤ difference ≤ 200, it is a Class II reservoir;

[0186] (3) If the difference is less than 100, it is a Class III reservoir.

[0187] This embodiment no longer uses the conventional rock pyrolysis method to complete the shale oil reservoir oil content evaluation and classification method, but introduces machine learning to construct a skeleton resistivity prediction model, and uses the skeleton resistivity prediction model to predict based on XRF element logging data to obtain a skeleton resistivity prediction value, and then reflects the reservoir oil content based on the difference between the skeleton resistivity prediction value and the skeleton resistivity measured value, thereby determining the corresponding reservoir classification result, thereby improving the accuracy and adaptability of the overall oil content evaluation of the shale oil reservoir without being affected by the pyrolysis parameters. In addition, the skeleton resistivity prediction model can be continuously learned and updated, dynamically adapting to the oil-based deep shale oil reservoirs in various regions, ensuring the accuracy of the prediction, thereby improving the accuracy of the reservoir classification.

[0188] Step S240, comprehensively analyzing the first reservoir classification result, the second reservoir classification result and the third reservoir classification result according to the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified.

[0189] In this embodiment, as shown in the attached Figure 6 As shown, step S240 specifically includes:

[0190] Step S241, scoring the first reservoir classification result, the second reservoir classification result and the third reservoir classification result according to the scoring rule, and summing the three scoring results.

[0191] In this step, the scoring rules can be constructed by selecting an adaptive method as needed, which may include but is not limited to determining the weights of the three categories, directly setting the scores according to the reservoirs, etc. In this embodiment, the scores of different categories of reservoirs can be directly set. For example, the scoring rule is that the first category of reservoirs gets 1 point, the second category of reservoirs gets 2 points, and the third category of reservoirs gets 3 points.

[0192] Step S242, bringing the summation result into the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified, wherein the comprehensive classification standard is composed of multiple summation threshold intervals and corresponding reservoir types.

[0193] The comprehensive classification standard in this step is composed of multiple summation threshold intervals and corresponding reservoir types, wherein the summation threshold intervals are set accordingly based on the analysis and summary of historical data, and can be, but not limited to, as follows:

[0194] (1) If the sum is ≤5, the comprehensive classification is one category;

[0195] (2) If 6≤sum result≤7, and the sum of any two of the scores of the first reservoir classification result, the second reservoir classification result, and the third reservoir classification result is less than 6, then the comprehensive classification is Class II;

[0196] (3) Other comprehensive classifications are divided into three categories.

[0197] Example 3: Taking the deep shale oil reservoir of 5430-5900m oil-based drilling in H7 well as an example, the classification method of deep shale oil reservoir disclosed in the present invention is used as an example. The specific classification process is as follows:

[0198] (1) Calculation of gas logging anomaly multiples and classification of gas logging evaluation

[0199] S1, calculate the total hydrocarbon value of the section 5430-5900m of Well H7, TG=C1+C2+C3+iC4+nC4+iC5+nC5;

[0200] S2, correcting the total hydrocarbon value based on drilling environment data;

[0201] TGx = drilling time × total hydrocarbon calculation / (3.14 × (0.216 × 0.216) / 4) × ((90 / drilling fluid temperature) × (0.1 / drilling fluid density × 0.00981 × well depth)) × 0.9 / 10000;

[0202] S3, calculate the overall trend line of total hydrocarbon value in the whole well section of shale oil reservoir, y = -4E-05*well depth + 1.3107;

[0203] S4, calculate the total hydrocarbon anomaly multiple TGy=TG / y;

[0204] S5. Analyze TGy according to the classification standard in Table 1 to obtain the first reservoir classification result.

[0205] (2) Calculation of element brittleness index and classification of rock brittleness evaluation

[0206] S1, collect the XRF element logging data of the cuttings in the H7 well section from 5430 to 5900 m, with a logging collection interval of 1 m, a total of 470 data samples (including element data, Bid data, and Bif data), and collect 35 elements including Si, Fe, Ca, and Mg. The obtained element data are normalized to establish an element analysis sample library. The normalization formula is as follows:

[0207]

[0208] Wherein, E is the measured value of the element; Emax is the maximum value of the element; Emin is the minimum value of the element; E* is the normalized value of the element;

[0209] S3, calculate the mineral brittleness index BId of 5430-5900m in H7 well, and use the Pearson correlation analysis method to screen out the elements with high correlation with the mineral brittleness index BId, which are elements Na, Al, Si, S, K, Ca, Fe, Ni, and Rb. Among them, the calculation formula of the mineral brittleness index BId is as follows:

[0210]

[0211] S4, perform multiple linear regression between the elements Na, Al, Si, S, K, Ca, Fe, Ni, Rb and the mineral brittleness index BId (element data is the independent variable, BId is the dependent variable) to obtain the element brittleness index, as follows:

[0212] BIf=-1.172+1.127×Na+1.67×Al+4.185×Si+0.628×S-1.025×K+4.257×Ca-1.159×Fe+0.422×Ni+0.251×Rb;

[0213] S5, analyzing the element brittleness index according to the classification standard in Table 1 to obtain the second reservoir classification result.

[0214] (3) Calculation of rock skeleton resistivity and classification of oil content in shale oil reservoirs

[0215] S1, based on the gain value, the seven most important features are determined as sensitive elements, namely S, Ca, Si, Al, K, Mg, and Na;

[0216] S2, using machine learning to build a skeleton resistivity prediction model;

[0217] S3, input the sensitive elements S, Ca, Si, Al, K, Mg, and Na of the section 5430-5900 m of the H7 well into the skeleton resistivity prediction model to obtain the skeleton resistivity prediction value;

[0218] S4, calculating the difference between the predicted skeleton resistivity value and the measured skeleton resistivity value, analyzing the difference according to the classification standard in Table 1, and obtaining the third reservoir classification result.

[0219] (4) According to the comprehensive classification standard, the first reservoir classification results, the second reservoir classification results and the third reservoir classification results are comprehensively analyzed to obtain the classification results of the deep shale oil reservoir to be classified. The classification results are shown in the attached figure. Figure 7 shown.

[0220] Table 1 Classification standard of deep shale oil reservoirs for oil-based drilling

[0221] .

[0222] Embodiment 4: As attached Figure 8 As shown, the embodiment of the present invention is a further optimization of the above embodiment, wherein the deep shale oil reservoirs to be classified are classified respectively by using the gas logging evaluation classification method, the cuttings pyrolysis logging evaluation classification method and the rock brittleness evaluation classification method, and all the classification results are comprehensively analyzed according to the comprehensive classification standard to obtain the classification results of the deep shale oil reservoirs to be classified, including:

[0223] Step S310, using gas logging data to determine the gas logging total hydrocarbon anomaly multiple, and combining the first classification standard to obtain a first reservoir classification result, wherein the first classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiple threshold intervals and corresponding reservoir types;

[0224] Step S320, determining the element brittleness index using the XRF element logging data, and obtaining a second reservoir classification result in combination with a second classification standard, wherein the second classification standard is composed of a plurality of element brittleness index threshold intervals and corresponding reservoir types;

[0225] Step S330, using the rock pyrolysis method to determine the pyrolysis index, and combining it with the fourth classification standard to obtain a fourth reservoir classification result, wherein the fourth classification standard is composed of a plurality of pyrolysis index threshold intervals and corresponding reservoir types;

[0226] Step S340, comprehensively analyzing the first reservoir classification result, the second reservoir classification result and the fourth reservoir classification result according to the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified.

[0227] In this embodiment, the specific implementation process of step S310, step S320 and step S340 is the same as that of embodiment 2 and will not be repeated herein.

[0228] In this embodiment, as shown in the attached Fig. 9 As shown, step S330 specifically includes:

[0229] Step S331, obtaining the pyrolysis parameters of the rock cuttings, and calculating the total pyrolysis amount and the light-to-heavy ratio based on the pyrolysis parameters of the rock cuttings; for example:

[0230] The rock cuttings pyrolysis parameters include S 0 , S 11 , S 21 , S 22 , S 23 , unit: mg / g;

[0231] Movable oil S 1 =S 11 +S 21 ×0.67;

[0232] Non-movable oil S 2 =S 21 ×0.33+S 22 +S 23 ;

[0233] Total pyrolysis amount S T =S 0 +S 11 +S 21 +S 22 +S 23 ;

[0234] Light-to-heavy ratio Qzhb=S 1 / S 2 .

[0235] Step S332, calculating the pyrolysis index based on the total pyrolysis amount and the light-to-heavy ratio;

[0236] Rzs=S T ×Q

[0237] Where Rzs is the thermal decomposition index; S Tis the total amount of pyrolysis; Qzhb is the light-heavy ratio;

[0238] Step S333, bringing the pyrolysis index into the fourth classification standard to determine the corresponding fourth reservoir classification result, wherein the fourth classification standard is composed of a plurality of pyrolysis index threshold intervals and corresponding reservoir types.

[0239] Example 5: Taking the classification of the deep shale oil reservoirs in the H7 well at 5430-5900 m water-based drilling as an example using the comprehensive classification of deep shale oil reservoirs disclosed in the present invention, the specific classification process is as follows:

[0240] (1) Calculation of gas logging anomaly multiples and classification of gas logging evaluation

[0241] S1, calculate the total hydrocarbon value of the section 5430-5900m of Well H7, TG=C1+C2+C3+iC4+nC4+iC5+nC5;

[0242] S2, correcting the total hydrocarbon value based on drilling environment data;

[0243] TGx = drilling time × total hydrocarbon calculation / (3.14 × (0.216 × 0.216) / 4) × ((90 / drilling fluid temperature) × (0.1 / drilling fluid density × 0.00981 × well depth) × 0.9) / 10000;

[0244] S3, calculate the overall trend line of total hydrocarbon value in the whole well section of shale oil reservoir, y = -4E-05*well depth + 1.3107;

[0245] S4, calculate the total hydrocarbon anomaly multiple TGy=TG / y;

[0246] S5, analyze TGy according to the classification standard in Table 2 to obtain the first reservoir classification result.

[0247] (2) Calculation of element brittleness index and classification of rock brittleness evaluation

[0248] S1, collect the XRF element logging data of the cuttings in the H7 well section from 5430 to 5900 m, with a logging collection interval of 1 m, a total of 470 data samples (including element data, Bid data, and Bif data), and collect 35 elements including Si, Fe, Ca, and Mg. The obtained element data are normalized to establish an element analysis sample library. The normalization formula is as follows:

[0249]

[0250] Wherein, E is the measured value of the element; Emax is the maximum value of the element; Emin is the minimum value of the element; E* is the normalized value of the element;

[0251] S3, calculate the mineral brittleness index BId of 5430-5900m in H7 well, and use the Pearson correlation analysis method to screen out the elements with high correlation with the mineral brittleness index BId, which are elements Na, Al, Si, S, K, Ca, Fe, Ni, and Rb. Among them, the calculation formula of the mineral brittleness index BId is as follows:

[0252]

[0253] S4, perform multiple linear regression between the elements Na, Al, Si, S, K, Ca, Fe, Ni, Rb and the mineral brittleness index BId (element data is the independent variable, BId is the dependent variable) to obtain the element brittleness index, as follows:

[0254] BIf=-1.172+1.127×Na+1.67×Al+4.185×Si+0.628×S-1.025×K+4.257×Ca-1.159×Fe+0.422×Ni+0.251×Rb;

[0255] S5, analyze the element brittleness index according to the classification standard in Table 2 to obtain the second reservoir classification result.

[0256] (3) Pyrolysis index calculation and classification of cuttings pyrolysis logging evaluation

[0257] S1, obtained cuttings pyrolysis parameters, S0, S11, S21, S22, S23, unit mg / g;

[0258] S2, movable oil S1 = S11 + S21 × 0.67;

[0259] S3, non-movable oil S2 = S21 × 0.33 + S22 + S23;

[0260] S4, total pyrolysis amount ST = S0 + S11 + S21 + S22 + S23;

[0261] S5, weight ratio Qzhb=S1 / S2.

[0262] S6, calculating the pyrolysis index based on the total pyrolysis amount and the light-to-heavy ratio;

[0263] Rzs=S T ×Q

[0264] Where Rzs is the thermal decomposition index; S T is the total amount of pyrolysis; Qzhb is the light-to-heavy ratio.

[0265] (4) According to the comprehensive classification standard, the first reservoir classification results, the second reservoir classification results and the fourth reservoir classification results are comprehensively analyzed to obtain the classification results of the deep shale oil reservoir to be classified. The classification results are shown in the attached figure. Fig.10 shown.

[0266] Table 2 Classification standard of deep shale oil reservoirs for water-based drilling

[0267] .

[0268] Embodiment 6: As attached Fig.11 As shown, the embodiment of the present invention discloses a deep shale oil reservoir comprehensive classification device, comprising:

[0269] A judgment unit, for judging the drilling fluid type of the deep shale oil reservoir to be classified;

[0270] The water-based classification unit, in response to the deep shale oil reservoir to be classified being a water-based drilling fluid shale oil reservoir, uses a gas logging evaluation classification method, a cuttings pyrolysis logging evaluation classification method, and a rock brittleness evaluation classification method to classify the deep shale oil reservoir to be classified, respectively, and conducts a comprehensive analysis of all classification results according to the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified, wherein the rock brittleness evaluation classification method uses XRF element logging data to determine the element brittleness index, and obtains the reservoir classification result based on the element brittleness index;

[0271] The oil-based classification unit responds to the deep shale oil reservoir to be classified as an oil-based drilling fluid shale oil reservoir, and the gas logging evaluation classification method, the rock brittleness evaluation classification method and the shale oil reservoir oil content evaluation classification method are used to classify the deep shale oil reservoir to be classified respectively, and all classification results are comprehensively analyzed according to the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified, wherein the rock brittleness evaluation classification method is to determine the element brittleness index by using the XRF element logging data, and obtain the reservoir classification result based on the element brittleness index, and the shale oil reservoir oil content evaluation classification method is to obtain the skeleton resistivity prediction value by using the skeleton resistivity prediction model, and obtain the reservoir classification result based on the difference between the skeleton resistivity prediction value and the skeleton resistivity measured value.

[0272] The specific implementation steps of each unit in this embodiment are as described in Embodiments 2 to 5, and will not be repeated here.

[0273] Example 7: An embodiment of the present invention discloses an electronic device, including a processor and a memory, wherein the memory stores a computer program, and the computer program is loaded and executed by the processor to implement a comprehensive classification method for deep shale oil reservoirs.

[0274] The processor may be a central processing unit (CPU), a general purpose processor, a digital signal processor (DSP), an ASIC, an FPGA or other programmable logic device, a transistor logic device, a hardware component or any combination thereof. It may implement or execute various exemplary logic blocks, modules and circuits described in conjunction with the disclosure of the present invention. It may also be a combination that implements computing functions, such as a combination of one or more microprocessors, a combination of a DSP and a microprocessor, and the like. The memory may include, but is not limited to, various media that can store computer programs, such as a USB flash drive, a read-only memory, a mobile hard disk, a magnetic disk or an optical disk.

[0275] It should be understood by those skilled in the art that the embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes. The solutions in the embodiments of the present invention may be implemented in various computer languages, for example, object-oriented programming language Java and interpreted scripting language JavaScript, etc.

[0276] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0277] These computer program instructions may also be stored in a computer readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0278] The above contents are only specific implementation methods of the present invention, which have strong adaptability and implementation effect, but the protection scope of the present invention is not limited to this. Any technician familiar with the technical field can easily think of changes or substitutions within the technical scope disclosed by the present invention, which should be covered within the protection scope of the present invention. Therefore, equivalent changes made according to the claims of the present invention are still within the scope covered by the present invention.

Claims

1. A comprehensive classification method for deep shale oil reservoirs, characterized in that: include: Determine the drilling fluid type of the deep shale oil reservoir to be classified; In response to the deep shale oil reservoir to be classified being a water-based drilling fluid shale oil reservoir, the deep shale oil reservoir to be classified is classified using the gas logging evaluation classification method, the cuttings pyrolysis logging evaluation classification method and the rock brittleness evaluation classification method, respectively, and all classification results are comprehensively analyzed according to the comprehensive classification standard to obtain the classification results of the deep shale oil reservoir to be classified, including: Determine the gas logging total hydrocarbon anomaly multiples by using gas logging well logging data, and obtain the first reservoir classification result by combining the first classification standard, wherein the first classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiples threshold intervals and corresponding reservoir types; The element brittleness index is determined by using XRF element logging data, and a second reservoir classification result is obtained by combining the second classification standard, wherein the second classification standard is composed of a plurality of element brittleness index threshold intervals and corresponding reservoir types; The pyrolysis index is determined by using the rock pyrolysis method, and the fourth reservoir classification result is obtained by combining the fourth classification standard, wherein the fourth classification standard is composed of a plurality of pyrolysis index threshold intervals and corresponding reservoir types; Comprehensively analyzing the first reservoir classification results, the second reservoir classification results and the fourth reservoir classification results according to the comprehensive classification standard, and obtaining the classification results of the deep shale oil reservoir to be classified; In response to the deep shale oil reservoir to be classified being an oil-based drilling fluid shale oil reservoir, the deep shale oil reservoir to be classified is classified using the gas logging evaluation classification method, the rock brittleness evaluation classification method and the shale oil reservoir oil content evaluation classification method, respectively, and all classification results are comprehensively analyzed according to the comprehensive classification standard to obtain the classification results of the deep shale oil reservoir to be classified, including: Determine the gas logging total hydrocarbon anomaly multiples by using gas logging well logging data, and obtain the first reservoir classification result by combining the first classification standard, wherein the first classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiples threshold intervals and corresponding reservoir types; The element brittleness index is determined by using XRF element logging data, and a second reservoir classification result is obtained by combining the second classification standard, wherein the second classification standard is composed of a plurality of element brittleness index threshold intervals and corresponding reservoir types; The predicted skeleton resistivity value is obtained by using XRF element logging data, and the corresponding third reservoir classification result is determined based on the difference between the predicted skeleton resistivity value and the measured skeleton resistivity value, wherein the measured skeleton resistivity value is obtained from the resistivity logging data while drilling; Comprehensively analyzing the first reservoir classification results, the second reservoir classification results, and the third reservoir classification results according to the comprehensive classification standard, and obtaining the classification results of the deep shale oil reservoir to be classified; The gas logging data is used to determine the gas logging total hydrocarbon anomaly multiple, and the first reservoir classification result is obtained in combination with the first classification standard, wherein the first classification standard is composed of multiple gas logging total hydrocarbon anomaly multiple threshold intervals and corresponding reservoir types, including: The total hydrocarbon value is determined using gas logging data and corrected based on drilling environment data; in, , The drilling fluid density collected in real time for logging. is the gravitational acceleration constant, is the well depth collected in real time by logging; ZS is the drilling time required to drill into 1m of the formation; TG is the total hydrocarbon value; is the bottom formation temperature, according to the well depth data of adjacent wells Temperature and regional geothermal gradient calculate, ; The outlet drilling fluid temperature collected in real time for logging; is the atmospheric pressure on the ground; The total hydrocarbon abnormality multiple is obtained by using the overall trend line of total hydrocarbon value in the whole well section of the shale oil reservoir and the corrected total hydrocarbon value; Among them, TGy is the total hydrocarbon anomaly multiple; TGx is the corrected total hydrocarbon value; y is the overall change trend line of the total hydrocarbon value of the entire well section of the shale oil reservoir; The total hydrocarbon anomaly multiple is brought into the first classification standard to determine the corresponding first reservoir classification result, wherein the first classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiple threshold intervals and corresponding reservoir types; Among them, the skeleton resistivity prediction value is obtained by using XRF element logging data, and the corresponding third reservoir classification result is determined based on the difference between the skeleton resistivity prediction value and the skeleton resistivity measured value, including: Select multiple sensitive elements in the XRF element logging data, input them into the skeleton resistivity prediction model, and obtain the skeleton resistivity prediction value, wherein the skeleton resistivity prediction model is obtained by training the XGBoost regression model using a number of samples, each of the several samples includes multiple sensitive elements and the identification of the skeleton resistivity measured value, and the multiple sensitive elements are obtained by obtaining the gain value of each element in the XRF element logging data and screening out the elements whose gain value is greater than the threshold value; The difference between the predicted skeleton resistivity value and the measured skeleton resistivity value is calculated, and is brought into the third classification standard to determine the corresponding third reservoir classification result, wherein the third classification standard is composed of multiple difference threshold intervals and corresponding reservoir types.

2. The comprehensive classification method for deep shale oil reservoirs according to claim 1, characterized in that: The process of comprehensively analyzing the first reservoir classification result, the second reservoir classification result and the third reservoir classification result according to the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified is the same as the process of comprehensively analyzing the first reservoir classification result, the second reservoir classification result and the fourth reservoir classification result according to the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified, wherein the comprehensive classification standard comprehensively analyzes the first reservoir classification result, the second reservoir classification result and the third reservoir classification result to obtain the classification result of the deep shale oil reservoir to be classified, including: Scoring the first reservoir classification result, the second reservoir classification result and the third reservoir classification result according to the scoring rule, and summing the three scoring results; The summation result is brought into the comprehensive classification standard to obtain the classification result of the deep shale oil reservoir to be classified, wherein the comprehensive classification standard is composed of multiple summation threshold intervals and corresponding reservoir types.

3. The comprehensive classification method for deep shale oil reservoirs according to claim 1 or 2, characterized in that: The method of determining the element brittleness index by using the XRF element logging data and obtaining the second reservoir classification result in combination with the second classification standard includes: The XRF element logging data was normalized, a mineral analysis sample library was established, and the mineral brittleness index was calculated based on the mineral analysis sample library. The formula for calculating the mineral brittleness index is as follows: The Pearson correlation analysis method is used to screen out multiple elements in the mineral analysis sample library whose correlation with the mineral brittleness index is greater than the preset value; All the screened elements were subjected to multiple linear regression with the mineral brittleness index to obtain the element brittleness index; The element brittleness index is introduced into the second classification standard to obtain the second reservoir classification result, wherein the second classification standard is composed of a plurality of element brittleness index threshold intervals and corresponding reservoir types.

4. The comprehensive classification method for deep shale oil reservoirs according to claim 1 or 2, characterized in that: The method of determining the pyrolysis index by using the rock pyrolysis method and obtaining the fourth reservoir classification result in combination with the fourth classification standard includes: Obtaining pyrolysis parameters of rock cuttings, and calculating the total pyrolysis amount and light-to-heavy ratio based on the pyrolysis parameters of rock cuttings; The pyrolysis index is calculated based on the total pyrolysis amount and the light-to-heavy ratio; Where Rzs is the thermal decomposition index; S T is the total amount of pyrolysis; Qzhb is the light-heavy ratio; The pyrolysis index is brought into the fourth classification standard to determine the corresponding fourth reservoir classification result, wherein the fourth classification standard is composed of a plurality of pyrolysis index threshold intervals and corresponding reservoir types.

5. A deep shale oil reservoir comprehensive classification device using the method as described in any one of claims 1 to 4, characterized in that: include: A judgment unit, for judging the drilling fluid type of the deep shale oil reservoir to be classified; The water-based classification unit, in response to the deep shale oil reservoir to be classified being a water-based drilling fluid shale oil reservoir, uses the gas logging evaluation classification method, the cuttings pyrolysis logging evaluation classification method and the rock brittleness evaluation classification method to classify the deep shale oil reservoir to be classified respectively, and comprehensively analyzes all classification results according to the comprehensive classification standard to obtain the classification results of the deep shale oil reservoir to be classified, including: Determine the gas logging total hydrocarbon anomaly multiples by using gas logging well logging data, and obtain the first reservoir classification result by combining the first classification standard, wherein the first classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiples threshold intervals and corresponding reservoir types; The element brittleness index is determined by using XRF element logging data, and a second reservoir classification result is obtained by combining the second classification standard, wherein the second classification standard is composed of a plurality of element brittleness index threshold intervals and corresponding reservoir types; The pyrolysis index is determined by using the rock pyrolysis method, and the fourth reservoir classification result is obtained by combining the fourth classification standard, wherein the fourth classification standard is composed of a plurality of pyrolysis index threshold intervals and corresponding reservoir types; Comprehensively analyzing the first reservoir classification results, the second reservoir classification results and the fourth reservoir classification results according to the comprehensive classification standard, and obtaining the classification results of the deep shale oil reservoir to be classified; The oil-based classification unit, in response to the deep shale oil reservoir to be classified being an oil-based drilling fluid shale oil reservoir, uses the gas logging evaluation classification method, the rock brittleness evaluation classification method, and the shale oil reservoir oil content evaluation classification method to classify the deep shale oil reservoir to be classified respectively, and comprehensively analyzes all classification results according to the comprehensive classification standard to obtain the classification results of the deep shale oil reservoir to be classified, including: Determine the gas logging total hydrocarbon anomaly multiples by using gas logging well logging data, and obtain the first reservoir classification result by combining the first classification standard, wherein the first classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiples threshold intervals and corresponding reservoir types; The element brittleness index is determined by using XRF element logging data, and a second reservoir classification result is obtained by combining the second classification standard, wherein the second classification standard is composed of a plurality of element brittleness index threshold intervals and corresponding reservoir types; The predicted skeleton resistivity value is obtained by using XRF element logging data, and the corresponding third reservoir classification result is determined based on the difference between the predicted skeleton resistivity value and the measured skeleton resistivity value, wherein the measured skeleton resistivity value is obtained from the resistivity logging data while drilling; Comprehensively analyzing the first reservoir classification results, the second reservoir classification results, and the third reservoir classification results according to the comprehensive classification standard, and obtaining the classification results of the deep shale oil reservoir to be classified; The gas logging data is used to determine the gas logging total hydrocarbon anomaly multiple, and the first reservoir classification result is obtained in combination with the first classification standard, wherein the first classification standard is composed of multiple gas logging total hydrocarbon anomaly multiple threshold intervals and corresponding reservoir types, including: The total hydrocarbon value is determined using gas logging data and corrected based on drilling environment data; in, , The drilling fluid density collected in real time for logging. is the gravitational acceleration constant, is the well depth collected in real time by logging; ZS is the drilling time required to drill into 1m of the formation; TG is the total hydrocarbon value; is the bottom formation temperature, according to the well depth data of adjacent wells Temperature and regional geothermal gradient calculate, ; The outlet drilling fluid temperature collected in real time for logging; is the atmospheric pressure on the ground; The total hydrocarbon abnormality multiple is obtained by using the overall trend line of total hydrocarbon value in the whole well section of the shale oil reservoir and the corrected total hydrocarbon value; Among them, TGy is the total hydrocarbon anomaly multiple; TGx is the corrected total hydrocarbon value; y is the overall change trend line of the total hydrocarbon value of the entire well section of the shale oil reservoir; The total hydrocarbon anomaly multiple is brought into the first classification standard to determine the corresponding first reservoir classification result, wherein the first classification standard is composed of a plurality of gas logging total hydrocarbon anomaly multiple threshold intervals and corresponding reservoir types; Among them, the skeleton resistivity prediction value is obtained by using XRF element logging data, and the corresponding third reservoir classification result is determined based on the difference between the skeleton resistivity prediction value and the skeleton resistivity measured value, including: Select multiple sensitive elements in the XRF element logging data, input them into the skeleton resistivity prediction model, and obtain the skeleton resistivity prediction value, wherein the skeleton resistivity prediction model is obtained by training the XGBoost regression model using a number of samples, each of the several samples includes multiple sensitive elements and the identification of the skeleton resistivity measured value, and the multiple sensitive elements are obtained by obtaining the gain value of each element in the XRF element logging data and screening out the elements whose gain value is greater than the threshold value; The difference between the predicted skeleton resistivity value and the measured skeleton resistivity value is calculated, and is brought into the third classification standard to determine the corresponding third reservoir classification result, wherein the third classification standard is composed of multiple difference threshold intervals and corresponding reservoir types.

6. An electronic device, characterized in that: The method comprises a processor and a memory, wherein a computer program is stored in the memory, and the computer program is loaded and executed by the processor to implement the steps in the method according to any one of claims 1 to 4.

Images