CN109145392A - Structure of downhole assembly Optimal Configuration Method and device - Google Patents
Structure of downhole assembly Optimal Configuration Method and device Download PDFInfo
- Publication number
- CN109145392A CN109145392A CN201810839432.XA CN201810839432A CN109145392A CN 109145392 A CN109145392 A CN 109145392A CN 201810839432 A CN201810839432 A CN 201810839432A CN 109145392 A CN109145392 A CN 109145392A
- Authority
- CN
- China
- Prior art keywords
- state vector
- drill assembly
- state
- micro unit
- under
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 239000011159 matrix material Substances 0.000 claims abstract description 45
- 238000012546 transfer Methods 0.000 claims abstract description 44
- 238000006073 displacement reaction Methods 0.000 claims abstract description 30
- 230000003068 static effect Effects 0.000 claims description 29
- 230000005284 excitation Effects 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000010008 shearing Methods 0.000 claims description 6
- 241001074085 Scophthalmus aquosus Species 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 5
- 230000005484 gravity Effects 0.000 claims description 5
- 230000010355 oscillation Effects 0.000 abstract description 7
- 238000005553 drilling Methods 0.000 description 20
- 238000010586 diagram Methods 0.000 description 16
- 238000005457 optimization Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000011158 quantitative evaluation Methods 0.000 description 3
- 238000000418 atomic force spectrum Methods 0.000 description 2
- 238000009472 formulation Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000000172 allergic effect Effects 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000037396 body weight Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Earth Drilling (AREA)
Abstract
The present invention provides a kind of structure of downhole assembly Optimal Configuration Method and devices, comprising: is based on Euler's Bernoulli's theory, establishes stress balance equation of the drill assembly micro unit under at least one status condition;By solving each stress balance equation, state vector transfer matrix of drill assembly micro unit under the conditions of corresponding state is obtained;According to state vector at drill bit and each state vector transfer matrix, state vector of drill assembly micro unit under the conditions of corresponding state is calculated;The element of state vector includes one or more of lateral displacement, deviation angle, moment of flexure and shear stress at the drill bit;Total state vector of the drill assembly micro unit is calculated according to the state vector under all status conditions, and distributes the structure of the affiliated drill assembly of drill assembly micro unit rationally according to the element of total state vector.The present invention can reduce the risk that oscillation crosswise occurs for drill assembly.
Description
Technical field
The present invention relates to exploration engineering field more particularly to a kind of structure of downhole assembly Optimal Configuration Methods and device.
Background technique
In drilling process, when drill string is rotated around own axes in well with certain revolving speed in the direction of the clock, pressure,
Under the synergy of pulling force, centrifugal force and torque, centralizer, drilling rod connecting are first-class will irregularly to jump out of the borehole wall in small annular space,
The frictional force that this way of contact generates can induce BHA (drill assembly) oscillation crosswise of certain displacement amplitude, i.e., more fulcrums are certainly
Swash lateral vibration.It can also revolve around hole axis while its main performance characteristic is BHA rotation, referred to as whirling motion (whirl);
With the increase of lateral vibration frequency, BHA can even be rotated around the axis of well counterclockwise, and reversion whirling motion (back occurs
Whirl), as shown in Fig. 1 a to Fig. 1 c, wherein Fig. 1 a be simulate drill assembly rotation schematic diagram, Fig. 1 b be by clockwise about
The dynamic schematic diagram of borehole axis line vortex, Fig. 1 c are by the schematic diagram moved counterclockwise around borehole axis line vortex.The high frequency width that return motion generates is handed over
Allergic effect Reeb transmits to well head along entire BHA and gradually decays, therefore earth's surface is difficult to discover, but each Component units meeting to BHA
Generate fatal destruction.
Summary of the invention
The present invention provides a kind of structure of downhole assembly Optimal Configuration Method and device, and transverse-vibration occurs to reduce drill assembly
Dynamic risk.
The embodiment of the present invention provides a kind of structure of downhole assembly Optimal Configuration Method, comprising: it is based on Euler's Bernoulli's theory,
Establish stress balance equation of the drill assembly micro unit under at least one status condition;By solving each stress balance side
Journey obtains state vector transfer matrix of drill assembly micro unit under the conditions of corresponding state;It is sweared according to state at drill bit
State arrow of drill assembly micro unit under the conditions of corresponding state is calculated in amount and each state vector transfer matrix
Amount;The element of state vector includes one or more of lateral displacement, deviation angle, moment of flexure and shear stress at the drill bit;
Total state vector of the drill assembly micro unit is calculated according to the state vector under all status conditions, and
Distribute the structure of the affiliated drill assembly of drill assembly micro unit rationally according to the element of total state vector.
The embodiment of the present invention also provides a kind of structure of downhole assembly and distributes device rationally, comprising: stress balance establishing equation
Unit is used for: being based on Euler's Bernoulli's theory, is established stress balance of the drill assembly micro unit under at least one status condition
Equation;State vector transfer matrix generation unit, is used for: by solving each stress balance equation, obtaining the drill set
Close state vector transfer matrix of micro unit under the conditions of corresponding state;State vector generation unit, is used for: according to shape at drill bit
State vector and each state vector transfer matrix, are calculated shape of drill assembly micro unit under the conditions of corresponding state
State vector;The element of state vector includes one or more in lateral displacement, deviation angle, moment of flexure and shear stress at the drill bit
It is a;Drill assembly optimizes unit, is used for: the drilling tool is calculated according to the state vector under all status conditions
Total state vector of micro unit is combined, and the drill assembly micro unit institute is distributed rationally according to the element of total state vector
Belong to the structure of drill assembly.
The embodiment of the present invention also provides a kind of computer readable storage medium, is stored thereon with computer program, the program
The step of above-described embodiment the method is realized when being executed by processor.
The embodiment of the present invention also provides a kind of computer equipment, including memory, processor and storage are on a memory simultaneously
The computer program that can be run on a processor, the processor realize above-described embodiment the method when executing described program
Step.
Structure of downhole assembly Optimal Configuration Method, the structure of downhole assembly of the embodiment of the present invention are distributed device rationally, are calculated
Machine readable storage medium storing program for executing and computer equipment establish stress balance equation by Euler's Bernoulli's theory, and flat according to the stress
Weighing apparatus equation obtains total state vector transfer matrix of drill assembly micro unit, can bored by total state vector transfer matrix
Force analysis is carried out to the micro unit distribution of the different configuration of drill assembly under the conditions of external excitation before well, and simulates calculating infinitesimal
Lateral displacement, deviation angle, moment of flexure and the shear stress of body, being capable of quantitative evaluation difference configuration drill assembly generation transverse-vibration with this
Dynamic risk, optimizes the configuration of drill assembly, is damaged to improve efficiency of breaking rock and reduce drilling tool fatigue.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
It obtains other drawings based on these drawings.In the accompanying drawings:
Fig. 1 a is the schematic diagram for simulating drill assembly rotation;
Fig. 1 b is by the schematic diagram moved clockwise about borehole axis line vortex;
Fig. 1 c is by the schematic diagram moved counterclockwise around borehole axis line vortex;
Fig. 2 is the flow diagram of the structure of downhole assembly Optimal Configuration Method of one embodiment of the invention;
Fig. 3 is the Calculation results curve graph of the different configuration BHA structures obtained using the method for the present invention;
Fig. 4 is that the structure of downhole assembly of one embodiment of the invention distributes the structural schematic diagram of device rationally.
Specific embodiment
Understand in order to make the object, technical scheme and advantages of the embodiment of the invention clearer, with reference to the accompanying drawing to this hair
Bright embodiment is described in further details.Here, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, but simultaneously
It is not as a limitation of the invention.
Fig. 2 is the flow diagram of the structure of downhole assembly Optimal Configuration Method of one embodiment of the invention.As shown in Fig. 2,
The structure of downhole assembly Optimal Configuration Method of the present embodiment, it may include:
Step S110: it is based on Euler's Bernoulli's theory, establishes drill assembly micro unit under at least one status condition
Stress balance equation.
The drill assembly micro unit can be a micro unit at any point in drill assembly, can be with drill bit at
Micro unit is as first micro unit.At least one status condition may include at least one in static state and dynamical state.
Stress balance equation under the conditions of different conditions can be different.Stress balance equation under various status conditions may include drill assembly
Parameter, for example, the axis that the moment of inertia of the Young's modulus of elasticity of drill assembly rod piece, drill assembly, the drill assembly body of rod are subject to
To stress etc..Stress balance equation under various status conditions may include the weight of drill assembly micro unit.Various status bars
Stress balance equation under part may include offset variable and drill assembly infinitesimal of the drill assembly micro unit relative to borehole axis line
Offset variable of the body relative to drill bit.Stress balance equation under the conditions of static state may includeIndicate drill assembly infinitesimal
Inclination angle of the body relative to borehole axis line.Stress balance equation under the conditions of dynamical state may include external excitation frequency.
Step S120: by solving each stress balance equation, the drill assembly micro unit is obtained in corresponding state
Under the conditions of state vector transfer matrix.
By solving the available drill assembly micro unit of stress balance equation under the conditions of static state
State vector transfer matrix under the conditions of static state;By solving stress balance equation under the conditions of dynamical state
State vector transfer matrix of available drill assembly micro unit under the conditions of dynamical state.It, can in specific embodiment
First to seek the characteristic root of stress balance equation, then according to the vector form of state vector at drill bit, stress balance equation and
The expression formula for the state vector that state vector indicates at drill assembly micro unit drill bit releases state vector transfer matrix.
Step S130: according to state vector at drill bit and each state vector transfer matrix, the drilling tool is calculated
Combine state vector of micro unit under the conditions of corresponding state;At the drill bit element of state vector include lateral displacement, partially
Move one or more of angle, moment of flexure and shear stress.
It can be calculated according to state vector at drill bit and state vector transfer matrix under the conditions of static state
To state vector of drill assembly micro unit under the conditions of static state;According to state vector at the drill bit and it is described
The drill assembly micro unit can be calculated in dynamical state condition in state vector transfer matrix under the conditions of dynamical state
Under state vector.In embodiment, phase can be arrived by the way that each state vector transfer matrix is multiplied with state vector at drill bit
The state vector for the drill assembly micro unit answered.State vector can be by lateral displacement, deviation angle, moment of flexure, cut at the drill bit
The vector that the five-element one that shearing stress and constant are constituted arrange.
Step S140: the drill assembly infinitesimal is calculated according to the state vector under all status conditions
Total state vector of body, and the affiliated drill set of drill assembly micro unit is distributed rationally according to the element of total state vector
The structure of conjunction.
It, can will be under the conditions of static state when status condition includes multiple, such as comprising static state and when dynamical state
State vector and state vector under the conditions of dynamical state carry out linear superposition and obtain the total of the drill assembly micro unit
State vector.The affiliated drill assembly of drill assembly micro unit can be distributed rationally according to the element of total state vector
Structure.When the element of state vector includes lateral displacement, deviation angle, moment of flexure and shear stress at drill bit, total state vector
Element also includes lateral displacement, deviation angle, moment of flexure and shear stress, and lateral displacement, deviation angle, the value of moment of flexure and shear stress are got over
It is small, it can be determined that structure of downhole assembly configuration more optimizes, so by adjusting structure of downhole assembly, for example, centralizer position,
Young's modulus, torque etc. may be implemented optimization and bore so that lateral displacement, deviation angle, the value of moment of flexure and shear stress are as small as possible
Has the configuration of composite structure.
In the present embodiment, stress balance equation is established by Euler's Bernoulli's theory, and obtain according to the stress balance equation
It, can be before drilling well externally by total state vector transfer matrix to total state vector transfer matrix of drill assembly micro unit
The micro unit distribution of different configuration of drill assembly under the conditions of excitation carries out force analysis, and simulates the transverse direction for calculating micro unit
Displacement, deviation angle, moment of flexure and shear stress are capable of the wind of quantitative evaluation difference configuration drill assembly generation oscillation crosswise with this
Danger, optimizes the configuration of drill assembly, is damaged to improve efficiency of breaking rock and reduce drilling tool fatigue.
In some embodiments, at least one status condition can include: static state and dynamical state.It is dynamic by considering
State status condition helps to obtain the state vector transfer matrix under the conditions of external excitation.By considering static state and dynamic simultaneously
Two kinds of state can more meet the actual conditions of drill assembly in use.
In some embodiments, above-mentioned steps S140, that is, calculated according to the state vector under all status conditions
Obtain total state vector of the drill assembly micro unit, it may include: according under the conditions of static state the state vector and
The state vector under the conditions of dynamical state carries out linear superposition and obtains total state vector of the drill assembly micro unit.Line
Property superposition when, the coefficient of the state vector of the state vector and dynamical state of static state can according to need setting.
In some embodiments, the formula for carrying out linear superposition can are as follows:
μ (x, t)=μs(x)+μd(x) sin (ω t),
Wherein, μ (x, t) indicates total state vector, and x indicates that drill assembly infinitesimal body position and t indicate time, μs(x) table
Show the state vector under the conditions of static state, μd(x) state vector under the conditions of dynamical state is indicated, ω indicates external excitation frequency
Rate.
State vector transfer matrix T in some embodiments, under the conditions of static stateSIt can are as follows:
Wherein, L indicates the length of drill assembly micro unit, λ1Indicate the spy of the stress balance equation under the conditions of static state
Root is levied,P indicates that pulling force suffered by drill assembly, W indicate the weight of drill assembly micro unit,Indicate drill assembly
The angle of micro unit and borehole axis, E indicate that the Young's modulus of elasticity of drill assembly rod piece, I indicate the moment of inertia of drill assembly.
State vector transfer matrix T in some embodiments, under the conditions of dynamical statedIt can are as follows:
Wherein, λ2The characteristic root of the stress balance equation under the conditions of dynamical state, characteristic root λ are indicated with κ2MeetωnIndicate that order is the external excitation frequency of n, n, which can use 1, g, indicates acceleration of gravity, special
Root κ is levied to meet
In some embodiments, the stress balance equation under the conditions of static state can are as follows:
Wherein, E indicates that the Young's modulus of elasticity of drill assembly rod piece, I indicate the moment of inertia of drill assembly, and y indicates drilling tool
Displacement of the micro unit relative to borehole axis line is combined, x indicates displacement of the drill assembly micro unit relative to drill bit, and p indicates drill set
The axial stress that the body of rod is subject to is closed, W indicates the weight of drill assembly micro unit,Indicate drill assembly micro unit and borehole axis
Angle.
In some embodiments, the stress balance equation under the conditions of dynamical state can are as follows:
Wherein, ω indicates external excitation frequency, and g indicates acceleration of gravity.
In order to configure preferably during Drilling Design to BHA structure before brill, configuration preferred process mainly includes number
According to three typing, BHA state vector analytical calculation and graphical representation aspects.In one embodiment, structure of downhole assembly optimization
Configuration method may include following three key steps:
Step 1 establishes BHA horizontal vibrating movement equation.
By linear decomposition, its state vector can be analyzed to static vector sum motion vector when oscillation crosswise occurs for BHA.Root
According to linear superposition theorem, state vector μ is by static vector μsWith motion vector μdComposition, respectively indicate BHA it is static with drilling two kinds
Motion state.When BHA is by external excitation forced vibration with frequency for ω, the state vector μ at BHA each point, which is represented by, to be appointed
Function when meaning position x and any time t:
μ (x, t)=μs(x)+μd(x)sin(ωt) (1)
Wherein, ω indicates external excitation frequency.
Micro unit state vector may be expressed as:
Step 2 establishes BHA micro unit state vector transfer matrix under stationary state, calculates the static vector μ of BHAs
(x)。
The expression formula of BHA micro unit can under stationary state are as follows:
Wherein, μ0For state vector at drill bit (being represented by the form of formula (2), can be known quantity),For static shape
The state vector of i-th of BHA micro unit, T under statesFor the transfer matrix of BHA micro unit state vector under stationary state.
Under the conditions of stationary state, the stress balance equation of BHA micro unit can are as follows:
Wherein, W is infinitesimal body weight, unit N in above formula;E indicates the Young's modulus of elasticity of rod piece, unit Pa;I is
The moment of inertia, unit m4;P is the body of rod by axial stress, unit N;Y indicates drill assembly micro unit relative to borehole axis line
Displacement, x indicate displacement of the drill assembly micro unit relative to drill bit,Indicate the folder of drill assembly micro unit and borehole axis
Angle.
It is solved according to formula (4) and obtains formulation character root:
Y=eβx, abbreviation obtains:
According to formulation character root, formula (4) general solution are as follows:
Y=ax2+bx+c+deβx+f-βx (5)
Wherein, a, b, c, d, f, β are parameter.
The displacement of drilling tool micro unit, inclination angle, moment of flexure and shear stress equation are obtained according to formula (4):
BHA micro unit transfer matrix T when simultaneous equations (2), (5) and (6) abbreviation obtains staticS:
Wherein, P is pulling force, unit N;L is infinitesimal body length;
Step 3 establishes BHA micro unit state vector transfer matrix under dynamical state, calculates the motion vector μ of BHAd
(x)。
The state vector of BHA micro unit under dynamical stateExpression formula can are as follows:
Wherein, μ0For state vector at drill bit,For the state vector of BHA micro unit under dynamical state, TdFor dynamic shape
The transfer matrix of BHA micro unit state vector under state.
Under the conditions of dynamical state, the stress balance equation of BHA micro unit can are as follows:
Wherein, y indicates displacement of the drill assembly micro unit relative to borehole axis line, x indicate drill assembly micro unit relative to
The displacement of drill bit, ω indicate external excitation frequency, and g indicates acceleration of gravity.
Solve to obtain formula (8), four characteristic roots of quadravalence homogeneous differential equation, respectivelyWherein k and
λ2Meet:
Wherein, ωnThe external excitation frequency for being 1 for order.
According to formula (8), solution obtains general solution:
Y=acosh (kx)+bsinh (kx)+ccos (λ2x)+dsin(λ2x) (10)
Wherein, a, b, c, d are parameter, k and λ2For the characteristic root of formula (8).
The displacement of drilling tool micro unit, inclination angle, moment of flexure and shear stress equation are obtained according to formula (8):
Simultaneous equations (2), BHA micro unit transfer matrix when (10) and (11) abbreviation obtains dynamic are as follows:
Wherein, k and λ2For the characteristic root of formula (8), meet above-mentioned formula (9).
The method of the present embodiment merges newton's equation of motion with Euler's Bernoulli equation, under the conditions of constructing external excitation
BHA force analysis distributed mass model.By state vector transfer matrix, different drill assembly can be configured and be carried out before brill
Simulation calculates, each micro unit state vector of comparative analysis drill assembly, preferably adjusts the position of centralizer, reduces drill assembly hair
The risk of raw oscillation crosswise, it is impaired with great importance to improving efficiency of breaking rock and reducing drilling tool fatigue.
The BHA method for optimizing configuration of the present embodiment can carry out Rationality Assessment to the BHA combination of design before brill, pass through
Lateral displacement, deviation angle, shearing force and the moment of flexure of the different configuration BHA micro units of analysis are calculated, and carries out parameter comparison evaluation BHA
The superiority and inferiority of configuration adjusts structure, can optimize BHA configuration, improve efficiency of breaking rock.This method can be used in Drilling Design process
In to BHA structure optimize, reduce occur underground transverse direction nuisance vibration risk, ensure downhole safety.
Fig. 3 is the Calculation results curve graph of the different configuration BHA structures obtained using the method for the present invention, such as Fig. 3 institute
Show, the different drill assembly A of design configurations, drill assembly B and drill assembly C, test three groups of drill assembly displacement curve,
Tilt curves, M curve and shearing force curve, contrast test result it is found that the displacement curve of drill assembly A, tilt curves,
M curve and the fluctuation of shearing force curve are minimum, so the configuration of drill assembly A is best.
The method for being used for BHA structure configuration optimization during Drilling Design provided in an embodiment of the present invention.Based on frequency domain
Under, have damping single-degree-of-freedom drilling tool oscillation crosswise prediction model, utilize drill assembly data, hole trajectory data, mud
Dynamic project parameter in energy data, geologic data and drilling course is sweared by the state that transfer matrix method describes BHA micro unit
Amount, including four lateral displacement, deviation angle, moment of flexure and shearing force parameters provide the excellent of different BHA structure configurations for technical staff
It is bad.
Based on inventive concept identical with structure of downhole assembly Optimal Configuration Method shown in Fig. 2, the embodiment of the present application is also
It provides a kind of structure of downhole assembly and distributes device rationally, as described in following example.Since structure of downhole assembly optimization is matched
It is similar to structure of downhole assembly Optimal Configuration Method to set the principle that device solves the problems, such as, therefore the structure of downhole assembly is distributed rationally
The implementation of device may refer to the implementation of structure of downhole assembly Optimal Configuration Method, and overlaps will not be repeated.
Fig. 4 is that the structure of downhole assembly of one embodiment of the invention distributes the structural schematic diagram of device rationally.As shown in figure 4,
The structure of downhole assembly of the present embodiment distributes device rationally, it may include: stress balance establishing equation unit 210, state vector pass
Matrix generation unit 220, state vector generation unit 230 and drill assembly optimization unit 240 are passed, above-mentioned each unit sequence connects
It connects.
Stress balance establishing equation unit 210, is used for: being based on Euler's Bernoulli's theory, establishes drill assembly micro unit and exist
Stress balance equation under at least one status condition;
State vector transfer matrix generation unit 220, is used for: by solving each stress balance equation, obtaining described
State vector transfer matrix of drill assembly micro unit under the conditions of corresponding state;
State vector generation unit 230, is used for: according to state vector at drill bit and each state vector transfer matrix,
State vector of drill assembly micro unit under the conditions of corresponding state is calculated;The element of state vector at the drill bit
Including one or more of lateral displacement, deviation angle, moment of flexure and shear stress;
Drill assembly optimizes unit 240, is used for: being calculated according to the state vector under all status conditions
Total state vector of the drill assembly micro unit, and the drill assembly is distributed rationally according to the element of total state vector
The structure of the affiliated drill assembly of micro unit.
In some embodiments, at least one status condition includes: static state and dynamical state.
In some embodiments, it includes data inputting module, BHA state vector that structure of downhole assembly, which distributes device rationally mainly,
Analytical calculation module and graphical representation module.BHA configuration optimization device can carry out different configuration of BHA structure before brill
Force analysis and optimization provide visualization tool for BHA optimization in Drilling Design.Merge Euler's Bernoulli equation and newton fortune
Dynamic equation, BHA infinitesimal weight are merged with Bending Deformation in a state vector transfer matrix function, utilize transfer matrix meter
Each micro unit state vector of entire BHA is calculated, and then assesses the superiority and inferiority of different BHA configurations.
The embodiment of the present invention also provides a kind of computer readable storage medium, is stored thereon with computer program, the program
The step of the various embodiments described above the method is realized when being executed by processor.
The embodiment of the present invention also provides a kind of computer equipment, including memory, processor and storage are on a memory simultaneously
The computer program that can be run on a processor, the processor realize the various embodiments described above the method when executing described program
The step of.
In conclusion structure of downhole assembly Optimal Configuration Method, the structure of downhole assembly of the embodiment of the present invention are distributed rationally
Device, computer readable storage medium and computer equipment establish stress balance equation by Euler's Bernoulli's theory, and according to
The stress balance equation obtains total state vector transfer matrix of drill assembly micro unit, passes through total state vector transfer matrix
Can before drilling well to external excitation under the conditions of the micro unit distribution of different configuration of drill assembly carry out force analysis, and simulate
Lateral displacement, deviation angle, moment of flexure and the shear stress of micro unit are calculated, it being capable of quantitative evaluation difference configuration drill assembly hair with this
The risk of raw oscillation crosswise, optimizes the configuration of drill assembly, is damaged to improve efficiency of breaking rock and reduce drilling tool fatigue.
In the description of this specification, reference term " one embodiment ", " specific embodiment ", " some implementations
Example ", " such as ", the description of " example ", " specific example " or " some examples " etc. mean it is described in conjunction with this embodiment or example
Particular features, structures, materials, or characteristics are included at least one embodiment or example of the invention.In the present specification,
Schematic expression of the above terms may not refer to the same embodiment or example.Moreover, the specific features of description, knot
Structure, material or feature can be combined in any suitable manner in any one or more of the embodiments or examples.Each embodiment
Involved in the step of sequence be used to schematically illustrate implementation of the invention, sequence of steps therein is not construed as limiting, can be as needed
It appropriately adjusts.
It should be understood by those skilled in the art that, the embodiment of the present invention can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the present invention
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the present invention, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
Particular embodiments described above has carried out further in detail the purpose of the present invention, technical scheme and beneficial effects
Describe in detail it is bright, it should be understood that the above is only a specific embodiment of the present invention, the guarantor being not intended to limit the present invention
Range is protected, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in this
Within the protection scope of invention.
Claims (10)
1. a kind of structure of downhole assembly Optimal Configuration Method characterized by comprising
Based on Euler's Bernoulli's theory, stress balance equation of the drill assembly micro unit under at least one status condition is established;
By solving each stress balance equation, state arrow of drill assembly micro unit under the conditions of corresponding state is obtained
Measure transfer matrix;
According to state vector at drill bit and each state vector transfer matrix, the drill assembly micro unit is calculated in phase
Answer the state vector under status condition;The element of state vector includes lateral displacement, deviation angle, moment of flexure and shearing at the drill bit
One or more of stress;
Total state arrow of the drill assembly micro unit is calculated according to the state vector under all status conditions
It measures, and distributes the structure of the affiliated drill assembly of drill assembly micro unit rationally according to the element of total state vector.
2. structure of downhole assembly Optimal Configuration Method as described in claim 1, which is characterized in that at least one status bar
Part includes: static state and dynamical state.
3. structure of downhole assembly Optimal Configuration Method as claimed in claim 2, which is characterized in that according to all status bars
Total state vector of the drill assembly micro unit is calculated in the state vector under part, comprising:
It is linearly folded according to the state vector under the conditions of static state and the state vector under the conditions of dynamical state
Add to obtain total state vector of the drill assembly micro unit.
4. structure of downhole assembly Optimal Configuration Method as claimed in claim 3, which is characterized in that carry out the formula of linear superposition
Are as follows:
μ (x, t)=μs(x)+μd(x) sin (ω t),
Wherein, μ (x, t) indicates total state vector, and x indicates that drill assembly infinitesimal body position and t indicate time, μs(x) indicate static
State vector under status condition, μd(x) state vector under the conditions of dynamical state is indicated, ω indicates external excitation frequency.
5. structure of downhole assembly Optimal Configuration Method as claimed in claim 2, which is characterized in that
State vector transfer matrix T under the conditions of static stateSAre as follows:
Wherein, L indicates the length of drill assembly micro unit, λ1Indicate the characteristic root of the stress balance equation under the conditions of static state,P indicates that pulling force suffered by drill assembly, W indicate the weight of drill assembly micro unit,Indicate drill assembly infinitesimal
The angle of body and borehole axis, E indicate that the Young's modulus of elasticity of drill assembly rod piece, I indicate the moment of inertia of drill assembly;
State vector transfer matrix T under the conditions of dynamical statedAre as follows:
Wherein, λ2The characteristic root of the stress balance equation under the conditions of dynamical state, characteristic root λ are indicated with κ2MeetωnIndicate external excitation frequency, g indicates that acceleration of gravity, characteristic root κ meet
6. structure of downhole assembly Optimal Configuration Method as claimed in claim 2, which is characterized in that
Stress balance equation under the conditions of static state are as follows:
Wherein, E indicates that the Young's modulus of elasticity of drill assembly rod piece, I indicate the moment of inertia of drill assembly, and y indicates drill assembly
Displacement of the micro unit relative to borehole axis line, x indicate displacement of the drill assembly micro unit relative to drill bit, and p indicates drill assembly bar
The axial stress that body is subject to, W indicate the weight of drill assembly micro unit,Indicate the folder of drill assembly micro unit and borehole axis
Angle;
Stress balance equation under the conditions of dynamical state are as follows:
Wherein, ω indicates external excitation frequency, and g indicates acceleration of gravity.
7. a kind of structure of downhole assembly distributes device rationally characterized by comprising
Stress balance establishing equation unit, is used for: being based on Euler's Bernoulli's theory, establishes drill assembly micro unit at least one
Stress balance equation under status condition;
State vector transfer matrix generation unit, is used for: by solving each stress balance equation, obtaining the drill assembly
State vector transfer matrix of micro unit under the conditions of corresponding state;
State vector generation unit, is used for: according to state vector at drill bit and each state vector transfer matrix, being calculated
State vector of drill assembly micro unit under the conditions of corresponding state;The element of state vector includes laterally at the drill bit
One or more of displacement, deviation angle, moment of flexure and shear stress;
Drill assembly optimizes unit, is used for: the brill is calculated according to the state vector under all status conditions
Total state vector of tool combination micro unit, and the drill assembly micro unit is distributed rationally according to the element of total state vector
The structure of affiliated drill assembly.
8. structure of downhole assembly as claimed in claim 7 distributes device rationally, which is characterized in that at least one status bar
Part includes: static state and dynamical state.
9. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the program is held by processor
The step of claim 1 to 6 the method is realized when row.
10. a kind of computer equipment including memory, processor and stores the meter that can be run on a memory and on a processor
Calculation machine program, which is characterized in that the step of processor realizes claim 1 to 6 the method when executing described program.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810839432.XA CN109145392A (en) | 2018-07-27 | 2018-07-27 | Structure of downhole assembly Optimal Configuration Method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810839432.XA CN109145392A (en) | 2018-07-27 | 2018-07-27 | Structure of downhole assembly Optimal Configuration Method and device |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109145392A true CN109145392A (en) | 2019-01-04 |
Family
ID=64799263
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810839432.XA Pending CN109145392A (en) | 2018-07-27 | 2018-07-27 | Structure of downhole assembly Optimal Configuration Method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109145392A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110866334A (en) * | 2019-11-05 | 2020-03-06 | 中国海洋石油集团有限公司 | Drilling tool state evaluation method and system based on digital twinning |
CN111460684A (en) * | 2020-04-21 | 2020-07-28 | 长江大学 | Mechanical analysis method of bottom hole assembly with bend angle |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100032165A1 (en) * | 2007-02-02 | 2010-02-11 | Bailey Jeffrey R | Modeling And Designing of Well Drilling System That Accounts For Vibrations |
US20110214878A1 (en) * | 2008-11-21 | 2011-09-08 | Bailey Jeffrey R | Methods and Systems For Modeling, Designing, and Conducting Drilling Operations That Consider Vibrations |
US20160186551A1 (en) * | 2013-10-11 | 2016-06-30 | Halliburton Energy Services Inc. | Optimal control of a drill path using path smoothing |
WO2018000211A1 (en) * | 2016-06-29 | 2018-01-04 | Schlumberger Technology Corporation | Drilling energy calculation based on transient dynamics simulation and its application to drilling optimization |
-
2018
- 2018-07-27 CN CN201810839432.XA patent/CN109145392A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100032165A1 (en) * | 2007-02-02 | 2010-02-11 | Bailey Jeffrey R | Modeling And Designing of Well Drilling System That Accounts For Vibrations |
US20110214878A1 (en) * | 2008-11-21 | 2011-09-08 | Bailey Jeffrey R | Methods and Systems For Modeling, Designing, and Conducting Drilling Operations That Consider Vibrations |
US20160186551A1 (en) * | 2013-10-11 | 2016-06-30 | Halliburton Energy Services Inc. | Optimal control of a drill path using path smoothing |
WO2018000211A1 (en) * | 2016-06-29 | 2018-01-04 | Schlumberger Technology Corporation | Drilling energy calculation based on transient dynamics simulation and its application to drilling optimization |
Non-Patent Citations (3)
Title |
---|
BAILEY J.R.等: "Development and Application of a BHA Vibrations Model" * |
CUI MENG等: "Development and Application of a BHA Vibrations Analysis Model" * |
韩加庚等: "基于地表数据的井下振动识别与控制" * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110866334A (en) * | 2019-11-05 | 2020-03-06 | 中国海洋石油集团有限公司 | Drilling tool state evaluation method and system based on digital twinning |
CN111460684A (en) * | 2020-04-21 | 2020-07-28 | 长江大学 | Mechanical analysis method of bottom hole assembly with bend angle |
CN111460684B (en) * | 2020-04-21 | 2022-04-12 | 长江大学 | Mechanical analysis method of bottom hole assembly with bend angle |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Petyt | Introduction to finite element vibration analysis | |
Gulyayev et al. | Free vibrations of drill strings in hyper deep vertical bore-wells | |
Esen | A new FEM procedure for transverse and longitudinal vibration analysis of thin rectangular plates subjected to a variable velocity moving load along an arbitrary trajectory | |
EA015308B1 (en) | Modeling and designing of well drilling system that accounts for vibrations | |
Servin et al. | Hybrid, multiresolution wires with massless frictional contacts | |
CN108820264A (en) | It is a kind of for removing the rope system towing system and method for space junk | |
Ondra et al. | Free vibration and stability analysis of a cantilever beam axially loaded by an intermittently attached tendon | |
Weiping et al. | Effect of prebent deflection on lateral vibration of stabilized drill collars | |
CN109145392A (en) | Structure of downhole assembly Optimal Configuration Method and device | |
Heydari et al. | Forced vibration analysis of a Timoshenko cracked beam using a continuous model for the crack | |
Peiris et al. | Soil–pile interaction of pile embedded in deep-layered marine sediment under seismic excitation | |
Zhao et al. | Characteristics of steady vibration in a rotating hub–beam system | |
Khajiyeva et al. | The effect of gas and fluid flows on nonlinear lateral vibrations of rotating drill strings | |
Silva et al. | Nonlinear flexural response of a slender cantilever beam of constant thickness and linearly-varying width to a primary resonance excitation | |
Perig et al. | Spherical pendulum small oscillations for slewing crane motion | |
Wang et al. | Frequency equation of flexural vibrating cantilever beam considering the rotary inertial moment of an attached mass | |
CN109165374A (en) | Down-hole equipment oscillation crosswise appraisal procedure and device | |
Hanson et al. | Dynamics modeling of PDC bits | |
Lin et al. | Analysis of rotating nonuniform pretwisted beams with an elastically restrained root and a tip mass | |
Lim et al. | Modal analysis of a rotating multi-packet blade system | |
Jia | On the bending coupled natural frequencies of a spinning, multispan Timoshenko shaft carrying elastic disks | |
Servin et al. | Massless Cable for Real‐time Simulation | |
Ondra et al. | Theoretical and experimental free vibration analysis of a beam-tendon system with an eccentrically placed tendon | |
Sergaliyev et al. | Flat flexural vibration of drill-string with an initial curvature | |
Wei et al. | Nonlinear dynamics of an electrorheological sandwich beam with rotary oscillation |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190104 |