In this paper, the iterative coupling approach is proposed for applications to solving multiphase flow equation systems in reservoir simulation, as it provides a more flexible time-stepping strategy than existing appr...In this paper, the iterative coupling approach is proposed for applications to solving multiphase flow equation systems in reservoir simulation, as it provides a more flexible time-stepping strategy than existing approaches. The iterative method decouples the whole equation systems into pressure and saturation/concentration equations, and then solves them in sequence, implicitly and semi-implicitly. At each time step, a series of iterations are computed, which involve solving linearized equations using specific tolerances that are iteration dependent. Following convergence of subproblems, material balance is checked. Convergence of time steps is based on material balance errors. Key components of the iterative method include phase scaling for deriving a pressure equation and use of several advanced numerical techniques. The iterative model is implemented for parallel computing platforms and shows high parallel efficiency and scalability.展开更多
Petroleum science has made remarkable progress in organic geochemistry and in the research into the theories of petroleum origin, its transport and accumulation. In estimating the oil-gas resources of a basin, the kno...Petroleum science has made remarkable progress in organic geochemistry and in the research into the theories of petroleum origin, its transport and accumulation. In estimating the oil-gas resources of a basin, the knowledge of its evolutionary history and especially the numerical computation of fluid flow and the history of its changes under heat is vital. The mathematical model call be described as a coupled system of nonlinear partial differentical equations with initial-boundary value problems. This thesis, from actual conditions such as the effect of fluid compressibility and the characteristic of large-scal science-engineering computalion. puts forward a kind of characteristic finite difference alternating-direction scheme. Optimal order estimates in L-2 norm are derived for the error in the approximate solutions.展开更多
This paper describes a way of solving the reservoir simulation pressure equation using mulligrid technique. The subroutine MG of four-grid method is presented. The result for 2-D two-phase problem is exactly the same ...This paper describes a way of solving the reservoir simulation pressure equation using mulligrid technique. The subroutine MG of four-grid method is presented. The result for 2-D two-phase problem is exactly the same as that of the SOR method and the CPU time is much less than that of the latter one.展开更多
This paper investigates the deposition of asphaltenes in the porous medium of the studied field in Russia and predicts production profiles based on uncertainty evaluation. This problem can be solved by dynamic modelin...This paper investigates the deposition of asphaltenes in the porous medium of the studied field in Russia and predicts production profiles based on uncertainty evaluation. This problem can be solved by dynamic modeling, during which production profiles are estimated in two scenarios: with and without the activation of the asphaltene option. Calculations are carried out for two development scenarios: field operation under natural depletion and water injection into the aquifer as a reservoir pressure maintenance system. A full-scale compositional reservoir simulation model of the Russian oilfield was created. Within a dynamic simulation, the asphaltene option was activated and the asphaltene behavior in oil and porous medium was tuned according to our own special laboratory experiments. The model was also matched to production historical data, and a pattern model was prepared using the full-scale simulation model. Technological and the asphaltene option parameters were used in sensitivity and an uncertainty evaluation. Furthermore, probable production profiles within a forecast period were estimated. The sensitivity analysis of the pattern model to input parameters of the asphaltene option allowed determining the following heavy-hitters on the objective function: the molar weight of dissolved asphaltenes as a function of pressure, the asphaltene dissociation rate, the asphaltene adsorption coefficient and the critical velocity of oil movement in the reservoir. Under the natural depletion scenario, our simulations show a significant decrease in reservoir pressure and the formation of drawdown cones leading to asphaltene deposition in the bottom-hole area of production wells, decreasing their productivity. Water injection generally allows us to significantly reduce the volume of asphaltene phase transitions and has a positive effect on cumulative oil production. Injecting water into aquifer can keep the formation pressure long above the pressure for asphaltene precipitation, preventing the asphaltene deposition resulted from interaction of oil and water, so this way has higher oil production.展开更多
Simulation study was applied in the development planning of East Unity oilfield, Sudan. A grid consisting of 2 000 cells was constructed. A major challenge of the study wasto evolve a full field development and future...Simulation study was applied in the development planning of East Unity oilfield, Sudan. A grid consisting of 2 000 cells was constructed. A major challenge of the study wasto evolve a full field development and future reservoir management strategy that would ensuremaximum recovery of oil based on well Un51. Simulation shows that Un51 as injection well inAradiebaC would yield better oil recovery than to be production well.展开更多
Krylov subspace projection methods are known to be highly efficient for solving large linear systems. Many different versions arise from different choices to the left and right subspaces. These methods were classified...Krylov subspace projection methods are known to be highly efficient for solving large linear systems. Many different versions arise from different choices to the left and right subspaces. These methods were classified into two groups in terms of the different forms of matrix H-m, the main properties in applications and the new versions of these two types of methods were briefly reviewed, then one of the most efficient versions, GMRES method was applied to oil reservoir simulation. The block Pseudo-Elimination method was used to generate the preconditioned matrix. Numerical results show much better performance of this preconditioned techniques and the GMRES method than that of preconditioned ORTHMIN method, which is now in use in oil reservoir simulation. Finally, some limitations of Krylov subspace methods and some potential improvements to this type of methods are further presented.展开更多
The exploitation of an oil field is a complex and multidisciplinary task, which demands a lot of prior knowledge, time, and money. A good reservoir characterization is deemed essential in the accomplishment of Enhance...The exploitation of an oil field is a complex and multidisciplinary task, which demands a lot of prior knowledge, time, and money. A good reservoir characterization is deemed essential in the accomplishment of Enhanced Oil Recovery (EOR) processes in order to estimate accurately the properties of the porous medium affecting the flow properties. Several techniques at a field scale are currently being used to determine these properties, which are time and money consuming. But these alone do not guarantee the success of the project. Reservoir simulation and numerical techniques were then included in the pre-development and follow-up studies as an effective tool to determine the productivity and future behavior of the oil field. As the computational power increased, more advanced and detailed models were developed, including different chemical and physical phenomena. But alongside this process, there was an active research in the area of reservoir simulation, improving the accuracy and efficiency of the numerical schemes used for the flow, transport, and energy equations. The aim of this review is to address the topics described. Firstly, the origin of an oil recovery process, the economic factors and field tests involved are introduced. Secondly, the oil and porous medium origin and characterization as well as an introduction to the fundamental concepts and equations are associated to reservoir simulation. Finally, a brief description and analysis of the techniques are used in reservoir simulation employing finite difference methods, their downsides and possible ways to overcome these problems.展开更多
In the current study,an artificial neural network(ANN)and a numerical reservoir simulation(NRS)technique are used to analyse reservoir performance under waterflooding in the ZH86 block of the Zhaozhouqiao oilfield,Chi...In the current study,an artificial neural network(ANN)and a numerical reservoir simulation(NRS)technique are used to analyse reservoir performance under waterflooding in the ZH86 block of the Zhaozhouqiao oilfield,China.Using five input datasets extracted from the history-matched NRS model,an NRS-ANN hybrid is trained using a trial-and-error approach.NRS-ANN hybrid model#46(which has 5,10,10,6,6,and 1 neurons in the input layer,four hidden layers,and output layer,respectively)is found to produce the minimal root mean square error on the test dataset.On the validation data,the prediction performance of the selected NRS-ANN hybrid model achieves a minimal root mean square error of 0.0274 m^(3)/day and maximal coefficient of determination and coefficient of correlation values of about 0.9999.The correlation between the block liquid production rate(BLPR,m^(3)/day),block water production rate(BWPR,m^(3)/day),block water cut(BWCT,%),block water injection rate(BWIR,m^(3)/day),and block reservoir pressure(BRP,bar)as input variables and the simulated oil production rate(SOPRH)as the output variable is investigated.There is a positive correlation between SOPRH and BLPR,BWIR,and BWCT,and a negative correlation between SOPRH and BRP and BWPR.Segment B of ZH86 block experiences a 3.8%increase in BLPR,while segments A and C show declines of 1.3%and 1.6%,respectively.These variations in the liquid production rate correspond to changes in SOPRH of 4.3%,1.9%,and 9.7%for segments A,B,and C,respectively.The prediction performance of the NRS-ANN hybrid model is compared with that of a simple NRS model.The accuracy of the NRS-ANN hybrid model in predicting oil production is found to be 1125 times that of the NRS model.Based on these results,it is concluded that the proposed NRS-ANN hybrid provides an accurate and useful tool for analysing reservoir performance under the waterflooding oil recovery technique.展开更多
Oilfields worldwide are increasingly grappling with challenges such as early water breakthrough and high water production,yet direct,targeted solutions remain elusive.In recent years,chemical flooding techniques desig...Oilfields worldwide are increasingly grappling with challenges such as early water breakthrough and high water production,yet direct,targeted solutions remain elusive.In recent years,chemical flooding techniques designed for tertiary oil recovery have garnered significant attention,with microgel flooding emerging as a particularly prominent area of research.Despite its promise,the complex mechanisms underlying microgel flooding have been rarely investigated numerically.This study aims to address these gaps by characterizing the distribution of microgel concentration and viscosity within different pore structures.To enhance the accuracy of these characterizations,the viscosity of microgels is adjusted to account for the shear effects induced by flow rate and the swelling effects driven by salinity variations.The absolute permeability of the rock and the relative permeability of both oil and microgel are also analyzed to elucidate the mechanisms of microgel flooding.Additionally,a connectivity model is employed to achieve a quantitative representation of fluid flow capacity.The proposed model is validated through conceptual examples and applied to real oilfield blocks,demonstrating its accuracy and practical applicability.展开更多
We study a class of preconditioners to solve large-scale linear systems arising from fully implicit reservoir simulation. These methods are discussed in the framework of the auxiliary space preconditioning method for ...We study a class of preconditioners to solve large-scale linear systems arising from fully implicit reservoir simulation. These methods are discussed in the framework of the auxiliary space preconditioning method for generality. Unlike in the case of classical algebraic preconditioning methods, we take several analytical and physical considerations into account. In addition, we choose appropriate auxiliary problems to design the robust solvers herein. More importantly, our methods are user-friendly and general enough to be easily ported to existing petroleum reservoir simulators. We test the efficiency and robustness of the proposed method by applying them to a couple of benchmark problems and real-world reservoir problems. The numerical results show that our methods are both efficient and robust for large reservoir models.展开更多
In this paper,we focus on graphical processing unit(GPU)and discuss how its architecture affects the choice of algorithm and implementation of fully-implicit petroleum reservoir simulation.In order to obtain satisfact...In this paper,we focus on graphical processing unit(GPU)and discuss how its architecture affects the choice of algorithm and implementation of fully-implicit petroleum reservoir simulation.In order to obtain satisfactory performance on new many-core architectures such as GPUs,the simulator developers must know a great deal on the specific hardware and spend a lot of time on fine tuning the code.Porting a large petroleum reservoir simulator to emerging hardware architectures is expensive and risky.We analyze major components of an in-house reservoir simulator and investigate how to port them to GPUs in a cost-effective way.Preliminary numerical experiments show that our GPU-based simulator is robust and effective.More importantly,these numerical results clearly identify the main bottlenecks to obtain ideal speedup on GPUs and possibly other many-core architectures.展开更多
The conventional approach for an EOR process is to compare the reservoir properties with those of successful worldwide projects.However,some proper cases may be neglected due to the lack of reliable data.A combination...The conventional approach for an EOR process is to compare the reservoir properties with those of successful worldwide projects.However,some proper cases may be neglected due to the lack of reliable data.A combination of experimental design and reservoir simulation is an alternative approach.In this work,the fractional factorial design suggests some numerical experiments which their results are analyzed by statistical inference.After determination of the main effects and interactions,the most important parameters of polymer flooding are studied by ANOVA method and Pareto and Tornado charts.Analysis of main effects shows that the oil viscosity,connate water saturation and the horizontal permeability are the 3 deciding factors in oil production.The proposed methodology can help to select the good candidate reservoirs for polymer flooding.展开更多
This present paper has proved the theorem of the Point Optimal Variable Successive Over Relaxation (OVSOR) method of the three-dimensional unsteady flow in the reservoir, and has put forward a formu- la for calculatin...This present paper has proved the theorem of the Point Optimal Variable Successive Over Relaxation (OVSOR) method of the three-dimensional unsteady flow in the reservoir, and has put forward a formu- la for calculating optimal parameters for OVSOR which vary with space points and time points. Using this method, internal memory of computer is the smallest, calculating work is the smallest, and calculating funds are the smallest. It is very easy to operate on microcomputers for three-dimensional res- ervoir simulation. The method is stable and convergent even if the time steps are taken to be large (for example, one year). The same applies for space steps. It is applicable both for homogeneous, isotropic porous mediums and for heterogeneous, anisotropic porous medium. On IBM microcomputers with internal memory of 512 thousand bytes, 8000 grid points may be cal- culated for three-dimensional simulation. It takes only two minutes to get convergence for one time step. It may be extended to three-dimensional heat conduction equation and three-dimensional simulation of the ground water flow. It looks much more advantageous for two-dimensional simulation.展开更多
The industry-standard constrained pressure residual(CPR)algorithm is often able to effectively improve the robustness behavior and the convergence speed of linear iterations for isothermal reservoir simulation.In this...The industry-standard constrained pressure residual(CPR)algorithm is often able to effectively improve the robustness behavior and the convergence speed of linear iterations for isothermal reservoir simulation.In this paper,we present and study an improved extension of CPR to the constrained pressure-temperature residual(CPTR)version for non-isothermal reservoir problems in heterogeneous porous media.In the proposed preconditioner,the corresponding approximations for the inverse of matrices are computed under a domain decomposition framework by using the restricted additive Schwarz(RAS)algorithm,to equally deal with the coupled thermalpressure-saturation reservoir system and highly exploit the parallelism of supercomputer platforms.Moreover,we introduce and develop a family of multilevel CPTR preconditioners with suitable coarse grid corrections,to further improve the applicability of this two-stage preconditioner for large-scale computation.Numerical results for strong heterogeneous flow problems show that the new approach can dramatically improve the convergence of linear iterations,and demonstrate the superiority of CPTR over the commonly used RAS preconditioners.The parallel scalability of the non-isothermal reservoir simulator is also studied versus a supercomputer with tens of thousands of processors.展开更多
Seismic quantitative reservoir simulations and characterizations have played a vital role in exploring stratigraphic traps,such as lateaggradational lowstands prograding wedge systems(LPWS)within lowstands systems tra...Seismic quantitative reservoir simulations and characterizations have played a vital role in exploring stratigraphic traps,such as lateaggradational lowstands prograding wedge systems(LPWS)within lowstands systems tracts(LST).However,seismic data acquisition operations are always dominated by exceptional seismic coherent noise events,e.g.,multiples,which reduce the signal strengths of the sourcegenerated incident seismic waves within vertically and laterally heterogeneous earth systems.Hence,these noise events create hurdles in predicting paleo-depositional impedance(PDI),paleo-thickness(PTS),paleo-dense fractured networks,erosional and depositional zones,faultcontrolled migrations,and types of seismic reflection configurations(SRFC),which are key elements in developing stratigraphic pinch-out traps.This research utilizes the state-of-the-art technologies of spectral wavelet-based instantaneous time-frequency analysis and seismic waveform frequency-controlled porosity-constrained static reservoir simulation(FDPVS)tools to quantify the LPWS inside the Onshore Basin,Pakistan.The use of conventional amplitude-based seismic attributes,such as the average energy,remained a better tool for deciphering the overall geological architecture of the LPWS.Conventional FDPVS realizations resolved a PDI of−1.391 gm./c.c.^(*)m/s to−0.97 gm./c.c.^(*)m/s for LPWS with PTS of 12 and 20 m,respectively.A 0.9 km lateral extent of paleo-dense fractured networks(PDFN)with a strong linear regression R^(2)=0.93 was also resolved.Average energy attribute-based instantaneous frequency FDPVS realizations enabled the imaging of parallel-toprograding SRFC with resolved magnitudes of−0.259 gm./c.c.^(*)m/s for PDI,20 m for PTS,and 0.73 km for PDFN with linear regression transforms at R^(2)=0.92,which indicates the deposition of onlap fill facies inside the LPWS during extensive sea-level fall.These realizations have also resolved frequency-controlled fault migrations on 27-Hz spectral waveform-based amplitude plots with 2.174 gm./c.c.^(*)m/s PDI for conduit fault systems and 27-Hz with 0.585 gm./c.c.^(*)m/s PDI for sealing fault systems.All these structural configurations are completely sealed up by transgressive seals of transgressive systems tracts and,hence,developed into pure stratigraphic-based oil and gas plays.This research has strong implications for side-tracking drilling locations and provides an analogue for basins with similar geology and stratigraphy worldwide.展开更多
A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite...A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite element method was used to determine the numerical solution and the accuracy of the model was verified. On this basis, the model was used to simulate productivity of multistage fractured horizontal wells in tight oil reservoirs. The results show that during the production of tight oil wells, the reservoir region close to artificial fractures deteriorated in physical properties significantly, e.g. the aperture and conductivity of artificial fractures dropped by 52.12% and 89.02% respectively. The simulations of 3000-day production of a horizontal well in tight oil reservoir showed that the predicted productivity by the uncoupled model had an error of 38.30% from that by the fully-coupled model. Apparently, ignoring the influence of fluid-solid interaction effect led to serious deviations of the productivity prediction results. The productivity of horizontal well in tight oil reservoir was most sensitive to the start-up pressure gradient, and second most sensitive to the opening of artificial fractures. Enhancing the initial conductivity of artificial fractures was helpful to improve the productivity of tight oil wells. The influence of conductivity, spacing, number and length of artificial fractures should be considered comprehensively in fracturing design. Increasing the number of artificial fractures unilaterally could not achieve the expected increase in production.展开更多
Unexpected noise in reservoir stochastic simulation realization may be too high to make the realization useful, especially when there is a lack of hard data. Through discussing the uncertainties, we present two ways t...Unexpected noise in reservoir stochastic simulation realization may be too high to make the realization useful, especially when there is a lack of hard data. Through discussing the uncertainties, we present two ways to control the uncertainty ratio that is brought by the algorithm of stochastic simulation. By reasonably reducing the random value of the stochastic simulation result, the unexpected values introduced by the residual that associates with random series can be controlled. Another way when the data disperse unevenly is to control the stochastic simulation order by grouping the points that need to be simulated to make those points which can be simulated by more neighborhood hard data calculated first. Both methods do not go against the core stochastic simulation algorithm.展开更多
Numerical simulation of groundwater in karst areas has long been restricted by the difficulty of generalizing the hydrogeological conditions of reservoirs and of determining the relevant parameters due to the anisotro...Numerical simulation of groundwater in karst areas has long been restricted by the difficulty of generalizing the hydrogeological conditions of reservoirs and of determining the relevant parameters due to the anisotropy and discontinuity of the karst water-bearing media in these areas. In this study, we used the Guang'an Longtan Coal mine in Sichuan as an example, and generalized the complex hydrogeological conditions in the reservoir area. A finite element numerical flow model was used to simulate current and future scenarios of roadway gushing at the bottom of the coal mine at pile number 1 + 700 m. The results show that the roadway section corresponding to valleys has a gushing quantity of 4323.8–4551.25 m^3/d before impoundment. Modeled water inflow after impoundment increased to 1.6 times the water inflow before impoundment, which threatens the impoundment as well as the roadway's normal operation. Therefore, roadway processing measures are needed to guarantee the safety of the impoundment and of the mining operation.展开更多
It is known that the pore media characteristics of glutenite reservoirs are different from those of conventional sandstone reservoirs.Low reservoir permeability and naturally developed microfractures make water inject...It is known that the pore media characteristics of glutenite reservoirs are different from those of conventional sandstone reservoirs.Low reservoir permeability and naturally developed microfractures make water injection in this kind of reservoir very difficult.In this study,new exploitation methods are explored.Using a real glutenite reservoir as a basis,a three-dimensional fine geological model is elaborated.Then,combining the model with reservoir performance information,and through a historical fitting analysis,the saturation abundance distribution of remaining oil in the reservoir is determined.It is shown that,using this information,predictions can be made about whether the considered reservoir is suitable for horizontal well fracturing or not.The direction,well length,well spacing and productivity of horizontal well are also obtained.展开更多
Due to the difficulties associated with preprocessing activities and poor grid convergence when simulating shale reservoirs in the context of traditional grid methods,in this study an innovative two-phase oil-water se...Due to the difficulties associated with preprocessing activities and poor grid convergence when simulating shale reservoirs in the context of traditional grid methods,in this study an innovative two-phase oil-water seepage model is elaborated.The modes is based on the radial basis meshless approach and is used to determine the pressure and water saturation in a sample reservoir.Two-dimensional examples demonstrate that,when compared to the finite difference method,the radial basis function method produces less errors and is more accurate in predicting daily oil production.The radial basis function and finite difference methods provide errors of 5.78 percent and 7.5 percent,respectively,when estimating the daily oil production data for a sample well.A sensitivity analysis of the key parameters that affect the radial basis function’s computation outcomes is also presented.展开更多
文摘In this paper, the iterative coupling approach is proposed for applications to solving multiphase flow equation systems in reservoir simulation, as it provides a more flexible time-stepping strategy than existing approaches. The iterative method decouples the whole equation systems into pressure and saturation/concentration equations, and then solves them in sequence, implicitly and semi-implicitly. At each time step, a series of iterations are computed, which involve solving linearized equations using specific tolerances that are iteration dependent. Following convergence of subproblems, material balance is checked. Convergence of time steps is based on material balance errors. Key components of the iterative method include phase scaling for deriving a pressure equation and use of several advanced numerical techniques. The iterative model is implemented for parallel computing platforms and shows high parallel efficiency and scalability.
文摘Petroleum science has made remarkable progress in organic geochemistry and in the research into the theories of petroleum origin, its transport and accumulation. In estimating the oil-gas resources of a basin, the knowledge of its evolutionary history and especially the numerical computation of fluid flow and the history of its changes under heat is vital. The mathematical model call be described as a coupled system of nonlinear partial differentical equations with initial-boundary value problems. This thesis, from actual conditions such as the effect of fluid compressibility and the characteristic of large-scal science-engineering computalion. puts forward a kind of characteristic finite difference alternating-direction scheme. Optimal order estimates in L-2 norm are derived for the error in the approximate solutions.
文摘This paper describes a way of solving the reservoir simulation pressure equation using mulligrid technique. The subroutine MG of four-grid method is presented. The result for 2-D two-phase problem is exactly the same as that of the SOR method and the CPU time is much less than that of the latter one.
文摘This paper investigates the deposition of asphaltenes in the porous medium of the studied field in Russia and predicts production profiles based on uncertainty evaluation. This problem can be solved by dynamic modeling, during which production profiles are estimated in two scenarios: with and without the activation of the asphaltene option. Calculations are carried out for two development scenarios: field operation under natural depletion and water injection into the aquifer as a reservoir pressure maintenance system. A full-scale compositional reservoir simulation model of the Russian oilfield was created. Within a dynamic simulation, the asphaltene option was activated and the asphaltene behavior in oil and porous medium was tuned according to our own special laboratory experiments. The model was also matched to production historical data, and a pattern model was prepared using the full-scale simulation model. Technological and the asphaltene option parameters were used in sensitivity and an uncertainty evaluation. Furthermore, probable production profiles within a forecast period were estimated. The sensitivity analysis of the pattern model to input parameters of the asphaltene option allowed determining the following heavy-hitters on the objective function: the molar weight of dissolved asphaltenes as a function of pressure, the asphaltene dissociation rate, the asphaltene adsorption coefficient and the critical velocity of oil movement in the reservoir. Under the natural depletion scenario, our simulations show a significant decrease in reservoir pressure and the formation of drawdown cones leading to asphaltene deposition in the bottom-hole area of production wells, decreasing their productivity. Water injection generally allows us to significantly reduce the volume of asphaltene phase transitions and has a positive effect on cumulative oil production. Injecting water into aquifer can keep the formation pressure long above the pressure for asphaltene precipitation, preventing the asphaltene deposition resulted from interaction of oil and water, so this way has higher oil production.
文摘Simulation study was applied in the development planning of East Unity oilfield, Sudan. A grid consisting of 2 000 cells was constructed. A major challenge of the study wasto evolve a full field development and future reservoir management strategy that would ensuremaximum recovery of oil based on well Un51. Simulation shows that Un51 as injection well inAradiebaC would yield better oil recovery than to be production well.
文摘Krylov subspace projection methods are known to be highly efficient for solving large linear systems. Many different versions arise from different choices to the left and right subspaces. These methods were classified into two groups in terms of the different forms of matrix H-m, the main properties in applications and the new versions of these two types of methods were briefly reviewed, then one of the most efficient versions, GMRES method was applied to oil reservoir simulation. The block Pseudo-Elimination method was used to generate the preconditioned matrix. Numerical results show much better performance of this preconditioned techniques and the GMRES method than that of preconditioned ORTHMIN method, which is now in use in oil reservoir simulation. Finally, some limitations of Krylov subspace methods and some potential improvements to this type of methods are further presented.
文摘The exploitation of an oil field is a complex and multidisciplinary task, which demands a lot of prior knowledge, time, and money. A good reservoir characterization is deemed essential in the accomplishment of Enhanced Oil Recovery (EOR) processes in order to estimate accurately the properties of the porous medium affecting the flow properties. Several techniques at a field scale are currently being used to determine these properties, which are time and money consuming. But these alone do not guarantee the success of the project. Reservoir simulation and numerical techniques were then included in the pre-development and follow-up studies as an effective tool to determine the productivity and future behavior of the oil field. As the computational power increased, more advanced and detailed models were developed, including different chemical and physical phenomena. But alongside this process, there was an active research in the area of reservoir simulation, improving the accuracy and efficiency of the numerical schemes used for the flow, transport, and energy equations. The aim of this review is to address the topics described. Firstly, the origin of an oil recovery process, the economic factors and field tests involved are introduced. Secondly, the oil and porous medium origin and characterization as well as an introduction to the fundamental concepts and equations are associated to reservoir simulation. Finally, a brief description and analysis of the techniques are used in reservoir simulation employing finite difference methods, their downsides and possible ways to overcome these problems.
基金Funding for this research was provided by the Department of Petroleum Engineering,China University of Geosciences,under the sponsorship of the China Scholarship Council (CSC)the National Natural Science Foundation of China through grant nos.41972326 and 51774258.
文摘In the current study,an artificial neural network(ANN)and a numerical reservoir simulation(NRS)technique are used to analyse reservoir performance under waterflooding in the ZH86 block of the Zhaozhouqiao oilfield,China.Using five input datasets extracted from the history-matched NRS model,an NRS-ANN hybrid is trained using a trial-and-error approach.NRS-ANN hybrid model#46(which has 5,10,10,6,6,and 1 neurons in the input layer,four hidden layers,and output layer,respectively)is found to produce the minimal root mean square error on the test dataset.On the validation data,the prediction performance of the selected NRS-ANN hybrid model achieves a minimal root mean square error of 0.0274 m^(3)/day and maximal coefficient of determination and coefficient of correlation values of about 0.9999.The correlation between the block liquid production rate(BLPR,m^(3)/day),block water production rate(BWPR,m^(3)/day),block water cut(BWCT,%),block water injection rate(BWIR,m^(3)/day),and block reservoir pressure(BRP,bar)as input variables and the simulated oil production rate(SOPRH)as the output variable is investigated.There is a positive correlation between SOPRH and BLPR,BWIR,and BWCT,and a negative correlation between SOPRH and BRP and BWPR.Segment B of ZH86 block experiences a 3.8%increase in BLPR,while segments A and C show declines of 1.3%and 1.6%,respectively.These variations in the liquid production rate correspond to changes in SOPRH of 4.3%,1.9%,and 9.7%for segments A,B,and C,respectively.The prediction performance of the NRS-ANN hybrid model is compared with that of a simple NRS model.The accuracy of the NRS-ANN hybrid model in predicting oil production is found to be 1125 times that of the NRS model.Based on these results,it is concluded that the proposed NRS-ANN hybrid provides an accurate and useful tool for analysing reservoir performance under the waterflooding oil recovery technique.
基金supported by the National Natural Science Foundation project“Micro-Scale Effect of Oil-Gas Flow and the Mechanism of Enhancing Shale Oil Recovery by Natural Gas Injection”(No.52074317)。
文摘Oilfields worldwide are increasingly grappling with challenges such as early water breakthrough and high water production,yet direct,targeted solutions remain elusive.In recent years,chemical flooding techniques designed for tertiary oil recovery have garnered significant attention,with microgel flooding emerging as a particularly prominent area of research.Despite its promise,the complex mechanisms underlying microgel flooding have been rarely investigated numerically.This study aims to address these gaps by characterizing the distribution of microgel concentration and viscosity within different pore structures.To enhance the accuracy of these characterizations,the viscosity of microgels is adjusted to account for the shear effects induced by flow rate and the swelling effects driven by salinity variations.The absolute permeability of the rock and the relative permeability of both oil and microgel are also analyzed to elucidate the mechanisms of microgel flooding.Additionally,a connectivity model is employed to achieve a quantitative representation of fluid flow capacity.The proposed model is validated through conceptual examples and applied to real oilfield blocks,demonstrating its accuracy and practical applicability.
基金supported by Petro-China Joint Research Funding(Grant No.12HT1050002654)National Science Foundation of USA(Grant No.DMS-1217142)+1 种基金the Dean’s Startup FundAcademy of Mathematics and System Sciences and the State High Tech Development Plan of China(863 Program)(GrantNo.2012AA01A309)
文摘We study a class of preconditioners to solve large-scale linear systems arising from fully implicit reservoir simulation. These methods are discussed in the framework of the auxiliary space preconditioning method for generality. Unlike in the case of classical algebraic preconditioning methods, we take several analytical and physical considerations into account. In addition, we choose appropriate auxiliary problems to design the robust solvers herein. More importantly, our methods are user-friendly and general enough to be easily ported to existing petroleum reservoir simulators. We test the efficiency and robustness of the proposed method by applying them to a couple of benchmark problems and real-world reservoir problems. The numerical results show that our methods are both efficient and robust for large reservoir models.
基金support from LSEC.The authors would like to thank RIPED,PetroChina,for providing data for the numerical tests and support through PetroChina New-generation Reservoir Simulation Software(No.2011A-1010)the Program of Research on Continental Sedimentary Oil Reservoir Simulation(No.z121100004912001)founded by Beijing Municipal Science&Technology Commission and PetroChina Joint Research Funding No.12HT1050002654.
文摘In this paper,we focus on graphical processing unit(GPU)and discuss how its architecture affects the choice of algorithm and implementation of fully-implicit petroleum reservoir simulation.In order to obtain satisfactory performance on new many-core architectures such as GPUs,the simulator developers must know a great deal on the specific hardware and spend a lot of time on fine tuning the code.Porting a large petroleum reservoir simulator to emerging hardware architectures is expensive and risky.We analyze major components of an in-house reservoir simulator and investigate how to port them to GPUs in a cost-effective way.Preliminary numerical experiments show that our GPU-based simulator is robust and effective.More importantly,these numerical results clearly identify the main bottlenecks to obtain ideal speedup on GPUs and possibly other many-core architectures.
文摘The conventional approach for an EOR process is to compare the reservoir properties with those of successful worldwide projects.However,some proper cases may be neglected due to the lack of reliable data.A combination of experimental design and reservoir simulation is an alternative approach.In this work,the fractional factorial design suggests some numerical experiments which their results are analyzed by statistical inference.After determination of the main effects and interactions,the most important parameters of polymer flooding are studied by ANOVA method and Pareto and Tornado charts.Analysis of main effects shows that the oil viscosity,connate water saturation and the horizontal permeability are the 3 deciding factors in oil production.The proposed methodology can help to select the good candidate reservoirs for polymer flooding.
文摘This present paper has proved the theorem of the Point Optimal Variable Successive Over Relaxation (OVSOR) method of the three-dimensional unsteady flow in the reservoir, and has put forward a formu- la for calculating optimal parameters for OVSOR which vary with space points and time points. Using this method, internal memory of computer is the smallest, calculating work is the smallest, and calculating funds are the smallest. It is very easy to operate on microcomputers for three-dimensional res- ervoir simulation. The method is stable and convergent even if the time steps are taken to be large (for example, one year). The same applies for space steps. It is applicable both for homogeneous, isotropic porous mediums and for heterogeneous, anisotropic porous medium. On IBM microcomputers with internal memory of 512 thousand bytes, 8000 grid points may be cal- culated for three-dimensional simulation. It takes only two minutes to get convergence for one time step. It may be extended to three-dimensional heat conduction equation and three-dimensional simulation of the ground water flow. It looks much more advantageous for two-dimensional simulation.
基金supported by the National Natural Science Foundation of China(No.12131002 and No.11971006)Shenzhen Science and Technology Program(No.JCYJ20210324130801003)+2 种基金Guangdong Basic and Applied Basic Research Foundation(No.2022A1515010147)Changsha science and technology bureau(No.kh2301001)The fourth author also greatly thanks for the support from King Abdullah University of Science and Technology(KAUST)through the grants BAS/1/1351-01 and URF/1/4074-01.
文摘The industry-standard constrained pressure residual(CPR)algorithm is often able to effectively improve the robustness behavior and the convergence speed of linear iterations for isothermal reservoir simulation.In this paper,we present and study an improved extension of CPR to the constrained pressure-temperature residual(CPTR)version for non-isothermal reservoir problems in heterogeneous porous media.In the proposed preconditioner,the corresponding approximations for the inverse of matrices are computed under a domain decomposition framework by using the restricted additive Schwarz(RAS)algorithm,to equally deal with the coupled thermalpressure-saturation reservoir system and highly exploit the parallelism of supercomputer platforms.Moreover,we introduce and develop a family of multilevel CPTR preconditioners with suitable coarse grid corrections,to further improve the applicability of this two-stage preconditioner for large-scale computation.Numerical results for strong heterogeneous flow problems show that the new approach can dramatically improve the convergence of linear iterations,and demonstrate the superiority of CPTR over the commonly used RAS preconditioners.The parallel scalability of the non-isothermal reservoir simulator is also studied versus a supercomputer with tens of thousands of processors.
文摘Seismic quantitative reservoir simulations and characterizations have played a vital role in exploring stratigraphic traps,such as lateaggradational lowstands prograding wedge systems(LPWS)within lowstands systems tracts(LST).However,seismic data acquisition operations are always dominated by exceptional seismic coherent noise events,e.g.,multiples,which reduce the signal strengths of the sourcegenerated incident seismic waves within vertically and laterally heterogeneous earth systems.Hence,these noise events create hurdles in predicting paleo-depositional impedance(PDI),paleo-thickness(PTS),paleo-dense fractured networks,erosional and depositional zones,faultcontrolled migrations,and types of seismic reflection configurations(SRFC),which are key elements in developing stratigraphic pinch-out traps.This research utilizes the state-of-the-art technologies of spectral wavelet-based instantaneous time-frequency analysis and seismic waveform frequency-controlled porosity-constrained static reservoir simulation(FDPVS)tools to quantify the LPWS inside the Onshore Basin,Pakistan.The use of conventional amplitude-based seismic attributes,such as the average energy,remained a better tool for deciphering the overall geological architecture of the LPWS.Conventional FDPVS realizations resolved a PDI of−1.391 gm./c.c.^(*)m/s to−0.97 gm./c.c.^(*)m/s for LPWS with PTS of 12 and 20 m,respectively.A 0.9 km lateral extent of paleo-dense fractured networks(PDFN)with a strong linear regression R^(2)=0.93 was also resolved.Average energy attribute-based instantaneous frequency FDPVS realizations enabled the imaging of parallel-toprograding SRFC with resolved magnitudes of−0.259 gm./c.c.^(*)m/s for PDI,20 m for PTS,and 0.73 km for PDFN with linear regression transforms at R^(2)=0.92,which indicates the deposition of onlap fill facies inside the LPWS during extensive sea-level fall.These realizations have also resolved frequency-controlled fault migrations on 27-Hz spectral waveform-based amplitude plots with 2.174 gm./c.c.^(*)m/s PDI for conduit fault systems and 27-Hz with 0.585 gm./c.c.^(*)m/s PDI for sealing fault systems.All these structural configurations are completely sealed up by transgressive seals of transgressive systems tracts and,hence,developed into pure stratigraphic-based oil and gas plays.This research has strong implications for side-tracking drilling locations and provides an analogue for basins with similar geology and stratigraphy worldwide.
基金Supported by the National Science and Technology Major Project (2017ZX05013-005)。
文摘A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite element method was used to determine the numerical solution and the accuracy of the model was verified. On this basis, the model was used to simulate productivity of multistage fractured horizontal wells in tight oil reservoirs. The results show that during the production of tight oil wells, the reservoir region close to artificial fractures deteriorated in physical properties significantly, e.g. the aperture and conductivity of artificial fractures dropped by 52.12% and 89.02% respectively. The simulations of 3000-day production of a horizontal well in tight oil reservoir showed that the predicted productivity by the uncoupled model had an error of 38.30% from that by the fully-coupled model. Apparently, ignoring the influence of fluid-solid interaction effect led to serious deviations of the productivity prediction results. The productivity of horizontal well in tight oil reservoir was most sensitive to the start-up pressure gradient, and second most sensitive to the opening of artificial fractures. Enhancing the initial conductivity of artificial fractures was helpful to improve the productivity of tight oil wells. The influence of conductivity, spacing, number and length of artificial fractures should be considered comprehensively in fracturing design. Increasing the number of artificial fractures unilaterally could not achieve the expected increase in production.
文摘Unexpected noise in reservoir stochastic simulation realization may be too high to make the realization useful, especially when there is a lack of hard data. Through discussing the uncertainties, we present two ways to control the uncertainty ratio that is brought by the algorithm of stochastic simulation. By reasonably reducing the random value of the stochastic simulation result, the unexpected values introduced by the residual that associates with random series can be controlled. Another way when the data disperse unevenly is to control the stochastic simulation order by grouping the points that need to be simulated to make those points which can be simulated by more neighborhood hard data calculated first. Both methods do not go against the core stochastic simulation algorithm.
基金supported by the National Natural Science Foundation of China (41272377)
文摘Numerical simulation of groundwater in karst areas has long been restricted by the difficulty of generalizing the hydrogeological conditions of reservoirs and of determining the relevant parameters due to the anisotropy and discontinuity of the karst water-bearing media in these areas. In this study, we used the Guang'an Longtan Coal mine in Sichuan as an example, and generalized the complex hydrogeological conditions in the reservoir area. A finite element numerical flow model was used to simulate current and future scenarios of roadway gushing at the bottom of the coal mine at pile number 1 + 700 m. The results show that the roadway section corresponding to valleys has a gushing quantity of 4323.8–4551.25 m^3/d before impoundment. Modeled water inflow after impoundment increased to 1.6 times the water inflow before impoundment, which threatens the impoundment as well as the roadway's normal operation. Therefore, roadway processing measures are needed to guarantee the safety of the impoundment and of the mining operation.
文摘It is known that the pore media characteristics of glutenite reservoirs are different from those of conventional sandstone reservoirs.Low reservoir permeability and naturally developed microfractures make water injection in this kind of reservoir very difficult.In this study,new exploitation methods are explored.Using a real glutenite reservoir as a basis,a three-dimensional fine geological model is elaborated.Then,combining the model with reservoir performance information,and through a historical fitting analysis,the saturation abundance distribution of remaining oil in the reservoir is determined.It is shown that,using this information,predictions can be made about whether the considered reservoir is suitable for horizontal well fracturing or not.The direction,well length,well spacing and productivity of horizontal well are also obtained.
基金supported by The China Postdoctoral Science Foundation(2021M702304)Natural Science Foundation of Shandong Province(ZR2021QE260).
文摘Due to the difficulties associated with preprocessing activities and poor grid convergence when simulating shale reservoirs in the context of traditional grid methods,in this study an innovative two-phase oil-water seepage model is elaborated.The modes is based on the radial basis meshless approach and is used to determine the pressure and water saturation in a sample reservoir.Two-dimensional examples demonstrate that,when compared to the finite difference method,the radial basis function method produces less errors and is more accurate in predicting daily oil production.The radial basis function and finite difference methods provide errors of 5.78 percent and 7.5 percent,respectively,when estimating the daily oil production data for a sample well.A sensitivity analysis of the key parameters that affect the radial basis function’s computation outcomes is also presented.
