-□ □-
Запропоновано аналтичн та алгоритмiчнi доповнення до матема-тичног моделi прогнозу роботи газо-конденсатного родовища в газовому режимi. Вдосконалено методику роз-рахун^в прогнозiв основних газоди-намiчних параметрiв та показни^в розробки продуктивних покладiв. Для перевiрки новог моделi використовува-лись фактичн дан родовищ Украгни. Створена програмна реалiзацiя роз-робленог технолого-математичног моделi впроваджена та апробована в газовидобувному виробництвi
Ключовi слова: газоконденсатне родовище, ттенсифшащя припливу, видобуток вуглеводтв, показникирозробки родовища, газодинамiчний роз-
рахунок
□-□
Предложены аналитические и алгоритмические дополнения к математической модели прогноза работы газоконденсатного месторождения в газовом режиме. Усовершенствована методика расчетов прогнозов основных газодинамических параметров продуктивных залежей. Для проверки новой модели использовались данные месторождений Украины. Созданная программная реализация разработанной технолого-математической модели внедрена и апробирована в газодобывающем производстве
Ключевые слова: газоконденсат-ное месторождение, интенсификация притока, добыча углеводородов, показатели разработки месторождения,
газодинамический расчет -□ □-
UDC 661.96:622.276.6
|DOI: 10.15587/1729-4061.2016.80073|
IMPROVEMENT OF TECHNOLOGICAL-MATHEMATICAL MODEL FOR THE MEDIUM-TERM PREDICTION OF THE WORK OF A GAS CONDENSATE FIELD
M. Kutia
Postgraduate Student Division of unconventional energy technologies* E-mail: kmm123@yandex.ru M. Fyk PhD**
E-mail: mfyk@ukr.net O. Kravchenko
PhD
Deputy Director for Science* E-mail: kravipmach@rambler.ru S. Palis
Professor
Department of Systems Management with distributed parameters Otto-von-Guericke-University Magdeburg Universitätsplatz, 2, Magdeburg, Deutschland, 39106 E-mail: stefan.palis@ovgu.de I. Fyk
Doctor of Technical Sciences, Professor** E-mail: mfyk@yandex.ua *A. N. Podgorny Institute for Mechanical Engineering Problems Dm. Pogarskogo str., 2/10, Kharkiv, Ukraine, 61046 **Department of "Oil, gas and condensate" National Technical University «Kharkiv Polytechnic Institute» Bagaliya str., 21, Kharkiv, Ukraine, 61002
1. Introduction
Predicting the gas condensate fields in the gas regime is one of the most important applications of mathematical modeling in the oil and gas sector [1]. This is connected, on the one hand, to the growing need for ergonomic gaseous fuel - natural gas - in the markets for natural resources. From the other hand, the gas regime in particular is the most common in the course of standard development of gas-condensate fields at depths accessible by contemporary wells [2].
The technological innovations of downhole gas production [3] undoubtedly require improvement in mathematical and software provision. In particular, it has always been important to predict gas production for a half century based on the classic gas laws, laws of conservation, motion and balance of substances in the underground layer-reservoir [4]. In the field of mathematical modeling of underground reservoirs, scientific paper [5] addresses many issues, as well
as the detailed practical guidelines [6]. In particular, article [6] demonstrates contemporary methods of mathematical modeling of the work of underground reservoirs.
For the short-term mathematical modeling of separate wells, the methods that are based on one of two principles are used: modeling the work of a field using the method of finite elements [6] and the method for compiling material balance of a gas field [1-3].
The problems, connected with the plecement of wells, geometric configuration of the well bottom, which may change after capital repairs and additional exploration of the field, require more detailed technological-mathematical modeling, optimization of arrangement of the points or horizontal lines of the well bottom [7], account of inaccuracy in geological models [8], multilayer development [9]. This extended multifaceted problematics requires further examination with the application of comprehensive technological-mathematical modeling. When the joint gas-dynamic
©
prediction of parameters of development of geometrically complex gas-condensate reservoir and multi-row multi-bottom wells over the short term is necessary, they pass to the multistage program study by special working groups. For developing each gas field, it is necessary to create a separate project of development and industrial operation, which contains geological-technological model, mechanical-technological model of wells, geological-field model of the bottom portion of the well, mechanical-technological model of wellhead equipment.
If the prediction over the medium-term period (5-10 years) is necessary for the fields, which are at the final stage of operation, it is possible to abandon a number of expensive design and scientific research in favor of accurate gas-dynamic modeling of reservoirs and wells according to the material balance method. If it is necessary to use comprehensive software program packages or operational-engineering programs, they pass to original software-mathematical designs of the comprehensive-systems type.
That is why it appears relevant to study the totality of gas-dynamic processes in the extractive exploration based on applying the methods of material-energy balances and using mathematical and integration software tools.
2. Literature review and problem statement
At present, modeling the work of a field is one of the contemporary technologies, applied in the development of productive zones and gas production, in the assessment of recoverable reserves. It is used for the comparison of different mechanisms of production, providing for a base of economic analysis of the potential scenarios of the field's work. The forecasts of mode parameters of gas-condensate fields obtained on this technological-mathematical base are exceptionally important for business planning in the gas extracting companies.
In practice, accurate medium-term prediction of parameters of development of gas-condensate fields proves difficult. Such prediction must consider a large number of special features and limitations of the technologies, planned to be applied, with a short and medium periods of implementation-realization. One of the simple but expensive methods to solve the questions and problems of such prediction is the acquisition of the software and hardware complexes from Halliburton and Schlumberger companies or their complexes-analogs. But when technological peculiarities emerge in the gas-extraction production, connected to the innovative methods of production intensification, the research institutes nevertheless prefer programming in the classic mathematical program packages. The existence in the mathematical packages of the integrated environment of development of the program code and computations with the use of scaled and relational databases, the tools of animated visualization, tying to the actual design of wells, is convenient [2, 3, 6, 9, 10].
Using the balance methods of predicting the work of a gas field, engineers have to accept a number of assumptions, which influence accuracy of the final result. Carrying out the analysis, it is possible to solve the problems, connected with the control over field development, to compare the design indices of development to actual data, to predict the production of gas and condensate, to track the most significant factors that influence the deviation of the planned and actual indicators of the field development [8].
Taking into account the complexity of geometry of the productive horizons, the ambiguity in interpretation of results of geophysical research into characteristics of the well bottom zone, as well as the change of these indices over time, the verification and estimation of relevance of these assumptions is possible only by way of retrospective analysis. This leads to the fact that a large number of variables that change over time are introduced into the model in the form of constant values, the result of which is the growing inaccuracy in the calculations: in addition to possible initial error, there appears the accumulated error [11]. Such factors do not make it possible to conduct precise estimation of the projected corrections of regime parameters of the field performance, which is especially important before conducting the intensification of gas inflow (with the injection of dry gas or using the exhaustion regime
[12]), when it is not possible, or undesirable, to stop wells for a long period for examining their work.
The described below mathematical and conceptual computer model with logical continuation of the ideas from [6, 8, 11-13, 15] partially solves similar problems owing to the creation of the new algorithm of calculation using the base of initial variables and their dynamic correction in the process of calculation. The correction of the key set of the gas-dynamic parameters made it possible to considerably increase the accuracy in predicting performance of gas fields. Here it is important to note the adaptive flexibility of the technological-mathematical model proposed not only to the averaged regime parameters of reservoir and physical-chemical properties of gases, but also to different technologies of the intensification used. Successful approaches to the solution of such a flexible program configuration are possible to find in the full integrated licensed packages Eclipse Petrel
[13], but the speed of modeling of the complete cycle of calculations is low due to the complicated interface and manual handling of the program's parameters. Furthermore, authors of the scientific papers [13, 15] failed to simulate in such a set of programs the operation of a two-row well or half thermally insulated column of PCP.
The technological peculiarities of intensification of the field's work require a more detailed account of regime parameters separately on the bottom and wellhead, change in the function of distribution of pressure and temperature not only along the layer, but also the column of a well. The coefficient of thermal conductivity from the water-bearing rocks to the flow of gas also influences the field's work to a large extent. The approach to solving the problem from the side of application of a separate program ANSIS (individual ther-mohydraulic calculation of a well) [14] and the subsequent iterative step by step calculations of flows and balances in the layers and wells - obtaining a significant error of prediction for the period longer than 5 years [15]. The programs of thermal simulation (similar or the analogs of the program ANSIS) or hydraulic modeling (similar or the analogs of the program WellPress) perform excellently, when used separately, but the transfer and copying of data to Eclipse Petrel or another analogous predicting program take too long.
In each of the software-simulating systems - prototypes to the designed system - we found different shortcomings in the course of the represented studies, since such specialized designs solved highly specialized tasks or did not consider a number of important initial data. The classic examples of such software products - ArcGIS, DV-SeisGeo - are explored in [16]. They are very relevant and necessary for the
solution of particular production-geological problems, but the modeling of gas-dynamic processes exceeds the scope of functionality.
Rather high accuracy was demonstrated by the calculations of comprehensive problems using the programs CGas-Flow, GasConDebit, WellPress, GasPro. The application of a set of such simplified programs of gas-dynamic and hydraulic simulation makes it possible to solve with engineering accuracy a fairly wide circle of issues when predicting gas-condensate fields. While comparing results of using the enumerated list of applied programs with the actual parameters of the verified fields, it was necessary after all to use several programs simultaneously, and then to connect the entire chain of predicting calculations all together. If only the medium-term forecast of the field's work is necessary, and the short term one is not relevant yet, then the time losses to the conceptual and preliminary engineering calculations are not acceptable.
And finally, the important component in the practice of taking account by the empirical dependencies of crystallization of chemical substances in the geological productive rocks [17] that is forgotten even by large corporations - developers of the appropriate specialized software. Authors of the present work propose to add the necessary formula of correction of viscosity and compressibility of the obtained gaseous fluids. In their opinion, this not only successfully complements the analytics of determining the long-term calculations of the predicting parameters [18], improving the accuracy of results of the medium-term forecasts of production, but it will lead to essential saving of fuel-energy resources as a whole [19].
Actually, in each particular case of developing or application of the software and hardware systems, specialized program products with the built-in mathematical model, the purpose must justify investments for the detailing of programming and testing of the adequacy of modeling the basic physical and physical-chemical processes [2-4, 9].
The combination of the growing relevance of specialized scientific-engineering tasks on the base of contemporary technologies of intensification [20], simulation of new geological [21] and hydrodynamic [22] specificity, development of the comprehensive programs of intensification of production at a late stage of operation of the gas-condensate fields [23], training-practical methodologies of prediction under conditions of contemporary technologies of the development of fields [2, 4, 18, 22] became the main stimulus for the development of improved specialized technological-mathematical model and its program realization. One of the relevant technological tasks of well control, for an example, energy-saving technologies, is the control over immersion and surface pumps and compressors based on the monitoring and predicting the regime parameters [26].
Thus, the combination of specialized scientific-engineering problems of the prediction of gas production and required by the account of possible innovation-technological intensification in the production of hydrocarbons leads to the need of creating a more universal mathematical model of the work of a gas field. Essential steps are taken in this aspect by the represented papers [24, 25]. But while [24] paid more attention to modeling the layer systems in dynamics, then [25] presents more accurate simulation of wells. An attempt at partial combination of positive results of the mentioned scientific trends into the uniform modeling space gave birth to a number of new functional properties. Import-
ant new property is the possibility to conduct rapid accurate correct calculations of the regime parameters from the reservoir to the comprehensive gas treatment unit (CGTU). This property must be included by its analytics into the temporal-spatial calculations of the predicting indices of the work of gas-condensate fields. Inclusion and full parametric integration of the more accurate gas-dynamic model layer-well - CGTU with a general mathematical model of temporal-spatial prediction will draw nearer the creation of the maximum-convenient integrated software-analytical environment of predicting the work of gas-condensate fields. Authors made an attempt at its creation and maximum adaptation for the operating fields in the application to the medium-term forecasts for the period of development of gas-condensate fields of up to 10 years.
3. Aim and tasks of the study
The purpose of the work is the modification and verification of the mathematical model of gas inflow and its lifting by the well taking into account a change in viscosity, gas compressibility, heat exchange processes of the system "gas-layer-column PCP" along the way of its motion, integrated into a general updated technological-mathematical model of gas-condensate production.
To accomplish the aim, the following main tasks were set:
- to develop and verify the improved particular mathematical models of gas inflow to the well bottom and its lifting on the example of parametric history of the main indices of development of real fields, history of the regime parameters of the operating wells;
- to integrate proposed applied empirical formulas and the computational-algorithmic techniques (logic-algorithmic bonds) for the standard applications of reservoir-engineering for the medium-term period;
- to develop the program realization of the improved general technological-mathematical model of development of gas-condensate field, which maximally adequately (by the actual history of the indices of development of test fields) describes the medium-term predictions of the main regime-technological parameters;
- to carry out verification of the program realization of the mathematical model in the operating field under conditions of specific measures for intensification of the well production.
4. Description of peculiarities of mathematical and algorithmic model
A general mathematical model of development of a gas field is based on the system of equations that describe the plane-radial inflow of gas to the well and vertical transport of gas along the pump and compressor pipes. The use of the equation of material balance provides for the continuity of calculation [1, 2, 7].
In the development of mathematical model, we used the following equations:
- equation of material balance of productive reservoir;
- equation of gas inflow to the well bottom;
- equation of gas transport by the lifting column of pump and compressor pipes (Adamov equation);
- equation of classic calculation and refinement of the filtration hydrodynamic coefficients of the layer;
- equation of refinement of dynamic viscosity of gas in the layer and the PCP column;
- equation of refinement of the coefficient of gas compressibility in the layer and the PCP column;
- equation of refinement of the coefficient of hydraulic resistance;
- the modified Adamov equation with the correction for nonisothermic lifting by the PCP column.
The special feature of mathematical model is taking account of chronological change in molar mass of the mixture in the modified nonisothermic Adamov equation [15] and the Starling-Ellington equation of viscosity.
A general form of the applied function of dynamic viscosity |(P, T, M) of Starling-Ellington should be demonstrated in more detail: it is substantially dependent not only on pressure and temperature in the layer but also on actual molar mass of the mixture:
, (9.41 + 0.02-M)-(1.8 ■ T)
u(P,T,M) = A-^—
v ; (209 +19 ■ M +1.8 ■ T)1C
T)107
xEXP
3.5 + 5478+0.01M I x 1.8 ■ T
P 103
Zskv(P,T) 83143 ■ T v 7 M
sj 2.4-0.2Î3.5+I547T+0.01M
(1)
Zskv (P,T) =
0.1P
0.4 x log(X)+0.73 T
(3)
At the calculated step, we take not the layer pressure but the mean dynamic pressure P between the starting value P1 and the final P2 in the productive zone along the axis of filtration in the following form:
p=2 3
Pr
Pl+P
(4)
27
Formula (4) is successfully applied in the practice of pipeline hydraulics, and, based on the hydraulic analogies at the distribution of pressures along the axis of the flow; this particular refining analytical construction is selected. In formula (2), one may see an example of application (4) for averaging the pressure in the section of the inflow between the layer pressure Ppl and the pressure on the bottom of the well P1. The latter novation made it possible to conduct additional correction of the parameters at each step.
The algorithmic construction of technological-mathematical model is in detail represented by Authors in [11], and the previous mathematical modification of the system of equations of nonisothermic lifting by the inclined well - in [15]. Let us consider the conceptual-computer part of modeling in more detail. Hierarchically, computer model consists of three main building blocks: regulating unit, the block of constructing the base of initial data and the calculation nucleus of the program (Fig. 1).
where ^ is the dynamic viscosity; P is the mean pressure; T is the mean temperature; M is the molar mass of the mixture; Pnk, Tnk are the pseudo critical pressure and temperature of gas; Zskv is the coefficient of gas compressibility (function from pressure and temperature).
It is evident from equation (1) that the formation of condensate in the layer and corresponding decrease in molar mass requires the correction of viscosity, which is included in the classic formula of calculation of linear coefficient of hydrodynamic resistance of productive layer at the planeradial filtration Af:
Af =
Zskv
P2
Pi +P
,T
2
Pi2
P +P
,T,M
2
P T
n k h Tt
xln
'R X
vRc
(2)
where | is the dynamic viscosity (function from three parameters); Pat, Tat are the atmospheric pressure and temperature; Ppl, Pj, P2 are the layer, bottom and well head pressures; k is the permeability of layer; h is the bulk of layer; Rk is the radius of power circuit; Rc is the radius of well at the bottom.
One can see two parameters in equation (2), which need correction at each step of calculations: dynamic viscosity and gas compressibility coefficient. When calculating the compressibility coefficient Zskv, Authors suggest introduction into the calculation formula of the Platonov- Gurevitch approximation that has the form:
Calculation nucleus
Regulating unit (adaptation)
Base of initial data
Fig. 1. Conceptual block diagram of the proposed computer technological-mathematical model of development of a gas-condensate field
The base of initial data (Fig. 1) is automatically formed after the line-by-line (for each well or a group of wells) introduction of the first layer of initial variables and point indication of the possible projected corrections, namely:
- the drained gas reserves; the gas, obtained prior to the beginning of calculation; the starting layer pressure;
- layer temperature; the mean temperature gradient; the component composition of gas;
- depth to the middle of the interval of perforation; PCP diameter; coefficients of filtration resistances A and B;
- the maximum permissible debit; the maximum permissible depression to the layer;
- the minimally permissible operating pressure (minimally permissible pressure on the well head);
- PCP roughness coefficient, well-performance coefficient.
As a result of the block's work, a three-dimensional matrix is formed, where initial data for the calculation of the j-th month are represented in each horizontal layer, and in the vertical - of the i-th well (a group of wells). Initial data of each well (a group of wells) are recorded in the line with the index i,j.
Regulating unit (Fig. 1) is responsible for the readout and line-by-line introduction of initial data from the formed base to the calculation nucleus of the program, closing the calculation on the equation of material balance, correction of initial data for each layer, according to the calculated
p
x
x
current indices, and formation of the resulting tables. The regulating and adaptation blocks jointly complement the classic algorithmic nucleus, described in the scientific literature [1, 13, 16].
In the calculation nucleus (Fig. 1) we perform initial data processing for each well (a group of wells), that are supplied from the regulating unit, and the five-stage calculation of necessary indices:
- at the first stage of calculation, the automated selection of regime (constant debit, constant depression to the layer, constant operating pressure) is conducted;
- at the second stage, the current layer pressure is specified;
- at the third stage, coefficients for calculating the vertical transport (at the unknown debit) are determined;
- at the fourth and fifth stages, by the iteration technique we determine the inflow of gas, depression to the layer and operating pressure on the well head, as well as the refinement of variables in the Adamov equation (at the known debit).
On finishing the work of calculation nucleus, the following indices return to the regulating unit: debit, depression to the layer, the well bottom pressure, operating pressure on the well head, static pressure on the well head, the refined value of layer pressure, the coefficient of PCP hydraulic resistance, the number of the applied mode of calculation.
After performing the calculations, the regulating unit forms i+1 tables with the predicting indicators of the work on the field as a whole and for each i-th well (a group of wells) separately.
control a change over time in the real radius of well or zone of drainage). The need for conducting correction by viscosity (1), (2) and coefficient of gas compressibility in the layer (3) is obvious, since they significantly influence the final result. Fig. 3 demonstrates the example of dependency of the coefficient of gas compressibility on the layer pressure, obtained when applying equation (3) for different values of layer pressure.
0 10 20 30 40 50 60
Flow rate (thousand m3/day)
Fig. 2. Example of the CHRPCP dependency on the well's debit
4. Correction of current initial data
As was noted above, mathematical model automatically corrects certain initial data, calculating the correction coefficients or recalculating the variable completely. As the first example we examine coefficient of hydraulic resistance of the pump and compressor pipes (CHRPCP). At the first stages of constructing the model, CHRPCP was considered as the constant magnitude or was calculated once. In the process of designing an algorithm for performing the calculation according to the Adamov formula, we made a decision to perform its recalculation for each well (a group of wells at each layer of the work of model. An example of the result of this recalculation is represented in Fig. 2. Graphic dependencies are built by the results of the program realization of mathematical model in the Matlab software.
Correct determining the coefficients of filtration resistance Af and Bf of the quadratic binomial equation of filtration, applied in the proposed model in the classic form of equation of gas inflow with debit Q to the bottom of the well [3, 4, 5, 18]
Fig. 3.
10 15
Layer pressure (MPa)
Example of dependency of coefficient of gas compressibility on the layer pressure
Also analogous is the implemented procedure of correction of coefficients of the filtration resistance of PCP column by multiplication by the correction coefficients (during each contact with the calculation nucleus according to the principle of the escalating multiplication):
Ppl2 -P,2 = Af • Q + Bf • Q2, (5)
where Ppl and P1 are the layer pressure and pressure on the well bottom, also important for conducting accurate predicting calculation. These coefficients describe the geometry of inflow, the physical properties of fluid and filtration characteristics of the layer. It is obvious that to tell that their values are constant and do not change over time is also impossible. To accurately determine some variables, which are part of these coefficients, is fairly problematic (it is difficult to
__K (6)
new prev.step cor' V /
where A,new is the new value of coefficient of hydraulic resistance of the PCP column; ^prevstep is the value of coefficient of hydraulic resistance of the PCP column at the preceding step of calculation; Kcor is the correction coefficient.
The correction coefficients for the coefficient of filtration resistance of the PCP column PCP also depend at each step of temporal layer on the calculated viscosity (1) and coefficient of gas compressibility (3).
5. Example of application of the model
The predicting calculation of parameters of development according to the initial data about the object of one of the fields of Poltava Region, Ukraine was taken as an example. The calculation was carried out by three variants:
1 - using the classic technique without connecting to the correction modules [18];
2 - with regard to the correction of initial data base relative to the current indices of calculation;
3 - analogous to the second variant but taking account of conducting intensification by comprehensive hydrogen thermobarochemical action on the well bottom zone [2 0], designed for the 25-27th months of the field operation relative to the start of calculation.
Fig. 4. Predicting operating pressure on the well head according to the first and second variants of calculation
Predicting indices by the first variant coincided fully with the results, obtained in the design organization: it is expected that the work of the field must be executed in the regime of constant depression to the layer with the transfer into the mode of constant operating pressure in 158 months of operation after the start of calculation. The second variant shows similar results, but the dynamics of indices somewhat differs. The most essential difference is the fact that the transfer into the mode of constant operating pressure must be executed somewhat later - in the 167-th month from the start of calculation (Fig. 4).
There are also differences in the dynamics of gas extraction (Fig. 5). The difference in the layer pressure (Fig. 7) according to the variants of calculation is due to the difference in the volumes of gas extraction.
The third variant of calculation involves performing the work on intensification. Using technique [11], we performed calculation of the new coefficient of permeability, and, accordingly, the new value of filtration resistance Af. Since the intensification of inflow will be carried out not simultaneously in all wells of the rerservoir, but in series in the course of three months, for the accuracy of calculation the one group of wells is divided into 3 parts (one well in each new group), but the initial calculation data are introduced into the base separately, according to the calendar plan of performing the work on intensification (Fig. 6).
Proposed model - Classic model
Fig. 5. Dynamics of change in the debit according to the variants
Fig. 6. Changing the layer pressure by the first and second variants
160
150
T3 140
S
T3 no
¡3
£ 120
a
O
v n 110
X!
1>
100
Q
90
A A
y
---- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
0
50
250
100 150 200 Months of calculation Without regard to intensification — With regard to intensification
Fig. 7. Daily gas extraction from the field
300
Fig. 7 displays calculated dependencies of daily gas extraction from the field taking into account and without taking into account of the work on intensification of the inflows. After introducing the appropriate changes to initial data, the program automatically adapts the algorithm of calculation for the work with three groups of wells. This action allowed us to more accurately calculate the extraction of gas during the period of performing the work on intensification, which leads to the more precise prediction of economic indices of the work of the field. The parameters of inflow, pressure on the well head, the bottom and other indices of work are calculated separately for each well without using so-called "average well", which is reflected in the increased part Fig. 7.
6. Discussion of results of the study of the modeling experiments and verification of applicability of analytics to different technologies of the intensification of wells
In addition to the block of examining the long-term outlook for gas extraction (Fig. 1-7), important are the studies of instantaneous change in the output when using technologies of the intensification of lifting part and the well bottom.
Fig. 8 demonstrates dependency of change in the debit of well on the change in different design and regime-technological
parameters of PCP. The figures display a change in the debit of well at a simultaneous change in two different parameters after conducting measures or measures for the intensification of gas inflow. Fig. 8, a shows the result of technological polishing, ceramic coating, coating with smooth films or with lubricant, polyurethane coating. Different applied thermal insulation introduces corrections to the indices of extraction. Fig. 8, b shows dependency at the discrete-assigned values of roughness of the PCP internal surface. Fig. 8, c analyzes a change in the debit of well with fluctuations of the wellhead pressure, the roughness is taken as discrete by the actual measurements with profilometer. It is known that the wellhead pressure at the finishing stage of development of a gas-condensate field often depends on the regime of operation of booster compressor station (BCS). The jumps in the consumption of gas or in the change of summary extraction continuously change the wellhead pressures. Fig. 8, d analyzes the influence of change in the diameter on the low values in the conditions of the discretely set temperature on the wellhead. The introduction of repairing or correction sleeves to PCP, the use of technologies for intensification with the introduction of pipe into pipe, the jammings of columns somewhat change hydraulically equivalent inside diameter of PCP. The latter fact leads to an abrupt change in the debit of well even at insignificant interferences in the the common flow section.
Dependency of additional debit on the roughness of mside diameter
• wellhead temperature 300 K wellhead temperature 310 K K
-wellheadtemperature 355 K
Dependency of additional debit on the change in the wellhead temperature
5*ic - l*:: 1.3x10
Change in roughness from the accepted (initial). M
Dependency of additional debit an the change in the wellhead pressure
--roughness - — - roughness — ■ — roughness i • • • roughness
0,05 mm ).l mm 0,25 mm
0 1 U: j:/ -n::
Change in pressure on the weUheadiom the accepted(initial with a throttle or with ECS I Pa
-roughness 0,03 mm - - - roughness 0,05 mm - ■ — roughness 0.1 mm • • • * roughness 0,25 mm
r--
Change in temperature on the wellheadlrom the accepted (initial). K b
Dependency of additional debit on the PCP restored diameter
• wellhead temperature 300 K ■ wellhead temperature 310 K
* wellhead temperature 330 K - wellhead temperature 355 K
0 5x10 n 1x10 J 1.5x10 J 2x10 J
Change in diameter from the accepted (initial, with scurf, film or jammed), M
c d
Fig. 8. Dependency of gas debit on the geometric and schematic-technological parameters of the PCP column: a — dependency of the additional debit of well on the change in roughness of the PCP internal surface (in comparison with the normative) at different temperatures on the wellhead; b — dependency of the additional debit of well on the change in temperature on the wellhead (in comparison with the actual average) at different roughness of the PCP internal surface; c — dependency of the additional debit of well on the change in pressure on the wellhead (in comparison with the actual average) at different roughness of the PCP internal surface; d — dependency of the additional debit of well on the change in the equivalent diameter (in comparison with the actual average) at different temperature on the wellhead
a
When they conduct innovation repair of the well and relative-operational measures for the intensification of gas extraction, there occurs a sharp jump of the regime parameters, which corresponds to initial data on the edges of the analyzed dependencies (Fig. 8, a-d). The jump of parameter of daily gas extraction is observed also in Fig. 7.
Further improvements in the model are planned to conduct in two directions. First, creation and usage of the base of highly specialized calculations for the non-standard cases of development (simultaneous and simultaneous-separate operation of multi-layer reservoirs, taking account of interlayer flows, the work of extraction by the intertube space). Second, to conduct the work for the creation of the real two-or three-dimensional corrected general model of a gas field, that is, to start withdrawal from the principles of complete linearization with a partial transition to the calculation, which considers real geology.
Let us draw some conclusions on the conducted research and analysis of graphs (Fig. 1-8):
- the improvement of the algorithms of general mathematical model of the performance of a gas field made it possible to considerably improve the predicting calculation of the indices of development of gas field at the gas regime (to 10-15 % in the longer term);
- it is demonstrated that the mathematical model automatically adapts the logic scheme of the calculation based on the entered initial data;
- the possibility of a particular well's calculation of the work of reservoir is shown;
- the possibility is shown of taking account of different technologies of intensification of the inflow, in particular, the application of comprehensive hydrogen thermalbaro-chemical action on the well's bottom zone, given by the example (Fig. 4-7), the application of PCP with the polymeric coatings (decreasing the roughness of pipes), the application of thermal-insulated PCP, reduction in the operating pressure in the wellhead.
The example of comparative calculation of the indices of work of a gas field according to the classic and proposed techniques revealed that taking account of the change over time of the physical properties of gas and parameters of lifting columns makes it possible to increase the accuracy of calculation of debit by 5-15 %.
7. Conclusions
1. As a result of this work, it was possible to establish possible changes in predicting the work of a gas-condensate field after introduction of the correction coefficients and analytical dependencies by 5-15 % by the total average annual extraction.
2. We created and fixed the program realization of the improved technological-mathematical model of development of a gas-condensate field. The scientific-applied peculiarities of the formulation of problem of research are taken into account, under conditions of contemporary technologies of intensification of gas inflow.
3. We developed the technique of comprehensive account of change in the physical properties of gas at a drop in the layer pressure and the speed of motion of gas by the PCP column. The results of calculations and comparison with production data testify to a considerable increase in accuracy of prediction as compared to the classic model. The automation of the process of selecting the optimal general mode of operation of productive layers and wells makes it possible to decrease effect of the human factor on the accuracy of prediction of natural gas extraction.
4. The technological-mathematical model provides the possibility to solve the extended circle of the applied problems of operational medium-term prediction of the work of gas-condensate fields with the tenfold savings of means for the software-hardware complexes of the short term prediction.
References
1. Mirzadzhanzade, A. H. Osnovy tehnologii dobychi nefti i gaza: Nauchnoe izdanie [Text] / A. H. Mirzadzhanzade, O. L. Kuznecov, K. S. Basniev, Z. S. Aliev. - Moscow: Nedra, 2003. - 880 p.
2. Grishanenko, V. P. Naukovi osnovi vdoskonalennja sistem rozrobki rodovish nafti i gazu [Text] / V. P. Grishanenko, Yu. O. Zarubin, V. M. Doroshenko et. al. - Kiiv, 2014. - 454 p.
3. Zakirov, S. N. Novye principy i tehnologii razrabotki mestorozhdenii nefti i gaza [Text] / S. N. Zakirov, Ye. S. Zakirov, I. S. Zakirov et. al. - Moscow: FGUP «PIK VINITI», 2004. - 521 p.
4. Deik, L. P. Osnovy razrabotki neftjanyh i gazovyh mestorozhdenii [Text] / L. P. Deik. - M.: OOO «Premium Inzhiniring», 2009. - 570 p.
5. Ertekin, T. Basic Applied Reservoir Simulation. Vol. 7 [Text]: monograph Series / T. Ertekin, J. H. Abou-Kassem, G. R. King. - Texas: SPE, Richardson, 2001. - 421 p.
6. Kanevskaja, R. D. Matematicheskoe modelirovanie gidrodinamicheskih processov razrabotki mestorozhdenii uglevodorodov [Text] / R. D. Kanevskaja. - Moscow-Izhevsk: Institut komp'yuternyh issledovanii, 2002. - 140 p.
7. Forouzanfar, F. Joint optimization of number of wells, well locations and controls using a gradient-based algorithm [Text] / F. Forouzanfar, A. C. Reynolds // Elsevier. Chemical engineering research and design. - 2014. - Vol. 92, Issue 7. - P. 1315-1328. doi: 10.1016/j.cherd.2013.11.006
8. Mirzadzhanzade, A. H. Modelirovanie processov neftegazodobychi. Nelineinost', neravnovesnost', neopredelennost'. [Text] / A. H. Mirzadzhanzade, M. M. Hasanov, R. N. Bahtizin. - Izhevsk: «Izhevskii institut komp'yuternyh issledovanii», 2005. - 368 p.
9. Fanchi, J. Principles of Applied Reservoir Simulation. 3rd Edition [Text] / J. Fanchi. - Elsevier GPP, 2006. - 532 p.
10. Nasrabadi, H. Well placement optimization: A survey with special focus on application for gas/gas-condensate reservoirs [Text] / H. Nasrabadi, A. Morales, D. Zhu // Journal of Natural Gas Science and Engineering. - 2012. - Vol. 5. - P. 6-16. doi: 10.1016/ j.jngse.2011.10.002
11. Kutja, M. M. Obshaja model' raschyota pokazatelei razrabotki gazokondensatnogo mestorozhdenija [Text] / M. M. Kutja, K. Yu. Skryl'nik, V. V. Shevchenko // Visnik Nacional'nogo tehnichnogo universitetu "KhPI". Serija: Matematichne modelyuvannja v tehnici ta tehnologijah. - 2014. - Vol. 39 - P. 91-97.
12. Kalugin, Yu. I. Mathematical Modeling and Optimization of Gas-Condensate Field Development [Text] / Yu. I. Kalugin, V. V. Yakovlev, A. Yu. Kalugin //Journal of Natural Gas Science and Engineering. - 2015. - Vol. 27. - P. 1195-1204. doi: 10.1016/j.jngse.2015.09.063
13. Fik, M. I. Osobennosti konceptual'no-tehnologicheskogo pohoda pri otkrytii vtorogo dyhanija neftegazokondensatnyh mestorozhdenii [Text]: konferencija / A. B. Sinyuk, M. I. Fik // Geopetrol' - 2014, 2014. - P. 565-570.
14. Vorkov, A. V. Teplovoi raschet rasteplenija gruntov v priust'evyh zonah skvazhin [Text] / A. V. Vorkov, S. A. Ovchinnikova, M. A. Fedotenko // Izvestija Samarskogo nauchnogo centra Rossiiskoi akademii nauk. - 2013. - Vol. 15, Issue 4-2. - P. 323-326.
15. Fik. M. I. Utochnennja rozrahunku efektivnosti roboti DKS v umovah faktichnih termogradientiv ta suchasnih pokrittiv NKT [Text] / M. I. Fik // Naftogazova promislovist' Ukraini. - 2014. - Vol. 1. - P. 25-28.
16. Aliev, Z. S. Rukovodstvo po proektirovaniyu razrabotki gazovyh i gazoneftjanyh mestorozhdenii [Text] / Z. S. Aliev, V. V. Bondarenko. -Pechora: Pechorskoe vremja, 2002. - 895 p.
17. Mangold, M. Nonlinear model reduction of a continuous fluidized bed crystallizer [Text] / M. Mangold, L. Feng, D. Khlopov, S. Palis, P. Benner, D. Binev, A. Seidel-Morgenstern // Journal of Computational and Applied Mathematics. - 2015. - Vol. 289. - P. 253-266. doi: 10.1016/j.cam.2015.01.028
18. Fyk, M. I. Razrabotka i yekspluatacija neftjanyh i gazovyh mestorozhdenii [Text]: uchebnik / M. I. Fyk, E. I. Hripko; I. M. Fyk (Ed.). -Kharkiv: Folio, 2015. - 399 p.
19. Soldatkov, S. G. Yenergosberegayushie tehnologii utilizacii plastovih poter' gaza pri yekspluatacii PHG [Text] / S. G. Soldatkov, I. V. Marushenko, S. A. Voronov // Gazovaja promyshlennost'. - 2010. - Vol. 9. - P. 63-65.
20. Kravchenko, O. V. An improved technology of a complex influence on productive layers of oil and gas wells [Text] / O. V. Kravchenko, D. A. Veligockii, A. N. Avramenko, R. A. Habibullin // Eastern-European Journal of Enterprise Technologies. - 2014. - Vol. 6, Issue 5 (72). - P. 4-9. doi: 10.15587/1729-4061.2014.29316
21. Chornii, M. I. Geologichni osnovi rozkrittja i viprobuvannja produktivnih plastiv [Text] / O. M. CHornii, I. M. Metoshop, I. M. Kuziv. - Ivano-Frankivs'k, 2013. - 306 p.
22. Gladkov, E. A. Geologicheskoe i gidrodinamicheskoe modelirovanie mestorozhdenii nefti i gaza [Text] / E. A. Gladkov. - Tomskii politehnicheskii universitet, 2012. - 99 p.
23. Fesenko, Y. L. Applied aspects of maintaining gas production in a gas condensate production field at a late stage of operation [Text] / Y. L. Fesenko, S. V. Kryvulia, B. B. Syniuk, M. I. Fyk // NAFTA-GAZ, ROK LXIX. - 2013. - Vol. 10. - P. 744-53.
24. Islam, M. R. Advanced Petroleum Reservoir Simulation. 2nd Edition [Text] / M. R. Islam, M. E. Hossain, H. Moussavizadegan, S. Mustafiz, J. H. Abou-Kassem. - Scrivener Publishing LLC, 2016. - 572 p. doi: 10.1002/9781119038573
25. Khasanov, М. A simple mechanistic model for void fraction and pressure-gradient prediction in vertical and inclined gas/liquid flow [Text] / М. Khasanov, R. Khabibullin, V. Krasnov, A. Pashali, V. Guk. - SPE Production & Operations, 2009. - P. 165-170.
26. Khamitov, R. N. Upravlenie pogruzhnyimi dvigatelyami ustanovok elektrotsentrobezhnyih nasosov po minimumu summarnyih poter [Text] / R. N. Khamitov, A. Yu. Kovalev // Promyishlennaya energetika. - 2011. - Vol. 1. - P. 42-46.