CN110688612A - Multi-producing-layer oil well yield prediction method based on temperature logging data - Google Patents
Multi-producing-layer oil well yield prediction method based on temperature logging data Download PDFInfo
- Publication number
- CN110688612A CN110688612A CN201910843143.1A CN201910843143A CN110688612A CN 110688612 A CN110688612 A CN 110688612A CN 201910843143 A CN201910843143 A CN 201910843143A CN 110688612 A CN110688612 A CN 110688612A
- Authority
- CN
- China
- Prior art keywords
- production
- oil well
- zone
- calculation
- well
- 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.)
- Granted
Links
- 239000003129 oil well Substances 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000004519 manufacturing process Methods 0.000 claims abstract description 31
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 238000012546 transfer Methods 0.000 claims abstract description 15
- 230000008878 coupling Effects 0.000 claims abstract description 5
- 238000010168 coupling process Methods 0.000 claims abstract description 5
- 238000005859 coupling reaction Methods 0.000 claims abstract description 5
- 239000012530 fluid Substances 0.000 claims description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 238000002347 injection Methods 0.000 claims description 5
- 239000007924 injection Substances 0.000 claims description 5
- 238000009826 distribution Methods 0.000 abstract description 16
- 238000012360 testing method Methods 0.000 abstract description 11
- 238000004088 simulation Methods 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000035699 permeability Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Business, Economics & Management (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Operations Research (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Development Economics (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- Entrepreneurship & Innovation (AREA)
- Agronomy & Crop Science (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a multi-producing zone oil well yield prediction method based on temperature logging data. Aiming at an oil well with a plurality of production layers, firstly, a comprehensive heat transfer model of the oil well with the plurality of production layers is established by analyzing the heat and mass transfer characteristics of a shaft, a stratum, a reservoir and coupling of the plurality of production layers and the shaft, then, on the basis of the established model, the difference of temperature change rules at different production layers is utilized, actual temperature logging data is substituted into a calculation model to calculate the flow distribution condition in the shaft by an inversion method, and the yield of each production layer is further calculated, and the problem of insufficient prediction precision of a single production layer caused by the adoption of an average distribution method when a traditional well testing method deals with the plurality of production layers is solved by utilizing the method; not only can the total production of the oil well be obtained, but also the production of a single production zone of a multi-production-zone oil well can be accurately predicted.
Description
Technical Field
The invention relates to the technical field of numerical reservoir simulation and temperature well testing, in particular to a method for predicting the yield of a multi-producing-layer oil well.
Background
The yield of oil reservoirs is the most important parameter in petroleum engineering, and the yield is closely related to the actual flowing process of fluids such as oil, water and the like in the production zone of an oil field during production, and most oil wells contain a plurality of production zones. Engineering often determines well production and other reservoir parameters by employing well testing methods. The well testing method is widely applied because the object is data convenient to measure in an oil well, and is not limited by technical conditions in data collection or data analysis, and a plurality of different well testing methods such as a pressure recovery well testing method, a pressure drop well testing method, a pulse well testing method and the like are developed at present.
However, most of oil wells contain a plurality of production zones, and most of the existing well testing methods generally adopt an average distribution method to predict the yield of the oil wells in the plurality of production zones, so that the overall yield can be predicted only, the yield prediction precision of a single production zone is obviously insufficient, and the comprehensive yield prediction effect of the oil wells is further influenced.
Disclosure of Invention
The invention provides a multi-producing zone oil well yield prediction method based on temperature logging data, aiming at the problem that the conventional well testing method is insufficient in single-producing zone productivity prediction precision of a multi-producing zone oil well.
The method comprises the steps of firstly analyzing the heat and mass transfer characteristics of a shaft, a stratum, a reservoir and a multi-zone-shaft coupling, and establishing a heat and mass transfer calculation model of a multi-zone oil well; and calculating the accurate yield of each production layer by using an inversion method based on the actual temperature logging data.
The method for predicting the oil well yield of the multiple producing zones based on the theory comprises the following main steps: (1) firstly, aiming at an oil well needing to be evaluated, collecting various parameters of a stratum, a reservoir, a shaft and injection conditions of the oil well, and providing basic parameters for an established numerical simulation calculation model; (2) collecting actual temperature logging data after the preparation of the calculation model is finished so as to provide a basis for inversion calculation; (3) and finally, giving a value range to the fluid flow of each production layer, and carrying out inversion calculation on the flow of each production layer according to a random approximation inversion method to obtain a prediction result.
The step (1) comprises the following steps: the required basic parameters include, in particular, formation thermal conductivity, specific heat capacity, geothermal gradient, formation density, reservoir permeability, reservoir porosity, well wall parameters, number and perforation positions of primary producing zones, well depth, injection temperature, pressure, and the like.
The step (2) comprises the following steps: the required actual temperature logging data are the temperature values of the fluid in the well and the corresponding depth values.
The step (3) comprises the following steps: the flow rate value range of each production zone can be very rough, then a group of flow rate distribution is randomly generated in the value range according to a random approximate inversion method, the flow rate distribution is substituted into a calculation model to obtain simulated borehole fluid temperature distribution, the simulated borehole fluid temperature distribution is compared with actual temperature logging data to calculate the root mean square error, then the random value of the flow rate is calculated in the next round, all the calculated root mean square error data are compared when the calculation is carried out to the specified round number, and the group with the minimum root mean square error is taken as the optimal prediction result.
The beneficial technical effects of the invention are as follows:
1. the invention solves the problem that the prediction precision of a single production zone is insufficient due to the adoption of an average distribution method when the traditional well testing method is used for dealing with a plurality of production zones.
2. The invention can not only obtain the total yield of the oil well, but also accurately predict the yield of a single production zone of a multi-production-zone oil well.
Drawings
FIG. 1 is a flow chart of the implementation steps of the stochastic approximation inversion method used in the present invention.
Fig. 2(a) is a comparison graph of the inversion result of the bulk temperature distribution value obtained in the embodiment and the actual temperature logging data.
Fig. 2(b) is a comparison graph of the inversion result of the temperature distribution value of the reservoir section obtained in the embodiment and the actual temperature logging data.
Fig. 3 is a reservoir water uptake distribution graph plotted against the flow distribution values of the various producing zones according to an embodiment of the present invention.
Detailed Description
The invention will now be further described by way of example with reference to the accompanying drawings.
The specific implementation mode is as follows: the embodiment provides the prediction of the oil well yield of the multi-producing zone based on the temperature logging data, and the numerical simulation calculation model of the oil well of the multi-producing zone established by the invention and the inversion of the temperature logging data predict the yield of each producing zone.
The numerical simulation calculation model of the multi-producing zone oil well, which is established by the invention, mainly comprises a relatively independent shaft model and a reservoir model which are connected through a flowing and thermal boundary condition, and is used as a stratum heat transfer similar to the boundary condition to calculate by utilizing a one-dimensional radial model and consider unsteady state:
where f (t) is the formation unsteady heat transfer function.
The wellbore model control equation is:
the momentum equation:
energy equation:
wherein a frictional resistance term tau is introduced in consideration of friction between the fluid and the borehole wallfAnd A is the cross section area of the injection well, the left side of the equation of the energy is the unit depth heat flow of the fluid in the well penetrating through the well wall, and the right side is an unsteady phase and z-direction convection term.
The reservoir model control equation is:
the continuous equation:
wherein SiSaturation of the water or oil phase, saturation of the water and oil phasesAnd the sum is zero,is the porosity of the porous medium.
The seepage velocity of the fluid in the reservoir during oil production is relatively slow, which can be considered as satisfying Darcy's law of seepage, then the momentum equation can be written as:
wherein KiFor each phase fluid permeability, the permeability of the fluid in an actual oilfield reservoir is affected by a number of factors, such as temperature, pressure, saturation, etc., which are set herein for simplicity as a function of the saturation of each phase:
Ki＝K0*Si 2
in the heat exchange process, the convective heat transfer between the fluids and the heat conduction between the fluids and the medium are considered, and then the energy equation is:
where the subscript por represents the reservoir porous media.
Aiming at a specific example in the method for predicting the yield of the multi-producing zone, a target oil well is firstly adapted to a numerical simulation calculation model, and parameters such as formation heat conductivity, specific heat capacity, injection temperature and pressure are selected from values given by an actually measured data source and some common engineering parameter values, namely lambdae＝2.7W/(m·K)；cp＝900J/(kg·K)；Tin＝273.15K；PinThe oil well depth is 3500m with 5 main production layers under 11MPa, and the flow rate of each production layer is set to be 0-30 (m) according to the actual measurement conditions3And d), then substituting the actual temperature logging data into the model to calculate according to the inversion calculation flow shown in the figure 1.
Table 1 inversion results of flow distribution
| Unit: (m)3/d) | Inversion result 1 | |
|
Mean value of |
| Reservoir 1 flow | 15.95 | 13.13 | 12.98 | 12.89 |
| Reservoir 2 flow | 11.97 | 12.41 | 12.57 | 12.61 |
| Reservoir 3 flow | 9.41 | 9.48 | 9.75 | 9.90 |
| Reservoir 4 flow | 3.10 | 3.56 | 3.47 | 3.58 |
| Reservoir 5 flow | 4.21 | 5.10 | 4.92 | 4.73 |
| Total flow rate | 44.64 | 43.68 | 43.69 | 43.71 |
Three groups of inversion results obtained according to the situation of the specific embodiment are shown in table 1, and it can be seen from table 1 that the inversion results converge on a flow distribution, and from comparing the temperature distribution inversion results shown in fig. 2 with the actual temperature logging data, it can be seen that the error between the temperature distribution calculation value corresponding to the inversion results and the actual logging data can be kept in a small range.
Claims (5)
1. A multi-producing zone oil well yield prediction method based on temperature logging data is characterized by comprising the following operation steps:
(1) analyzing the heat and mass transfer characteristics of a shaft, a stratum, a reservoir and a multi-zone-shaft coupling, and establishing a heat and mass transfer calculation model of a multi-zone oil well;
(2) and calculating the accurate yield of each production layer by using an inversion method based on the actual temperature logging data.
2. The method of claim 1 for predicting well production from a multi-zone well based on temperature log data, wherein: the heat and mass transfer calculation model of the multi-zone oil well comprehensively considers the respective calculation models of a shaft, a stratum and a reservoir and the flow-thermal coupling boundary conditions among the calculation models; the shaft is considered to be axial one-dimensional heat transfer, flow and one-dimensional heat transfer between the radial direction and the stratum, and the reservoir is considered to be two-dimensional porous medium two-phase seepage and heat transfer.
3. The method of claim 1 for predicting well production from a multi-zone well based on temperature log data, wherein: the parameters required by the heat and mass transfer calculation model of the multi-zone oil well comprise: the system comprises a shaft injection parameter, a well wall structure parameter, a formation thermophysical parameter, a reservoir seepage characteristic parameter and a thermophysical parameter.
4. The method of claim 2, wherein the method comprises the steps of: in the heat and mass transfer calculation model of the multi-production-layer oil well, a co-located grid semi-implicit algorithm is adopted for partial flow thermal coupling calculation of the reservoir.
5. The method of claim 2, wherein the method comprises the steps of: the inversion calculation method for the flow of each production zone adopts a random approximation method, and comprises the following specific steps:
(1) firstly, determining a fuzzy flow value range according to actual working conditions;
(2) randomly taking values of the flow of each layer in the range;
(3) substituting the random value-taking result into a calculation model to calculate the corresponding wellbore fluid temperature-depth relation;
(4) comparing the actual temperature logging data with the calculation data to calculate the root mean square error between the actual temperature logging data and the calculation data and recording the root mean square error;
(5) after the steps (2), (3) and (4) are repeated for a certain number of times, comparing the root mean square error calculated for each time;
(6) and taking a group of flow values with the minimum root mean square error as an inversion calculation result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910843143.1A CN110688612B (en) | 2019-09-06 | 2019-09-06 | Multi-producing-layer oil well yield prediction method based on temperature logging data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910843143.1A CN110688612B (en) | 2019-09-06 | 2019-09-06 | Multi-producing-layer oil well yield prediction method based on temperature logging data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110688612A true CN110688612A (en) | 2020-01-14 |
| CN110688612B CN110688612B (en) | 2021-07-06 |
Family
ID=69107863
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910843143.1A Active CN110688612B (en) | 2019-09-06 | 2019-09-06 | Multi-producing-layer oil well yield prediction method based on temperature logging data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110688612B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111648764A (en) * | 2020-07-20 | 2020-09-11 | 西南石油大学 | Interpretation and evaluation method for underground distributed temperature monitoring output profile of multilayer gas reservoir |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101421489A (en) * | 2006-04-07 | 2009-04-29 | 国际壳牌研究有限公司 | Method for production metering of oil wells |
| CN101906966A (en) * | 2010-07-16 | 2010-12-08 | 中国石油天然气股份有限公司 | Method and device for forecasting reservoir yield |
| CN105574318A (en) * | 2014-11-07 | 2016-05-11 | 中国石油化工股份有限公司 | Oil well production automatic splitting device and method |
| KR101706245B1 (en) * | 2015-09-14 | 2017-02-14 | 동아대학교 산학협력단 | Method for controlling production rate using artificial neural network in digital oil field |
| CN106869914A (en) * | 2017-03-09 | 2017-06-20 | 长江大学 | The PRODUCTION FORECASTING METHODS that seepage flow is coupled with flowing in pit shaft in a kind of oil reservoir |
| CN107563899A (en) * | 2016-06-30 | 2018-01-09 | 中国石油天然气股份有限公司 | Oil ＆ Gas Productivity Forecasting Methodology and device |
| CN108547610A (en) * | 2018-02-07 | 2018-09-18 | 中国石油天然气股份有限公司 | The determination method and apparatus of horizontal productivity under volume fracturing |
| WO2018212674A1 (en) * | 2017-05-18 | 2018-11-22 | Владимир Георгиевич КИРЯЧЕК | Method of deriving hydrocarbons from oil-prone kerogen-rich formations and technological complex. |
-
2019
- 2019-09-06 CN CN201910843143.1A patent/CN110688612B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101421489A (en) * | 2006-04-07 | 2009-04-29 | 国际壳牌研究有限公司 | Method for production metering of oil wells |
| CN101906966A (en) * | 2010-07-16 | 2010-12-08 | 中国石油天然气股份有限公司 | Method and device for forecasting reservoir yield |
| CN105574318A (en) * | 2014-11-07 | 2016-05-11 | 中国石油化工股份有限公司 | Oil well production automatic splitting device and method |
| KR101706245B1 (en) * | 2015-09-14 | 2017-02-14 | 동아대학교 산학협력단 | Method for controlling production rate using artificial neural network in digital oil field |
| CN107563899A (en) * | 2016-06-30 | 2018-01-09 | 中国石油天然气股份有限公司 | Oil ＆ Gas Productivity Forecasting Methodology and device |
| CN106869914A (en) * | 2017-03-09 | 2017-06-20 | 长江大学 | The PRODUCTION FORECASTING METHODS that seepage flow is coupled with flowing in pit shaft in a kind of oil reservoir |
| WO2018212674A1 (en) * | 2017-05-18 | 2018-11-22 | Владимир Георгиевич КИРЯЧЕК | Method of deriving hydrocarbons from oil-prone kerogen-rich formations and technological complex. |
| CN108547610A (en) * | 2018-02-07 | 2018-09-18 | 中国石油天然气股份有限公司 | The determination method and apparatus of horizontal productivity under volume fracturing |
Non-Patent Citations (3)
| Title |
|---|
| LIANG HUI-ZHEN 等: "The BP Network Study of the Time Series Overrolling Model for Forecasting the Oilfield Output", 《2009 IITA INTERNATIONAL CONFERENCE ON SERVICES SCIENCE, MANAGEMENT AND ENGINEERING》 * |
| 年永乐: "稠油热采传热和基于测温数据的油井参数预测分析", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
| 杨顺辉 等: "多层合采智能井井筒温度场预测模型及应用", 《石油钻探技术》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111648764A (en) * | 2020-07-20 | 2020-09-11 | 西南石油大学 | Interpretation and evaluation method for underground distributed temperature monitoring output profile of multilayer gas reservoir |
| CN111648764B (en) * | 2020-07-20 | 2021-03-19 | 西南石油大学 | Interpretation and evaluation method for underground distributed temperature monitoring output profile of multilayer gas reservoir |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110688612B (en) | 2021-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kujawa et al. | Utilization of existing deep geological wells for acquisitions of geothermal energy | |
| Saeid et al. | A prototype design model for deep low-enthalpy hydrothermal systems | |
| García et al. | Estimation of temperatures in geothermal wells during circulation and shut-in in the presence of lost circulation | |
| BR112015005136B1 (en) | METHOD AND APPARATUS TO PROVIDE PRESSURE WINDOW FOR DRILLING DRILLING HOLE | |
| WO2005035944A1 (en) | System and method for determining a flow profile in a deviated injection well | |
| He et al. | Pressure-transient behavior of multisegment horizontal wells with nonuniform production: theory and case study | |
| Shook et al. | Use of tracers and temperature to estimate fracture surface area for EGS reservoirs | |
| Beier et al. | Thermal response tests on deep boreholes through multiple ground layers | |
| Wołoszyn et al. | Sensitivity analysis of efficiency thermal energy storage on selected rock mass and grout parameters using design of experiment method | |
| CN102682195B (en) | Semisubmersible platform transient state bored shaft temperature computation method | |
| Buchardt et al. | Estimating the basal melt rate at NorthGRIP using a Monte Carlo technique | |
| CN111814100A (en) | Method for dynamically simulating circulating temperature change of marine subsea pump lifting drilling system without marine riser | |
| JP6230025B2 (en) | Thermal response test analysis method and analysis program | |
| CN104847314A (en) | High-temperature high-pressure oil and gas vertical-well single-phase flow perforation well completion parameter optimization method | |
| CN110688612B (en) | Multi-producing-layer oil well yield prediction method based on temperature logging data | |
| Vaferi et al. | Simulation of dynamic pressure response of finite gas reservoirs experiencing time varying flux in the external boundary | |
| CN114580100B (en) | Method and device for calculating full wellbore pressure of fractured horizontal well and computer readable storage medium | |
| Roux et al. | An improved approach to estimating true reservoir temperature from transient temperature data | |
| CN108319771A (en) | A kind of low-permeability oil deposit temperature in wellbore computational methods | |
| CN106934075A (en) | Drilling fluid density determines that method and static equal yield density determine method | |
| Chevarunotai | Analytical Models for Flowing-Fluid Temperature Distribution in Single-Phase Oil Reservoirs Accounting for Joule-Thomson Effect | |
| US9891343B2 (en) | Method and system for mapping a three-dimensional structure using motes | |
| Prytula et al. | REVEALING PATTERNS IN THE INFLUENCE OF VARIABLE PERMEABILITY OF WELL BOTTOMHOLE ZONES ON THE OPERATIONAL MODES OF UNDERGROUND GAS STORAGE FACILITIES. | |
| Akhlaghi et al. | Effect of Underground Convective Flows on the Performance of Ground Source Heat Exchanger Systems | |
| Akhlaghi et al. | Modeling the Temperature Behavior of Ground Source Heat Pump Systems |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |