CN117436319B - Oil pumping well production gas-oil ratio calculation method based on ground indicator diagram - Google Patents
Oil pumping well production gas-oil ratio calculation method based on ground indicator diagram Download PDFInfo
- Publication number
- CN117436319B CN117436319B CN202311761535.6A CN202311761535A CN117436319B CN 117436319 B CN117436319 B CN 117436319B CN 202311761535 A CN202311761535 A CN 202311761535A CN 117436319 B CN117436319 B CN 117436319B
- Authority
- CN
- China
- Prior art keywords
- pump
- oil
- gas
- sucker rod
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 59
- 238000010586 diagram Methods 0.000 title claims abstract description 57
- 238000004364 calculation method Methods 0.000 title claims abstract description 38
- 238000005086 pumping Methods 0.000 title claims abstract description 20
- 238000000034 method Methods 0.000 claims abstract description 55
- 239000003129 oil well Substances 0.000 claims abstract description 22
- 238000006073 displacement reaction Methods 0.000 claims abstract description 13
- 239000007788 liquid Substances 0.000 claims abstract description 13
- 230000008859 change Effects 0.000 claims abstract description 11
- 238000013016 damping Methods 0.000 claims abstract description 9
- 239000003921 oil Substances 0.000 claims description 70
- 230000008569 process Effects 0.000 claims description 22
- 239000007789 gas Substances 0.000 claims description 16
- 239000010779 crude oil Substances 0.000 claims description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 14
- 239000000243 solution Substances 0.000 claims description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 4
- 238000007872 degassing Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 239000011259 mixed solution Substances 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 239000000463 material Substances 0.000 claims description 2
- 239000003345 natural gas Substances 0.000 claims description 2
- 229920006395 saturated elastomer Polymers 0.000 claims description 2
- 239000003208 petroleum Substances 0.000 abstract description 2
- 238000012544 monitoring process Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 206010063385 Intellectualisation Diseases 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000035699 permeability Effects 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/008—Monitoring of down-hole pump systems, e.g. for the detection of "pumped-off" conditions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Theoretical Computer Science (AREA)
- Geology (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Agronomy & Crop Science (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention belongs to the technical field of petroleum exploitation, and particularly relates to a calculation method of gas-oil ratio of oil pumping well production based on a ground indicator diagram, which adopts ground displacement and load as boundary conditions, adopts a finite difference solution to solve a one-dimensional wave equation with damping, and calculates to obtain a downhole pump indicator diagram; establishing an underground pump diagram curvature solving model, and solving the position of an opening and closing point of a pump valve; and (3) establishing a calculation model of pump efficiency loss caused by gas and underfill and a calculation model of production gas-oil ratio, and solving the production gas-oil ratio by combining calculation of pump efficiency loss caused by actual pump efficiency, stroke loss, pump leakage and liquid volume change of an oil well. The method can calculate the production gas-oil ratio only by acquiring the actually measured ground indicator diagram and the production data, has high precision and quick dynamic response, can effectively improve the data acquisition capacity and informatization management level of the oil pumping well, and provides support for intelligent exploitation, cost reduction and synergy of oil gas.
Description
Technical Field
The invention belongs to the technical field of petroleum exploitation, and particularly relates to a calculation method of a gas-oil ratio of oil pumping well production based on a ground indicator diagram.
Background
The gas-oil ratio refers to the ratio of daily gas production to daily oil production of an oil well, and represents the amount of gas produced with each ton of crude oil produced. In the process of exploiting crude oil from stratum to ground, the pressure of fluid is continuously reduced, when the pressure is lower than the saturation pressure of crude oil, light components in crude oil are separated to form gas phase, if the crude oil is degassed in stratum, the viscosity of underground crude oil is increased, meanwhile, the relative permeability of oil phase is reduced due to the existence of gas phase, and the flow resistance is increased; if degassing occurs in the well bore, the pump efficiency of the pump, the oil well production, the efficiency of the production system, etc. are affected. Therefore, the production gas-oil ratio is an important parameter essential for optimizing and adjusting the oil gas development and exploitation work.
Before the digitalized and informationized construction of the oil-gas field, the oil-gas-water three-phase separator and the gas-liquid flowmeter are generally adopted on site to measure the oil production, the water production and the gas production, and then the gas-oil ratio or the gas-liquid ratio is calculated, so that the equipment cost and the maintenance cost are high, and the problems of potential safety hazard and the like exist. At present, digitization, informatization and intellectualization of oil and gas fields rapidly develop, the digitization coverage rate of a production site reaches more than 90%, the labor organization mode and the production operation mode are greatly changed, a metering device at the wellhead of an oil well is basically canceled, the measurement of the liquid production amount is realized, the gas-oil ratio is difficult to measure, and the defect of the important parameter is caused.
Therefore, by means of mechanism analysis, a calculation method of the production gas-oil ratio of the oil pumping unit well based on the ground indicator diagram is established, real-time automatic monitoring of the production gas-oil ratio of the oil pumping unit well is achieved, the problems of high parameter acquisition difficulty and high investment are solved, the intelligent and informationized management level of an oil well is improved, and theoretical and technical support is provided for real-time optimization of oil well production.
Disclosure of Invention
The invention aims to provide a calculation method of the gas-oil ratio of the production of the rod-pumped well based on a ground indicator diagram so as to realize real-time automatic monitoring of the gas-oil ratio of the production of the rod-pumped well.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
a calculation method of the gas-oil ratio of the oil pumping well production based on the ground indicator diagram comprises the following steps:
s1, solving a one-dimensional band damping wave equation by using ground displacement and load (ground indicator diagram) as boundary conditions and adopting a finite difference solution to obtain a downhole pump indicator diagram.
The propagation of stress waves in a sucker rod string can be described by a one-dimensional damped wave equation:
;
in the formula, u is%xT) is any section of the sucker rod stringxAt the position oftTime displacement, namely, downward positive and m;cthe propagation speed of stress wave in the sucker rod is m/s;vis damping coefficient, 1/s;
for a single stage sucker rod, the finite difference of the model is decomposed into:
;
the boundary conditions are:
;
in the method, in the process of the invention,E r the elastic modulus of the sucker rod is MPa;A r is the cross-sectional area of the sucker rod, m 2 ;
Convergence conditions:;
for a multi-stage sucker rod, solving according to continuous conditions comprises the following steps:
;
s2, establishing a curvature solving model of the underground pump work diagram, and solving the position of the point with the largest curvature change in the pump work diagram, namely the position of the opening and closing points of the valve.
The curvature K of each discrete point of the pump diagram can be utilized by the discrete pointsX i (x i ,y i )The previous discrete pointX i-1 (x i-1 ,y i-1 )And the latter discrete pointX i+1 (x i+1 ,y i+1 )The geometric relation among three points is calculated, and the mathematical model of the curvature of a certain discrete point is as follows:
;
the specific steps for calculating the position of the valve opening and closing point are as follows:
(1) calculating the coordinate average value of each point by using a five-point average method;
(2) respectively solving the maximum value of the abscissa of the discrete points of the pump work diagramX max And minimum valueX min Maximum value of ordinateY max And minimum valueY min ;
(3) The discrete points are normalized, and the normalization formula is as follows:
;
(4) expanding the normalized pump work diagram along the plunger stroke, and changing the pump work diagram from a closed curve to a single-value curve;
(5) calculating the curvature value of each discrete point according to the curvature calculation model of any discrete pointK i ;
(6) According toFinding the curvature of arbitrary discrete pointsK i Curvature from its subsequent discrete pointK i+1 Variation of (2)δ i ;
(7) Average load introducing normalized loaddfaAssume thatdfu=dfa+0.1,dfd=dfa-0.1, then during the upstroke, at normalized load greater thandfuDetermining the opening and closing points of the fixed valve according to two points with the largest curvature variation in the range; during the downstroke, at normalized load less thandfdThe maximum change in curvatureThe two points of the traveling valve are determined as the opening and closing points.
S3, calculating pump efficiency loss caused by stroke loss, pump leakage and liquid volume change;
drawing a pump work diagram in a plane rectangular coordinate system, marking the position of an opening and closing point of the valve obtained by solving in the pump work diagram,P SO (X SO ,Y SO )is a fixed valve opening point;P SC (X SC ,Y SC )is a fixed valve closing point;P TO (X TO ,Y TO )is the trip valve opening point;P TC (X TC ,Y TC )is a traveling valve closing point, whereinX L 、X R The abscissa of the left and right end points of the pump indicator diagram are respectively shown.
S31 pump loss due to pump leakage
(1) The pump efficiency loss caused by the leakage of the traveling valve is as follows:
;
(2) the pump efficiency loss caused by the leakage of the fixed valve is as follows:
as can be seen from the figures of the drawing,S TO the pump stroke loss caused by the leakage of the fixed valve during the downstroke can be considered as the pump stroke loss caused by the leakage of the fixed valve during the downstrokeS TOl Equal toS TC It is possible to obtain:
;
the pump efficiency loss caused by the S32 stroke loss is:
;
wherein eta is λ To consider the ratio of plunger stroke to polish rod stroke after the sucker rod string and tubing string are elastically stretched; s is(s) p M is used for considering the plunger stroke after the sucker rod string and the tubing string are elastically stretched; s is the stroke of the polish rod, m;
the pump efficiency loss caused by the change of the volume of the S33 liquid is as follows:
;
in the method, in the process of the invention,n w is the standard condition%P st ,T st ) The water content of the lower mixed solution,%;B O is the volume coefficient of crude oil under stratum condition, m 3 / m 3 ;B w Is the volume coefficient of water under stratum condition, m 3 / m 3 ;
S4, calculating actual pump efficiency of the oil well according to the actual yield and the theoretical displacement of the pump;
;
in the method, in the process of the invention,ηthe actual pump efficiency of the oil well is achieved;Q actual practice is that of For actual liquid production of oil well, m 3 /d;Q Theory of Is the theoretical displacement of oil well, m 3 /d。
S5, calculating a pump filling coefficient, establishing a calculation model of the production gas-oil ratio, and solving the production gas-oil ratio.
The pump fullness coefficient is:
;
in the method, in the process of the invention,βis the filling coefficient of the pump, decimal;
according toβAnd the relation of the production gas-oil ratio to obtain a production gas-oil ratio calculation model:
;
in the method, in the process of the invention,R p to produce the gas-oil ratio, m 3 / m 3 ;R bs For the ratio of dissolved gas to oil in the pump, m 3 / m 3 ;p i Is sinking pressure, MPa;f w water content,%.
From the above, the calculation method of the gas-oil ratio of the oil pumping well production based on the ground indicator diagram only needs to realize real-time automatic monitoring of the gas-oil ratio of the oil pumping well production according to the ground actually measured indicator diagram and production data, and has higher precision; the real-time data acquisition and informatization and intelligent management level of the oil well can be improved; and the labor intensity is reduced, the labor amount is reduced, the production efficiency is improved, and the production cost is saved.
Drawings
For a clearer description of one or more embodiments of the present description or of the solutions of the prior art, the drawings that are necessary for the description of the embodiments or of the prior art will be briefly described, it being apparent that the drawings in the description below are only one or more embodiments of the present description, from which other drawings can be obtained, without inventive effort, for a person skilled in the art.
FIG. 1 is a flow chart of a method for calculating the gas-oil ratio of the production of an oil pumping unit well based on a ground indicator diagram;
FIG. 2 is a schematic diagram of a discrete point curvature solution model;
FIG. 3 is a typical pumping diagram for a traveling valve, a fixed valve leak, and pump underfill and gas influencing conditions;
FIG. 4 is a graph comparing the results of the calculation of the gas-oil ratio of 17 wells.
Detailed Description
The present disclosure will be further described in detail with reference to specific embodiments in order to make the objects, technical solutions and advantages of the present disclosure more apparent.
It is noted that unless otherwise defined, technical or scientific terms used in one or more embodiments of the present disclosure should be taken in a general sense as understood by one of ordinary skill in the art to which the present disclosure pertains. The use of the terms "first," "second," and the like in one or more embodiments of the present description does not denote any order, quantity, or importance, but rather the terms "first," "second," and the like are used to distinguish one element from another. The word "comprising" or "comprises", and the like, means that elements or items preceding the word are included in the element or item listed after the word and equivalents thereof, but does not exclude other elements or items. The terms "connected" or "connected," and the like, are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "upper", "lower", "left", "right", etc. are used merely to indicate relative positional relationships, which may also be changed when the absolute position of the object to be described is changed.
The application provides a calculation example based on 17 wells on site, and provides a calculation method (figure 1) of the gas-oil ratio of the oil pumping unit well production based on a ground indicator diagram, which comprises the following steps:
s1, solving a one-dimensional band damping wave equation by using ground displacement and load (ground indicator diagram) as boundary conditions and adopting a finite difference solution to obtain a downhole pump indicator diagram.
The propagation of stress waves in a sucker rod string can be described by a one-dimensional damped wave equation:
;
in the formula, u is%xT) is any section of the sucker rod stringxAt the position oftTime displacement, namely, downward positive and m;cthe propagation speed of stress wave in the sucker rod is m/s;vis a damping coefficient, 1/s.
Wherein,cthe calculation method of (2) is as follows:
;
in the method, in the process of the invention,E r the elastic modulus of the sucker rod is MPa;ρ r density of sucker rod, kg/m 3 。
The damping coefficient calculating method comprises the following steps:
;
in the method, in the process of the invention,μis the dynamic viscosity of the oil well fluid,;ρ r for sucker rod material density->;A r Is the cross-sectional area of the sucker rod->;D r The diameter of the sucker rod, m;D C the outer diameter of the sucker rod coupling is m;D ti the inner diameter of the oil pipe is m;L r length of single sucker rod, m.
For a single stage sucker rod, the finite difference of the model is decomposed into:
;
the boundary conditions are:
;
in the method, in the process of the invention,E r the elastic modulus of the sucker rod is MPa;A r is the cross-sectional area of the sucker rod, m 2 ;
Convergence conditions:;
for a multi-stage sucker rod, solving according to continuous conditions comprises the following steps:
;
s2, establishing a curvature solving model (figure 2) of the underground pump work diagram, and solving the position of the point with the largest curvature change in the pump work diagram, namely the position of the opening and closing points of the valve.
The curvature K of each discrete point of the pump diagram can be utilized by the discrete pointsX i (x i ,y i )The previous discrete pointX i-1 (x i-1 ,y i-1 )And the latter discrete pointX i+1 (x i+1 ,y i+1 )The geometric relation among three points is calculated, and the mathematical model of the curvature of a certain discrete point is as follows:
;
wherein S is Δ The calculation method comprises the following steps:
;
the specific steps for calculating the position of the valve opening and closing point are as follows:
(1) calculating the coordinate average value of each point by using a five-point average method;
(2) respectively solving the maximum value of the abscissa of the discrete points of the pump work diagramX max And minimum valueX min Maximum value of ordinateY max And minimum valueY min ;
(3) The discrete points are normalized, and the normalization formula is as follows:
;
(4) expanding the normalized pump work diagram along the plunger stroke, and changing the pump work diagram from a closed curve to a single-value curve;
(5) calculating the curvature value of each discrete point according to the curvature calculation model of any discrete pointK i ;
(6) According toFinding the curvature of arbitrary discrete pointsK i Curvature from its subsequent discrete pointK i+1 Variation of (2)δ i ;
(7) Average load introducing normalized loaddfaAssume thatdfu=dfa+0.1,dfd=dfa-0.1, then during the upstroke, at normalized load greater thandfuDetermining the opening and closing points of the fixed valve according to two points with the largest curvature variation in the range; during the downstroke, at normalized load less thandfdAnd determining the opening and closing points of the traveling valve according to the two points with the largest curvature variation in the range.
S3, calculating pump efficiency loss caused by stroke loss, pump leakage and liquid volume change;
drawing a pump work diagram in a plane rectangular coordinate system, marking the position of an opening and closing point of the valve obtained by solving in the pump work diagram (figure 3),P SO (X SO ,Y SO )is a fixed valve opening point;P SC (X SC ,Y SC )is a fixed valve closing point;P TO (X TO ,Y TO )is the trip valve opening point;P TC (X TC ,Y TC )is a traveling valve closing point, whereinX L 、X R The abscissa of the left and right end points of the pump indicator diagram are respectively shown.
S31 pump loss due to pump leakage
(1) The pump efficiency loss caused by the leakage of the traveling valve is as follows:
;
(2) the pump efficiency loss caused by the leakage of the fixed valve is
As can be seen from the figures of the drawing,S TO including pump gas and pump stroke losses due to loss of the fixed valve during downstroke due to underfill effects. At this time, it can be considered that the pump stroke is lost due to the leakage amount of the down-stroke fixed valveS TOl Equal toS TC It is possible to obtain:
;
the pump efficiency loss caused by the S32 stroke loss is:
;
wherein eta is λ To consider the ratio of plunger stroke to polish rod stroke after the sucker rod string and tubing string are elastically stretched; s is(s) p M is used for considering the plunger stroke after the sucker rod string and the tubing string are elastically stretched; s is the stroke of the polish rod, m;
the pump efficiency loss caused by the change of the volume of the S33 liquid is as follows:
;
in the method, in the process of the invention,n w is the standard condition%P st ,T st ) The water content of the lower mixed solution,%;B O is the volume coefficient of crude oil under stratum condition, m 3 / m 3 ;B w Is the volume coefficient of water under stratum condition, m 3 / m 3 。
Wherein,B O the calculation method comprises the following steps:
;
wherein R is s To dissolve the gas-oil ratio, m 3 / m 3 ;r o The relative density of the crude oil is the ground degassing, and the dimensionless is achieved;r g for the relative density of the separator gas (taking air with a density of 1.0kg/m 3 ) Dimensionless;Ttemperature, K; ρ 0 Is the density of crude oil, kg/m 3 ;P i Is the pressure and MPa.
B w The calculation method comprises the following steps:
;
in the method, in the process of the invention,D 0 、D 1 、D 2 as the coefficient of the light-emitting diode,,θ t is the formation temperature, F; coefficients ofa 0 、a 1 、a 2 The values of (2) are calculated according to the saturated water of the natural gas, and are specifically shown in table 1.
S4, calculating actual pump efficiency of the oil well according to the actual yield and the theoretical displacement of the pump;
;
in the method, in the process of the invention,ηthe actual pump efficiency of the oil well is achieved;Q actual practice is that of For actual liquid production of oil well, m 3 /d;Q Theory of Is the theoretical displacement of oil well, m 3 /d。
S5, calculating a pump filling coefficient, establishing a calculation model of the production gas-oil ratio, and solving the production gas-oil ratio.
The pump fullness coefficient is:
;
in the method, in the process of the invention,βis the filling coefficient of the pump, decimal;
according toβAnd the relation of the production gas-oil ratio to obtain a production gas-oil ratio calculation model:
;
in the method, in the process of the invention,R p to produce the gas-oil ratio, m 3 / m 3 ;R bs For the ratio of dissolved gas to oil in the pump, m 3 / m 3 ;p i Is sinking pressure, MPa;f w water content,%.
To verify the effect of the above gasoline ratio calculation method, it is known from the calculation and analysis of the on-site 17-well data (see table 2, fig. 4): the maximum relative error calculated by the gas-oil ratio of the oil well is 9.75%, the minimum relative error is 1.96%, and the average relative error is 6.46%. The results show that: the production gas-oil ratio calculation model has better precision, and the prediction error is in the engineering allowable range, which indicates that the method realizes the real-time automatic monitoring of the production gas-oil ratio of the pumping unit well.
The calculation method of the production gas-oil ratio of the oil pumping unit well based on the ground indicator diagram can realize real-time and accurate calculation of the production gas-oil ratio of a single well according to the ground indicator diagram and daily production data, and solves the difficulty of real-time monitoring of the production gas-oil ratio.
Those of ordinary skill in the art will appreciate that: the discussion of any of the embodiments above is merely exemplary and is not intended to suggest that the scope of the disclosure, including the claims, is limited to these examples; combinations of features of the above embodiments or in different embodiments are also possible within the spirit of the present disclosure, steps may be implemented in any order, and there are many other variations of the different aspects of one or more embodiments described above which are not provided in detail for the sake of brevity.
The present disclosure is intended to embrace all such alternatives, modifications and variances which fall within the broad scope of the appended claims. Any omissions, modifications, equivalents, improvements, and the like, which are within the spirit and principles of the one or more embodiments of the disclosure, are therefore intended to be included within the scope of the disclosure.
Claims (6)
1. The calculation method of the gas-oil ratio of the oil pumping well production based on the ground indicator diagram is characterized by comprising the following steps:
s1, solving a one-dimensional band damping wave equation by using ground displacement and load as boundary conditions and adopting a finite difference solution to obtain an underground pump diagram:
the propagation of stress waves in a sucker rod string can be described by a one-dimensional damped wave equation:
;
in the formula, u is%xT) is any section of the sucker rod stringxAt the position oftTime displacement, namely, downward positive and m;cthe propagation speed of stress wave in the sucker rod is m/s;vis damping coefficient, 1/s;
for a single stage sucker rod, the finite difference of the model is decomposed into:
;
the boundary conditions are:
;
in the method, in the process of the invention,E r the elastic modulus of the sucker rod is MPa;A r is the cross-sectional area of the sucker rod, m 2 ;
Convergence conditions:;
for a multi-stage sucker rod, solving according to continuous conditions comprises the following steps:
;
s2, establishing a curvature solving model of the underground pump work diagram, and solving the position of the point with the largest curvature change in the pump work diagram, namely the position of the opening and closing points of the valve:
the curvature K of each discrete point of the pump diagram can be utilized by the discrete pointsX i (x i ,y i )The previous discrete pointX i-1 (x i-1 ,y i-1 )And the latter discrete pointX i+1 (x i+1 ,y i+1 )The geometric relation among three points is calculated, and the mathematical model of the curvature of a certain discrete point is as follows:
;
the specific steps for calculating the position of the valve opening and closing point are as follows:
(1) calculating the coordinate average value of each point by using a five-point average method;
(2) respectively solving the maximum value of the abscissa of the discrete points of the pump work diagramX max And minimum valueX min Maximum value of ordinateY max And minimum valueY min ;
(3) The discrete points are normalized, and the normalization formula is as follows:
;
(4) expanding the normalized pump work diagram along the plunger stroke, and changing the pump work diagram from a closed curve to a single-value curve;
(5) calculating the curvature value of each discrete point according to the curvature calculation model of any discrete pointK i ;
(6) According toFinding the curvature of arbitrary discrete pointsK i Curvature from its subsequent discrete pointK i+1 Variation of (2)δ i ;
(7) Average load introducing normalized loaddfaAssume thatdfu=dfa+0.1,dfd=dfa-0.1, then during the upstroke, at normalized load greater thandfuDetermining the opening and closing points of the fixed valve according to two points with the largest curvature variation in the range; during the downstroke, at normalized load less thandfdDetermining the opening and closing points of the traveling valve according to two points with the largest curvature variation in the range;
s3, calculating pump efficiency loss caused by stroke loss, pump leakage and liquid volume change;
drawing a pump work diagram in a plane rectangular coordinate system, marking the position of an opening and closing point of the valve obtained by solving in the pump work diagram,P SO (X SO ,Y SO )is a fixed valve opening point;P SC (X SC ,Y SC )is a fixed valve closing point;P TO (X TO ,Y TO )is the trip valve opening point;P TC (X TC ,Y TC )is a traveling valve closing point, whereinX L 、X R The abscissa of the left end point and the right end point of the pump indicator diagram are respectively;
pump efficiency loss due to S31 pump loss:
(1) the pump efficiency loss caused by the leakage of the traveling valve is as follows:
;
(2) the pump efficiency loss caused by the leakage of the fixed valve is as follows:
ST in the pump diagram O The pump stroke loss due to the leakage of the fixed valve during the downstroke including the influence of the pump gas and the underfill can be considered as the pump stroke loss due to the leakage of the fixed valve during the downstrokeS TOl Equal toS TC It is possible to obtain:
;
the pump efficiency loss caused by the S32 stroke loss is:
;
wherein eta is λ To consider the ratio of plunger stroke to polish rod stroke after the sucker rod string and tubing string are elastically stretched; s is(s) p M is used for considering the plunger stroke after the sucker rod string and the tubing string are elastically stretched; s is the stroke of the polish rod, m;
the pump efficiency loss caused by the change of the volume of the S33 liquid is as follows:
;
in the method, in the process of the invention,n w is the standard condition%P st ,T st ) The water content of the lower mixed solution,%;B O is the volume coefficient of crude oil under stratum condition, m 3 / m 3 ;B w Is the volume coefficient of water under stratum condition, m 3 / m 3 ;
S4, calculating actual pump efficiency of the oil well according to the actual yield and the theoretical displacement of the pump;
;
in the method, in the process of the invention,ηthe actual pump efficiency of the oil well is achieved;Q actual practice is that of For actual liquid production of oil well, m 3 /d;Q Theory of Is the theoretical displacement of oil well, m 3 /d;
S5, calculating a pump filling coefficient, establishing a calculation model of a production gas-oil ratio, and solving the production gas-oil ratio;
the pump fullness coefficient is:
;
in the method, in the process of the invention,βis the filling coefficient of the pump, decimal;
according toβAnd the relation of the production gas-oil ratio to obtain a production gas-oil ratio calculation model:
;
in the method, in the process of the invention,R p to produce the gas-oil ratio, m 3 / m 3 ;R bs For the ratio of dissolved gas to oil in the pump, m 3 / m 3 ;p i Is sinking pressure, MPa;f w water content,%.
2. The method for calculating the gas-oil ratio of the production of the rod-pumped well based on the ground indicator diagram according to claim 1, wherein in the step S1, the propagation speed of the stress wave in the sucker rodcThe calculation method of (2) is as follows:
;
in the method, in the process of the invention,E r the elastic modulus of the sucker rod is MPa;ρ r density of sucker rod, kg/m 3 。
3. The method for calculating the gas-oil ratio of the oil pumping well production based on the ground indicator diagram according to claim 1, wherein in the step S1, the damping coefficient isvThe calculation method comprises the following steps:
;
in the method, in the process of the invention,μis the dynamic viscosity of the oil well fluid,;ρ r for sucker rod material density->;A r Is the cross-sectional area of the sucker rod->;D r The diameter of the sucker rod, m;D C the outer diameter of the sucker rod coupling is m;D ti the inner diameter of the oil pipe is m;L r length of single sucker rod, m.
4. The method for calculating the gas-oil ratio of the oil pumping well production based on the ground indicator diagram according to claim 1, wherein in the step S2, S Δ The calculation method comprises the following steps:
。
5. the method for calculating the gas-oil ratio of the oil pumping well production based on the ground indicator diagram according to claim 1, wherein in the step S3, the volume coefficient of the crude oil under the stratum conditionB O The calculation method comprises the following steps:
;
wherein R is s To dissolve the gas-oil ratio, m 3 / m 3 ;r o The relative density of the crude oil is the ground degassing, and the dimensionless is achieved;r g the relative density of the separator gas is dimensionless;Ttemperature, K;ρ 0 is the density of crude oil, kg/m 3 ;P i Is the pressure and MPa.
6. The method for calculating the gas-oil ratio of the oil pumping well production based on the ground indicator diagram according to claim 1, wherein in the step S3, the volume coefficient of water under the stratum conditionB w The calculation method comprises the following steps:
;
in the method, in the process of the invention,D 0 、D 1 、D 2 as the coefficient of the light-emitting diode,,θ t is the formation temperature, F; coefficients ofa 0 、a 1 、a 2 The value of (2) is calculated according to the saturated water of the natural gas.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311761535.6A CN117436319B (en) | 2023-12-20 | 2023-12-20 | Oil pumping well production gas-oil ratio calculation method based on ground indicator diagram |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311761535.6A CN117436319B (en) | 2023-12-20 | 2023-12-20 | Oil pumping well production gas-oil ratio calculation method based on ground indicator diagram |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117436319A CN117436319A (en) | 2024-01-23 |
CN117436319B true CN117436319B (en) | 2024-03-19 |
Family
ID=89553962
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311761535.6A Active CN117436319B (en) | 2023-12-20 | 2023-12-20 | Oil pumping well production gas-oil ratio calculation method based on ground indicator diagram |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117436319B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108222897A (en) * | 2016-12-22 | 2018-06-29 | 中国科学院沈阳自动化研究所 | It is a kind of that the method for drawing Dynamic Control Chart is combined with statistics based on calculating |
CN108223343A (en) * | 2016-12-22 | 2018-06-29 | 中国科学院沈阳自动化研究所 | A kind of pumpingh well method for controlling frequency conversion based on Dynamic Control Chart and indicator card |
CN110206536A (en) * | 2019-07-05 | 2019-09-06 | 海默潘多拉数据科技(深圳)有限公司 | A kind of well head Liquid output acquisition method based on pump dynamometers |
CN114790885A (en) * | 2021-01-08 | 2022-07-26 | 中国石油天然气股份有限公司 | Method and device for measuring output of oil pumping well |
-
2023
- 2023-12-20 CN CN202311761535.6A patent/CN117436319B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108222897A (en) * | 2016-12-22 | 2018-06-29 | 中国科学院沈阳自动化研究所 | It is a kind of that the method for drawing Dynamic Control Chart is combined with statistics based on calculating |
CN108223343A (en) * | 2016-12-22 | 2018-06-29 | 中国科学院沈阳自动化研究所 | A kind of pumpingh well method for controlling frequency conversion based on Dynamic Control Chart and indicator card |
CN110206536A (en) * | 2019-07-05 | 2019-09-06 | 海默潘多拉数据科技(深圳)有限公司 | A kind of well head Liquid output acquisition method based on pump dynamometers |
CN114790885A (en) * | 2021-01-08 | 2022-07-26 | 中国石油天然气股份有限公司 | Method and device for measuring output of oil pumping well |
Non-Patent Citations (3)
Title |
---|
Experimental study on the performance of foamy oil flow under different solution gas-oil ratios;Songyan Li 等;《RSC Advances》;20151231;第66797-66806页 * |
基于功图的油井动液面计算模型及应用;陈德春;张瑞超;孟红霞;谢文献;王欣辉;;科学技术与工程;20151118(32) * |
电潜柱塞泵优化设计与诊断;褚浩元;《中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑)》;20210215;正文第2-4章 * |
Also Published As
Publication number | Publication date |
---|---|
CN117436319A (en) | 2024-01-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103498663B (en) | A kind of sucker rod pump hoisting system pumping technological parameter defining method and device | |
CN107608940B (en) | Method for determining oil well interval pumping period | |
CN103590824A (en) | Productivity calculation method for tight gas horizontal wells modified by multi-stage fracturing | |
Baranov et al. | Constructing and calculating of multistage sucker rod string according to reduced stress | |
CN106437682B (en) | A method of prediction oil well indicator card | |
CN110206536A (en) | A kind of well head Liquid output acquisition method based on pump dynamometers | |
CN111810101B (en) | Dynamic analysis method and device for water-drive reservoir | |
CN105464644A (en) | Wellhead pressure monitoring method and control system | |
CN103726815B (en) | A kind of CO 2drive produced well pit shaft fluidised form is determined and parameter optimization method | |
CN106547930A (en) | Consider the gas drainage radius computational methods of tight gas reservoir seepage flow mechanism | |
CN105672997A (en) | Monitoring method for formation leakage of drilling fluid | |
CN111963151A (en) | Method for determining formation pressure through suspension point static load of oil pumping unit | |
CN104153982B (en) | A kind of acquisition characteristic method and device of rod-pumped well downhole system | |
CN103498661B (en) | A kind of method determining oil reservoir physical data under high pressure | |
CN110761743A (en) | Thick oil and thin oil mixed gas lift process pipe column and implementation method | |
CN106761680A (en) | A kind of chemical viscosity reduction auxiliary threaded rod pump lifts the determination methods of viscous crude technique | |
CN102562024B (en) | Fracturing design method for optimizing uniform spreading concentration | |
Jiang et al. | An efficient downhole oil/water-separation system with sucker-rod pump | |
CN117436319B (en) | Oil pumping well production gas-oil ratio calculation method based on ground indicator diagram | |
CN111489090A (en) | Method and device for analyzing influence of wellhead back pressure on oil well yield and storage medium | |
CN104405364B (en) | A kind of evaluation method and device for oil well production characteristics | |
RU2700738C1 (en) | Method of improving reliability of water cut monitoring of products of oil producing wells equipped with sucker-rod bottom pumps | |
CN104453796B (en) | Intelligent water drainage and gas production device for gas well | |
CN106682790B (en) | Method for prejudging lifting of shallow heavy oil reservoir single-phase flow screw pump | |
US20210270124A1 (en) | Method for distinguishing authenticity of high-pressure physical property parameters of oil reservoirs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |