Integral-Balance Solution to Nonlinear Subdiffusion Equation

Author(s):  
Jordan Hristov
2015 ◽  
Vol 23 (17) ◽  
pp. 2795-2818 ◽  
Author(s):  
Jordan Hristov

An approximate integral-balance solution of the fractional subdiffusion equation by a double-integration technique has been conceived. The time-fractional linear subdiffusion equation with Dirichlet boundary condition (and zero initial condition) has been chosen as a test example. Approximations of time-fractional Riemann-Liouville and Caputo derivatives when the distribution is assumed as a parabolic profile with unspecified exponent have been developed. Problems pertinent to determination of the optimal exponent of the parabolic profile and approximations of the time-fractional derivative of by different approaches have been formulated. Solved and unresolved problems in determination of the optimal exponents have been demonstrated. Examples with predetermined quadratic and cubic assumed profiles are analyzed, too. Comparative numerical studies with exact solutions expressed by the Mainardi function in terms of a similarity variable have been performed.


2002 ◽  
Vol 5 (01) ◽  
pp. 49-59 ◽  
Author(s):  
J.L. Pletcher

Summary Experience with material-balance data sets from the field and from simulation has revealed some procedures that can be used to improve analysis of both oil and gas reservoirs:Failure to account for a weak waterdrive can result in significant material-balance errors.The assertion of previous authors that weak waterdrive exhibits a negative slope on the Cole (gas) and Campbell (oil) plots has been confirmed. A weak waterdrive is much more unambiguous on these plots than on commonly used plots, such as the p/z plot for gas.A modified version of the Cole plot is proposed to account for formation compressibility.The reservoir drive indices are a useful tool for determining the correctness of the material-balance solution because they must sum to unity. The drive indices should never be normalized to sum to unity because this obscures their usefulness and leads to a false sense of security.A modified version of the Roach plot (for gas) is proposed that improves interpretation in some waterdrive situations.Material balance has not been replaced by reservoir simulation; rather, it is complementary to simulation and can provide valuable insights to reservoir performance that cannot be obtained by simulation. Introduction Classical material balance is one of the fundamental tools of reservoir engineering. Many authors have addressed the difficult problem of solving the material balance in the presence of a waterdrive (Refs. 1 through 5 are just a few of the more significant ones). The emphasis in the literature has been on strong and moderate waterdrives. In this paper, examples of weak waterdrives are shown in which the effects on the material balance are significant. All aquifers studied here are of the "pot aquifer" type, which is time-independent. In gas reservoirs, the plot of p/z vs. cumulative gas production, Gp, is a widely accepted method for solving the gas material balance1 under depletion-drive conditions. Extrapolation of the plot to atmospheric pressure provides a reliable estimate of original gas in place (OGIP). If a waterdrive is present, the plot often appears to be linear, but the extrapolation will give an erroneously high value for OGIP. Many authors have addressed this problem (including those in Refs. 2 and 5 through 8), especially in cases of strong or moderate waterdrives. The p/z plot is actually more ambiguous in weak waterdrives than in strong or moderate ones. The Cole plot7,9 has proven to be a valuable diagnostic tool for distinguishing between depletion-drive gas reservoirs and those that are producing under a waterdrive. The analogous plot for oil reservoirs is the Campbell plot.10 The literature has emphasized strong and moderate waterdrives, the signature shapes of which are a positive slope and a hump-shaped curve, respectively, on these plots. Previous authors have recognized that weak waterdrives can produce negative slopes on these two diagnostic plots, but this author is not aware of any example plots in the literature. This paper shows examples, using simulation and actual field data, wherein a negative slope clearly reveals a weak waterdrive. These plots are much more diagnostic than the p/z plot. Once a weak waterdrive has been diagnosed, the appropriate steps can be taken in the material-balance equations to yield more accurate results. The Cole plot assumes that formation compressibility can be neglected, which is frequently the case with gas. However, in those reservoirs in which formation compressibility is significant, a modification to the Cole plot is presented that incorporates formation compressibility and gives more accurate results. The reservoir drive indices have been used to quantify the relative magnitude of the various energy sources active in a reservoir. It is shown here that the drive indices are also a useful diagnostic tool for determining the correctness of a material balance solution because they must sum to unity. If they do not sum to unity, a correct solution has not been obtained. In some commercial material-balance software, the drive indices are automatically normalized to sum to unity, which not only obscures their usefulness but also leads to the false impression of having achieved a correct solution. The Roach plot has been presented11 as a tool for solving the gas material balance when formation compressibility is unknown, with or without the presence of waterdrive. This paper shows that for waterdrives that fit the small pot aquifer model, incorporating cumulative water production into the x-axis plotting term improves the linearity of the Roach plot and gives more accurate values for OGIP. Finally, it is argued that even in those reservoirs for which a simulation study is performed, classical material-balance evaluation should be performed on a stand-alone basis. Simulation should not be viewed as a replacement for material balance because the latter can yield valuable insights that can be obscured during simulation. Performing a separate material balance study usually will improve overall reservoir understanding and enhance any subsequent simulation study. Material balance should be viewed as a complement to simulation, not as a competing approach. In this paper, formation compressibility, cf, is assumed to be constant and unchanging over the reservoir life under investigation. References are given for recommended methods to be used in those cases in which cf is variable.


2022 ◽  
pp. 138-148
Author(s):  
A. N. Pyshkin

Approaches to the methodology for calculating key performance indicators for development institutions as part of their contribution to national development goals achievement has been reviewed. It has been shown that, depending on the methodology used, the value of this indicator can vary by up to 5 units, with each approach being relevant for a particular task. Using the example of an indicator in one of the areas of the national project “Small and Medium-Sized Business” (training activities) as an example, an integrated approach has been proposed to balance solution of the objectives of optimal spending of the federal budget and increasing the public outreach of its results. 


2021 ◽  
Vol 57 (2) ◽  
pp. 297-309
Author(s):  
J. I. Castorena ◽  
A. C. Raga ◽  
A. Esquivel ◽  
A. Rodríguez-González ◽  
L. Hernández-Martínez ◽  
...  

We study the problem of a Herbig-Haro jet with a uniformly accelerating ejection velocity, travelling into a uniform environment. For the ejection density we consider two cases: a time-independent density, and a time-independent mass loss rate. For these two cases, we obtain analytic solutions for the motion of the jet head using a ram-pressure balance and a center of mass equation of motion. We also compute axisymmetric numerical simulations of the same flow, and compare the time-dependent positions of the leading working surface shocks with the predictions of the two analytic models. We find that if the jet is over-dense and over-pressured (with respect to the environment) during its evolution, a good agreement is obtained with the analytic models, with the flow initially following the center of mass analytic solution, and (for the constant ejection density case) at later times approaching the ram-pressure balance solution.


Sign in / Sign up

Export Citation Format

Share Document