Sandstone sticking method, device, medium, and electronic device

Abstract

The application discloses a sandstone layer sticking method, device, medium and electronic equipment. The method comprises the following steps: predicting the top surface depth of the target well to be stuck in the sandstone layer; reducing the engineering parameter collection interval, and acquiring the logging engineering parameter at the current depth; maintaining the constant drilling fluid performance and the stable logging engineering parameter; calculating the normalized coefficient of the rotary table revolution, the normalized coefficient of the standpipe pressure, the drilling time correction parameter and the torque correction parameter; calculating the drilling time correction parameter change rate, the drilling pressure change rate, the torque correction parameter change rate and the total engineering parameter change rate, and establishing the sandstone layer sticking standard; and performing the target sandstone layer sticking and the rapid identification of the sandstone layer. The technical scheme provided by the application can improve the accuracy of the sandstone layer sticking and the identification ability of the permeable thin sandstone layer, effectively reduces the geological cycle times of the layer sticking operation, and reduces the drilling cost.

Description

Technical Field

[0001] This invention application belongs to the field of oil drilling and exploration technology, and particularly relates to a method, device, medium and electronic equipment for sandstone layer trapping. Background Technology

[0002] Sandstone is the most common type of stratigraphic rock in oilfields both domestically and internationally. Some sandstone formations can serve as oil and gas reservoirs, as well as as marker beds for stratigraphic deposition. Evaluating sandstone reservoirs can clarify their petrological characteristics, diagenetic features, and physical properties. Analyzing sandstone marker beds can reveal their development and allow for the study of stratigraphic depositional characteristics using methods such as stratigraphic correlation. Therefore, during drilling, the timely adjustment of drilling parameters based on sandstone core samples is essential for accurately obtaining sandstone cores.

[0003] Chinese patent CN113495135A discloses a method and apparatus for determining the salt bottom layer of a massive gypsum-salt rock. The method involves sampling rock fragments at a first location in the massive gypsum-salt rock to be analyzed, detecting the content of chemical elements in the rock fragments, calculating a first salt bottom index, and combining this with the content of magnesium and chloride ions to determine the salt bottom layer of the massive gypsum-salt rock. This method improves the accuracy of determining the salt bottom layer of massive gypsum-salt rocks. However, the detection cycle for magnesium and chloride ion content in a single sample is long, resulting in poor timeliness of operational guidance when used for on-site layer determination.

[0004] Chinese patent CN104133250A discloses a geological layer-locking method for gypsum-salt formations. This method predicts the top and bottom depths and thicknesses of gypsum-salt layers through seismic data interpretation and combines this with changes in well logging, gas logging, and drilling fluid performance parameters to achieve better layer-locking of gypsum-salt formations. However, this method does not consider the error in predicted depth caused by the similarity of seismic data between some salt rock layers and gypsum-salt layers. Furthermore, this method is not effective in locking thin layers of gypsum-salt or mud-bearing gypsum-salt.

[0005] Chinese patent CN111305750A discloses a method for salt bottom detection based on drilling time logging. It utilizes the correlation between drilling time and lithology to determine the lithology of the formation to be drilled. Based on the drilling time data, real-time lithology, and the lithology of the drilled and exposed formations, it determines whether the salt bottom has been reached. This invention uses information such as drilling time logging to determine whether the salt bottom has been reached, ensuring the reliability of the judgment results and achieving salt bottom detection. Furthermore, the basic data is easy to obtain, the steps are simple, highly targeted, and extremely easy to implement. This invention requires maintaining constant logging torque, pump pressure, and drilling pressure parameters during drilling. In actual operation, it is difficult to artificially maintain the stability of these parameters. At the same time, there will inevitably be changes in torque, drilling pressure, and other parameters caused by geological factors at different lithological interfaces or transition layers. Relying solely on drilling time characteristics is also insufficient to accurately identify muddy and dense salt bottoms.

[0006] Chinese patent CN112049627A discloses a method for rapidly identifying the lithology of Upper Paleozoic strata during drilling. This method primarily utilizes gas logging characteristic parameters to establish a linear relationship diagram for identifying gas-bearing sandstone / coal seams, thus determining the interpretation layer of the lithology to be identified. This method is simple, rapid, and accurate, providing a reliable basis for accurately pinpointing completed drilling layers, minimizing ineffective drilling footage, and selecting gas testing layers post-drilling. However, this method uses gas logging as the primary technique for identifying gas-bearing sandstone / coal seams. It does not consider that changes in the gas logging component content after the gas-bearing sandstone and coal seam become water-bearing will alter the relationship diagram of the gas logging characteristic parameters, thereby affecting the accurate identification of gas-bearing sandstone / coal seams.

[0007] Chinese patent CN112049627A discloses a method for predicting stuck layers in high-quality shale core sampling. This method utilizes seismic data combined with elemental analysis, geochemical logging, cuttings color variations, and seismic prediction to identify early warning points for core sampling. Simultaneously, it closely integrates with firsthand field data to analyze the location of these warning points in real time, predicting the depth of the core sampling location in a timely manner. This guides the core sampling work, ensuring accurate core placement and achieving cost savings, shorter construction periods, and reduced project risks. However, this invention primarily focuses on lithological analysis and identification after cuttings return, resulting in poor timeliness for on-site operational guidance. Furthermore, it does not describe how to assess formation lithology and control drilling during the drilling process.

[0008] Currently, published patents mainly focus on lithological analysis and layer identification of gypsum-salt layers, salt rock layers, and coal seams. Published research papers primarily utilize micro-drilling time analysis (MDTA) methods for identifying sandstone and mudstone formations. This typically involves analyzing the drilling time changes corresponding to a 0.1m-0.2m logging interval near the sandstone layer to be identified, confirming whether the sandstone layer has been reached. However, in practice, relying solely on MDTA methods cannot accurately identify thin, uneven sandstone layers, sandstone layers without significant drilling time differences from the upper layers, or layers with poor inter-well correlation. Furthermore, current patents and research papers do not present a technical method that can effectively identify lithologies with good permeability without requiring cuttings and drilling fluid to return to the wellhead.

[0009] In view of the above problems and practical needs, in order to accurately capture thin sandstone layers, sandstone layers that do not have obvious drilling time differences with the upper mudstone, and sandstone layers with poor inter-well contrast, and to quickly identify sandstone layers with good permeability during drilling, effectively reduce the number of geological cycles for layer capture operations, reduce drilling costs, and improve the accuracy of sandstone layer capture and the ability to identify thin permeable sandstone layers, this invention provides a sandstone layer capture method to achieve the above technical requirements. Summary of the Invention

[0010] This invention provides a sandstone layer-locking method, apparatus, medium, and electronic equipment. The sandstone layer-locking method provided by this invention can improve the accuracy of sandstone layer-locking and the ability to identify thin, permeable sand layers. It can also effectively reduce the number of geological cycles required for layer-locking operations, thereby reducing drilling costs.

[0011] Other features and advantages of this invention will become apparent from the following detailed description, or may be learned in part by practice of this invention.

[0012] According to a first aspect of the present invention, a method for sandstone layering is provided, characterized in that the method comprises:

[0013] S1. Based on the actual drilling formation data of adjacent wells and the seismic attribute data of the target well, predict the top surface depth of the sandstone layer to be blocked in the target well;

[0014] S2. At the set depth of mudstone above the top surface of the sandstone layer to be blocked, reduce the interval of engineering parameter acquisition and obtain the logging engineering parameters at the current depth. The logging engineering parameters include drilling time, drilling pressure, torque, standpipe pressure and rotary table revolutions.

[0015] S3. Maintain constant drilling fluid properties and keep logging parameters stable;

[0016] S4. Calculate the rotary table speed normalization coefficient and use the rotary table speed normalization coefficient to calculate the drilling time correction parameters;

[0017] S5. Calculate the riser pressure normalization coefficient and use the riser pressure normalization coefficient to calculate the torque correction parameters.

[0018] S6. Calculate the rate of change of drilling correction parameters, the rate of change of drilling pressure, and the rate of change of torque correction parameters;

[0019] S7. Calculate the total rate of change of engineering parameters;

[0020] S8. Establish a sandstone strata standard based on drilling time correction parameters, torque correction parameters, drilling pressure, rate of change of correction parameters, rate of change of drilling pressure, rate of change of torque correction parameters, and total rate of change of engineering parameters.

[0021] S9. Based on the sandstone layer standard, perform target sandstone layer identification and rapid sandstone layer identification.

[0022] In some embodiments of this invention application, based on the foregoing scheme, the calculation of the turntable rotation speed normalization coefficient and the riser pressure normalization coefficient includes:

[0023] The normalization coefficient for the turntable rotation speed is t. RPM The formula for calculating the normalization coefficient of the turntable rotation speed is as follows:

[0024]

[0025] Among them, t RPMi Let RPM be the normalization factor for the turntable revolutions per minute (RPM) at depth i. i Let be the number of revolutions per minute (RPM) of the turntable at depth i. This represents the average number of revolutions per minute (RPM) of the turntable between depths from point 1 to point i-1.

[0026] The normalization coefficient for the riser pressure is t. SPP The formula for calculating the normalization coefficient of the riser pressure is as follows:

[0027]

[0028] Among them, t SPPi Let SPP be the normalized coefficient for riser pressure (SPP) at depth i. i Let be the riser pressure (SPP) value at depth i. The average riser pressure (SPP) between depths from point 1 to point i-1.

[0029] In some embodiments of this invention application, based on the foregoing scheme, the calculation of the drilling time correction parameters and the torque correction parameters includes:

[0030] The drilling time correction parameter is X. ROP The calculation formula for the drilling time correction parameter is as follows:

[0031]

[0032] Among them, X ROPi Let ROP be the drilling time (ROP) correction value at the depth of point i. i Let t be the drilling time (ROP) value at the depth of point i. RPMi is the normalization factor for the number of revolutions per minute (RPM) of the turntable at depth i.

[0033] The torque correction parameter is X. TOR The calculation formula for the torque correction parameter is as follows:

[0034]

[0035] Among them, X TORi The torque (TOR) correction value at depth i is given by TOR. i Let t be the torque (TOR) value at depth i. SPPi is the normalization coefficient for riser pressure (SPP) at depth i.

[0036] In some embodiments of this invention application, based on the foregoing scheme, calculating the rate of change of the drilling time correction parameter, the rate of change of the drilling pressure, and the rate of change of the torque correction parameter includes:

[0037] The rate of change of the drilling time correction parameter is α ROP The formula for calculating the rate of change of the drilling correction parameters is as follows:

[0038]

[0039] Where, α ROPi X is the depth of point i. ROP Relative rate of change, X ROPi Let be the drilling time (ROP) correction value at the depth of point i. X is the depth between point i and point i-5. ROP The average value;

[0040] The rate of change of drilling pressure is α WOB The formula for calculating the rate of change of drilling pressure is as follows:

[0041]

[0042] Where, α WOBi Let WOB be the relative rate of change of depth at point i. i Let be the drill weight (WOB) value at the depth of point i. This represents the average drill pressure (WOB) value between depths from point i to depths i-5.

[0043] The rate of change of the torque correction parameter is α TOR The formula for calculating the rate of change of the torque correction parameter is as follows:

[0044]

[0045] Where, α TORi X is the depth of point i. TOR Relative rate of change, X TORi The torque (TOR) correction value at depth i is given. X is the depth between point i and point i-5. TOR The average value.

[0046] In some embodiments of this invention application, the calculation of the total rate of change of the engineering parameters based on the foregoing scheme includes:

[0047] The total rate of change of the engineering parameters is β, and the formula for calculating the total rate of change of the engineering parameters is as follows:

[0048] β i =α ROPi +α WOBi+α TORi

[0049] Where, β i Let α be the total rate of change of engineering parameters at depth i. ROPi X is the depth of point i. ROP relative rate of change, α WOBi Let α be the relative rate of change of WOB at the depth of point i. TORi X is the depth of point i. TOR Relative rate of change.

[0050] In some embodiments of this invention application, based on the aforementioned scheme, the target sandstone layering is performed, including:

[0051] N1. Perform steps S1 to S7 of the sandstone layer trapping method, and calculate the drilling time correction parameter, torque correction parameter, drilling pressure, correction parameter change rate, drilling pressure change rate, torque correction parameter change rate and total engineering parameter change rate.

[0052] N2. Based on the sandstone layering standard in step S8 of the sandstone layering method, determine the lithology of the current stratum;

[0053] N3. If the drilling time correction parameters, torque correction parameters, drilling pressure, correction parameter change rate, drilling pressure change rate, torque correction parameter change rate, and total change rate of engineering parameters of the current formation to be analyzed all meet the sandstone characteristics in the sandstone stratigraphic standard described in claim 1, then the current formation is determined to be the target sandstone layer, and geological circulation is performed. The lithology of the rock cuttings returned from the wellhead is determined to be the target layer sandstone. If the returned rock cuttings are the target layer sandstone, then core drilling is performed. If the returned rock cuttings are not the target layer sandstone, then drilling continues, and steps N1 to N2 are repeated.

[0054] In some embodiments of this invention application, based on the aforementioned scheme, rapid identification of sandstone layers includes:

[0055] L1. Perform steps S1 to S7 of the sandstone layer trapping method, and calculate drilling time correction parameters, torque correction parameters, drilling pressure, correction parameter change rate, drilling pressure change rate, torque correction parameter change rate, and total engineering parameter change rate.

[0056] L2. Based on the sandstone layering standard in step S8 of the sandstone layering method, determine the lithology of the current stratum;

[0057] L3. If the drilling time correction parameters, torque correction parameters, drilling pressure, correction parameter change rate, drilling pressure change rate, torque correction parameter change rate, and total change rate of engineering parameters of the current formation to be analyzed all meet the characteristics of sandstone or argillaceous sandstone in the sandstone stratum standard described in claim 1, then record the top and bottom depths of the sandstone layer or argillaceous sandstone layer, continue drilling operations, and re-execute steps L1 to L2.

[0058] According to a second aspect of the present invention, a sandstone layer trapping device is provided, characterized in that the device comprises: a prediction unit, configured to predict the top depth of the sandstone layer to be trapped in the target well based on actual drilling formation data from adjacent wells and seismic attribute data from the target well; an acquisition unit, configured to acquire logging engineering parameters at a predetermined depth above the top depth of the sandstone layer to be trapped, by reducing the engineering parameter acquisition interval, the logging engineering parameters including drilling time, drilling pressure, torque, standpipe pressure, and rotary table revolutions; and a calculation unit, configured to calculate drilling time correction parameters, torque correction parameters, drilling time correction parameter change rate, drilling pressure change rate, and torque correction parameters based on the acquired logging engineering parameters. The unit is configured to determine the lithology of the strata at the current depth based on the drilling time correction parameters, the torque correction parameters, the drilling pressure, the rate of change of the drilling time correction parameters, the rate of change of the drilling pressure, the rate of change of the torque correction parameters, the rate of change of the total change of the engineering parameters, and the sandstone layer standard. It is also configured to perform target sandstone layer identification and rapid sandstone layer recognition based on the sandstone layer standard. The sandstone layer standard is used to characterize the qualitative features and quantitative numerical distribution ranges of the drilling time correction parameters, torque correction parameters, drilling pressure, rate of change of the drilling time correction parameters, rate of change of the drilling pressure, rate of change of the torque correction parameters, and the rate of change of the total change of the engineering parameters corresponding to different strata.

[0059] According to a third aspect of the present invention, a computer-readable storage medium is provided, characterized in that the computer-readable storage medium stores at least one piece of program code, the at least one piece of program code being loaded and executed by a processor to implement the operations performed by the method.

[0060] According to a fourth aspect of the present invention, an electronic device is provided, characterized in that the electronic device includes one or more processors and one or more memories, wherein the one or more memories store at least one piece of program code, the at least one piece of program code being loaded and executed by the one or more processors to implement the operation performed by the method.

[0061] Compared with the prior art, the present invention has at least the following beneficial effects:

[0062] Before predicting the top depth of the sandstone layer to be blocked, this invention establishes a sandstone layer blocking standard by conducting actual drilling operations on adjacent wells. The sandstone layer blocking standard is used to characterize the drilling time correction parameters, torque correction parameters, drilling pressure, drilling time correction parameter change rate, drilling pressure change rate, torque correction parameter change rate, and total change rate of engineering parameters corresponding to different formations.

[0063] After establishing the sandstone layer standard, the top depth of the sandstone layer to be trapped is predicted based on the actual drilling data of adjacent wells and the seismic attribute data of the target well. Then, at a predetermined depth of mudstone above the top depth of the sandstone layer to be trapped, the interval of the logging instrument's engineering parameter acquisition is reduced, and the logging engineering parameters at the current depth are obtained. It should be noted that since the top depth of the sandstone layer to be trapped has been predicted before drilling the target well, reducing the interval of the logging instrument's engineering parameter acquisition can increase the number of logging engineering parameter acquisitions, making the obtained logging engineering parameters more valuable for reference, thereby improving the accuracy of the target sandstone layer trapping in this invention.

[0064] Specifically, during the drilling process of the target well, engineers need to maintain constant drilling fluid properties to ensure the stability of logging parameters. These logging parameters include drilling time, drilling pressure, torque, standpipe pressure, and rotary table revolutions.

[0065] Based on the acquired logging parameters, the drilling time correction parameters, torque correction parameters, drilling time correction parameter change rate, drilling pressure change rate, torque correction parameter change rate, and total change rate of engineering parameters at the current depth are calculated. By comparing the acquired logging parameters and the calculated drilling time correction parameters, torque correction parameters, drilling time correction parameter change rate, drilling pressure change rate, torque correction parameter change rate, and total change rate of engineering parameters with the sandstone layer standard corresponding to the target sandstone layer, it can be determined whether the formation at the current depth is the target sandstone layer. Then, based on the sandstone layer standard, target sandstone layer identification and rapid sandstone layer identification are performed.

[0066] Based on this, the sandstone layer-catching method provided by the present invention can improve the accuracy of sandstone layer-catching and the ability to identify thin, permeable sandstone layers. It can also accurately catch sandstone layers that do not show significant differences in drilling time from the upper layers and sandstone layers with poor inter-well correlation.

[0067] At the same time, since the present invention determines the lithology of the formation by using drilling time correction parameters, torque correction parameters, drilling time correction parameter change rate, drilling pressure change rate, torque correction parameter change rate and total change rate of engineering parameters, it can effectively reduce the number of geological cycles for stuck formation operations, thereby reducing drilling costs.

[0068] It should be understood that the above general description and the following detailed description are exemplary and explanatory only, and do not limit the scope of this application. Attached Figure Description

[0069] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present invention and, together with the description, serve to explain the principles of the present invention. Obviously, the drawings described below are merely some embodiments of the present invention, and those skilled in the art can obtain other drawings based on these drawings without any inventive effort. In the drawings:

[0070] Figure 1 The steps of the sandstone layering method in the embodiments of this application are shown;

[0071] Figure 2 A flowchart of the target sandstone layer in an embodiment of this invention is shown;

[0072] Figure 3 A flowchart illustrating the rapid identification of sandstone layers in an embodiment of this invention is shown.

[0073] Figure 4 An example diagram of a sandstone stratum in an embodiment of this invention is shown;

[0074] Figure 5 A schematic diagram of the sandstone layering device in an embodiment of this invention is shown;

[0075] Figure 6 A schematic diagram of the structure of an electronic device according to an embodiment of the present invention is shown. Detailed Implementation

[0076] The technical solutions of the embodiments of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this invention, and not all of them. Based on the embodiments of this invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this invention.

[0077] Furthermore, the described features, structures, or characteristics can be combined in any suitable manner in one or more embodiments. In the following description, numerous specific details are provided to give a thorough understanding of embodiments of this application. However, those skilled in the art will recognize that the technical solutions of this application can be practiced without one or more of the specific details, or other methods, components, apparatuses, steps, etc., can be employed. In other instances, well-known methods, apparatuses, implementations, or operations are not shown or described in detail to avoid obscuring various aspects of this application.

[0078] The block diagrams shown in the accompanying drawings are merely functional entities and do not necessarily correspond to physically independent entities. That is, these functional entities can be implemented in software, in one or more hardware modules or integrated circuits, or in different network and / or processor devices and / or microcontroller devices.

[0079] The flowcharts shown in the accompanying drawings are merely illustrative and do not necessarily include all content and operations / steps, nor do they necessarily have to be performed in the described order. For example, some operations / steps can be broken down, while others can be combined or partially combined; therefore, the actual execution order may change depending on the specific circumstances.

[0080] The following section will elaborate on this invention application:

[0081] Figure 1 The steps of the sandstone layer-locking method in the embodiments of this application are shown below:

[0082] Step S1: Based on the actual drilling formation data of adjacent wells and the seismic attribute data of the target well, predict the top depth of the sandstone layer to be blocked in the target well.

[0083] The prediction of the top depth of the sandstone layer to be blocked mainly employs two techniques. First, by comparing and analyzing the actual drilling formation data and sand body data from the target well and adjacent wells, the top depth of the sandstone layer to be blocked in the target well is predicted. Second, by analyzing the distribution of seismic attribute parameters such as the root mean square amplitude and wave impedance of the sandstone layer to be blocked in the geological space, the top depth of the sandstone layer to be blocked in the target well is predicted.

[0084] It should be noted that the two techniques mentioned above can be used in combination to further improve the accuracy of predicting the depth of the top surface of the sandstone layer to be sealed through multi-layer comparison, layer-by-layer approximation, and seismic attribute optimization.

[0085] Step S2: At the set depth of mudstone above the top surface of the sandstone layer to be logged, reduce the interval of engineering parameter acquisition and obtain the logging engineering parameters at the current depth. The logging engineering parameters include drilling time, drilling pressure, torque, standpipe pressure and rotary table revolutions.

[0086] By reducing the logging tool's parameter acquisition interval at a set depth above the top surface of the sandstone formation to be captured, more intensive acquisition of drilling time (ROP), bit pressure (WOB), torque (TOR), standpipe pressure (SPP), and rotary table revolutions (RPM) data can be achieved. This allows for better characterization of the detailed changes in formation lithology across these logging parameters, thus enabling precise capture of the top surface of the sandstone formation. Specifically, acquiring the logging parameters at the current depth refers to acquiring the drilling time, bit pressure, torque, standpipe pressure, and rotary table revolutions corresponding to the current formation depth each time the drilling equipment advances by one interval of the parameter acquisition interval.

[0087] Specifically, the set depth can be a different distance, such as 10m, 20m, or 30m above the sand layer to be trapped. The reduced engineering parameter acquisition interval can be a different interval, such as 0.1m, 0.2m, or 0.5m.

[0088] It should be noted that, since the top depth of the sandstone layer to be trapped was predicted before drilling, and given that the top depth of the sandstone layer is typically quite large, and the logging parameters for a significant portion of the formation above the top depth of the sandstone layer have limited reference value for achieving the target sandstone layer capture, the drilling process can utilize normal parameter acquisition intervals to improve work efficiency.

[0089] Step S3: Maintain constant drilling fluid properties and keep logging parameters stable.

[0090] The term "constant drilling fluid performance" refers to the fact that, during the drilling process to the sandstone where the drilling fluid is stuck, engineers do not arbitrarily change the parameters such as drilling fluid density, viscosity, and shear force when measuring these parameters, while maintaining minimal fluctuations in these parameters during the sandstone stuck process.

[0091] The stability of logging parameters refers to the fact that, during drilling into the stuck sandstone, engineers cannot manually or arbitrarily change the values ​​of logging parameters such as drilling time (ROP), weight on bit (WOB), torque (TOR), standpipe pressure (SPP), and rotary table revolutions (RPM). However, variations in these logging parameters due to different formation rock properties are permissible. Specifically, directional drilling operations cannot be carried out during the stuck sandstone drilling process; only conventional drilling operations with the rotary table are allowed.

[0092] Step S4: Calculate the rotary table rotation normalization coefficient and use the rotary table rotation normalization coefficient to calculate the drilling time correction parameters.

[0093] The normalization coefficient for the turntable rotation speed is t.RPM The formula for calculating the normalization coefficient of the turntable rotation speed is as follows:

[0094]

[0095] Among them, t RPMi Let RPM be the normalization factor for the number of revolutions per minute (RPM) at depth i, which is dimensionless; RPM i Let be the number of revolutions per minute (RPM) of the turntable at the depth of point i, in r / min; The average rotary table revolutions (RPM) between the depths of point 1 and point i-1 is given in r / min; the depth of point 1 is the drilling depth corresponding to the reduced spacing of the engineering parameter acquisition in step S2 above.

[0096] The drilling time correction parameter is X. ROP The calculation formula for the drilling time correction parameter is as follows:

[0097]

[0098] Among them, X ROPi The ROP (Recovery Time Per Minute) correction value for the depth at point i, in min / m; ROP i t represents the drilling time (ROP) at the depth of point i, in min / m. RPMi is the normalization coefficient for the number of revolutions per minute (RPM) of the turntable at depth i, which is dimensionless.

[0099] During drilling, with a fixed weight on bit (WOB), drilling time (ROP) is significantly affected by rotary table revolutions (RPM). As the RPM increases, the rate at which the drill bit breaks rocks per unit time at the bottom of the well increases, and the drilling time (ROP) decreases accordingly. Therefore, by correcting the drilling time (ROP) using the normalization factor of the rotary table revolutions (RPM), the influence of changes in the RPM value on the drilling time (ROP) value can be eliminated.

[0100] Step S5: Calculate the riser pressure normalization coefficient and use the riser pressure normalization coefficient to calculate the torque correction parameters.

[0101] The normalization coefficient for the riser pressure is t. SPP The formula for calculating the normalization coefficient of the riser pressure is as follows:

[0102]

[0103] Among them, t SPPi SPP is the normalized coefficient for riser pressure (SPP) at depth i, dimensionless; i Let S be the riser pressure (SPP) value at depth i, in MPa; The standpipe pressure (SPP) is the average value in MPa between depths from point 1 to depth i-1. The depth of point 1 is the drilling depth corresponding to the reduced interval of engineering parameter acquisition in step S2 above.

[0104] The torque correction parameter is X. TOR The calculation formula for the torque correction parameter is as follows:

[0105]

[0106] Among them, X TORi Let be the torque (TOR) correction value at depth i, in kN.m; TORi be the torque (TOR) value at depth i, in kN.m; t SPPi is the normalization coefficient for riser pressure (SPP) at depth i, which is dimensionless.

[0107] During drilling, with a constant wobble pressure (WOB), torque (TOR) is significantly affected by standpipe pressure (SPP). As the standpipe pressure (SPP) increases, the drilling fluid flow rate in the wellbore and annulus increases, the cuttings return rate at the bottom of the well increases, the annulus becomes relatively cleaner, the rotational resistance of the drill string decreases, and the torque (TOR) decreases. Therefore, by correcting the torque (TOR) using the standpipe pressure (SPP) normalization coefficient, the influence of changes in the standpipe pressure (SPP) value on the torque (TOR) value can be eliminated.

[0108] Step S6: Calculate the rate of change of drilling correction parameters, the rate of change of drilling pressure, and the rate of change of torque correction parameters.

[0109] The rate of change of the drilling time correction parameter is α ROP The formula for calculating the rate of change of the drilling correction parameters is as follows:

[0110]

[0111] Where, α ROPi For the depth of point i XROP Relative rate of change, %; X ROPi is the drilling time (ROP) correction value at the depth of point i, in r / min; X is the depth between point i and point i-5. ROP The average value, r / min. In the α ROP In the calculation formula, the minimum value of i is 6.

[0112] Specifically, based on the engineering parameter acquisition interval settings of the logging tool in step S2 and personal preference, the average value of the drilling time (ROP) correction value can also be set to the average value between 4 or 10 points above the current depth. ROP The rate of change α at the current depth ROPCharacterize X at the current depth ROP Value relative to the current depth before X ROP The magnitude of the average value reflects the current X ROP The degree of drastic increase or decrease in value.

[0113] The rate of change of drilling pressure is α WOB The formula for calculating the rate of change of drilling pressure is as follows:

[0114]

[0115] Where, α WOBi is the relative change rate of WOB at depth i, %; WOBi is the drill weight (WOB) value at depth i, t; t is the average drill pressure (WOB) value between depths from point i to depth i-5. In the α... WOB In the calculation formula, the minimum value of i is 6.

[0116] Specifically, based on the engineering parameter acquisition interval settings of the logging instrument in step S2 and personal preference, the average value of the drilling pressure (WOB) value can also be set to the average value of 4 or 10 points above the current depth.

[0117] The rate of change of the torque correction parameter is α TOR The formula for calculating the rate of change of the torque correction parameter is as follows:

[0118]

[0119] Where, α TORi X is the depth of point i. TOR Relative rate of change, %; X TORi Here is the torque (TOR) correction value at depth i, in kN.m; X is the depth between point i and point i-5. TOR The average value, KN.m. In the α TOR In the calculation formula, the minimum value of i is 6.

[0120] Specifically, based on the engineering parameter acquisition interval settings of the logging instrument in step S2 and personal preference, the average value of the torque (TOR) correction value can also be set to the average value between 4 or 10 points above the current depth.

[0121] Specifically, in the α TOR In the calculation formula, X ROP Average, average drill bit weight (WOB) and X TOR The depth spacing used in the average value calculation must be consistent. For example, X ROPThe average value is calculated using the values ​​from the five depths preceding the current depth. Therefore, the average drill weight (WOB) and X... TOR The average value should also be calculated using the values ​​of the five depth points preceding the current depth point.

[0122] The α TOR The principle behind the calculation formula is that during the drilling process from mudstone strata to relatively loose sandstone strata, the loose nature of sandstone makes it prone to breakage. Therefore, the drilling time per unit thickness of strata is relatively reduced, and the drilling time correction value X... ROP Decrease, X ROP Rate of change α at the current depth ROP The drilling footage per unit time increases, the weight load on the drill string borne by the hook load on the derrick increases, resulting in a decrease in the weight on the drill bit (WOB), and the rate of change of WOB at the current depth α. WOB Increase. As mentioned above, during the drilling process from mudstone to relatively loose sandstone formations, the drilling time (ROP) decreases, the drill bit footage per unit time increases, and the amount of rock cuttings broken by the drill bit relatively increases. With a constant standpipe pressure (SPP) and rotary table revolutions per minute (RPM), the resistance to the drill string rotation within the wellbore increases, and the corresponding torque correction value X increases. TOR Increase, X TOR Rate of change α at the current depth TOR Increase.

[0123] Step S7: Calculate the total rate of change of engineering parameters.

[0124] The total rate of change of the engineering parameters is β, and the formula for calculating the total rate of change of the engineering parameters is as follows:

[0125] β i =α ROPi +α WOBi +α TORi

[0126] Where, β i Let α be the total rate of change of engineering parameters at depth i, in %; ROPi X is the depth of point i. ROP Relative rate of change, %; α WOBi Let α be the relative change rate of WOB at the depth of point i, in %; TORi X is the depth of point i. TOR Relative rate of change, %.

[0127] As described in step S7, the general rule during the drilling process from relatively dense mudstone to relatively loose sandstone is α. ROP α WOB and α TORAll of these factors relatively increase, but different sandstone compositions, structures, diagenetic processes, and tectonic features, as well as variations in drilling fluid properties during the drilling process, can all cause α to increase. ROP α WOB and α TOR Different numerical change responses. That is, among the aforementioned relative rate of change parameters, some parameters increase significantly, while others increase only slightly. The total rate of change β of the engineering parameter effectively characterizes the α value during drilling from mudstone to sandstone. ROP α WOB and α TOR The overall change is the impact of drilling from mudstone to sandstone on α ROP α WOB and α TOR The overall changes are relatively magnified. Therefore, by applying the total change rate β parameter of the engineering parameter, engineers can more accurately determine the top depth of sandstone development and promptly extract sandstone. At the same time, during non-stretching operations, the total change rate β parameter of the engineering parameter can also be used to analyze the hardness of the formation rocks in real time during drilling.

[0128] Step S8: Establish a sandstone strata standard based on drilling time correction parameters, torque correction parameters, drilling pressure, rate of change of correction parameters, rate of change of drilling pressure, rate of change of torque correction parameters, and total rate of change of engineering parameters.

[0129] Based on the measured logging parameters of the adjacent well, the measured formation lithology of the adjacent well, and the calculation formulas in steps S4 to S7, the X corresponding to different formations in the adjacent well can be calculated. ROP X TOR WOB, α ROP α TOR α WOB And the β parameter. Data calculated based on steps S4 to S7 indicates that when the drill bit reaches a non-permeable, relatively dense mudstone layer, X ROP WOB has high parameter values, while X TOR α ROP α TOR α WOB And β have low parameter values; when the drill bit drills into a relatively loose sandstone layer, X ROP WOB has lower parameter values, while X TOR α ROP α TOR α WOB And β have high parameter values. Therefore, X values ​​of mudstone, argillaceous sandstone, and sandstone that have been verified by drilling core samples were selected. ROP X TOR WOB, α ROP α TOR α WOB β data, analyze the XROP X TOR WOB, α ROP α TOR α WOB The distribution of the β parameter in mudstone, argillaceous sandstone, and sandstone lithology was analyzed, and a model based on the X parameter was established. ROP X TOR WOB, α ROP α TOR α WOB The sandstone layer standard with β parameter can be used to accurately capture the target sandstone layer.

[0130] It should also be noted that the numerical variation ranges for different lithologies defined in the sandstone stratigraphic standard will vary in different regions and under different drilling techniques. This invention only provides one example. Referring to Table 1, one sandstone stratigraphic standard from an embodiment of this invention is shown.

[0131] Table 1

[0132]

[0133] In particular, mudstone, argillaceous sandstone, and sandstone at different strata in different regions possess varying stress fields, rock structures and textures, and diagenetic processes, while also being influenced by drilling engineers' preferences for controlling drilling parameters. Therefore, it is difficult to form specific X-type stratigraphic patterns within the mudstone, argillaceous sandstone, and sandstone categories mentioned in the sandstone stratigraphic standard. ROP X TOR The WOB value range can be specified, but the X value of the mudstone above the adjacent sandstone layer to be tested can be included. ROP X TOR Using WOB values ​​as a baseline, the lower lithology X was analyzed. ROP X TOR The relative variation and trend of WOB data values ​​can be used to determine whether sandstone formations have been encountered during drilling.

[0134] α ROP α TOR α WOB Mainly used to characterize X ROP X TOR The relative rate of change of the WOB parameter curves exhibits similar numerical variation ranges in mudstone, argillaceous sandstone, and sandstone that have undergone core sampling and lithological verification in different regions and at different stratigraphic levels. Therefore, sandstone stratigraphic standards can be established based on specific numerical variation ranges. Similarly, β is characterized by α. ROP α TOR and α WOB The overall rate of change can also be used to establish sandstone stratigraphic standards by measuring the numerical variation ranges on mudstone, argillaceous sandstone, and sandstone.

[0135] Step S9: Based on the sandstone layer standard, perform target sandstone layer identification and rapid sandstone layer recognition.

[0136] When locating the target sandstone layer, the following steps can be performed:

[0137] N1. Perform steps S1 to S7 of the sandstone layer trapping method to calculate drilling time correction parameters, torque correction parameters, drilling pressure, correction parameter change rate, drilling pressure change rate, torque correction parameter change rate, and total change rate of engineering parameters.

[0138] N2. Based on the sandstone layering standard in step S8 of the sandstone layering method, determine the lithology of the current stratum.

[0139] N3. If the drilling time correction parameters, torque correction parameters, drilling pressure, correction parameter change rate, drilling pressure change rate, torque correction parameter change rate, and total change rate of engineering parameters of the current formation to be analyzed all meet the sandstone characteristics in the sandstone stratigraphic standard described in claim 1, then the current formation is determined to be the target sandstone layer, and geological circulation is performed. The lithology of the rock cuttings returned from the wellhead is determined to be the target layer sandstone. If the returned rock cuttings are the target layer sandstone, then core drilling is performed. If the returned rock cuttings are not the target layer sandstone, then drilling continues, and steps N1 to N2 are repeated.

[0140] Specifically, refer to Figure 2 The diagram illustrates the process of identifying the target sandstone layer in an embodiment of this invention. The main workflow involves calculating X at a predetermined depth in the mudstone above the top surface of the sandstone layer to be identified. ROP X TOR And obtain α ROP α TOR α WOB And β value. Using the sandstone layer standard in step S8, determine whether the current depth has entered a sandstone layer. If the target sandstone layer has been entered, stop drilling, perform geological circulation, and analyze the cuttings returned from the wellhead. If the cuttings lithology is the expected sandstone, proceed with subsequent drilling and coring operations. If the cuttings lithology is not the expected sandstone, continue normal drilling operations, and recalculate X based on the acquisition interval according to the engineering parameters. ROP X TOR and regain α ROP α TOR α WOB And β value, and at the same time use the sandstone strata standard in step S8 to determine whether the subsequent strata enter the sandstone layer.

[0141] When quickly identifying sandstone layers, the following steps can be performed:

[0142] L1. Perform steps S1 to S7 of the sandstone layer trapping method to calculate drilling time correction parameters, torque correction parameters, drilling pressure, correction parameter change rate, drilling pressure change rate, torque correction parameter change rate, and total change rate of engineering parameters.

[0143] L2. Based on the sandstone layering standard in step S8 of the sandstone layering method, determine the lithology of the current stratum.

[0144] L3. If the drilling time correction parameters, torque correction parameters, drilling pressure, correction parameter change rate, drilling pressure change rate, torque correction parameter change rate, and total change rate of engineering parameters of the current formation to be analyzed all meet the characteristics of sandstone or argillaceous sandstone in the sandstone stratum standard described in claim 1, then record the top and bottom depths of the sandstone layer or argillaceous sandstone layer, continue drilling operations, and re-execute steps L1 to L2.

[0145] Specifically, refer to Figure 3 The diagram illustrates a flowchart of a rapid sandstone layer identification process according to an embodiment of this invention. The main workflow involves calculating X at a predetermined depth above the top surface of the sandstone layer to be identified, within which mudstone is positioned. ROP X TOR And obtain α ROP α TOR α WOB And β value. Using the sandstone layer standard in step S8, determine whether the current depth has entered the sandstone layer. If the target sandstone layer has been entered, record the top depth of the sandstone layer, and simultaneously determine the X value. ROP X TOR WOB, α ROP α TOR α WOB Based on the changes in the β parameter value, and using the sandstone layer standard in step S8, the bottom depth of the sandstone layer is recorded, and normal drilling continues. The next target sandstone layer is also determined using the sandstone layer standard in step S8. Specifically, during this process, the rock cuttings returned from the wellhead can be analyzed to further determine the lithology of the current sandstone layer.

[0146] Based on the sandstone layer trapping method described in this invention application, the coring operation process of the target well will be described below with a specific embodiment.

[0147] Specifically, refer to Figure 4 The diagram illustrates an example of a sandstone strata in an embodiment of this invention. The normalization coefficient t of the rotary disk revolutions per minute (RPM) is used. RPM Calculate the drilling time (ROP) correction parameter X ROP Using the riser pressure (SPP) normalization coefficient t SPP Calculate the torque (TOR) correction parameter X TOR At the same time, X was calculated. ROP rate of change αROP Rate of change of drilling pressure (WOB) α WOB and X TOR rate of change α TOR And the total rate of change of engineering parameters β.

[0148] Using sandstone stratigraphic standards, the lithology at the current depth is determined. In the target well (2501.1m–2501.7m), the reduction in drilling time (ROP) is not significant, the decrease in drilling pressure (WOB) is relatively small, and the increase in torque (TOR) is significant. WOB and α TOR The increase is significant, β increases significantly, and according to the sandstone stratigraphic standard in this method, α WOB α TOR The β parameter value is within the sandstone parameter range, so the 2501.1m to 2501.7m range is identified as sandstone. The logging parameters of this layer with sandstone properties are not high, and the duration is short, meaning the formation thickness is thin. Therefore, drilling continues.

[0149] 2502.9m~2503.4m strata α ROP No increase, α WOB The increase was not significant, only α TOR The formation exhibits increasing characteristics. Based on the sandstone layer standard in this method, the area from 1502.9m to 1503.4m is determined to be argillaceous sandstone, not the target layer sandstone. Drilling should continue.

[0150] For target wells deeper than 2505.5m, drilling time (ROP) and drilling pressure (WOB) gradually decrease, while torque (TOR) gradually increases, and correspondingly α ROP α WOB α TOR α and β also gradually increase. At depths above 2506.6 m, drilling time (ROP) and drilling pressure (WOB) decrease significantly, indicating a significant increase in formation drillability, while torque (TOR) increases significantly, indicating a significant increase in cuttings content within the wellbore. Correspondingly, α... ROP α WOB α TOR The values ​​of α and β are significantly increased. Applying the sandstone stratigraphic standard in this method, all logging parameters at depths above 2506.6m conform to sandstone characteristics. Furthermore, when the drill bit reached 2507.0m, the characteristics of the aforementioned logging parameters remained relatively stable. Simultaneously, the curves corresponding to each logging parameter showed a further trend of sandstone response. Therefore, 2505.6m was determined to be the top depth of the target sandstone layer.

[0151] After determining 2505.6m as the top depth of the target sandstone layer, geological circulation was carried out. The rock cuttings returned to the wellhead from 2505.6m to 2507.0m were sandstone, which is consistent with the characteristics of the sandstone layer to be blocked. Therefore, drilling and coring operations can be carried out.

[0152] Following core drilling, significant sandstone was encountered at the top of the core sample at 2507.0m. Wireline logging after the target well was completed also revealed significant sandstone in the natural gamma ray (GR) and spontaneous potential (SP) curves between 2505.6m and 2507m. This also confirmed the development of a thin sandstone layer between 2501.1m and 2501.7m in the target well. Using this method, accurate sandstone extraction was achieved with only one geological cycle, improving sandstone layer extraction efficiency and reducing drilling costs.

[0153] Based on the same inventive concept, this application also provides a sandstone layering device, referring to... Figure 5 This illustration shows a structural schematic diagram of a sandstone layer-trapping device according to an embodiment of the present invention. The sandstone layer-trapping device 500 includes: a prediction unit 501, used to predict the top depth of the sandstone layer to be trapped in the target well based on actual drilling formation data from adjacent wells and seismic attribute data from the target well; an acquisition unit 502, used to acquire logging engineering parameters at a set depth of mudstone above the top depth of the sandstone layer to be trapped, by reducing the engineering parameter acquisition interval, the logging engineering parameters including drilling time, drilling pressure, torque, standpipe pressure, and rotary table revolutions; and a calculation unit 503, used to calculate drilling time correction parameters, torque correction parameters, drilling time correction parameter change rate, drilling pressure change rate, torque correction parameter change rate, and engineering parameters based on the acquired logging engineering parameters. Total rate of change; Judgment unit 504 is used to determine the lithology of the strata at the current depth based on the drilling time correction parameters, the torque correction parameters, the drilling pressure, the rate of change of the drilling time correction parameters, the rate of change of the drilling pressure, the rate of change of the torque correction parameters, the total rate of change of the engineering parameters, and the sandstone layer standard, and to perform target sandstone layer identification and rapid identification of sandstone layers based on the sandstone layer standard. The sandstone layer standard is used to characterize the qualitative characteristics and quantitative numerical distribution range of the drilling time correction parameters, torque correction parameters, drilling pressure, rate of change of the drilling time correction parameters, rate of change of the drilling pressure, rate of change of the torque correction parameters, and the total rate of change of the engineering parameters corresponding to different strata.

[0154] For details not disclosed in the device embodiments of this invention, please refer to the embodiments of the method described above in this invention.

[0155] Based on the same inventive concept, this application also provides a computer-readable storage medium, characterized in that the computer-readable storage medium stores at least one piece of program code, which is loaded and executed by a processor to implement the operation performed by the method.

[0156] Based on the same inventive concept, this application also provides an electronic device, referring to... Figure 6 , Figure 6A schematic diagram of the structure of an electronic device according to an embodiment of the present invention is shown.

[0157] The electronic device includes one or more memories 604, one or more processors 602, and at least one computer program (program code) stored in the memory 604 and executable on the processor 602, wherein the processor 602 executes the computer program to implement the method described above.

[0158] Among them, Figure 6 In this document, a bus architecture (represented by bus 600) is used. Bus 600 may include any number of interconnected buses and bridges, linking various circuits including one or more processors represented by processor 602 and memory represented by memory 604. Bus 600 may also link various other circuits such as peripheral devices, voltage regulators, and power management circuits, which are well known in the art and therefore will not be described further herein. Bus interface 605 provides an interface between bus 600 and receiver 601 and transmitter 603. Receiver 601 and transmitter 603 may be the same element, i.e., a transceiver, providing a unit for communicating with various other devices over a transmission medium. Processor 602 is responsible for managing bus 600 and general processing, while memory 604 can be used to store data used by processor 602 during operation.

[0159] The functions described herein may be implemented in hardware, software executed by a processor, firmware, or any combination thereof. If implemented in software executed by a processor, the functions may be stored as one or more instructions or codes on or transmitted via a computer-readable medium. Other examples and embodiments are within the scope and spirit of this application and the appended claims. For example, due to the nature of software, the functions described above may be implemented using software executed by a processor, hardware, firmware, hardwired, or any combination thereof. Furthermore, the functional units may be integrated into a single processing unit, or each unit may exist physically separately, or two or more units may be integrated into a single unit.

[0160] In the several embodiments provided in this application, it should be understood that the disclosed technical content can be implemented in other ways. The device embodiments described above are merely illustrative; for example, the division of units can be a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the displayed or discussed mutual coupling, direct coupling, or communication connection may be through some interfaces; the indirect coupling or communication connection of units or modules may be electrical or other forms.

[0161] The units described as separate components may or may not be physically separate. Similarly, the components of the control device may or may not be physical units; they may be located in one place or distributed across multiple units. Some or all of the units can be selected to achieve the purpose of this embodiment, depending on actual needs.

[0162] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this invention application, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this invention application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, read-only memory (ROM), random access memory (RAM), portable hard drives, magnetic disks, or optical disks.

[0163] The above description is merely an embodiment of the present invention and is not intended to limit the present invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of the claims of the present invention.

Claims

1. A sandstone layering method, characterized by, The method steps include: S1, predicting the top surface depth of the target well according to the adjacent well real drilled formation data and the target well seismic attribute data; S2, reducing the engineering parameter acquisition interval at the depth of the mudstone above the top surface depth of the target sandstone layer, and obtaining the logging engineering parameters at the current depth, the logging engineering parameters including drilling time, drilling pressure, torque, standpipe pressure and rotary table revolutions; S3, maintaining the drilling fluid performance constant and keeping the logging engineering parameters stable; S4, calculating the rotary table revolution normalization coefficient and using the rotary table revolution normalization coefficient to calculate the drilling time correction parameter; S5, calculating the standpipe pressure normalization coefficient and using the standpipe pressure normalization coefficient to calculate the torque correction parameter; S6, calculating the drilling time correction parameter change rate, the drilling pressure change rate and the torque correction parameter change rate; Wherein, the drilling time correction parameter change rate, the drilling pressure change rate and the torque correction parameter change rate are calculated, including: The rate of change of the drilling time correction parameter is α ROP The formula for calculating the rate of change of the drilling time correction parameter is as follows: wherein, α ROPi is the X ROP relative change rate, X ROPi is the ROP correction value of the i-th point depth, is the average value of X ROP from the i-th point depth to the i-5th point depth; The WOB change rate is α WOB The WOB change rate is calculated according to the following formula: wherein, α WOBi is the WOB relative change rate at the i th point depth, WOB i is the WOB value at the i th point depth, is the average value of the WOB values between the i th point depth and the i-5 th point depth. The torque correction parameter change rate is α TOR The formula for calculating the torque correction parameter change rate is as follows: wherein α TORi is the X TOR relative change rate, X TORi is the torque (TOR) correction value of the i-th point depth, is the average value of X TOR from the i-th point depth to the i-5th point depth; S7, calculating the total change rate of the engineering parameter; Wherein, the total engineering parameter change rate is calculated, including: The total engineering parameter change rate is β, and the calculation formula of the total engineering parameter change rate is as follows: wherein, β i is the total variation rate of engineering parameters at the depth of the i th point, α ROPi is the X ROP relative variation rate at the depth of the i th point, α WOBi is the WOB relative variation rate at the depth of the i th point, α TORi is the X TOR relative variation rate at the depth of the i th point; S8, establishing a sandstone sticking layer standard based on the drilling time correction parameter, the torque correction parameter, the WOB, the correction parameter variation rate, the WOB variation rate, the torque correction parameter variation rate and the total variation rate of engineering parameters; S9, according to the sandstone sticking layer standard, carrying out target sandstone sticking layer and fast identification of the sandstone layer.

2. The method of claim 1, wherein, The rotary table revolution normalization coefficient and the standpipe pressure normalization coefficient are calculated, including: The rotation number normalization coefficient of the turntable is t RPM The calculation formula of the rotation number normalization coefficient of the turntable is as follows: wherein, t RPMi is the normalized coefficient of the number of revolutions per minute (RPM) of the turntable at the depth of the i-th point, RPM i is the value of the number of revolutions per minute (RPM) of the turntable at the depth of the i-th point, is the average value of the number of revolutions per minute (RPM) of the turntable between the depths of the 1st point and the i-1st point; The riser pressure normalization factor is t SPP The riser pressure normalization factor is calculated as follows: where t SPPi is the riser pressure (SPP) normalization factor for the depth of the ith point, SPP i is the riser pressure (SPP) value for the depth of the ith point, is the average riser pressure (SPP) value between the depths of the 1stpoint and the ith-1 point.

3. The method of claim 1, wherein, The drilling time correction parameter and the torque correction parameter are calculated, including: The drilling time correction parameter is X ROP The calculation formula of the drilling time correction parameter is as follows: wherein X ROPi is the ROP correction value for the depth of the i-th point, ROP i is the ROP value for the depth of the i-th point, t RPMi is the RPM normalization coefficient for the depth of the i-th point; The torque correction parameter is X TOR The calculation formula of the torque correction parameter is as follows: where X TORi is the torque (TOR) correction value for the i-th point depth, TOR i is the torque (TOR) value for the i-th point depth, t SPPi is the riser pressure (SPP) normalization coefficient for the i-th point depth.

4. The method of claim 1, wherein, The target sandstone layer is formed, including: N1, performing steps S1 to S7 of claim 1 to calculate the drilling time correction parameter, the torque correction parameter, the drilling pressure, the correction parameter change rate, the drilling pressure change rate, the torque correction parameter change rate and the total engineering parameter change rate parameter; N2, judging the lithology of the current formation according to the sandstone layer standard in step S8 of claim 1; N3, if the drilling time correction parameter, the torque correction parameter, the drilling pressure, the correction parameter change rate, the drilling pressure change rate, the torque correction parameter change rate and the total engineering parameter change rate parameter of the current formation to be analyzed all meet the sandstone characteristics in the sandstone layer standard of claim 1, it is determined that the current formation is the target sandstone layer, and the geological cycle is performed, and it is judged whether the lithology of the wellhead returned cuttings is the target layer sandstone, if the returned cuttings are the target layer sandstone, drilling coring operation is performed; if the returned cuttings are non-target layer sandstone, drilling is continued, and steps N1 to N2 are re-executed.

5. The method of claim 1, wherein, The sandstone layer is quickly identified, including: L1, performing steps S1 to S7 of claim 1 to calculate the drilling time correction parameter, the torque correction parameter, the drilling pressure, the correction parameter change rate, the drilling pressure change rate, the torque correction parameter change rate and the total engineering parameter change rate parameter; L2, judging the lithology of the current formation according to the sandstone layer standard in step S8 of claim 1; L3, if the drilling time correction parameter, the torque correction parameter, the drilling pressure, the correction parameter change rate, the drilling pressure change rate, the torque correction parameter change rate and the total engineering parameter change rate parameter of the current formation to be analyzed all meet the sandstone or argillaceous sandstone characteristics in the sandstone layer standard of claim 1, the top depth and the bottom depth of the sandstone layer or the argillaceous sandstone layer are recorded, and drilling operation is continued, and steps L1 to L2 are re-executed.

6. A sandstone layering device characterized by, The device includes: The prediction unit is configured to predict the top surface depth of the to-be-stuck sandstone layer of the target well according to the adjacent well real drilled formation data and the target well seismic attribute data. The acquisition unit is configured to acquire the logging engineering parameters at the current depth by reducing the engineering parameter acquisition interval at the set depth of the mudstone above the top surface depth of the to-be-stuck sandstone layer, wherein the logging engineering parameters include drilling time, drilling pressure, torque, standpipe pressure and rotary table revolutions. The calculation unit is configured to calculate the drilling time correction parameter, the torque correction parameter, the drilling time correction parameter change rate, the drilling pressure change rate, the torque correction parameter change rate and the engineering parameter total change rate according to the acquired logging engineering parameters. The calculation of the drilling time correction parameter change rate, the drilling pressure change rate and the torque correction parameter change rate includes: The rate of change of the drilling time correction parameter is α ROP The formula for calculating the rate of change of the drilling time correction parameter is as follows: wherein, α ROPi is the X ROP relative change rate, X ROPi is the ROP correction value of the i-th point depth, is the average value of X ROP between the i-th point depth and the i-5th point depth; The WOB change rate is α WOB The WOB change rate is calculated according to the following formula: wherein, α WOBi is the WOB relative change rate at the i th point depth, WOB i is the WOB value at the i th point depth, is the average of the WOB values between the i th point depth and the i-5 th point depth; The torque correction parameter change rate is α TOR The formula for calculating the torque correction parameter change rate is as follows: wherein α TORi is the X TOR relative change rate, X TORi is the torque (TOR) correction value of the i-th point depth, is the average value of X TOR between the i-th point depth and the i-5th point depth. The calculation of the engineering parameter total change rate includes: The engineering parameter total change rate is β, and the calculation formula of the engineering parameter total change rate is as follows: wherein β i is the total rate of change of the engineering parameter at the depth of the ith point, α ROPi is the relative rate of change of X ROP at the depth of the ith point, α WOBi is the relative rate of change of WOB at the depth of the ith point, α TORi is the relative rate of change of X TOR at the depth of the ith point. The judgment unit is configured to judge the lithology of the formation at the current depth according to the drilling time correction parameter, the torque correction parameter, the drilling pressure, the drilling time correction parameter change rate, the drilling pressure change rate, the torque correction parameter change rate, the engineering parameter total change rate and the sandstone sticking layer standard, and to perform target sandstone sticking layer and rapid identification of sandstone layer according to the sandstone sticking layer standard, wherein the sandstone sticking layer standard is used to represent the qualitative characteristics and quantitative numerical distribution interval of the drilling time correction parameter, the torque correction parameter, the drilling pressure, the drilling time correction parameter change rate, the drilling pressure change rate, the torque correction parameter change rate and the engineering parameter total change rate corresponding to different formations.

7. A computer-readable storage medium, characterized in that, The computer readable storage medium stores at least one program code, and the at least one program code is loaded and executed by the processor to implement the operations performed by the method of any one of claims 1 to 5.

8. An electronic device, comprising: The electronic device includes one or more processors and one or more memories, and the one or more memories store at least one program code, and the at least one program code is loaded and executed by the one or more processors to implement the operations performed by the method of any one of claims 1 to 5.

Images