Underground Gas Storage Process Optimisation with Respect to Reservoir Parameters and Production Equipment

: The storage of natural gas in geological structures such as depleted ﬁelds, aquifers and salt caverns plays an important role in a gas supply system as it balances the ﬂuctuation of gas demand and price. Hydraulic loss due to ﬂuid ﬂow through gas storage production equipment and an interfering effect from nonequal productivity index of storage wells may have an important inﬂuence on gas storage performance. An integrated mathematical model is developed based on underground gas storage facility production data. Using this model, the hydraulic loss is determined. A real test case that consists of a gas storage reservoir linked to the surface facility is analysed. The mathematical model uses an experimentally determined pressure drop coefﬁcient in chokes. The base case scenario created using real gas storage facility data enables the achievement of a good history match with the given parameters of the gas storage reservoir. Using the history match simulation case as an initial scenario (a base case), two different scenarios are created to determine the injection and withdrawal performance of the gas storage ﬁeld. The results indicate that the pressure drop in chokes, when fully open as a constraints in an underground gas storage facility, has a signiﬁcant impact on gas storage operations and deliverability.


Introduction
Underground gas storage facilities (UGS) represent an important segment of gas infrastructure. Their primary function is to ensure a reliable and efficient supply of natural gas to private and public consumers. The storage of natural gas is a critical component of the natural gas supply chain that must be optimised to balance the gas supply and demand [1]. Natural gas is stored within underground geological formations (reservoirs), and in most cases, a former producing field is used and transformed into a storage facility [2][3][4][5]. In some cases, storage facilities have been developed in aquifers that did not previously hold gas [6] or salt caverns [7]. Withdrawal and injection capacity are important technical features of every gas storage facility, and they are a function of reservoir pressure, physical properties of the geological formation and the number of wells. Maximum gas injection capacity is usually sustainable throughout the whole storage cycle, but maximum withdrawal storage capacity is only sustainable to a certain point of reservoir pressure, which is continuously reduced by gas depletion. Since the UGS facility is likely to work intensively during the winter months, it is desirable that the UGS maximum withdrawal capacity is sustainable at lower values of reservoir pressure. This is important because UGS contributes to the safety of gas supply in winter periods during peak consumption loads. In this paper, we present findings from a study conducted to evaluate the impact of technological factors affecting the performance of the UGS storage facility and the sustainability of gas withdrawal capacity. This is generally related to hydraulic losses that occur during the flow of gas through production equipment. The influence of hydraulic losses is mainly emphasised with a high gas velocity that occurs at low reservoir pressure at the end of the gas withdrawal cycle. Parts of production equipment that have been examined and tested are wellhead chokes and production strings. The mentioned problems and their impact on underground gas storage facility performance have not been studied in a systematic way until now. In order to examine this issue, a full-scale mathematical model of a gas storage facility was developed, which enables the continuous monitoring of physical parameters during flow simulation from the reservoir to processing parts of the facility. Mathematical models are a powerful and heuristic tool for investigating complex systems. Once a model has been developed, it incorporates a large bulk of data and knowledge, and it can be used to make predictions regarding how a system (UGS facility) would respond under various conditions. With a model, it was possible to examine the impact of the production equipment's internal diameter on the sustainability of maximum gas storage capacity. Discharge coefficient values for wellhead and reduction station chokes at 100% opening were determined experimentally, and the obtained values are incorporated into the mathematical model.

Literature Review
Earlier research regarding the development of mathematical models of underground gas storage applied individual models of the production system. Process optimisation using a single model is based on the parameter's analysis of a single production process (i.e., gas flow through reservoir rock). Alternatively, the interdependence between the reservoir, the surface pipeline network, and the process facility is observed simultaneously with an Integrated UGS production (mathematical) model. The integrated model helps to determine some specific problems that are undetectable using stand-alone model simulation (e.g., a flow-assurance problem during CO 2 injection [8]). Khodri et al. (1997) also used the same methodology to model gas storage reservoirs and address the problems of improving the peak flow rate, maximising the performance and deliverability and decision on drilling new wells [9]. Brown et al. (1999) used a homogeneous volumetric tank model to solve various storage problems, such as determining the minimum withdrawal time (or the peak withdrawal rate), and determining the minimum injection time at peak injection rate [10]. McVay and Spivey (2001) used a reservoir simulation model to examine a practical range of compression horsepower, the number of wells and amount of cushion gas and determined the combination that minimises the objective function, defined as the initial development cost [11]. These studies did not take into account the overall integration of the gas storage system components. Kuncir et al. (2003) used a tank model with a single-phase gas flow coupled with a surface flow system to evaluate the development and expansion options for a gas storage reservoir. They developed an optimisation code to evaluate different horsepower and well count configurations to meet a demand schedule [12]. Bagci and Ozturk (2007) provided a simulation study using both 3-D full-field black oil reservoir models to evaluate and demonstrate the feasibility of creating UGS in a depleted gas field. Their effort was to fit the wellhead pressure, average reservoir pressure and cumulative production generated by the numerical model with actual field data gathered from 1998 to 2002 [13]. Moradi (2009) simulated the effect of natural gas injection and production on variations in gas relative permeability and condensate saturation within near and far away areas of the wellbore by using the commercial simulator Eclipse 300. It has been found that the injection of a higher volume of gas in the first injection plan could prevent condensate from forming near the wellbore area [14]. Sun et al. (2017) provided a geomechanical simulation method for assessing the risk of caprock integrity failure and fault leakage of large UGS in a produced gas field during long-term operations. The large produced gas field located in China was converted to a UGS facility in 2013. In this study, the authors built a static regional-scale 3D geomechanical model through the integrated analysis of geologic, seismic, well drilling, logging and experimental rock data and established a plan for the design of UGS in complex, faulted, depleted gas field [15]. The objective of this work is to propose an integrated mathematical model as a solution and investigation tool for specific UGS facility problems.

Problem Statement
During gas flow through process equipment, friction losses are a function of the system geometry, fluid properties and flow rate. Their influence on gas storage operations is particularly emphasised at low reservoir pressure (usually at the end of the gas withdrawal cycle) due to high gas velocity in pipelines. The parts of the production equipment where most of the total friction loss occurs are the wellhead and the reduction station chokes. The reason for their application in UGS is the control of fluid flow from the reservoir by reducing its pressure. Flow through the chokes could be critical or subcritical depending on the gas velocity through the choke. The "critical flow" occurs when the fluid velocity at restriction (choke) is equal to the velocity of sound in that medium. When the velocity is less than the velocity of sound, it is supercritical flow. The basic difference between critical and subcritical flows is how the flow rate through the restriction is affected by a pressure drop across the restriction. During "subcritical flow", the flow rate is related to the pressure drop across the restriction. At "critical flow", the gas rate is only related to the upstream pressure [16]. At the beginning of the gas storage withdrawal cycle, typically, critical flow occurs at chokes during high-pressure differential between the reservoir pressure and the gas transport pressure. As reservoir pressure declines, critical flow leans towards subcritical flow. Eventually, with a further decrease in reservoir pressure, chokes no longer have a function in regulation and become an unnecessary flow restriction that has an influence on gas storage physical characteristics even in a fully open state. This paper presents a method to investigate the impact of major production equipment hydraulic losses (chokes and production string diameter) on gas storage operation and performance.

Methodology and Experiments
Within this paper, an integrated mathematical model of the UGS facility was applied. With a model, it is possible to simulate gas flow from porous rock through a production string and surface equipment in real-time. A model consists of several single mathematical models (a reservoir model, a well model, a choke model, and a surface process model) integrated into one unit. Using a model, by varying certain parameters within technically acceptable limits, it is possible to examine the effect of hydraulic losses on gas storage withdrawal capacity duration. Field experiments were also carried out to investigate pressure drop in chokes assuming isothermal, compressible and one-dimensional gas flow. A particular focus was placed on the discharge coefficient (Cd) number determination at a fully open cross-section area of the choke valve. The experiments were performed for various gas inlet pressures and temperatures, and the obtained Cd value was integrated into the mathematical model.

Mathematical Model of Underground Gas Storage
A mathematical model of the UGS facility was built and validated based on field and reservoir data from the Croatia UGS facility, collected during production and the injection storage working cycle. Once a model has been validated, it is ready to be used both as a predictive and as a surveillance tool that can honour a variety of system constraints and contract scenarios. Data uncertainty is tested through sensitivity analysis. The mentioned UGS facility in Croatia was developed with 23 working wells and has working storage volume of 553 × 10 6 m 3 [17]. The model consists of several interconnected separate models that mathematically describe fluid flow of multiple production segments. A reservoir model was built with Eclipse software. IPM software was used to model fluid flow in wells and process pipelines; interconnection and integration of individual models was archived using IPM RESOLVE program (see Figure 1). 3. A model of the surface gathering system is matched on real dynamic production data (pressure and temperature drops through pipes and restrictions).
Using this model, it is possible to monitor a pressure drop from well perforations to the main gathering system gas separator as the last node of the storage system. With this level of integration, it is possible to monitor the influence of any flow restriction on gas storage operation simultaneously. Withdrawal cycle starts at the beginning of October and ends at the beginning of April. The flow-line network (surface model) connecting storage wells to the first-stage gas separator is represented by the steady-state thermodynamic model, where input data vary with time because they are determined by the reservoir model. In this model, GAP is used to calculate pressure loss in the pipeline. The pressure drop is solved through back-calculation using the Weymouth equation [18]. Surface-pipeline model assumes the following assumptions: The model's adjustment to the measurement results are as follows [17]:

1.
A reservoir 3D model in Eclipse is calibrated on production history data and injection/withdrawal cycles.

2.
A storage well model is calibrated with data obtained from modified isochronal tests performed on gas storage wells during injection and withdrawal cycles. 3.
A model of the surface gathering system is matched on real dynamic production data (pressure and temperature drops through pipes and restrictions).
Using this model, it is possible to monitor a pressure drop from well perforations to the main gathering system gas separator as the last node of the storage system. With this level of integration, it is possible to monitor the influence of any flow restriction on gas storage operation simultaneously.
Input settings and restrictions applied within the RESOLVE model are as follows: The flow-line network (surface model) connecting storage wells to the first-stage gas separator is represented by the steady-state thermodynamic model, where input data vary with time because they are determined by the reservoir model. In this model, GAP is used to calculate pressure loss in the pipeline. The pressure drop is solved through back-calculation using the Weymouth equation [18]. Surface-pipeline model assumes the following assumptions: -Heat transfer to the ground is assumed to be at steady state and the same material is assumed for all pipes; - The outlet pressure and temperature are calculated based on incremental energy and mass balances. Figures 2 and 3 show a surface model built in GAP software, both for gas withdrawal and the injection process. Chokes (SAP#1, SAP#2 and SAP#3) are present only in the withdrawal process model. In Table 1, pipeline data are indicated in more detail. − Heat transfer to the ground is assumed to be at steady state and the same material is assumed for all pipes; − The outlet pressure and temperature are calculated based on incremental energy and mass balances.  Table 1, pipeline data are indicated in more detail.

Experimental Determination of Choke Discharge Coefficient
Discharge coefficient was determined analytically using measured data of pressure drop across fully opened chokes during gas flow. In order to eliminate any background noise that could cause data error, measured data were acquired after flow stabilisation (see Figure 4). Gas volume flow was measured in position after the choke with two flow measuring devices (ultrasonic meter and orifice) according to Figure 4. The accuracy of the measurement was around ±0.5%. Measurement data incorporate a full range of dynamic pressures and volume flows that may occur during the gas storage operation.
The discharge coefficient for wellhead and process chokes is calculated using Equation (1) and the measurement results of the pressure drop across chokes during the gas storage withdrawal cycle [17]. where: C d = discharge coefficient; Z = compressibility factor; T 1 = gas temperature before nozzle, K; M a = molar mass of air, 28.966 × 10 −3 kg/mol; R = general gas constant, 8.3145 J K −1 mol −1 ; γ g = relative gas density; T sc = standard temperature, K; p sc = standard pressure, Pa; D = diameter of the nozzle opening, m; κ = adiabatic exponent, defined by the specific heat ratio c p /c v ; p 1 = gas pressure before nozzle, Pa; p 2 = gas pressure after nozzle, Pa.
Specific heat values (c p , c v ) for natural gas and z factor at referent PT conditions are calculated using Peng and Robinson Equation of state (EOS). noise that could cause data error, measured data were acquired after flow stabilisation (see Figure 4). Gas volume flow was measured in position after the choke with two flow measuring devices (ultrasonic meter and orifice) according to Figure 4. The accuracy of the measurement was around ±0.5%. Measurement data incorporate a full range of dynamic pressures and volume flows that may occur during the gas storage operation. The discharge coefficient for wellhead and process chokes is calculated using Equation (1) and the measurement results of the pressure drop across chokes during the gas storage withdrawal cycle [17].

Model Simulation Scenarios
The influence of restriction diameters on gas storage parameters was investigated with two scenarios whose simulation results were compared with the basic scenario (matched on real gas storage data). The base case scenario was designed on the basis of gas storage production/injection data from all existing wells. Twenty-two wells are working wells (producer/injection), and two of them are water injectors. The base case scenario will be run until 2022 and used for comparison with all other cases. The purpose of the base case model simulation scenario is to simulate the actual flow parameters of the UGS facility. To prove the hypothesis, within the mathematical model, two test scenarios were created. The scenarios' results were then analysed and compared with the base case scenario. The optimisation of the injector/producers rates, wellhead chokes, process chokes and production string diameters are the main issues that have been investigated. Maximum permissible reservoir pressure drawdown on all storage wells of 3 MPa was implemented within all the simulation scenarios. In all the scenarios, the end part of the withdrawal simulation model is input-output separator with a fixed pressure of 4 MPa. Additionally, in the injection cycle model, the input/output gas water separator represents the starting point of the model as compressor fixed suction pressure. The two other production scenarios analyse the impact of production equipment sizing on UGS working capabilities. Production equipment diameters for the applied scenarios are shown in Table 2. In all simulation scenarios, production equipment diameters for other wells not mentioned in Table 2 have the same tubing diameter. This also applies to the diameter of the wellhead and the process chokes. All these scenarios will be run until 2022. In scenario 1 (model), the production string (tubing) internal diameter of the selected wells ( Table 2) was set to 0.1016 m. All other gas storage wells have a "base case" internal tubing diameters set to 0.076 m. The wellhead and process reduction station chokes diameters were set to 0.0254 m and 0.054 m. This was completed to investigate the impact of a tubing's internal diameter on gas storage withdrawal capabilities.

Case 2
In scenario 2, the production string (tubing) internal diameter for all gas storage wells was set to 0.076 m. The diameter of the wellhead and the process reduction station chokes were set to 0.05 m and 0.072 m. This was completed to investigate the impact of a choke's internal diameter on gas storage withdrawal capabilities.

Data Presentation and Analysis of Results
There are two sets of results shown in this chapter. First are the results related to the discharge coefficient calculation (when chokes are fully opened). Those results are obtained during pressure drop measurement across a choke. Second are the simulation scenarios' results, and they are shown later in this paper.

Discharge Coefficient Calculation and Results
Combining measurement data and Equation (1), the discharge coefficient for UGS chokes was calculated (Table 3). After obtaining the universal value of Cd, further calculation of the pressure drop on chokes was performed in order to select proper Choke ID for the case 2 model scenario. Pressure drop calculation with the obtained Cd value for different choke diameters is shown later in the text. The average value of discharge coefficient for wellhead and process chokes was 0.76, and that value was implemented in the simulation model.

Influence of Choke Diameter on Head Loss
Using Equation (1) 6 show the interdependence of hydraulic loss and gas flow for a single set of input data. How these losses affect the performance of the gas storage facility could  6 show the interdependence of hydraulic loss and gas flow for a single set of input data. How these losses affect the performance of the gas storage facility could only be determined with an integrated dynamic (mathematical) model. For that purpose, the mentioned scenarios were selected, and the specific model results are presented in this chapter. Gas storage withdrawal capacity strongly depends on reservoir parameters, and the number of working wells and facility construction. At some point during the production cycle, as the reservoir pressure is constantly decreasing, the designed maximum gas rate is no longer attainable. The scenario results show the dependency of maximum withdrawal storage capacity on reservoir pressure for every case scenario. The evident change in the reservoir pressure value at which it is still possible to operate the storage with maximum withdrawal capacity is a direct indicator of changes made in scenarios 1 and 2 (in relation to the baseline scenario). Other relevant data that indirectly affect the withdrawal capacity are also shown in the figures below (pressure drop across chokes trends, drawdown trends, etc.).

Figures 5 and
Simulation results (pressure dependence vs. withdrawal capacity, other related data) for the base case scenario are shown in Figures 7-9. Corresponding simulation results for Scenario 1 are shown in Figures 10 and 11 and for Scenario 2 in Figures 12-15. relation to the baseline scenario). Other relevant data that indirectly affect the withdrawal capacity are also shown in the figures below (pressure drop across chokes trends, drawdown trends, etc.).
Simulation results (pressure dependence vs. withdrawal capacity, other related data) for the base case scenario are shown in Figures 7-9. Corresponding simulation results for Scenario 1 are shown in Figures 10 and 11 and for Scenario 2 in Figures 12-15.     In the base case scenario, the sustainability of the maximum withdrawal capacity decreases with a value of reservoir pressure below 116 bar. This value matches the real gas storage performance data. For the settings of scenario 1, the maximum withdrawal capacity is sustainable at a lower reservoir pressure (regarding the base case scenario) at 112 bar. This is a direct indicator of changes introduced in Scenario 1, and their influence on storage withdrawal capacity performance.  For the settings of scenario 1, the maximum withdrawal capacity is sustainable at a lower reservoir pressure (regarding the base case scenario) at 112 bar. This is a direct indicator of changes introduced in Scenario 1, and their influence on storage withdrawal capacity performance. For production wells that were subject to a change in tubing diameter (Scenario 1), the difference between the formation pressure and the bottom hole hydrostatic pressure (Drawdown) is shown. Due to the lower hydraulic losses occurring on the production string (comparing with the base case scenario shown in Figure 8), a higher average drawdown per well was achieved during the withdrawal storage cycle.  Due to the demand for accuracy in the calculation of the pressure drop across the choke, determination of the discharge coefficient (Cd) for a specific wellhead and process chokes (with a circular cross-section) was necessary. By incorporating these values in a mathematical (simulation) model, it was possible to investigate the influence of hydraulic losses at chokes on gas storage technical capabilities. Scenario 2 simulation results are shown in Figures 12-15. The minimum reservoir pressure value at which maximum withdrawal capacity is possible is given in Figure 12. Figure 13 shows the working wells' drawdown during a simulation run.
It should be emphasised that a longer sustainability of UGS withdrawal capacity was achieved with a low-pressure drop through the chokes, as seen in Figures 14 and 15. For the settings of scenario 2, the maximum withdrawal capacity is sustainable at lower reservoir pressure (regarding the base case scenario and Scenario 1) at 102 bar. From this, it can be concluded that scenario 2 has the most impact on gas storage withdrawal capacity performance, and hydraulic losses occurring on restriction (chokes) are not negligible. For production wells that were subject to a change in wellhead choke diameter (Scenario 2), the difference between the formation pressure and the bottom hole hydrostatic pressure (drawdown) is shown in Figure 13. Due to the lower hydraulic losses occurring on wellhead and facility chokes (Figures 14 and 15), a higher average drawdown per well was achieved during the withdrawal storage cycle.   For production wells that were subject to a change in wellhead choke diameter (Scenario 2), the difference between the formation pressure and the bottom hole hydrostatic pressure (drawdown) is shown in Figure 13. Due to the lower hydraulic losses occurring on wellhead and facility chokes (Figures 14 and 15), a higher average drawdown per well was achieved during the withdrawal storage cycle. In Figure 14, the pressure drop values across the wellhead chokes are presented for Scenario 2 chokes during the withdrawal cycle. The pressure drop values are negligible, and thus have a direct influence on the wells' productivity performance.  Figure 15 shows pressure drop values on the reduction station nozzles during the gas storage withdrawal cycle. This pressure drop affects drawdown on all gas storage wells during the production cycle. Pressure drop is negligible, which has a direct influence on the well's productivity performance. A lover pressure drop on the nozzles increases the gas storage deliverability pattern due to lower hydraulic losses.

Discussion
The conducted experimental research combined with mathematical modelling confirmed the impact of hydraulic losses on gas storage capacity sustainability. It has been shown that hydraulic resistance through the underground gas storage chokes in the state of their complete openness significantly affects the gas withdrawal capacity. This is particularly evident at low reservoir pressure during the end of the withdrawal cycle. The investigation results may influence the choice of technical production equipment in underground gas storage as this problem has not yet been systematically investigated.
An additional contribution of this research represents the experimentally determined value of the discharge coefficient for wellhead and process chokes with a circular crosssection area. The established value of Cd for subcritical flow completes the existing results of the two-phase fluid flow research. This paper provides the methodology for an integrated mathematical model of UGS development with real data obtained from an underground gas storage facility in Croatia. The presented concept with the necessary modifications could be applied to all types of underground gas storage.

Conclusions
Based on the conducted research, the following can be concluded:

•
For the subcritical flow of natural gas through the wellhead chokes (type Needle and Seat choke valve), the average value for the discharge coefficient of 0.76 was determined. The obtained coefficient was applied in the mathematical model for pressuredrop calculation on UGS chokes. With a developed mathematical model of the UGS facility, the impact of hydraulic losses (head loss) on withdrawal capacity was observed. In the base case scenario, the sustainability of the maximum withdrawal capacity decreases with a value of reservoir pressure below 116 bar. This value matches the real gas storage performance data.
For the settings of scenario 1, the maximum withdrawal capacity is sustainable at a lower reservoir pressure (regarding the base case scenario) at 112 bar. This is a direct indicator of changes introduced in Scenario 1, and their influence on storage withdrawal capacity performance.
For production wells that were subject to a change in tubing diameter (Scenario 1), the difference between the formation pressure and the bottom hole hydrostatic pressure (Drawdown) is shown. Due to the lower hydraulic losses occurring on the production string (comparing with the base case scenario shown in Figure 8), a higher average drawdown per well was achieved during the withdrawal storage cycle.
Due to the demand for accuracy in the calculation of the pressure drop across the choke, determination of the discharge coefficient (Cd) for a specific wellhead and process chokes (with a circular cross-section) was necessary. By incorporating these values in a mathematical (simulation) model, it was possible to investigate the influence of hydraulic losses at chokes on gas storage technical capabilities. Scenario 2 simulation results are shown in Figures 12-15. The minimum reservoir pressure value at which maximum withdrawal capacity is possible is given in Figure 12. Figure 13 shows the working wells' drawdown during a simulation run.
It should be emphasised that a longer sustainability of UGS withdrawal capacity was achieved with a low-pressure drop through the chokes, as seen in Figures 14 and 15.
For the settings of scenario 2, the maximum withdrawal capacity is sustainable at lower reservoir pressure (regarding the base case scenario and Scenario 1) at 102 bar. From this, it can be concluded that scenario 2 has the most impact on gas storage withdrawal capacity performance, and hydraulic losses occurring on restriction (chokes) are not negligible.
For production wells that were subject to a change in wellhead choke diameter (Scenario 2), the difference between the formation pressure and the bottom hole hydrostatic pressure (drawdown) is shown in Figure 13. Due to the lower hydraulic losses occurring on wellhead and facility chokes (Figures 14 and 15), a higher average drawdown per well was achieved during the withdrawal storage cycle.
In Figure 14, the pressure drop values across the wellhead chokes are presented for Scenario 2 chokes during the withdrawal cycle. The pressure drop values are negligible, and thus have a direct influence on the wells' productivity performance. Figure 15 shows pressure drop values on the reduction station nozzles during the gas storage withdrawal cycle. This pressure drop affects drawdown on all gas storage wells during the production cycle. Pressure drop is negligible, which has a direct influence on the well's productivity performance. A lover pressure drop on the nozzles increases the gas storage deliverability pattern due to lower hydraulic losses.

Discussion
The conducted experimental research combined with mathematical modelling confirmed the impact of hydraulic losses on gas storage capacity sustainability. It has been shown that hydraulic resistance through the underground gas storage chokes in the state of their complete openness significantly affects the gas withdrawal capacity. This is particularly evident at low reservoir pressure during the end of the withdrawal cycle. The investigation results may influence the choice of technical production equipment in underground gas storage as this problem has not yet been systematically investigated.
An additional contribution of this research represents the experimentally determined value of the discharge coefficient for wellhead and process chokes with a circular crosssection area. The established value of Cd for subcritical flow completes the existing results of the two-phase fluid flow research. This paper provides the methodology for an integrated mathematical model of UGS development with real data obtained from an underground gas storage facility in Croatia. The presented concept with the necessary modifications could be applied to all types of underground gas storage.

Conclusions
Based on the conducted research, the following can be concluded:

•
For the subcritical flow of natural gas through the wellhead chokes (type Needle and Seat choke valve), the average value for the discharge coefficient of 0.76 was determined. The obtained coefficient was applied in the mathematical model for pressure-drop calculation on UGS chokes. With a developed mathematical model of the UGS facility, the impact of hydraulic losses (head loss) on withdrawal capacity was observed. • Larger string diameter results in lower friction between gas particles and the tubing (pipe) wall, enabling higher reservoir pressure drawdown on well perforation level and gas withdrawal capacity extension by 10% related to the initial state (base case model).

•
Hydraulic losses in the chokes during the sub-critical flow of gas in the stage of full choke openness also significantly affect the UGS withdrawal capacity. This is particularly evident at low reservoir pressure due to an increase in gas flow velocity through production equipment. The implementation of a double sized choke diameter (regarding the base case scenario) increases the gas withdrawal capacity extension by 18% for a specific case scenario (Scenario 2). • For the first time, a dynamic mathematical model was used to valorise the impact of hydraulic losses of production equipment on its working capabilities. Without the use of an integrated mathematical model, it would not be possible to systematically examine the above except by directly changing the equipment and measuring, which is an unprofitable and inappropriate procedure.