CN109635448B - Method for calculating critical whirling speed of slewing device - Google Patents

Method for calculating critical whirling speed of slewing device Download PDF

Info

Publication number
CN109635448B
CN109635448B CN201811536616.5A CN201811536616A CN109635448B CN 109635448 B CN109635448 B CN 109635448B CN 201811536616 A CN201811536616 A CN 201811536616A CN 109635448 B CN109635448 B CN 109635448B
Authority
CN
China
Prior art keywords
critical
boundary
steps
following
condition
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811536616.5A
Other languages
Chinese (zh)
Other versions
CN109635448A (en
Inventor
曹银萍
窦益华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Shiyou University
Original Assignee
Xian Shiyou University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Shiyou University filed Critical Xian Shiyou University
Priority to CN201811536616.5A priority Critical patent/CN109635448B/en
Publication of CN109635448A publication Critical patent/CN109635448A/en
Application granted granted Critical
Publication of CN109635448B publication Critical patent/CN109635448B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13Differential equations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD

Abstract

A method for calculating the critical whirling speed of rotator features that the centrifugal force (mu) per unit length is used to calculate the dynamic effect of rotator 1 ω 2 y represents that the critical vortex motion speed under the condition of the simple support-simple support boundary is calculated, and the critical vortex motion speed under the condition of the simple support-simple support boundary is obtained by calculating the shape curve of the device, the bending strain energy, the work done by the external force P, the centrifugal force of unit length and the total energy of the system; critical whirling velocity under the boundary condition of clamped-clamped; the boundary conditions of solid branch-solid branch are: the invention can deduce the critical whirling speed formula according to the constraint conditions at two ends of the revolving body structure, and substitute the structure parameters into the critical rotational speed formula of the sleeve whirling, thus the critical whirling speed can be calculated quantitatively, the rotational speed of the sleeve during drilling can be guided according to the calculation result, and the damage of the whirling phenomenon can be avoided.

Description

Method for calculating critical whirling speed of slewing device
Technical Field
The invention relates to the technical field of oil and gas well exploitation, in particular to a method for calculating a critical whirl velocity of a slewing device.
Background
In petroleum drilling, in order to reduce the non-production time and reduce the loss of drilling fluid, a casing drilling method is proposed. Due to the restriction of the borehole, the casing is quite complex to move during drilling, and the phenomenon of whirl of the drill string occurs similarly to conventional drilling. At present, the research on the whirl phenomenon mostly focuses on the research on the whirl mechanism and law of the drill string. Since the sixty-seven decades of the last century, drill string dynamics research and practice have confirmed that the lower drill string has a common whirl phenomenon in the borehole (meaning that the drill string revolves around the axis of the borehole in addition to rotating); the lower drill string is usually in a reverse whirl (counter-rotation direction) state, while the forward whirl (forward rotation direction) is relatively less; reverse whirl of a lower drill string easily induces fatigue damage of a drilling tool, early damage of a drill bit and the like. For the reasons, the research on the whirl mechanism and the law of the lower drill string is very important at home and abroad. With the development of casing drilling technology, the damage caused by the casing whirl phenomenon is increasingly prominent, and how to accurately calculate the critical whirl speed when the casing rotates, so as to guide the rotating speed selection during construction and avoid the whirl phenomenon when the casing is drilled is more important.
At present, the research on the whirl phenomenon mostly focuses on the research on the whirl mechanism and the whirl rule of the drill string, and the method is also mostly limited to indoor test simulation. On one hand, the method is high in cost and time consumption, on the other hand, the method is too strong in pertinence, only suitable for a specific structure and a specific boundary, and the universality is lacked.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention aims to provide a method for calculating the critical whirling speed of a slewing device. The rotating speed of the casing during drilling can be guided according to the calculation result, and the damage of the vortex phenomenon is avoided.
In order to achieve the purpose, the invention adopts the technical scheme that:
a method for calculating the critical whirling speed of a slewing device comprises the following steps;
based on the d 'Alembert's principle, the dynamic effect of the rotating device is realized by the centrifugal force mu with a unit length 1 ω 2 y represents; the method comprises the following steps of calculating the critical whirling speed under the condition of a simple support-simple support boundary and calculating the critical whirling speed under the condition of a fixed support-fixed support boundary, wherein:
the calculation steps of the critical whirling speed under the condition of the simple support-simple support boundary are as follows:
(1) For simple-simple support, the boundary conditions are as follows:
y = y "= 0 (at x =0 and x = l)
(2) The shape curve of the device is:
Figure BDA0001907043870000021
wherein L is the distance between the two support points, and the differential is taken from the formula (1) to obtain: .
Figure BDA0001907043870000022
Figure BDA0001907043870000023
(3) Bending strain energy U b Comprises the following steps:
Figure BDA0001907043870000024
(4) Work W done by external force P 1 Comprises the following steps:
Figure BDA0001907043870000025
(5) The centrifugal force per unit length is:
dT=μ 1 dxω 2 y (6)
work W done by centrifugal force 2 Comprises the following steps:
Figure BDA0001907043870000031
(6) The total energy of the system is as follows:
Figure BDA0001907043870000032
(7) When the total potential energy is minimum at the critical vortex speed, the first order differential of the total energy of the system should be 0, that is:
Figure BDA0001907043870000033
further obtaining:
Figure BDA0001907043870000034
after finishing to obtain
Figure BDA0001907043870000035
Wherein the content of the first and second substances,
Figure BDA0001907043870000036
is euler buckling load at the simple-simple boundary,
Figure BDA0001907043870000037
then
Figure BDA0001907043870000038
Namely the critical whirling speed under the condition of a simple support-simple support boundary;
the calculation steps of the critical whirling speed under the boundary condition of the fixed branch and the fixed branch are as follows:
(1) The boundary conditions of solid branch-solid branch are:
y = y' =0 (at x =0 and x = l)
(2) The shape curve of the device is:
Figure BDA0001907043870000039
(3) Bending strain energy U b Comprises the following steps:
Figure BDA0001907043870000041
(4) Work W done by external force P 1 Comprises the following steps:
Figure BDA0001907043870000042
(5) Work W done by centrifugal force 2 Comprises the following steps:
Figure BDA0001907043870000043
and (3) obtaining the following conditions after sorting according to the analogy simple support-simple support boundary conditions:
Figure BDA0001907043870000044
wherein the content of the first and second substances,
Figure BDA0001907043870000045
is euler bending load of clamped-clamped
Figure BDA0001907043870000046
I.e. the critical whirling speed in the case of a solidus-solidus boundary.
The invention has the beneficial effects that:
the invention can make up the defects of the experimental method adopted by the conventional whirling law research institute, and has the characteristics of time saving, cost saving and strong applicability. According to the constraint conditions at two ends of the revolving body structure, the actual structure parameters are substituted into the critical rotational speed formula of the sleeve whirl, so that the critical whirl speed of the sleeve during rotation can be accurately calculated, the rotational speed of the sleeve during drilling is further guided, and the damage of the whirl phenomenon is avoided.
Drawings
Fig. 1 is a schematic view of a rotating device rotating at an angular velocity ω under a pressure P.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings.
A method for calculating the critical whirling speed of a slewing device comprises the following steps;
based on the d 'Alembert's principle, the dynamics problem can be converted into dynamics problem by using the centrifugal force mu with unit length for the dynamics effect of the rotating device 1 ω 2 y, from which the schematic shown in fig. 1 can be built to describe this problem.
Critical whirl velocity under simple-simple boundary condition
(1) For simple-simple support, the boundary conditions are as follows:
y = y "= 0 (at x =0 and x = l)
(2) The shape curve of the device is:
Figure BDA0001907043870000051
wherein L is the distance between the two support points, and the differential is taken from the formula (1) to obtain:
Figure BDA0001907043870000052
Figure BDA0001907043870000053
(3) Bending strain energy U b Comprises the following steps:
Figure BDA0001907043870000054
(4) Work W done by external force P 1 Comprises the following steps:
Figure BDA0001907043870000055
(5) The centrifugal force per unit length is:
dT=μ 1 dxω 2 y (6)
work W done by centrifugal force 2 Comprises the following steps:
Figure BDA0001907043870000056
(6) The total energy of the system is as follows:
Figure BDA0001907043870000057
(7) When the total potential energy is minimum at the critical vortex speed, the first order differential of the total energy of the system should be 0, that is:
Figure BDA0001907043870000061
further obtaining:
Figure BDA0001907043870000062
after finishing to obtain
Figure BDA0001907043870000063
Wherein the content of the first and second substances,
Figure BDA0001907043870000064
euler buckling load at the simple-simple boundary,
Figure BDA0001907043870000065
then the
Figure BDA0001907043870000066
Namely the critical whirling speed under the condition of a simple support-simple support boundary;
critical vortex velocity at clamped-clamped boundary condition
(1) The boundary conditions of solid branch-solid branch are:
y = y' =0 (at x =0 and x = l)
(2) The shape curve of the device is:
Figure BDA0001907043870000067
(3) Bending strain energy U b Comprises the following steps:
Figure BDA0001907043870000068
(4) Work W done by external force P 1 Comprises the following steps:
Figure BDA0001907043870000069
(5) Work W done by centrifugal force 2 Comprises the following steps:
Figure BDA00019070438700000610
analogy simple-simple boundary conditions, after sorting, we obtain:
Figure BDA0001907043870000071
wherein the content of the first and second substances,
Figure BDA0001907043870000072
is euler bending load of clamped-clamped
Figure BDA0001907043870000073
I.e. the critical whirling speed in the case of a solidus-solidus boundary.
According to the constraint conditions at two ends of the revolving body structure, the structural parameters are substituted into a casing whirl critical rotating speed formula, so that the critical whirl speed can be quantitatively calculated, the rotating speed of the casing during drilling can be guided according to the calculation result, and the harm of the whirl phenomenon is avoided.

Claims (1)

1. A method for calculating the critical whirling speed of a rotating device is characterized in that,
based on the d 'Alembert's principle, the dynamic effect of the rotating device is realized by the centrifugal force mu with a unit length 1 ω 2 y represents; the method comprises the following steps of calculating the critical whirling speed under the condition of a simple support-simple support boundary and calculating the critical whirling speed under the condition of a fixed support-fixed support boundary, wherein:
the calculation steps of the critical whirling speed under the condition of the simple support-simple support boundary are as follows:
(1) For simple-simple support, the boundary conditions are as follows:
y = y "= 0 (at x =0 and x = l)
(2) The shape curve of the device is:
Figure FDA0001907043860000011
wherein L is the distance between the two support points, and the differential is taken from the formula (1) to obtain: .
Figure FDA0001907043860000012
Figure FDA0001907043860000013
(3) Bending strain energy U b Comprises the following steps:
Figure FDA0001907043860000014
(4) Work W done by external force P 1 Comprises the following steps:
Figure FDA0001907043860000015
(5) The centrifugal force per unit length is:
dT=μ 1 dxω 2 y (6)
work W done by centrifugal force 2 Comprises the following steps:
Figure FDA0001907043860000021
(6) The total energy of the system is as follows:
Figure FDA0001907043860000022
(7) When the total potential energy is minimum at the critical vortex speed, the first order differential of the total energy of the system should be 0, that is:
Figure FDA0001907043860000023
further obtaining:
Figure FDA0001907043860000024
after finishing to obtain
Figure FDA0001907043860000025
Wherein, the first and the second end of the pipe are connected with each other,
Figure FDA0001907043860000026
is euler buckling load at the simple-simple boundary,
Figure FDA0001907043860000027
then the
Figure FDA0001907043860000028
I.e. the critical whirl velocity in the case of a simple-simple boundary.
The calculation steps of the critical whirling speed under the condition of the clamped-clamped boundary are as follows:
(1) The boundary conditions of solid branch-solid branch are:
y = y' =0 (at x =0 and x = l)
(2) The shape curve of the device is:
Figure FDA0001907043860000029
(3) Bending strain energy U b Comprises the following steps:
Figure FDA00019070438600000210
(4) Work W done by external force P 1 Comprises the following steps:
Figure FDA0001907043860000031
(5) Work W done by centrifugal force 2 Comprises the following steps:
Figure FDA0001907043860000032
analogy simple-simple boundary conditions, after sorting, we obtain:
Figure FDA0001907043860000033
wherein the content of the first and second substances,
Figure FDA0001907043860000034
is euler buckling load of clamped-clamped
Figure FDA0001907043860000035
I.e. the critical whirling speed in the case of a solidus-solidus boundary.
CN201811536616.5A 2018-12-14 2018-12-14 Method for calculating critical whirling speed of slewing device Active CN109635448B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811536616.5A CN109635448B (en) 2018-12-14 2018-12-14 Method for calculating critical whirling speed of slewing device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811536616.5A CN109635448B (en) 2018-12-14 2018-12-14 Method for calculating critical whirling speed of slewing device

Publications (2)

Publication Number Publication Date
CN109635448A CN109635448A (en) 2019-04-16
CN109635448B true CN109635448B (en) 2023-02-28

Family

ID=66074224

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811536616.5A Active CN109635448B (en) 2018-12-14 2018-12-14 Method for calculating critical whirling speed of slewing device

Country Status (1)

Country Link
CN (1) CN109635448B (en)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6059533B2 (en) * 2009-11-09 2017-01-11 ユニバーシティ オブ ワシントン センター フォー コマーシャライゼーション Functionalized chromogenic polymer dots and bioconjugates thereof
US9841329B2 (en) * 2015-06-08 2017-12-12 Pioner Engineering Company Strain gage based system and method for failure detection of a fluid film bearing
CN106503399B (en) * 2016-11-19 2017-09-15 东北石油大学 Peupendicular hole hangs the determination method of tubing string Helical Buckling Critical Load
CN107506529B (en) * 2017-08-01 2020-09-18 中国航空工业集团公司西安飞机设计研究所 Method for calculating axial compression stability of composite material reinforced wall plate

Also Published As

Publication number Publication date
CN109635448A (en) 2019-04-16

Similar Documents

Publication Publication Date Title
CN106761382B (en) A kind of deep-well coring device and its operating method
CN101260783B (en) Prebending kinetics deviation control and fast drilling method
CN104895517A (en) Worm-driven double-valve hydraulic oscillator
CN108060896A (en) The mechanical underground instrument of vertical drilling is realized using pin cam
CN205638230U (en) Well drilling speed -raising device
CN109635448B (en) Method for calculating critical whirling speed of slewing device
CN108547569B (en) Drill bit speed increaser and drilling device
CN1804622B (en) Boring and coring process for roller compacted concrete
CN203284783U (en) High-pressure grouting plug
CN108678674A (en) Shale gas develops double acting and couples coiled tubing resistance absorber
CN207485369U (en) A kind of extended reach well antifriction subtracts torsion combination drill string
CN102654035A (en) Bottom drilling tool combination structure for coiled tubing drilling or drilling
CN107448192A (en) The actual bottom hole WOB Forecasting Methodology of static pushing type rotary steering drilling tool
CN205778581U (en) Alleviate the drill string assembly of directional well slipping drilling backing pressure
CN204703815U (en) A kind of turbine drives bivalve hydroscillator
CN210977391U (en) Vortex generator
CN204532329U (en) Borehole accuracy precise control device
CN107560999A (en) A kind of full-scale aluminium drill pipe Dynamic Corrosion experimental bench
CN201661250U (en) Drill collar capable of relieving stick-slip phenomenon
CN206513313U (en) A kind of diameter is telescopic to use screw terminal with pipe
CN206458344U (en) The anti-whirlpool supercharging instrument in underground
CN106735376B (en) A kind of drill clamp device
CN205714182U (en) Without magnetic motor
CN101117884B (en) Boring tool drag reduction device
CN2937475Y (en) Landing centering guide of radio inclinometer during drilling

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant