A real-time efficiency estimation method for the transmission belt of a beam pumping unit

Abstract

The invention discloses a real-time efficiency estimation method for a walking beam oil pump transmission belt. The method is designed to solve the problem that a traditional belt transmission efficiency calculation method needs to establish a dynamic model of torsional vibration according to the torsional vibration characteristics of a belt transmission device, and has the disadvantages of complex modeling and low data accuracy. The method proposed by the invention simulates the belt loss power, and brings the simulated loss power into the calculation of the dead point line slope and the loading line slope under different losses, and determines the angle between the dead point line and the loading line by the slope. The simulated angle is compared with the theoretical angle to select the simulated angle closest to the theoretical angle, and the corresponding simulated loss is the belt loss power. The belt transmission efficiency at that moment can be obtained according to the belt loss power. The real-time belt transmission efficiency calculated by the method has high accuracy, does not need to establish a complex model, and is simple, saving additional hardware costs.

Description

Technical Field

[0001] The invention relates to the field of oil production engineering, and in particular to a real-time efficiency estimation method for a driving belt of a beam pumping unit. Background Art

[0002] Pumping units are the main force in oilfield production. Compared with other types of pumping units, beam pumping units have the characteristics of strong applicability, reliable operation, and low environmental requirements. They have been widely used in the oil production industry. The power transmission system of the beam pumping unit mainly includes motors, belts and reduction boxes. The motor converts electrical energy into mechanical energy. The generated power is transmitted to the reduction box after being decelerated by the belt drive. The reduction box drives the cranks on both sides to rotate, and then drives the crank-connecting rod mechanism to rotate at a low speed to extract underground oil. The oil well indicator diagram is an important representation of the working status of the oil well. It can fully reflect the operating status of the pumping unit and the production of crude oil. In order to accurately calculate the indicator diagram, the mechanical loss of power on the belt must be considered. The transmission efficiency of the belt in this process affects the transmission mechanical loss in the transmission process. Therefore, we need to provide a method for the real-time efficiency of the transmission belt.

[0003] At present, the dynamic model of the torsional vibration of the belt drive device of the walking beam pump is generally established, the torsional vibration characteristics of the belt drive device are analyzed, and a simulation model of the transient transmission efficiency of the belt drive is established. This method is used to calculate the real-time efficiency of the transmission belt, but it has the disadvantages of complex model establishment and low data accuracy. Summary of the invention

[0004] The present invention proposes a method for estimating the real-time efficiency of a transmission belt of a walking beam oil pump. The method simulates the belt loss power, brings the simulated loss power into the calculation of the dead point line slope and the loading line slope under different losses, determines the angle between the dead point line and the loading line by the slope, compares the simulated angle with the theoretical angle, and selects the simulated angle closest to the theoretical angle. The corresponding simulated loss is the belt loss power. The real-time transmission efficiency of the belt calculated by the method has high accuracy, does not require the establishment of a complex model, is simple, and saves additional hardware costs.

[0005] To achieve the above object, the present invention provides a technical solution: a method for estimating the real-time efficiency of a beam pumping unit transmission belt, comprising the following steps:

[0006] Step S1, according to the gravity of the sucker rod in the oil Calculate the bottom dead point of the theoretical power diagram , the calculation method is as follows:

[0007] ,

[0008] ,

[0009] in, is the displacement when the suspension point of the theoretical power diagram reaches the bottom dead center, is the load when the suspension point of the theoretical power diagram reaches the bottom dead center, is the weight of the sucker rod in the oil.

[0010] Step S2, according to the stroke S, the gravity of the liquid column above the pump , the gravity of the sucker rod in the oil Calculate the top dead point of the theoretical power diagram , which is calculated as follows:

[0011] ,

[0012] ,

[0013] in, is the displacement of the theoretical power diagram suspension point when it reaches the top dead center, S For the stroke, is the load when the suspension point of the theoretical power diagram reaches the top dead center, is the weight of the liquid column above the pump, is the weight of the sucker rod in the oil.

[0014] Step S3, according to the bottom dead point of the theoretical power diagram 、Theoretical power diagram top dead point The slope of the dead point line of the theoretical power diagram can be obtained .

[0015] Step S4, according to the gravity of the liquid column above the pump , pump depth H, elastic modulus E of the sucker rod material, cross-sectional area of ​​the sucker rod , cross-sectional area of ​​the pumping pipe , calculate the loading point of the theoretical work diagram .

[0016] Step S5: according to the bottom dead point of the theoretical power diagram 、Theoretical power diagram loading point The slope of the loading line of the theoretical power diagram can be obtained .

[0017] Step S6, calculating the angle between the loading line of the theoretical power diagram and the dead point line of the theoretical power diagram .

[0018] Step S7, record the initial power emitted from the motor at a certain moment P The simulated belt power loss gradually increases from 0% to n %( ), so the transmission power after belt loss is , nEqual to 0, 1, 2, 3... n , m , n is a constant, n+1 It must include a current moment actual belt transmission power , Belongs to n+1 one of the.

[0019] Step S8, according to the crank angular velocity , Transmitted power after simulated belt loss Calculate the crank torque for n+1 different simulated power diagrams ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0020] Step S9, calculate n+1 simulated power diagram bottom dead points , n+1 simulated power diagram top dead point , according to the simulation diagram bottom dead center , Simulation diagram top dead point Calculate the slope of the dead point line of n+1 simulated power diagrams , ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0021] Step S10, calculating n+1 simulated work diagram loading points , calculate the slope of n+1 simulated power diagram loading lines according to the simulated power diagram bottom dead center and simulated power diagram loading point , ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0022] Step S11, according to the slope of the dead point line of n+1 simulated power diagrams and the slope of the loading line of n+1 simulated power diagram The angle between n+1 simulated power diagram loading lines and the simulated power diagram dead point line can be calculated , ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0023] Step S12: The angle between the n+1 simulated power diagram loading lines and the simulated power diagram dead point line is The angle between the loading line of the theoretical power diagram and the dead point line of the theoretical power diagram In comparison, choose equal or closest of , Belongs to n+1 As one of the angles between the actual load line and the dead point line, The corresponding simulation diagram bottom dead point , Belongs to n+1 One of the following, simulates the top dead center of the power diagram , Belongs to n+1 One of the simulation diagram loading points , Belongs to n+1 As the actual bottom dead point, actual dead point, and actual loading point at the current moment, its corresponding is the actual belt transmission power at the current moment, Belongs to n+1 one of the,( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0024] Step S13, according to the actual belt transmission power at the current moment Calculate and output the belt transmission efficiency at that moment .

[0025] Further, in step S3, the slope of the dead point line of the theoretical power diagram The calculation method is as follows:

[0026] ,

[0027] in, is the slope of the dead point line of the theoretical power diagram, is the load when the suspension point of the theoretical power diagram reaches the top dead center, is the load when the suspension point of the theoretical power diagram reaches the bottom dead center, is the displacement of the theoretical power diagram suspension point when it reaches the top dead center, is the displacement when the suspension point of the theoretical power diagram reaches the bottom dead center, is the weight of the liquid column above the pump, S For stroke.

[0028] Further, in step S4, the theoretical work diagram loading point The calculation method is as follows:

[0029] ,

[0030] ,

[0031] in, is the displacement when the suspension point of the theoretical work diagram reaches the loading point, is the gravity of the liquid column above the pump, H is the depth of the pump, E is the elastic modulus of the sucker rod material, is the cross-sectional area of ​​the sucker rod, is the cross-sectional area of ​​the pumping pipe, is the load when the suspension point of the theoretical power diagram reaches the loading point, It is the load when the suspension point of the theoretical power diagram reaches the top dead center.

[0032] Further, in step S5, the slope of the loading line of the theoretical power diagram The calculation method is as follows:

[0033] ,

[0034] in, is the slope of the loading line of the theoretical power diagram, is the load when the suspension point of the work diagram reaches the loading point, is the load when the suspension point of the theoretical power diagram reaches the bottom dead center, is the displacement when the suspension point of the theoretical work diagram reaches the loading point, It is the displacement when the suspension point of the theoretical power diagram reaches the bottom dead center.

[0035] Further, in step S6, the angle between the loading line of the theoretical power diagram and the dead point line of the theoretical power diagram is The calculation formula is as follows:

[0036] ,

[0037] in, is the angle between the loading line of the theoretical power diagram and the dead point line of the theoretical power diagram, arctan is the inverse tangent function, is the slope of the loading line of the theoretical power diagram, is the slope of the dead point line of the theoretical power diagram.

[0038] Further, in step S8, the crank torque of the n+1 simulated power diagram The calculation formula is as follows:

[0039] ,

[0040] in, To simulate the crank torque of the power diagram, is the transmission power after belt loss, is the crank angular velocity, ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0041] Further, in step S9, n+1 simulated power diagram bottom dead points In and The calculation formula is as follows:

[0042] ,

[0043] ,

[0044] in, To simulate the load when the suspension point of the power diagram reaches the bottom dead center, To simulate the crank torque when the suspension point of the power diagram reaches the bottom dead center, is the weight of the crank and counterweight, sin is the sine function, To simulate the crank angle when the suspension point of the power diagram reaches the bottom dead center, b is the length of the rear wall of the walking beam, To simulate the angle between the beam and the connecting rod when the suspension point of the power diagram reaches the bottom dead center, r is the crank radius, is the length of the forearm of the walking beam, To simulate the angle between the crank and the connecting rod when the suspension point of the power diagram reaches the bottom dead center, To simulate the displacement of the power diagram suspension point when it reaches the bottom dead center, S For the stroke, is the position angle of the walking beam at the highest point, To simulate the position angle of the beam when the suspension point of the power diagram reaches the bottom dead center, is the position angle of the walking beam at the lowest point, ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0045] n+1 simulated power diagram top dead point In and The calculation formula is as follows:

[0046] ,

[0047] ,

[0048] in, To simulate the load when the suspension point of the power diagram reaches the top dead center, To simulate the crank torque when the suspension point of the power diagram reaches the top dead center, is the weight of the crank and counterweight, sin is the sine function, To simulate the crank angle when the suspension point of the power diagram reaches the top dead center, b is the length of the rear wall of the walking beam, To simulate the angle between the beam and the connecting rod when the suspension point of the power diagram reaches the top dead center, r is the crank radius, is the length of the forearm of the walking beam, To simulate the angle between the crank and the connecting rod when the suspension point of the power diagram reaches the top dead center, To simulate the displacement of the power diagram suspension point when it reaches the top dead center, S For the stroke, is the position angle of the walking beam at the highest point, To simulate the position angle of the beam when the suspension point of the power diagram reaches the top dead center, is the position angle of the walking beam at the lowest point, ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0049] Slope of dead point line of analog power diagram , ( ) is calculated as follows:

[0050] ,

[0051] in, is the slope of the dead point line of the simulated power diagram, To simulate the load when the suspension point of the power diagram reaches the top dead center, To simulate the load when the suspension point of the power diagram reaches the bottom dead center, To simulate the displacement of the power diagram suspension point when it reaches the top dead center, To simulate the displacement of the power diagram suspension point when it reaches the bottom dead center, ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0052] Further, in step S10, n+1 simulated work diagram loading points In and The calculation formula is as follows:

[0053] ,

[0054] ,

[0055] in, To simulate the load when the suspension point of the power diagram reaches the loading point, To simulate the crank torque when the suspension point reaches the loading point, is the weight of the crank and counterweight, sin is the sine function, is the circumference of a circle, is the angle between the walking beam and the connecting rod when the suspension point reaches the loading point, is the angle between the crank and the connecting rod when the suspension point reaches the loading point, is the position angle of the walking beam when the suspension point reaches the loading point, is the angle between the line from the center of the oil beam to the center of the crank and the vertical direction, b is the length of the rear wall of the walking beam, r is the crank radius, is the length of the forearm of the walking beam, To simulate the displacement of the suspension point when it reaches the loading point, S For the stroke, is the position angle of the walking beam at the highest point, is the position angle of the walking beam at the lowest point, ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0056] Slope of loading line of simulated power diagram The calculation formula is as follows:

[0057] ,

[0058] in, is the slope of the simulated power diagram loading line, To simulate the load when the suspension point of the power diagram reaches the loading point, To simulate the load when the suspension point of the power diagram reaches the bottom dead center, To simulate the displacement of the suspension point when it reaches the loading point, To simulate the displacement of the power diagram suspension point when it reaches the bottom dead center, ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0059] Furthermore, in step S11, the angle between the n+1 simulated power diagram loading lines and the simulated power diagram dead point line is The calculation formula is as follows:

[0060] ,

[0061] in, is the angle between the simulated power diagram loading line and the simulated power diagram dead point line, arctan is the inverse tangent trigonometric function, is the slope of the simulated power diagram loading line, is the slope of the dead point line of the simulated power diagram, ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0062] Further, in step S13, the belt transmission efficiency The calculation formula is as follows:

[0063] ,

[0064] in, is the belt transmission efficiency, is the actual transferred power and P is the initial power.

[0065] The technical effects of the present invention are:

[0066] 1. No need to build a complex model, the method is simple and accurate;

[0067] 2. Determine the belt transmission efficiency by software algorithm, saving additional hardware costs;

[0068] 3. Compared with the traditional method of calculating the actual power diagram by using fixed belt transmission efficiency, the real-time belt transmission efficiency calculation method using this method is more accurate and reliable. BRIEF DESCRIPTION OF THE DRAWINGS

[0069] Figure 1 A flowchart of the method for estimating the belt real-time efficiency provided by the present invention;

[0070] Figure 2 A schematic diagram of a theoretical work diagram provided by the present invention;

[0071] Figure 3 An example of a belt real-time efficiency estimation method provided by the present invention n +1 schematic diagram of the simulation work diagram;

[0072] Figure 4 The present invention provides n + 1 actual power diagram selected from the simulated power diagram;

[0073] Figure 5A diagram showing the variation of belt efficiency within one day provided by the present invention. DETAILED DESCRIPTION

[0074] The following is only a preferred embodiment of the present invention. The present invention is further described in detail below in conjunction with the accompanying drawings and specific embodiments, so that those skilled in the art can understand the present invention. It should be pointed out that for those of ordinary skill in the art, as long as various changes are within the spirit and scope of the present invention defined and determined by the attached claims without departing from the principle of the present invention, all inventions and creations using the concept of the present invention are protected.

[0075] Figure 1 The process step diagram of the belt real-time efficiency estimation method provided by the present invention is as follows: Figure 1 As shown, a method for estimating the real-time efficiency of a beam pumping unit transmission belt comprises the following steps:

[0076] Step S1, according to the gravity of the sucker rod in the oil Calculate the bottom dead point of the theoretical power diagram , the calculation method is as follows:

[0077] ,

[0078] ,

[0079] in, is the displacement when the suspension point of the theoretical power diagram reaches the bottom dead center, is the load when the suspension point of the theoretical power diagram reaches the bottom dead center, is the weight of the sucker rod in the oil.

[0080] Step S2, according to the stroke S, the gravity of the liquid column above the pump , the gravity of the sucker rod in the oil Calculate the top dead point of the theoretical power diagram , which is calculated as follows:

[0081] ,

[0082] ,

[0083] in, is the displacement of the theoretical power diagram suspension point when it reaches the top dead center, S For the stroke, is the load when the suspension point of the theoretical power diagram reaches the top dead center, is the weight of the liquid column above the pump, is the weight of the sucker rod in the oil.

[0084] Step S3, according to the bottom dead point of the theoretical power diagram 、Theoretical power diagram top dead point The slope of the dead point line of the theoretical power diagram can be obtained ; Among them, the slope of the dead point line of the theoretical power diagram is The calculation method is as follows:

[0085] ,

[0086] in, is the slope of the dead point line of the theoretical power diagram, is the load when the suspension point of the theoretical power diagram reaches the top dead center, is the load when the suspension point of the theoretical power diagram reaches the bottom dead center, is the displacement of the theoretical power diagram suspension point when it reaches the top dead center, is the displacement when the suspension point of the theoretical power diagram reaches the bottom dead center, is the weight of the liquid column above the pump, S For stroke.

[0087] Step S4, according to the gravity of the liquid column above the pump , pump depth H, elastic modulus E of the sucker rod material, cross-sectional area of ​​the sucker rod , cross-sectional area of ​​the pumping pipe , calculate the loading point of the theoretical work diagram ; Among them, the theoretical power diagram loading point The calculation method is as follows:

[0088] ,

[0089] ,

[0090] in, is the displacement when the suspension point of the theoretical work diagram reaches the loading point, is the gravity of the liquid column above the pump, H is the depth of the pump, E is the elastic modulus of the sucker rod material, is the cross-sectional area of ​​the sucker rod, is the cross-sectional area of ​​the pumping pipe, is the load when the suspension point of the theoretical power diagram reaches the loading point, It is the load when the suspension point of the theoretical power diagram reaches the top dead center.

[0091] Step S5: according to the bottom dead point of the theoretical power diagram 、Theoretical power diagram loading point The slope of the loading line of the theoretical power diagram can be obtained ; Among them, the slope of the loading line of the theoretical power diagram is The calculation method is as follows:

[0092] ,

[0093] in, is the slope of the loading line of the theoretical power diagram, is the load when the suspension point of the work diagram reaches the loading point, is the load when the suspension point of the theoretical power diagram reaches the bottom dead center, is the displacement when the suspension point of the theoretical work diagram reaches the loading point, It is the displacement when the suspension point of the theoretical power diagram reaches the bottom dead center.

[0094] Step S6, calculating the angle between the loading line of the theoretical power diagram and the dead point line of the theoretical power diagram ; Among them, the angle between the loading line of the theoretical power diagram and the dead point line of the theoretical power diagram is The calculation formula is as follows:

[0095] ,

[0096] in, is the angle between the loading line of the theoretical power diagram and the dead point line of the theoretical power diagram, arctan is the inverse tangent function, is the slope of the loading line of the theoretical power diagram, is the slope of the dead point line of the theoretical power diagram.

[0097] Step S7, record the initial power emitted from the motor at a certain moment P The simulated belt power loss gradually increases from 0% to n %( ), so the transmission power after belt loss is , n Equal to 0, 1, 2, 3... n , m , n is a constant, n+1 It must include a current moment actual belt transmission power , Belongs to n+1 one of the.

[0098] Step S8, according to the crank angular velocity , Transmitted power after simulated belt loss Calculate the crank torque for n+1 different simulated power diagrams ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant; among them, the crank torque of n+1 simulated power diagram The calculation formula is as follows:

[0099] ,

[0100] in, To simulate the crank torque of the power diagram, is the transmission power after belt loss, is the crank angular velocity, ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0101] Step S9, calculate n+1 simulation power diagram bottom dead points , n+1 simulated power diagram top dead point , according to the simulation diagram bottom dead center , Simulation diagram top dead point Calculate the slope of the dead point line of n+1 simulated power diagrams , ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant; among them, n+1 simulation power diagram bottom dead point In and The calculation formula is as follows:

[0102] ,

[0103] ,

[0104] in, To simulate the load when the suspension point of the power diagram reaches the bottom dead center, To simulate the crank torque when the suspension point of the power diagram reaches the bottom dead center, is the weight of the crank and counterweight, sin is the sine function, To simulate the crank angle when the suspension point of the power diagram reaches the bottom dead center, b is the length of the rear wall of the walking beam, To simulate the angle between the beam and the connecting rod when the suspension point of the power diagram reaches the bottom dead center, r is the crank radius, is the length of the forearm of the walking beam, To simulate the angle between the crank and the connecting rod when the suspension point of the power diagram reaches the bottom dead center, To simulate the displacement of the power diagram suspension point when it reaches the bottom dead center, S For the stroke, is the position angle of the walking beam at the highest point, To simulate the position angle of the beam when the suspension point of the power diagram reaches the bottom dead center, is the position angle of the walking beam at the lowest point, ( ), n Equal to 0, 1, 2, 3...n , m , n is a constant.

[0105] Among them, n+1 simulated power diagram top dead point In and The calculation formula is as follows:

[0106] ,

[0107] ,

[0108] in, To simulate the load when the suspension point of the power diagram reaches the top dead center, To simulate the crank torque when the suspension point of the power diagram reaches the top dead center, is the weight of the crank and counterweight, sin is the sine function, To simulate the crank angle when the suspension point of the power diagram reaches the top dead center, b is the length of the rear wall of the walking beam, To simulate the angle between the beam and the connecting rod when the suspension point of the power diagram reaches the top dead center, r is the crank radius, is the length of the forearm of the walking beam, To simulate the angle between the crank and the connecting rod when the suspension point of the power diagram reaches the top dead center, To simulate the displacement of the power diagram suspension point when it reaches the top dead center, S For the stroke, is the position angle of the walking beam at the highest point, To simulate the position angle of the beam when the suspension point of the power diagram reaches the top dead center, is the position angle of the walking beam at the lowest point, ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0109] Among them, the slope of the dead point line of the simulated power diagram is , ( ) is calculated as follows:

[0110] ,

[0111] in, is the slope of the dead point line of the simulated power diagram, To simulate the load when the suspension point of the power diagram reaches the top dead center, To simulate the load when the suspension point of the power diagram reaches the bottom dead center, To simulate the displacement of the power diagram suspension point when it reaches the top dead center, To simulate the displacement of the power diagram suspension point when it reaches the bottom dead center, ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0112] Step S10, calculating n+1 simulated work diagram loading points , calculate the slope of n+1 simulated power diagram loading lines according to the simulated power diagram bottom dead center and simulated power diagram loading point , ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant; among them, n+1 simulated power diagram loading points In and The calculation formula is as follows:

[0113] ,

[0114] ,

[0115] in, To simulate the load when the suspension point of the power diagram reaches the loading point, To simulate the crank torque when the suspension point reaches the loading point, is the weight of the crank and counterweight, sin is the sine function, is the circumference of a circle, is the angle between the walking beam and the connecting rod when the suspension point reaches the loading point, is the angle between the crank and the connecting rod when the suspension point reaches the loading point, is the position angle of the walking beam when the suspension point reaches the loading point, is the angle between the line from the center of the oil beam to the center of the crank and the vertical direction, b is the length of the rear wall of the walking beam, r is the crank radius, is the length of the forearm of the walking beam, To simulate the displacement of the suspension point when it reaches the loading point, S For the stroke, is the position angle of the walking beam at the highest point, is the position angle of the walking beam at the lowest point, ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0116] Among them, the slope of the simulated power diagram loading line The calculation formula is as follows:

[0117] ,

[0118] in, is the slope of the simulated power diagram loading line, To simulate the load when the suspension point of the power diagram reaches the loading point, To simulate the load when the suspension point of the power diagram reaches the bottom dead center, To simulate the displacement of the suspension point when it reaches the loading point, To simulate the displacement of the power diagram suspension point when it reaches the bottom dead center, ( ), n Equal to 0, 1, 2, 3... n , m , n Step S11, according to the slope of the dead point line of the n+1 simulated power diagram and the slope of the loading line of n+1 simulated power diagram The angle between n+1 simulated power diagram loading lines and the simulated power diagram dead point line can be calculated , ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0119] Among them, the angle between the n+1 simulated power diagram loading line and the simulated power diagram dead point line is The calculation formula is as follows:

[0120] ,

[0121] in, is the angle between the simulated power diagram loading line and the simulated power diagram dead point line, arctan is the inverse tangent trigonometric function, is the slope of the simulated power diagram loading line, is the slope of the dead point line of the simulated power diagram, ( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0122] Step S12: The angle between the n+1 simulated power diagram loading lines and the simulated power diagram dead point line is The angle between the loading line of the theoretical power diagram and the dead point line of the theoretical power diagram In comparison, choose equal or closest of , Belongs to n+1 As one of the angles between the actual load line and the dead point line, The corresponding simulation diagram bottom dead point , Belongs to n+1 One of the following, simulates the top dead center of the power diagram , Belongs to n+1 One of the simulation diagram loading points , Belongs to n+1 As the actual bottom dead point, actual dead point, and actual loading point at the current moment, its corresponding is the actual belt transmission power at the current moment, Belongs to n+1 one of the,( ), n Equal to 0, 1, 2, 3... n , m , n is a constant.

[0123] Step S13, according to the actual belt transmission power at the current moment Calculate and output the belt transmission efficiency at that moment ; Among them, the belt transmission efficiency The calculation formula is as follows:

[0124] ,

[0125] in, is the belt transmission efficiency, is the actual transferred power and P is the initial power.

[0126] Figure 2 The schematic diagram of the theoretical work diagram provided by the present invention is composed of Figure 2 It can be obtained that the theoretical power diagram is a standard parallelogram, with the top dead center C With the loading point B The vertical axis is parallel.

[0127] Figure 3 An example of a belt real-time efficiency estimation method provided by the present invention n +1 simulation diagram, by Figure 3 It can be obtained that n+1 different simulated power diagrams can be obtained by simulating different belt loss powers, among which there is a set of simulated power diagrams that are closest to the theoretical power diagram.

[0128] Figure 4 The present invention provides n+ 1 actual power diagram selected from the simulated power diagram, by Figure 4 Shown from n +1 actual power diagram selected from the simulated power diagram, the belt power loss represented by the power diagram is the actual belt power loss at the current moment.

[0129] Figure 5 The belt efficiency variation diagram provided by the present invention within one day is as follows: Figure 5 As shown, the real-time belt efficiency changes irregularly within a day and is not a fixed value given by the manufacturer. The dynamometer diagram calculated based on the real-time belt efficiency is more accurate and reliable.

[0130] Although the specific implementation of the invention is described in detail in conjunction with the drawings, it should not be understood as limiting the scope of protection of this patent. Within the scope described in the claims, various modifications and variations that can be made by those skilled in the art without creative work still fall within the scope of protection of this patent.

Claims

1. A method for estimating the real-time efficiency of a beam pumping unit drive belt, characterized in that: The following steps are involved: Step S1, according to the gravity of the sucker rod in the oil Calculate the bottom dead point of the theoretical power diagram , the calculation method is as follows: ， ， in, is the displacement when the suspension point of the theoretical power diagram reaches the bottom dead center, is the load when the suspension point of the theoretical power diagram reaches the bottom dead center, is the weight of the sucker rod in the oil; Step S2, according to the stroke S, the gravity of the liquid column above the pump , the gravity of the sucker rod in the oil Calculate the top dead point of the theoretical power diagram , which is calculated as follows: ， ， in, is the displacement of the theoretical power diagram suspension point when it reaches the top dead center, S For the stroke, is the load when the suspension point of the theoretical power diagram reaches the top dead center, is the weight of the liquid column above the pump, is the weight of the sucker rod in the oil; Step S3, according to the bottom dead point of the theoretical power diagram 、Theoretical power diagram top dead point The slope of the dead point line of the theoretical power diagram can be obtained ; Step S4, according to the gravity of the liquid column above the pump , pump depth H, elastic modulus E of the sucker rod material, cross-sectional area of ​​the sucker rod , cross-sectional area of ​​the pumping pipe , calculate the theoretical power diagram loading point , is the displacement when the suspension point of the theoretical work diagram reaches the loading point, is the load when the suspension point of the theoretical power diagram reaches the loading point; Step S5: according to the bottom dead point of the theoretical power diagram 、Theoretical power diagram loading point The slope of the loading line of the theoretical power diagram can be obtained ; Step S6, calculating the angle θ between the loading line of the theoretical power diagram and the dead point line of the theoretical power diagram; Step S7, record the initial power emitted from the motor at a certain moment P The simulated belt power loss gradually increases from 0% to n %, , so the transmission power after belt loss is , , n Equal to 0, 1, 2, 3... n , m , n is a constant, n+1 It must include a current moment actual belt transmission power , Belongs to n+1 one of the; Step S8, according to the crank angular velocity , Transmitted power after simulated belt loss Calculate the crank torque for n+1 different simulated power diagrams , , n Equal to 0, 1, 2, 3... n , m , n is a constant; Step S9, calculate n+1 simulated power diagram bottom dead points , n+1 simulated power diagram top dead point , according to the simulation diagram bottom dead center , Simulation diagram top dead point Calculate the slope of the dead point line of n+1 simulated power diagrams , , n Equal to 0, 1, 2, 3... n , m , n is a constant; Step S10, calculating n+1 simulated work diagram loading points , calculate the slope of n+1 simulated power diagram loading lines according to the simulated power diagram bottom dead center and simulated power diagram loading point , , n Equal to 0, 1, 2, 3... n , m , n is a constant; Step S11, according to the slope of the dead point line of n+1 simulated power diagrams and the slope of the loading line of n+1 simulated power diagram The angle between n+1 simulated power diagram loading lines and the simulated power diagram dead point line can be calculated , , n Equal to 0, 1, 2, 3... n , m , n is a constant; Step S12: The angle between the n+1 simulated power diagram loading lines and the simulated power diagram dead point line is Compared with the angle θ between the loading line of the theoretical power diagram and the dead point line of the theoretical power diagram, select the angle equal to θ or closest to θ. , Belongs to n+1 As one of the angles between the actual load line and the dead point line, The corresponding simulation diagram bottom dead point , Belongs to n+1 One of the following, simulates the top dead center of the power diagram , Belongs to n+1 One of the simulation diagram loading points , Belongs to n+1 As the actual bottom dead point, actual dead point, and actual loading point at the current moment, its corresponding is the actual belt transmission power at the current moment, Belongs to n+1 one of the, , n Equal to 0, 1, 2, 3... n , m , n is a constant; Step S13, according to the actual belt transmission power at the current moment Calculate and output the belt transmission efficiency at that moment .

2. A method for estimating the real-time efficiency of a beam pumping unit drive belt according to claim 1, characterized in that: In step S3, the slope of the dead point line of the theoretical power diagram The calculation method is as follows: ， in, is the slope of the dead point line of the theoretical power diagram, is the load when the suspension point of the theoretical power diagram reaches the top dead center, is the load when the suspension point of the theoretical power diagram reaches the bottom dead center, is the displacement of the theoretical power diagram suspension point when it reaches the top dead center, is the displacement when the suspension point of the theoretical power diagram reaches the bottom dead center, is the weight of the liquid column above the pump, S For stroke.

3. The method for estimating the real-time efficiency of a beam pumping unit drive belt according to claim 1, characterized in that: In step S4, the theoretical power diagram loading point The calculation method is as follows: ， ， in, is the displacement when the suspension point of the theoretical work diagram reaches the loading point, is the gravity of the liquid column above the pump, H is the depth of the pump, E is the elastic modulus of the sucker rod material, is the cross-sectional area of ​​the sucker rod, is the cross-sectional area of ​​the pumping pipe, is the load when the suspension point of the theoretical power diagram reaches the loading point, It is the load when the suspension point of the theoretical power diagram reaches the top dead center.

4. A method for estimating the real-time efficiency of a beam pumping unit drive belt according to claim 1, characterized in that: In step S5, the slope of the loading line of the theoretical power diagram The calculation method is as follows: ， in, is the slope of the loading line of the theoretical power diagram, is the load when the suspension point of the theoretical power diagram reaches the loading point, is the load when the suspension point of the theoretical power diagram reaches the bottom dead center, is the displacement when the suspension point of the theoretical work diagram reaches the loading point, It is the displacement when the suspension point of the theoretical power diagram reaches the bottom dead center.

5. The method for estimating the real-time efficiency of a beam pumping unit drive belt according to claim 1, characterized in that: In step S6, the angle θ between the loading line of the theoretical power diagram and the dead point line of the theoretical power diagram is calculated as follows: ， Where θ is the angle between the loading line of the theoretical power diagram and the dead point line of the theoretical power diagram, arctan is the inverse tangent function, is the slope of the loading line of the theoretical power diagram, is the slope of the dead point line of the theoretical power diagram.

6. A method for estimating the real-time efficiency of a beam pumping unit drive belt according to claim 1, characterized in that: In step S8, the n+1 simulated power diagram crank torque The calculation formula is as follows: ， in, To simulate the crank torque of the power diagram, is the transmission power after belt loss, is the crank angular velocity, , n Equal to 0, 1, 2, 3... n , m , n is a constant.

7. A method for estimating the real-time efficiency of a beam pumping unit drive belt according to claim 1, characterized in that: In step S9, n+1 simulated power diagram bottom dead points In and The calculation formula is as follows: ， ， in, To simulate the load when the suspension point of the power diagram reaches the bottom dead center, To simulate the crank torque when the suspension point of the power diagram reaches the bottom dead center, is the weight of the crank and counterweight, sin is the sine function, To simulate the crank angle when the suspension point of the power diagram reaches the bottom dead center, b is the length of the rear wall of the walking beam, To simulate the angle between the beam and the connecting rod when the suspension point of the power diagram reaches the bottom dead center, r is the crank radius, is the length of the forearm of the walking beam, To simulate the angle between the crank and the connecting rod when the suspension point of the power diagram reaches the bottom dead center, To simulate the displacement of the power diagram suspension point when it reaches the bottom dead center, S For the stroke, is the position angle of the walking beam at the highest point, To simulate the position angle of the beam when the suspension point of the power diagram reaches the bottom dead center, is the position angle of the walking beam at the lowest point, , n Equal to 0, 1, 2, 3... n , m , n is a constant; n+1 simulated power diagram top dead point In and The calculation formula is as follows: ， ， in, To simulate the load when the suspension point of the power diagram reaches the top dead center, To simulate the crank torque when the suspension point of the power diagram reaches the top dead center, is the weight of the crank and counterweight, sin is the sine function, To simulate the crank angle when the suspension point of the power diagram reaches the top dead center, b is the length of the rear wall of the walking beam, To simulate the angle between the beam and the connecting rod when the suspension point of the power diagram reaches the top dead center, r is the crank radius, is the length of the forearm of the walking beam, To simulate the angle between the crank and the connecting rod when the suspension point of the power diagram reaches the top dead center, To simulate the displacement of the power diagram suspension point when it reaches the top dead center, S For the stroke, is the position angle of the walking beam at the highest point, To simulate the position angle of the beam when the suspension point of the power diagram reaches the top dead center, is the position angle of the walking beam at the lowest point, , n Equal to 0, 1, 2, 3... n , m , n is a constant; Slope of dead point line of analog power diagram , The calculation formula is as follows: ， in, is the slope of the dead point line of the simulated power diagram, To simulate the load when the suspension point of the power diagram reaches the top dead center, To simulate the load when the suspension point of the power diagram reaches the bottom dead center, To simulate the displacement of the power diagram suspension point when it reaches the top dead center, To simulate the displacement of the power diagram suspension point when it reaches the bottom dead center, , n Equal to 0, 1, 2, 3... n , m , n is a constant.

8. A method for estimating the real-time efficiency of a beam pumping unit drive belt according to claim 1, characterized in that: In step S10, n+1 simulated work diagram loading points In and The calculation formula is as follows: ， ， in, To simulate the load when the suspension point of the power diagram reaches the loading point, To simulate the crank torque when the suspension point reaches the loading point, is the weight of the crank and counterweight, sin is the sine function, is the circumference of a circle, is the angle between the walking beam and the connecting rod when the suspension point reaches the loading point, is the angle between the crank and the connecting rod when the suspension point reaches the loading point, is the position angle of the walking beam when the suspension point reaches the loading point, is the angle between the line from the center of the oil beam to the center of the crank and the vertical direction, b is the length of the rear wall of the walking beam, r is the crank radius, is the length of the forearm of the walking beam, To simulate the displacement of the suspension point when it reaches the loading point, S For the stroke, is the position angle of the walking beam at the highest point, is the position angle of the walking beam at the lowest point, , n Equal to 0, 1, 2, 3... n , m , n is a constant; Slope of loading line of simulated power diagram The calculation formula is as follows: ， in, is the slope of the simulated power diagram loading line, To simulate the load when the suspension point of the power diagram reaches the loading point, To simulate the load when the suspension point of the power diagram reaches the bottom dead center, To simulate the displacement of the suspension point when it reaches the loading point, To simulate the displacement of the power diagram suspension point when it reaches the bottom dead center, , n Equal to 0, 1, 2, 3... n , m , n is a constant.

9. A method for estimating the real-time efficiency of a beam pumping unit drive belt according to claim 1, characterized in that: In step S11, the angle between the n+1 simulated power diagram loading lines and the simulated power diagram dead point line The calculation formula is as follows: ， in, is the angle between the simulated power diagram loading line and the simulated power diagram dead point line, arctan is the inverse tangent trigonometric function, is the slope of the simulated power diagram loading line, is the slope of the dead point line of the simulated power diagram, , n Equal to 0, 1, 2, 3... n , m , n is a constant.

10. The method for estimating the real-time efficiency of a beam pumping unit drive belt according to claim 1, characterized in that: In step S13, the belt transmission efficiency The calculation formula is as follows: ， in, is the belt transmission efficiency, is the actual transmitted power and P is the initial power.

Images