Real-time drilling parameter optimization method and system based on extremum optimizing
Technical Field
The invention relates to the technical field of oil and gas drilling, in particular to a real-time drilling parameter optimization method and system based on extremum optimizing.
Background
With the deep development of oil and gas exploration, the quality of resources is worsened, and unconventional and ultra-deep oil and gas resources become a major discovery matrix for oil and gas in China. The complex drilling and accelerating technology is affected by multiple factors such as complex drilling environment, engineering technology level and the like, the exploration and development difficulty is increased continuously, and the complex drilling and accelerating technology is widely applied. However, the well drilling parameter optimization belongs to systematic engineering, and has the problems of strong artificial dependence of effective cognition of technological and equipment limitations, immediate countermeasures and the like, especially in deep difficult-to-drill stratum, the difficulty of well drilling parameter regulation and control is increased due to the heterogeneity of stratum and the concealment of underground working conditions, if the well drilling parameters are improperly applied, the low-efficiency rock breaking is easily caused, the well drilling period is increased, the underground abnormal vibration is easily caused, the damage of a drill bit is accelerated, and the failure risk of a drill string is increased. Therefore, the optimization of drilling parameters is an important technical means for realizing safe and efficient drilling in complex drilling environments.
On one hand, the traditional drilling parameter optimization method based on the regional drilling rate model depends on regression statistical analysis of regional drilling engineering data, model prediction accuracy is low, and is limited to pre-drilling scheme optimization, so that the correlation between the underground screw output characteristics and drilling parameters under the condition of a composite drilling process is ignored, the drilling parameter optimization in the real drilling process depends on empirical regulation, no basis exists for quantitative regulation and control basis under the influence of geological engineering fusion factors such as working conditions, formation lithology changes and the like, and whether the method is an optimal scheme cannot be judged. On the other hand, the mechanical drilling speed is improved as an optimization target, so that abnormal influences of abnormal drilling conditions such as drill bit stick slip, whirling and balling on drilling efficiency in the real drilling process are ignored, and negative influences such as early damage of the drill bit and damage of a downhole tool are easily caused by applying drilling parameters exceeding a rock breaking unstable threshold. In addition, the current real-time optimization of drilling parameters is in a starting stage, and due to the nonlinear correlation of a large amount of monitoring data in the real-time drilling process, the traditional optimization algorithm of the drilling parameters has the problem of excessive fitting, and the model is easy to fall into a local optimal solution, so that the speed-up and synergy effects of drilling are affected. Accordingly, in view of the above, there is a need for a method and system for optimizing real-time drilling parameters that overcomes the above-mentioned drawbacks and is simple and practical.
Disclosure of Invention
The invention aims to provide a real-time drilling parameter optimization method and system based on extremum optimizing.
In order to achieve the above purpose, the invention is implemented according to the following technical scheme:
the invention comprises the following steps:
S1, acquiring logging engineering parameters under a real-time working condition, wherein the logging engineering parameters comprise drilling pressure, rotating speed of a rotating disc, torque, drilling fluid discharge capacity and drilling fluid density.
S2, acquiring a rock breaking strength baseline in the optimal drilling state of a target well based on a rock Confined Compressive Strength (CCS) calculation model, and enabling the rock breaking strength baseline to be equivalent to theoretical minimum breaking energy (VCS min), wherein the concrete process of acquiring the rock Confined Compressive Strength (CCS) is as follows:
s21, acquiring gamma, density and acoustic logging data of an adjacent well;
S22, acquiring adjacent well Shi Shuju, wherein the method specifically comprises real drilling stratum lithology and stratum sequence data;
S23, obtaining predicted formation lithology of a target well, predicted formation pressure, predicted formation sequence and designed borehole orbit data;
S24, according to the logging information of the adjacent wells, eliminating cycle jumping points through a rock Uniaxial Compressive Strength (UCS) calculation model, and performing smooth pretreatment to construct a single-well longitudinal uniaxial compressive strength logging response curve. The rock Uniaxial Compressive Strength (UCS) model calculation method is that corresponding empirical calculation formulas are respectively selected according to stratum lithology, sandstone, mudstone and carbonate stratum.
S25, mapping the rock uniaxial compressive strength logging response curve of the adjacent well to the corresponding well depth of the target well by adopting a spatial interpolation extrapolation method according to the target well data and the adjacent well Shi Ziliao;
S26, calculating real-time bottom hole drilling fluid Equivalent Circulating Density (ECD) according to the real-time logging engineering data of the target well;
s27, predicting formation pressure according to the target well, and calculating real-time downhole drilling fluid Equivalent Circulating Density (ECD) to obtain effective confining pressure of the well bottom;
and S28, constructing a rock Confined Compressive Strength (CCS) calculation model according to the uniaxial compressive strength and the shaft bottom confining pressure value of the target well, calculating to obtain a longitudinal upper confined rock compressive strength response curve of the target well through the rock confined compressive strength calculation model, and equivalently using the confined rock compressive strength response curve as a rock breaking strength baseline of the target well in the optimal drilling state. The calculation formula of the lateral limit compressive strength of the rock is as follows:
;
Wherein CCS is rock confined compressive strength, MPa, UCS is rock uniaxial compressive strength, MPa, D P is effective confining pressure at the bottom of the well, MPa; Is the internal friction angle of rock
S3, constructing a crushed rock energy model under a composite working condition considering the output loss of the screw, and calculating virtual rock crushing strength (VCS) under the current engineering parameter state according to the crushed rock energy model;
s31, considering the friction resistance between a rotor and a stator and the leakage between a sealing cavity of a screw drilling tool under the actual composite drilling working condition, and the mechanical loss and the hydraulic loss influence of a transmission shaft and a bearing joint, and determining the actual output torque and the pressure difference of the screw according to the volumetric efficiency under the leakage effect and the mechanical efficiency under the torque loss;
s32, correcting the ground weight and torque to the weight and torque received at the drill bit at the bottom of the well based on friction stress analysis of different wells Duan Zuanzhu of the directional well and the well wall and sliding friction analysis of the drill bit;
S33, taking a drill bit and rock interface as an energy input/output interface, comprehensively considering ground energy input and screw energy input and energy loss caused by drill string transmission, friction resistance and downhole vibration, and obtaining a broken rock energy model under a composite working condition considering screw output loss;
The virtual rock breaking strength calculation model is as follows:
;
the VCS virtual rock breaking strength is MPa, the work done by breaking rock in unit time is W, the rock volume is broken in unit time is V, the mechanical efficiency of an E m new drill bit is V, the weight on bit is measured on the WOB ground, the MPa, the friction coefficient of a mu s drill string, the well bottom well inclination angle of gamma b, the degree, the cross section area of the A b drill bit, m 2, the output efficiency of an eta screw, the pressure drop of the delta P m screw, the MPa, the displacement of Q drilling fluid, L/s and the rate of ROP mechanical drilling, m/h.
S34, based on the theory of rock breaking and a mole-Coulomb (Mohr-Coulomb) rock strength breaking mechanism, the rock breaking energy which needs to be overcome when the PDC drill bit cutting teeth bite into the stratum is equivalent to the lateral compressive strength of the rock, so that the actual broken rock energy is equivalent to the virtual rock breaking strength.
And S4, defining rock breaking energy transmission efficiency (E) according to the rock breaking strength base line and the real-time virtual rock breaking strength, judging the rock breaking efficiency and the underground abnormal working condition in real time, assisting in searching an unstable drilling parameter threshold under the composite drilling working condition, and comparing the deviation between the virtual rock breaking strength and the rock breaking strength base line in real time while drilling, and defining the rock breaking energy transmission efficiency to feed back the real-time rock breaking efficiency. The rock breaking energy transfer efficiency (E) is specifically as follows:
;
Wherein E is the rock breaking energy transmission efficiency, dimensionless, VCS min is the theoretical minimum rock breaking energy, MPa
The step S5 includes the steps of:
S5, taking the broken rock strength baseline as a minimum virtual rock breaking strength adjustment target, comprehensively considering the unstable drilling parameter threshold value, and acquiring a recommended drilling parameter combination under the minimum virtual rock breaking strength by adopting a multivariate extremum optimizing algorithm.
S51, taking the drilling pressure and the rotating speed which are obtained in real time as input values;
S52, obtaining gradient of drilling pressure and rotating speed change by adopting different filters, and then integrating variable values to obtain control quantity of an algorithm after synthesis and demodulation;
s53, judging whether the optimal fitness is the optimal fitness and updating the optimal fitness value;
s54, updating input parameters, and repeating the steps S51-S53;
and S55, searching for the optimal drilling parameter combination by continuously adjusting the drilling weight and the drilling speed.
The system for executing the real-time drilling parameter optimization method based on extremum optimizing comprises the steps of collecting on-site logging data, drilling daily report and drilling fluid report in real time, screening logging engineering data to remove noise points, inputting neural network calculation, automatically converting comprehensive logging second data into input parameters of engineering modeling calculation, selecting drilling tool combinations, inputting drilling fluid performance and inclinometry data, and finally realizing engineering calculation modeling and real-time dynamic analysis of the real-time second data through an engineering calculation engine to realize real-time scientific optimization decision of drilling parameters.
The beneficial effects of the invention are as follows:
Compared with the prior art, the real-time drilling parameter optimization method and system based on extremum optimization can realize starting from a geological engineering double-factor fusion analysis method, and based on single-well rock compressive strength calculation and fine geological data, take limiting factors such as drilling tool combination output characteristics, output loss and the like under a drilling working condition into consideration, and obtain virtual rock crushing strength under a composite drilling working condition. The key parameters after denoising are input into a neural network for calculation by means of real-time acquisition and processing of logging engineering parameters in the drilling process, rock breaking efficiency of a drill bit under the well is comprehensively analyzed by drilling data, then a real-time global optimal solution is calculated by a multivariate extremum optimizing algorithm, a virtual rock breaking strength rapid calculation and mechanical drilling speed real-time optimizing model based on a machine learning method and an optimizing algorithm is built, finally the model is embedded into a parameter optimizing system for on-site operation guidance, mechanical energy loss and negative effects caused by drill bit stick-slip vibration and drilling tool buckling vibration are reduced by taking the minimum virtual rock breaking strength as an index, and the maximum cutting efficiency is obtained in a reinforced parameter interval to realize drilling acceleration. The data blind area optimized by the conventional drilling parameters is broken, the rock breaking efficiency in the whole drilling process is ensured to be maximized, and the limit lifting of the mechanical drilling speed is realized. Has stronger application value, and provides support for saving drilling cost and improving economic benefit in drilling engineering.
Drawings
FIG. 1 is a schematic diagram of a real-time drilling parameter optimization method based on extremum seeking in the present invention;
FIG. 2 is a schematic diagram of a process for establishing a broken rock strength baseline in accordance with the present invention;
FIG. 3 is an exemplary diagram of a method for determining rock breaking efficiency and abnormal downhole conditions in real time according to the present invention;
FIG. 4 is a schematic diagram of a method of optimizing minimum virtual rock breaking strength according to the present invention;
FIG. 5 is a schematic flow chart of a minimum virtual rock fracture strength optimizing technique of the present invention;
FIG. 6 is a schematic diagram of the real-time drilling parameter optimization system based on extremum seeking.
Detailed Description
The invention will be further described with reference to the accompanying drawings and specific embodiments, wherein the exemplary embodiments and descriptions of the invention are for purposes of illustration, but are not intended to be limiting.
As shown in FIG. 1, the real-time drilling parameter optimization method based on extremum optimizing comprises the following steps:
S1, acquiring logging engineering parameters under a real-time working condition, wherein the logging engineering parameters comprise drilling pressure, rotating speed of a rotating disc, torque, drilling fluid discharge capacity and drilling fluid density.
As shown in FIG. 2, S2, based on a rock Confined Compressive Strength (CCS) calculation model, acquiring a rock breaking strength baseline of a target well in an optimal drilling state, and equivalent the rock breaking strength baseline as theoretical minimum rock breaking energy (VCS min), wherein the specific process of acquiring the rock Confined Compressive Strength (CCS) is as follows:
s21, acquiring gamma, density and acoustic logging data of an adjacent well;
S22, acquiring adjacent well Shi Shuju, wherein the method specifically comprises real drilling stratum lithology and stratum sequence data;
S23, obtaining predicted formation lithology of a target well, predicted formation pressure, predicted formation sequence and designed borehole orbit data;
S24, according to the logging information of the adjacent wells, eliminating cycle jumping points through a rock Uniaxial Compressive Strength (UCS) calculation model, and performing smooth pretreatment to construct a single-well longitudinal uniaxial compressive strength logging response curve. The rock Uniaxial Compressive Strength (UCS) model calculation method is that corresponding empirical calculation formulas are respectively selected according to stratum lithology, sandstone, mudstone and carbonate stratum.
S25, mapping the rock uniaxial compressive strength logging response curve of the adjacent well to the corresponding well depth of the target well by adopting a spatial interpolation extrapolation method according to the target well data and the adjacent well Shi Ziliao;
S26, calculating real-time bottom hole drilling fluid Equivalent Circulating Density (ECD) according to the real-time logging engineering data of the target well;
S27, predicting formation pressure according to the target well, and calculating real-time downhole drilling fluid Equivalent Circulating Density (ECD) to obtain effective confining pressure of the well bottom;
and S28, constructing a rock Confined Compressive Strength (CCS) calculation model according to the uniaxial compressive strength and the bottom hole confining pressure value of the target well. Calculating to obtain a longitudinal upper lateral rock compressive strength response curve of the target well through a rock lateral compressive strength calculation model, and equivalently using the lateral rock compressive strength response curve as a rock breaking strength baseline of the target well in an optimal drilling state;
the calculation formula of the lateral limit compressive strength of the rock is as follows:
;
Wherein CCS is rock confined compressive strength, MPa, UCS is rock uniaxial compressive strength, MPa, DP is effective confining pressure at the bottom of the well, and MPa; is the internal friction angle of the rock (°).
S3, constructing a crushed rock energy model under a composite working condition considering the output loss of the screw, and calculating virtual rock crushing strength (VCS) under the current engineering parameter state according to the crushed rock energy model;
s31, considering the friction resistance between a rotor and a stator and the leakage between a sealing cavity of a screw drilling tool under the actual composite drilling working condition, and the mechanical loss and the hydraulic loss influence of a transmission shaft and a bearing joint, and determining the actual output torque and the pressure difference of the screw according to the volumetric efficiency under the leakage effect and the mechanical efficiency under the torque loss;
s32, correcting the ground weight and torque to the weight and torque received at the drill bit at the bottom of the well based on friction stress analysis of different wells Duan Zuanzhu of the directional well and the well wall and sliding friction analysis of the drill bit;
S33, taking a drill bit and rock interface as an energy input/output interface, comprehensively considering ground energy input and screw energy input and energy loss caused by drill string transmission, friction resistance and downhole vibration, and obtaining a broken rock energy model under a composite working condition considering screw output loss;
The virtual rock breaking strength calculation model is as follows:
;
The VCS virtual rock breaking strength is MPa, the work done by breaking rock in unit time is W, the rock volume is broken in unit time is V, the mechanical efficiency of a new drill bit is E m, the weight on bit is measured on the ground of WOB, the MPa, the friction coefficient of a mu s drill string, the well bottom well inclination angle of gamma b, the degree, the cross section area of the drill bit A b, m 2, the output efficiency of an eta screw, the pressure drop of the delta Pm screw, the pressure drop of the MPa, the displacement of drilling fluid Q, L/s and the rate of ROP mechanical drilling, m/h.
S34, based on the theory of rock breaking and a mole-Coulomb (Mohr-Coulomb) rock strength breaking mechanism, the rock breaking energy which needs to be overcome when the PDC cutting teeth eat the stratum is equivalent to the lateral limit compressive strength of the rock, so that the actual broken rock energy is equivalent to the virtual rock breaking strength.
And S4, defining rock breaking energy transmission efficiency (E) according to the rock breaking strength base line and the real-time virtual rock breaking strength, judging the rock breaking efficiency and the underground abnormal working condition in real time, assisting in searching an unstable drilling parameter threshold under the composite drilling working condition, and comparing the deviation between the virtual rock breaking strength and the rock breaking strength base line in real time while drilling, and feeding back the real-time rock breaking efficiency by adopting the rock breaking energy transmission efficiency.
In an ideal state, the rock breaking energy transmission efficiency tends to be 1, and the lower the rock breaking energy transmission efficiency value deviates from 1, the lower the rock breaking efficiency is, which indicates that the underground drill bit is in an unstable rock breaking state.
In the method, rock breaking effects are specifically divided into 2 categories, the rock breaking capacity transfer efficiency is 0.8-1.2, the efficient drilling is indicated, and the rock breaking energy transfer efficiency lower than 0.8 indicates the inefficient drilling.
In a specific embodiment, as shown in fig. 3, the efficiency of rock breaking energy transmission is calculated in real time to be lower than 0.8, and the virtual rock breaking strength is in a spike-shaped area and shows a large fluctuation, and meanwhile, the acquired logging engineering parameters show low frequency, high amplitude and periodic fluctuation, which indicate that the threshold value of unstable weight on bit is exceeded at the moment, the weight on bit is required to be reduced step by step, and the specific adjustment amplitude is controlled to be 10% until the efficiency of rock breaking energy transmission is recovered to 0.8-1.2.
In a specific embodiment, the energy transmission efficiency of breaking rock is calculated in real time to be lower than 0.8, if the lithology is unchanged, the breaking strength of the virtual rock is steadily increased, the fluctuation range is about 5%, the condition that the fluctuation range is lower than the lowest threshold value of hydraulic energy at the moment is indicated, and the drill bit Shui Mali is required to be increased until the energy transmission efficiency of breaking rock is restored to 0.8-1.2.
In a specific embodiment, the real-time calculation of the rock breaking energy transfer efficiency is lower than 0.8, if the lithology is unchanged, the virtual rock breaking strength is continuously increased, and no improvement is achieved through drilling parameter adjustment, and the condition that the rock breaking energy transfer efficiency is lower than the threshold of the passivation degree of the drill bit is indicated at the moment, and the drill bit needs to be started to be replaced.
The rock breaking energy transfer efficiency is specifically as follows:
;
wherein E is the rock breaking energy transmission efficiency, dimensionless, and VCS min is the theoretical minimum rock breaking energy, MPa.
And S5, taking the broken rock strength baseline as a minimum virtual rock breaking strength adjustment target, comprehensively considering the unstable drilling parameter threshold, and acquiring a recommended drilling parameter combination under the minimum virtual rock breaking strength by adopting a multivariate extremum optimizing algorithm, as shown in figures 4 and 5.
When the rock breaking capacity transmission efficiency is between 0.8 and 1.2, a multivariate extremum searching algorithm is considered, and the drilling operation is considered as a multi-stage decision process (the number of stages (the optimal drilling parameter combination) is unknown). In the variable extremum searching algorithm, the controllable variables for realizing the optimized drilling are the weight on bit and the rotating speed. The lateral confined compressive strength of the rock is taken as an objective function, the minimum virtual rock breaking strength is realized by taking the combination of the quantized bit pressure and the rotating speed as an objective, and the optimization of drilling parameters is realized.
The algorithm uses a discrete time representation as measurements of the drilling parameters and commands provided to the control system on the rig are performed at regular intervals. It is assumed that the correlation measurements are performed at time intervals of Dt seconds, and that the top drive and automatic drilling machine may receive updated set points for target weight-on-bit and rotational speed every Dt seconds. The current time step is denoted by t and the command for the iteration time step is denoted by the symbol t + Dt.
S51, taking the drilling pressure and the rotating speed which are obtained in real time as input values;
S52, obtaining gradient of drilling pressure and rotating speed change by adopting different filters, and then integrating variable values to obtain control quantity of an algorithm after synthesis and demodulation;
s53, judging whether the optimal fitness is the optimal fitness and updating the optimal fitness value;
s54, updating input parameters, and repeating the steps S51-S53;
and S55, searching for the optimal drilling parameter combination by continuously adjusting the drilling weight and the drilling speed.
As shown in fig. 6, the system for executing the real-time drilling parameter optimization method based on extreme value optimizing is characterized in that through collecting on-site logging data, drilling daily report and drilling fluid report in real time, logging engineering data are screened for noise points and input into a neural network for calculation, comprehensive logging second data are automatically converted into input parameters of engineering modeling calculation, drilling tool combinations are selected, drilling fluid performance and inclinometry data are input, engineering calculation modeling and real-time dynamic analysis of the real-time second data are finally realized through an engineering calculation engine, and real-time scientific optimization decision of the drilling parameters is realized.
The technical scheme of the invention is not limited to the specific embodiment, and all technical modifications made according to the technical scheme of the invention fall within the protection scope of the invention.