CN108763657A - Horizontal gas well continuously takes liquid critical flow and determines method and device - Google Patents
Horizontal gas well continuously takes liquid critical flow and determines method and device Download PDFInfo
- Publication number
- CN108763657A CN108763657A CN201810419162.7A CN201810419162A CN108763657A CN 108763657 A CN108763657 A CN 108763657A CN 201810419162 A CN201810419162 A CN 201810419162A CN 108763657 A CN108763657 A CN 108763657A
- Authority
- CN
- China
- Prior art keywords
- gas well
- horizontal gas
- critical
- liquid
- horizontal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000007788 liquid Substances 0.000 title claims abstract description 197
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000004519 manufacturing process Methods 0.000 claims abstract description 17
- 239000007791 liquid phase Substances 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 238000010008 shearing Methods 0.000 claims description 11
- 230000005484 gravity Effects 0.000 claims description 8
- 238000011958 production data acquisition Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 6
- 238000011084 recovery Methods 0.000 abstract description 6
- 238000011065 in-situ storage Methods 0.000 abstract description 5
- 239000010408 film Substances 0.000 description 62
- 238000004364 calculation method Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 9
- 241000521257 Hydrops Species 0.000 description 7
- 238000004590 computer program Methods 0.000 description 7
- 239000010409 thin film Substances 0.000 description 7
- 206010030113 Oedema Diseases 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000009795 derivation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- -1 pressure difference Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000005211 surface analysis Methods 0.000 description 1
- 230000005514 two-phase flow 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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Geochemistry & Mineralogy (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A kind of horizontal gas well of present invention offer continuously takes liquid critical flow and determines method and device, wherein this approach includes the following steps:Establish horizontal gas well critical liquid carrying flow rate computation model;Obtain horizontal gas well actual production data;Based on horizontal gas well critical liquid carrying flow rate computation model, according to the horizontal gas well creation data, determine that horizontal gas well continuously takes liquid critical flow.Since the program considers influence of the factors such as pressure difference and gas turbulence effect to liquid film when establishing horizontal gas well critical liquid carrying flow rate computation model, error can be reduced when calculating critical liquid carrying flow rate, the ultimate recovery for improving in-situ horizontal gas well has certain directive function.
Description
Technical field
The present invention relates to conventional water outlet gas well production technique field, more particularly to a kind of horizontal gas well continuously takes liquid critical flow
Amount determines method and device.
Background technology
With the continuous increase of normal gas pools Exploitation degree, strata pressure is gradually reduced with gas production, causes taking for gas
Liquid energy power declines;And side water is constantly invaded with low water, shaft bottom hydrops amount gradually increases, and is significantly dropped so as to cause gas well yield
It is low, or even stop production.Particularly with horizontal gas well, with the continuous increase of its at home and abroad each Gas Fields application scale, once go out
Existing hydrops phenomenon leads to the gas well underproduction, can greatly increase the cost of winning of gas well, influence economic benefit.Therefore accurately predict water
Flat gas well critical liquid carrying flow rate improves its ultimate recovery and is of great significance for reducing shaft bottom hydrops.
Many models for calculating critical liquid carrying flow rate are proposed in the prior art, and one kind is drop theory model:Recognize
It is to be carried and migrate with drops for liquid phase, if maximum gauge drop does not slip in pit shaft, shaft bottom is not in
Hydrops phenomenon;Another kind of is liquid film theoretical model:Think that liquid phase is to be carried out pit shaft with form of liquid film.By theoretical and experiment
Prove, the casing programme of horizontal gas well is complex, in pit shaft the biphase gas and liquid flow characteristics of motion and take liquid mechanism with straight well and calmly
It is larger to well difference;Drop can not exist steadily in the long term in horizontal gas well, mainly be carried out pit shaft by form of liquid film,
Therefore to take liquid mechanism more accurate for continuous based on liquid film theoretical model.But the liquid film model proposed at present is primarily directed to straight well
And directional well, difference of the horizontal well with vertical well and directional well in liquid film distribution and force-bearing situation is not considered, has ignored pressure
The influence of the factors to liquid film model such as difference and gas turbulence effect, error is larger, or even causes mistake.
Invention content
It continuously takes liquid critical flow an embodiment of the present invention provides a kind of horizontal gas well and determines method and device, it is contemplated that pressure
The influence of the factors to liquid film such as difference and gas turbulence effect, establishes horizontal gas well critical liquid carrying flow rate computation model, is calculating
Error can be reduced when critical liquid carrying flow rate, there is the ultimate recovery for improving in-situ horizontal gas well certain guidance to make
With.
The horizontal gas well continuously takes liquid critical flow and determines that method includes:
Establish horizontal gas well critical liquid carrying flow rate computation model;
Obtain horizontal gas well actual production data;
Horizontal gas is determined according to the horizontal gas well creation data based on horizontal gas well critical liquid carrying flow rate computation model
Well continuously takes liquid critical flow.
The horizontal gas well continuously takes liquid critical flow determining device:
Horizontal gas well critical liquid carrying flow rate computation model establishes module, is calculated for establishing horizontal gas well critical liquid carrying flow rate
Model;
Horizontal gas well actual production data acquisition module, for obtaining horizontal gas well actual production data;
Horizontal gas well continuously takes liquid critical flow determining module, for calculating mould based on horizontal gas well critical liquid carrying flow rate
Type determines that horizontal gas well continuously takes liquid critical flow according to the horizontal gas well creation data.
In embodiments of the present invention, horizontal gas well critical liquid carrying flow rate computation model is established, with existing liquid film model
It compares with directional well primarily directed to straight well, force analysis has been carried out to liquid film in horizontal gas well, it is contemplated that pressure difference and gas
The influence of the factors such as turbulence effect calculates critical liquid carrying flow rate by the horizontal gas well critical liquid carrying flow rate computation model, can be with
Reduce error, the ultimate recovery for improving in-situ horizontal gas well has certain directive function.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
Obtain other attached drawings according to these attached drawings.
Fig. 1 is that a kind of horizontal gas well provided in an embodiment of the present invention continuously takes liquid critical flow and determines method flow diagram;
Fig. 2 is a kind of horizontal well thin film drainage model schematic provided in an embodiment of the present invention;
Fig. 3 is the relation curve of modifying factor and the lifting factor under a kind of different pipe apertures provided in an embodiment of the present invention
Figure;
Fig. 4 is the graph of relation of modifying factor and viscosity under a kind of different pipe apertures provided in an embodiment of the present invention;
Fig. 5 is that a kind of horizontal gas well provided in an embodiment of the present invention continuously takes liquid critical flow determining device structure diagram.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation describes, it is clear that described embodiment is only a part of the embodiment of the present invention, instead of all the embodiments.Based on this
Embodiment in invention, every other reality obtained by those of ordinary skill in the art without making creative efforts
Example is applied, shall fall within the protection scope of the present invention.
In embodiments of the present invention, it provides a kind of horizontal gas well and continuously takes liquid critical flow and determine method, such as Fig. 1 institutes
Show, this method includes:
Step 101:Establish horizontal gas well critical liquid carrying flow rate computation model;
Step 102:Obtain horizontal gas well actual production data;
Step 103:Based on horizontal gas well critical liquid carrying flow rate computation model, according to the horizontal gas well creation data, really
Determine horizontal gas well and continuously takes liquid critical flow.
When it is implemented, step 101 is to establish horizontal gas well critical liquid carrying flow rate computation model as follows:
(1) horizontal gas well liquid film physical model is established, the parameter of horizontal gas well liquid film physical model is obtained.
Specifically, being experimentally confirmed, in horizontal gas well recovery process, liquid phase is existed in the form of liquid film.And liquid
Force-bearing situation and vertical well and directional well of the film in horizontal wellbore are completely different, therefore are pushed away according to straight well and directional well liquid film theory
The critical liquid carrying flow rate computation model led is no longer desirable for horizontal well.In horizontal wellbore, under the action of central gas stream shearing force,
Liquid film is carried out pit shaft;When the shearing stress between central gas stream and liquid film, pressure difference, liquid film and the pipe on thin film drainage direction
When shearing stress between wall tends to balance, the mobile of liquid film is stagnated, and causes shaft bottom hydrops.Therefore, in horizontal gas well
Liquid film physical model is as shown in Figure 2.
In Fig. 2, the parameter of horizontal gas well liquid film physical model includes:VsgIndicate the conversion flow velocity of central gas stream, m/s;Vl
Indicate the flow velocity of liquid film, m/s;D indicates tubing string internal diameter, m;τgIndicate the shearing stress between central gas stream and liquid film, MPa;τwIndicate liquid
Shearing stress between film and tube wall, MPa;DgIndicate the equivalent diameter of central gas stream, m, δcrIndicate the critical thickness of air-flow top liquid film
Degree, m.Wherein, the x-axis in figure indicates that the direction of thin film drainage, y-axis are indicated perpendicular to the direction of thin film drainage.
(2) according to the parameter of the horizontal well liquid film physical model of above-mentioned acquisition, establish that horizontal gas well is critical to take flow velocity meter
Calculate model.
Specifically, when liquid phase is in the minimal flow condition for stablizing migration, there are critical thickness for the liquid film in horizontal wellbore
Degree;And the influence due to gravity is shown according to Related Experimental Study, air-flow top liquid film critical thickness in horizontal wellbore
Less than lower part liquid film critical thickness.It discusses below using air-flow top liquid film as research object.Perpendicular to thin film drainage
Direction, for liquid film there are stress balance relationship, central gas stream is equal to liquid film lifting force vertically upward and buoyancy the gravity of liquid film,
I.e.
Fl+Fb=Fg (1)
In formula:FlIndicate lifting force of the central gas stream to liquid film vertically upward, N;FbIndicate buoyancy, N;FgIndicate liquid film
Gravity, N.
Lifting force is the power that turbulent gas applies liquid film in vertical direction.It is lifted according to suffered by known spherical bubbles
The computation model of power is deduced the general formulae that lifting force suffered by liquid phase is calculated in gas well:
In formula:ClIndicate the lifting factor, dimensionless;AlIndicate that lifting force acts on the sectional area on liquid film, m2;ρgIt indicates
The density of gas, kg/m3;A indicates the voidage of biphase gas and liquid flow, dimensionless.
It is analyzed by imfinitesimal method, its own gravity and buoyancy suffered by the liquid film of air-flow top are represented by:
In formula:δ indicates thickness of liquid film, m;ρlIndicate the density of liquid, kg/m3;G indicates acceleration of gravity, m/s2。
Since central gas stream sectional area existsTherefore the equivalent diameter of central gas stream is represented by:
Formula (2), formula (3) and formula (4) are substituted into formula (1), the expression formula of central gas stream conversion flow velocity can be obtained:
It is discharged gas well for conventional, liquid film can be considered incompressible Newtonian fluid, it is assumed that the thin film drainage in horizontal wellbore
For stable state laminar flow, then its governing equation is represented by:
In formula:P indicates to act on the pressure on liquid film, MPa.
For central gas stream, between liquid film and central gas stream, there are stress balance relationships, central gas stream to be applied by liquid film
Shearing stress and the resultant force of thin film drainage pressure drop are 0:
Formula (7) is arranged to and substituted into formula (6), can be obtained:
The boundary condition of the governing equation is as follows:
As y=0, Vl=0;As y=δ,
According to above-mentioned Boundary Condition for Solving formula (8), the expression formula of liquid phase flow rate can be obtained:
In formula:μlFor the viscosity of liquid phase, mPas.
Formula (9) is integrated by the thickness of liquid film to air-flow top, can be obtained under the liquid film unit perimeter of air-flow top
Flow:
In formula:QlFor the flow under the liquid film unit perimeter of air-flow top, m2/s。
Formula (10) can be obtained by arrangement:
When liquid film, which is in, is stablized the minimal flow condition carried, there are critical thickness for the top liquid film in horizontal wellbore
Degree, therefore formula (11) answers existence and unique solution, need to meet the following conditions:
Under the conditions of gas liquid two-phase flow, thickness of liquid film, voidage and tubing string internal diameter these three parameters relationship can indicate
For:
Simultaneous formula (5), formula (12) and formula (13), by arranging, the critical flow velocity of taking that horizontal gas well can be obtained calculates public affairs
Formula:
According to hydrodynamics correlative study it is found that Bond number is due to effect of surface tension and a dimensionless of determination
Number, i.e.,:
In formula:NBFor Bond number, dimensionless;σ is gas-liquid surface tension, N/m.
It is compared for convenience with classical critical liquid carrying flow rate calculation formula Turner models, by formula (14) and formula (15)
It is arranged, can be obtained:
For the ease of research, defining k-factor is:
Formula (17) is substituted into formula (16), can be obtained:
In formula, K is modifying factor;AndFor the result of calculation of Turner models.
From formula (17) it will be evident that modifying factor is mainly by liquid phase viscosity (μl), Bond number (NB) and the lifting factor (Cl)
It determines.Oil Field is often 50.53mm, 62.00mm, 76.00mm and 100.53mm with pipe aperture.Lower surface analysis is different once
The relationship of modifying factor and modifying factor and viscosity under the relationship of the lifting factor, different pipe apertures under pipe aperture.
Wherein, Fig. 3 is the graph of relation of modifying factor and the lifting factor under different pipe apertures.As can be seen from Figure 3:With
The continuous increase of pipe aperture, modifying factor gradually increases, therefore critical liquid carrying flow rate also accordingly increases, and horizontal gas well takes liquid
Ability gradually weakens;Lift the factor it is smaller when, modifying factor is affected by pipe aperture, with lifting the factor it is continuous
Increase, modifying factor is influenced to be gradually reduced by pipe aperture.With the increase of the lifting factor, modifying factor constantly reduces, and
The amplitude of reduction is gradually reduced with the increase of the lifting factor, therefore critical liquid carrying flow rate also reduces therewith, and horizontal gas well is taken
Liquid energy power gradually increases.
Fig. 4 is the graph of relation of modifying factor and viscosity under different pipe apertures.As can be seen from Figure 4:With in oil pipe
The continuous increase of diameter, modifying factor gradually increase, and when liquid phase viscosity is smaller, modifying factor is influenced very by pipe aperture
Small, with the continuous increase of viscosity, the influence of pipe aperture constantly increases, therefore critical liquid carrying flow rate also accordingly increases, horizontal gas
The fluid-carrying capability of well gradually weakens.With the continuous increase of liquid phase viscosity, modifying factor gradually increases, and increase amplitude with
The increase of viscosity and gradually increase, therefore critical liquid carrying flow rate also increases therewith, and the fluid-carrying capability of horizontal gas well gradually weakens.
Table 1 is modifying factor zoom table, in actual use, can be selected according to the actual conditions of in-situ horizontal gas well corresponding
Coefficient be modified.
1 modifying factor sublist of table
(3) flow velocity computation model is taken according to the horizontal gas well is critical, establishes the calculating of horizontal gas well critical liquid carrying flow rate
Model.
Specifically, the liquid critical flow velocity of taking that the model that above-mentioned (14) or (18) determine is calculated to gained is changed into the status of criterion
Under gas production, can show that horizontal gas well continuously takes liquid critical flow:
In formula:QscIndicate that horizontal gas well continuously takes liquid critical flow, 104m3/d;A indicates oil pipe sectional area, m2;pwIt indicates
Pressure, MPa;Z indicates the compressed coefficient of gas under certain temperature and pressure, dimensionless;T indicates temperature, K.
Embodiment
There are two mouthfuls of representative horizontal well A wells, B wells in certain gas field, and model is carried out using the daily production report material of this two mouthfuls of wells
The Fitting Calculation and precision analysis.Wherein, A wells are continuing to produce due to not hydrops, and B wells due to wellbore effusion it is more serious and
Closing well.The Field Production Data of A wells and B wells is as shown in table 2 (B wells are the creation data before closing well), critical to take the calculating of liquid model
The results are shown in Table 3.
2 mouthfuls of certain gas field of table 2 well Field Production Data
Pound sign | A wells | B wells |
Maximum hole angle/(°) | 88.9 | 89.2 |
Wellhead temperature/DEG C | 20 | 18 |
Pipe aperture/mm | 50.3 | 62 |
Well head pressure/MPa | 8.02 | 7.28 |
Gas production/104m3·d-1 | 3.01 | 0.57 |
Water yield/m3·d-1 | 0.01 | 7.13 |
Production status | Not hydrops | Hydrops |
Table 3 is critical to take the comparison of liquid the model calculation
Result of calculation by comparing 4 kinds of classical critical liquid carrying flow rate models is can be found that:Turner models and Li Minmo
Type relative error is larger, is above 36%, highest error is to the greatest extent up to 70% or so, the main reason is that model is to be based on drop theory
The derivation of progress, and by theoretical and it is demonstrated experimentally that horizontal gas well liquid phase mainly exists and is carried in the form of liquid film;And
And Turner models and Li Min models are to be based on vertical well condition, and cut in force-bearing situation and vertical well of the liquid phase in horizontal well
It is so different;Therefore, the result of calculation deviation of above-mentioned two model is larger.And the result of calculation relative error of Xiao Gaomian models also compared with
Greatly, be higher than 20%, be primarily due to the liquid film model be based on directional well under the conditions of derivation, while not considering the front and back pressure of liquid film
Influence of the difference to critical liquid carrying flow rate, and pressure difference and horizontal gas well critical liquid carrying flow rate are negatively correlated before and after liquid film, are affected,
Therefore when for calculated level gas well critical liquid carrying flow rate, result of calculation is bigger than normal.And the result of calculation of Belfroid models is opposite accidentally
Difference is above 18%, the main reason is that the model be also based on directional well under the conditions of derivation, and the model thinks that drop is
It moves up along wellbore centre line, is not inconsistent with the liquid situation of actually taking of horizontal gas well.And model of the present invention calculates two mouthfuls of horizontal gas
The relative error for the critical liquid carrying flow rate that well obtains is respectively 5.4% and 7.6%, is below 10%, with several common calculating moulds
Type is compared, and precision improves 10.5%~62.8%, and the on-site actual situations goodness of fit better than the accuracy of several classical models
It is higher, and then the reliability of new model is demonstrated, the critical liquid carrying flow rate of accurately prediction level gas well is can be used to, is effectively instructed
The lasting Efficient Development at scene.
Based on same inventive concept, a kind of horizontal gas well is additionally provided in the embodiment of the present invention, and continuously to take liquid critical flow true
Device is determined, as described in the following examples.Since horizontal gas well continuously takes the principle that liquid critical flow determining device solves the problems, such as
Liquid critical flow is continuously taken to horizontal gas well and determines that method is similar, therefore horizontal gas well continuously takes liquid critical flow determining device
Implementation may refer to horizontal gas well and continuously take the implementation that liquid critical flow determines method, and overlaps will not be repeated.It is following to be made
, the combination of the software and/or hardware of predetermined function may be implemented in term " unit " or " module ".Although following embodiment
Described device preferably realized with software, but the realization of the combination of hardware or software and hardware be also may be simultaneously
It is contemplated.
Fig. 5 is that the horizontal gas well of the embodiment of the present invention continuously takes a kind of structure diagram of liquid critical flow determining device, such as
Shown in Fig. 5, including:
Horizontal gas well critical liquid carrying flow rate computation model establishes module 501, for establishing horizontal gas well critical liquid carrying flow rate
Computation model;
Horizontal gas well actual production data acquisition module 502, for obtaining horizontal gas well actual production data;
Horizontal gas well continuously takes liquid critical flow determining module 503, for being calculated based on horizontal gas well critical liquid carrying flow rate
Model determines that horizontal gas well continuously takes liquid critical flow according to the horizontal gas well creation data.
The structure is illustrated below.
It is specifically used for when it is implemented, the horizontal gas well critical liquid carrying flow rate computation model establishes module 501:
Horizontal gas well critical liquid carrying flow rate computation model is established as follows:
Horizontal gas well liquid film physical model is established, the parameter of horizontal gas well liquid film physical model is obtained;
According to the parameter of the horizontal well liquid film physical model, establish that horizontal gas well is critical to take flow velocity computation model;
Flow velocity computation model is taken according to the horizontal gas well is critical, horizontal gas well critical liquid carrying flow rate is established and calculates mould
Type.
When it is implemented, the parameter of the horizontal gas well liquid film physical model includes:Conversion flow velocity, the liquid film of central gas stream
Flow velocity, the tubing string internal diameter of horizontal gas well, the shearing stress between central gas stream and liquid film, cutting between liquid film and the tube wall of horizontal gas well
The critical thickness of stress, the equivalent diameter of central gas stream and air-flow top liquid film.
It is specifically used for when it is implemented, the horizontal gas well critical liquid carrying flow rate computation model establishes module 501:
Establish that horizontal gas well is critical to take flow velocity computation model according to following formula:
Wherein, VsgIndicate the conversion flow velocity of central gas stream, m/s;μlFor the viscosity of liquid phase, mPas;ρlIndicate liquid
Density, kg/m3;D indicates tubing string internal diameter, m;ClIndicate the lifting factor, dimensionless;ρgIndicate the density of gas, kg/m3;G is indicated
Acceleration of gravity, m/s2;K is modifying factor,NBFor Bond number, dimensionless,
σ is gas-liquid surface tension, N/m.
It is specifically used for when it is implemented, the horizontal gas well critical liquid carrying flow rate computation model establishes module 501:
Horizontal gas well critical liquid carrying flow rate computation model is established according to following formula:
Wherein, QscIndicate that horizontal gas well continuously takes liquid critical flow, 104m3/d;A indicates oil pipe sectional area, m2;pwIt indicates
Pressure, MPa;VsgIndicate the conversion flow velocity of central gas stream, m/s;Z indicates the compressed coefficient of gas under certain temperature and pressure, nothing
Dimension;T indicates temperature, K.
In conclusion the present invention carries out force analysis to liquid film in horizontal gas well, pressure difference and gas turbulence effect etc. are considered
The influence of factor can be used for accurate calculated level gas well and connect to establish horizontal gas well critical liquid carrying flow rate computation model
It is continuous to take liquid critical flow, solve the problems, such as to calculate in existing method error it is big this, for improving the final of in-situ horizontal gas well
Recovery ratio has certain directive function.
It should be understood by those skilled in the art that, the embodiment of the present invention can be provided as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the present invention
Apply the form of example.Moreover, the present invention can be used in one or more wherein include computer usable program code computer
The computer program production implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.)
The form of product.
The present invention be with reference to according to the method for the embodiment of the present invention, the flow of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that can be realized by computer program instructions every first-class in flowchart and/or the block diagram
The combination of flow and/or box in journey and/or box and flowchart and/or the block diagram.These computer programs can be provided
Instruct the processor of all-purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine so that the instruction executed by computer or the processor of other programmable data processing devices is generated for real
The device for the function of being specified in present one flow of flow chart or one box of multiple flows and/or block diagram or multiple boxes.
These computer program instructions, which may also be stored in, can guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works so that instruction generation stored in the computer readable memory includes referring to
Enable the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one box of block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device so that count
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, in computer or
The instruction executed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one
The step of function of being specified in a box or multiple boxes.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the embodiment of the present invention can have various modifications and variations.All within the spirits and principles of the present invention, made by
Any modification, equivalent substitution, improvement and etc. should all be included in the protection scope of the present invention.
Claims (10)
1. a kind of horizontal gas well continuously takes liquid critical flow and determines method, which is characterized in that including:
Establish horizontal gas well critical liquid carrying flow rate computation model;
Obtain horizontal gas well actual production data;
Based on horizontal gas well critical liquid carrying flow rate computation model, according to the horizontal gas well creation data, determine that horizontal gas well connects
It is continuous to take liquid critical flow.
2. horizontal gas well as described in claim 1 continuously takes liquid critical flow and determines method, which is characterized in that according to such as lower section
Formula establishes horizontal gas well critical liquid carrying flow rate computation model:
Horizontal gas well liquid film physical model is established, the parameter of horizontal gas well liquid film physical model is obtained;
According to the parameter of the horizontal well liquid film physical model, establish that horizontal gas well is critical to take flow velocity computation model;
Flow velocity computation model is taken according to the horizontal gas well is critical, establishes horizontal gas well critical liquid carrying flow rate computation model.
3. horizontal gas well as claimed in claim 2 continuously takes liquid critical flow and determines method, which is characterized in that the horizontal gas
The parameter of well liquid film physical model includes:Central gas stream conversion flow velocity, the flow velocity of liquid film, horizontal gas well tubing string internal diameter, in
The equivalent diameter of shearing stress, central gas stream between shearing stress, liquid film between motive stream and liquid film and the tube wall of horizontal gas well is gentle
Flow the critical thickness of top liquid film.
4. horizontal gas well as claimed in claim 3 continuously takes liquid critical flow and determines method, which is characterized in that according to following public affairs
Formula establishes that horizontal gas well is critical to take flow velocity computation model:
Wherein, VsgIndicate the conversion flow velocity of central gas stream, m/s;μlFor the viscosity of liquid phase, mPas;ρlIndicate the density of liquid,
kg/m3;D indicates tubing string internal diameter, m;ClIndicate the lifting factor, dimensionless;ρgIndicate the density of gas, kg/m3;G indicates that gravity adds
Speed, m/s2;K is modifying factor,NBFor Bond number, dimensionless,σ is gas
Liquid surface tension, N/m.
5. horizontal gas well as claimed in claim 4 continuously takes liquid critical flow and determines method, which is characterized in that according to following public affairs
Formula establishes horizontal gas well critical liquid carrying flow rate computation model:
Wherein, QscIndicate that horizontal gas well continuously takes liquid critical flow, 104m3/d;A indicates oil pipe sectional area, m2;pwIndicate pressure,
MPa;VsgIndicate the conversion flow velocity of central gas stream, m/s;Z indicates the compressed coefficient of gas under certain temperature and pressure, dimensionless;T
Indicate temperature, K.
6. a kind of horizontal gas well continuously takes liquid critical flow determining device, which is characterized in that including:
Horizontal gas well critical liquid carrying flow rate computation model establishes module, and mould is calculated for establishing horizontal gas well critical liquid carrying flow rate
Type;
Horizontal gas well actual production data acquisition module, for obtaining horizontal gas well actual production data;
Horizontal gas well continuously takes liquid critical flow determining module, for being based on horizontal gas well critical liquid carrying flow rate computation model, root
According to the horizontal gas well creation data, determine that horizontal gas well continuously takes liquid critical flow.
7. horizontal gas well as claimed in claim 6 continuously takes liquid critical flow determining device, which is characterized in that the horizontal gas
Well critical liquid carrying flow rate computation model is established module and is specifically used for:
Horizontal gas well critical liquid carrying flow rate computation model is established as follows:
Horizontal gas well liquid film physical model is established, the parameter of horizontal gas well liquid film physical model is obtained;
According to the parameter of the horizontal well liquid film physical model, establish that horizontal gas well is critical to take flow velocity computation model;
Flow velocity computation model is taken according to the horizontal gas well is critical, establishes horizontal gas well critical liquid carrying flow rate computation model.
8. horizontal gas well as claimed in claim 7 continuously takes liquid critical flow determining device, which is characterized in that the horizontal gas
The parameter of well liquid film physical model includes:Central gas stream conversion flow velocity, the flow velocity of liquid film, horizontal gas well tubing string internal diameter, in
The equivalent diameter of shearing stress, central gas stream between shearing stress, liquid film between motive stream and liquid film and the tube wall of horizontal gas well is gentle
Flow the critical thickness of top liquid film.
9. horizontal gas well as claimed in claim 8 continuously takes liquid critical flow determining device, which is characterized in that the horizontal gas
Well critical liquid carrying flow rate computation model is established module and is specifically used for:
Establish that horizontal gas well is critical to take flow velocity computation model according to following formula:
Wherein, VsgIndicate the conversion flow velocity of central gas stream, m/s;μlFor the viscosity of liquid phase, mPas;ρlIndicate the density of liquid,
kg/m3;D indicates tubing string internal diameter, m;ClIndicate the lifting factor, dimensionless;ρgIndicate the density of gas, kg/m3;G indicates that gravity adds
Speed, m/s2;K is modifying factor,NBFor Bond number, dimensionless,σ is gas
Liquid surface tension, N/m.
10. horizontal gas well as claimed in claim 9 continuously takes liquid critical flow determining device, which is characterized in that the level
Gas well critical liquid carrying flow rate computation model is established module and is specifically used for:
Horizontal gas well critical liquid carrying flow rate computation model is established according to following formula:
Wherein, QscIndicate that horizontal gas well continuously takes liquid critical flow, 104m3/d;A indicates oil pipe sectional area, m2;pwIndicate pressure,
MPa;VsgIndicate the conversion flow velocity of central gas stream, m/s;Z indicates the compressed coefficient of gas under certain temperature and pressure, dimensionless;T
Indicate temperature, K.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810419162.7A CN108763657A (en) | 2018-05-04 | 2018-05-04 | Horizontal gas well continuously takes liquid critical flow and determines method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810419162.7A CN108763657A (en) | 2018-05-04 | 2018-05-04 | Horizontal gas well continuously takes liquid critical flow and determines method and device |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108763657A true CN108763657A (en) | 2018-11-06 |
Family
ID=64010011
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810419162.7A Pending CN108763657A (en) | 2018-05-04 | 2018-05-04 | Horizontal gas well continuously takes liquid critical flow and determines method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108763657A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110580656A (en) * | 2019-10-10 | 2019-12-17 | 中国石油集团川庆钻探工程有限公司 | method for predicting underground limit liquid carrying flow of horizontal well cylinder |
CN111723531A (en) * | 2020-05-25 | 2020-09-29 | 长江大学 | Shale gas reservoir horizontal well critical liquid-carrying gas production rate prediction method and device |
CN112145153A (en) * | 2020-11-26 | 2020-12-29 | 西南石油大学 | Shale gas well gas production pipe column size optimization method |
CN113468826A (en) * | 2021-06-17 | 2021-10-01 | 西南石油大学 | Shale gas horizontal well critical liquid carrying flow prediction method based on real liquid film distribution |
CN113931621A (en) * | 2020-07-14 | 2022-01-14 | 中国石油天然气股份有限公司 | Method and device for determining gas well accumulated liquid information and storage medium |
CN114233272A (en) * | 2021-12-17 | 2022-03-25 | 西安安森智能仪器股份有限公司 | Intelligent exploitation control method and device for natural gas well |
CN114263450A (en) * | 2021-12-31 | 2022-04-01 | 中国石油天然气股份有限公司 | Method, device and terminal for optimizing flowback system of tight gas well |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5983164A (en) * | 1997-02-25 | 1999-11-09 | Stella, Llc | Method and apparatus for measuring and controlling the flow of natural gas from gas wells |
US6332111B1 (en) * | 1997-09-24 | 2001-12-18 | Bechtel Bwxt Idaho, Llc | Method and system for measuring multiphase flow using multiple pressure differentials |
US20020069022A1 (en) * | 1997-09-24 | 2002-06-06 | Fincke James R. | Multiphase flow calculation software |
CN103670352A (en) * | 2012-09-18 | 2014-03-26 | 中国石油天然气股份有限公司 | Automatic control method for removal of gas well deposition liquid |
CN103867184A (en) * | 2014-02-10 | 2014-06-18 | 中国石油天然气股份有限公司 | Gas well critical liquid carrying flow rate determining method and device |
US20160265331A1 (en) * | 2011-11-04 | 2016-09-15 | Schlumberger Technology Corporation | Modeling of interaction of hydraulic fractures in complex fracture networks |
CN106570273A (en) * | 2016-11-04 | 2017-04-19 | 中国石油天然气股份有限公司 | Building method and application of three parameter-gas well drainage gas recovery technology optimizing model |
CN106600442A (en) * | 2016-12-06 | 2017-04-26 | 中国海洋石油总公司 | Critical liquid carrying condition optimization method in deep water gas well testing |
CN107975363A (en) * | 2017-10-20 | 2018-05-01 | 中国石油天然气集团公司 | Condensate gas well takes liquid critical flow Forecasting Methodology and device |
-
2018
- 2018-05-04 CN CN201810419162.7A patent/CN108763657A/en active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5983164A (en) * | 1997-02-25 | 1999-11-09 | Stella, Llc | Method and apparatus for measuring and controlling the flow of natural gas from gas wells |
US6332111B1 (en) * | 1997-09-24 | 2001-12-18 | Bechtel Bwxt Idaho, Llc | Method and system for measuring multiphase flow using multiple pressure differentials |
US20020069022A1 (en) * | 1997-09-24 | 2002-06-06 | Fincke James R. | Multiphase flow calculation software |
US20160265331A1 (en) * | 2011-11-04 | 2016-09-15 | Schlumberger Technology Corporation | Modeling of interaction of hydraulic fractures in complex fracture networks |
CN103670352A (en) * | 2012-09-18 | 2014-03-26 | 中国石油天然气股份有限公司 | Automatic control method for removal of gas well deposition liquid |
CN103867184A (en) * | 2014-02-10 | 2014-06-18 | 中国石油天然气股份有限公司 | Gas well critical liquid carrying flow rate determining method and device |
CN106570273A (en) * | 2016-11-04 | 2017-04-19 | 中国石油天然气股份有限公司 | Building method and application of three parameter-gas well drainage gas recovery technology optimizing model |
CN106600442A (en) * | 2016-12-06 | 2017-04-26 | 中国海洋石油总公司 | Critical liquid carrying condition optimization method in deep water gas well testing |
CN107975363A (en) * | 2017-10-20 | 2018-05-01 | 中国石油天然气集团公司 | Condensate gas well takes liquid critical flow Forecasting Methodology and device |
Non-Patent Citations (1)
Title |
---|
陈德春: "定向气井临界携液流量预测新模型", 《天然气工业》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110580656A (en) * | 2019-10-10 | 2019-12-17 | 中国石油集团川庆钻探工程有限公司 | method for predicting underground limit liquid carrying flow of horizontal well cylinder |
CN110580656B (en) * | 2019-10-10 | 2022-11-22 | 中国石油集团川庆钻探工程有限公司 | Horizontal well cylinder underground limit liquid carrying flow prediction method |
CN111723531A (en) * | 2020-05-25 | 2020-09-29 | 长江大学 | Shale gas reservoir horizontal well critical liquid-carrying gas production rate prediction method and device |
CN113931621A (en) * | 2020-07-14 | 2022-01-14 | 中国石油天然气股份有限公司 | Method and device for determining gas well accumulated liquid information and storage medium |
CN113931621B (en) * | 2020-07-14 | 2023-08-22 | 中国石油天然气股份有限公司 | Method, device and storage medium for determining gas well hydrops information |
CN112145153A (en) * | 2020-11-26 | 2020-12-29 | 西南石油大学 | Shale gas well gas production pipe column size optimization method |
CN113468826A (en) * | 2021-06-17 | 2021-10-01 | 西南石油大学 | Shale gas horizontal well critical liquid carrying flow prediction method based on real liquid film distribution |
CN113468826B (en) * | 2021-06-17 | 2024-05-28 | 西南石油大学 | Shale gas horizontal well critical carrying flow prediction method based on real liquid film distribution |
CN114233272A (en) * | 2021-12-17 | 2022-03-25 | 西安安森智能仪器股份有限公司 | Intelligent exploitation control method and device for natural gas well |
CN114233272B (en) * | 2021-12-17 | 2023-09-22 | 西安安森智能仪器股份有限公司 | Intelligent exploitation control method and device for natural gas well |
CN114263450A (en) * | 2021-12-31 | 2022-04-01 | 中国石油天然气股份有限公司 | Method, device and terminal for optimizing flowback system of tight gas well |
CN114263450B (en) * | 2021-12-31 | 2024-05-28 | 中国石油天然气股份有限公司 | Flowback system optimization method, flowback system optimization device and flowback system optimization terminal applicable to tight gas well |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108763657A (en) | Horizontal gas well continuously takes liquid critical flow and determines method and device | |
CN107045671A (en) | Water-producing gas well hydrops Risk Forecast Method | |
CN105626036B (en) | A kind of reasonable Liquid output reservoir engineering calculation method of determining oil reservoir | |
US11105187B2 (en) | Apparatus and method for simulating and/or controlling fluid injection | |
CN110866321B (en) | Accumulated liquid judging method suitable for oil-gas-water co-production horizontal well | |
CN104504604B (en) | A kind of method of qualitative Wellbore of Gas Wells hydrops | |
WO2019242296A1 (en) | Method for obtaining relationship model between water injection well and oil production well, and method for determining production rate and water injection amount | |
US9200511B2 (en) | Enhanced dynamic well model for reservoir pressure determination | |
CN105160071B (en) | A kind of suitable gas-liquid is the same as the method for discrimination of production horizontal well underground working | |
CN106156439A (en) | A kind of heavy crude reservoir dissolved gas drive method for numerical simulation with foam oil phenomenon | |
CN109033518A (en) | The water breakthrough time prediction technique and device of bottom water gas condensate reservoir | |
Nikoofard et al. | Design and comparison of constrained mpc with pid controller for heave disturbance attenuation in offshore managed pressure drilling systems | |
CN108830410A (en) | The water breakthrough time prediction technique and device of hyposmosis bottomwater gas field | |
CN104809363B (en) | Ternary composite driving optimization method based on the approximate Dynamic Programmings of BK | |
CN111506865B (en) | Shale gas well production fitting and predicting method | |
WO2019206463A1 (en) | Improved flow measurement | |
GB2611428A (en) | A method of modelling a production well | |
CN112836438A (en) | Method and device for predicting water breakthrough time of gas well of heterogeneous edge water-gas reservoir | |
Pedersen et al. | Model predictive control of flow and pressure in underbalanced drilling | |
Saeid et al. | Numerical investigations of drilling mud flow characteristics in vertical well | |
CN113236203B (en) | Water invasion dynamic production allocation method for carbonate rock with water-gas reservoir | |
CN115217467A (en) | Method and device for determining breakthrough time of bottom water of multi-medium bottom water gas reservoir horizontal well | |
CN114169204A (en) | Sand prevention opportunity determination method for offshore oil and gas field development and production | |
Sutton et al. | Investigation of gas carryover with a downward liquid flow | |
CN111206919A (en) | Long-well-section high-yield gas well reservoir section wellbore pressure calculation method |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181106 |
|
RJ01 | Rejection of invention patent application after publication |