The pressure buildup test is performed by shutting in a producing well at time $t=t_p$, after which a smooth rise of the well head pressure can be observed. For a geothermal reservoir, the buildup test result is interpreted using a Horner plot in order to evaluate the reservoir permeability or transmissivity. In this benchmark, observation data from a geothermal well is employed to parameterize the numerical model. In the model, a time dependent nodal source term was set up to represent the shutin operation. The simulated pressure profile is then verified against the analytical solution.
This benchmark represents a scenario in which the well had been producing geothermal brine for $118\ \mathrm{h}$ at a rate of $78\ t/\mathrm{h}$ and was then shut in for a buildup test. At the given reservoir temperature (260$^\circ$C) and pressure ($47\ \mathrm{bar}$), the density and viscosity of the water and steam mixture were calculated as $\rho=78.68\ \mathrm{kg/m^3}$ and $\mu=1.295\times10^{4}\ \mathrm{Pa\ s}$. The compressibility of the mixture was estimated as $nc_t=0.0805\ \mathrm{bar^{1}}$ ($n$ refers to the porosity of the reservoir), which yields a specific storage coefficient of $S=\rho gnc_t=6.21\times 10^{4}\ \mathrm{m^{1}}$. The observed pressure readings during the buildup test are cited from Chapter 6 of the book Geothermal Power Generation [1], and the data is archived in the Appendix.
The permeability of the reservoir can be estimated by means of a Horner plot, in which the pressure $p$ is plotted against $(t_p+\Delta t)/\Delta t$, i.e. Horner time on a semilogarithmic scale (cf. Figure 1). In the Horner plot, the data points form a straight line in the latetime period of the test. Note that the time increases in the opposite direction of the Xaxis. Therefore, the linear section appears at the left side of the diagram.
\centering
Figure 1: Horner plot ($p$ vs $(t_p+\Delta t)/\Delta t$) for buildup test showing the inferred Horner straight line
The slope $m$ of the Horner straight line is expressed as: $$m=0.1832\frac{Q\mu}{\kappa b}$$ in which $Q\ \mathrm{[L^3/T]}\ (Q>0)$ is the production rate of the well before shutin, $\kappa\ \mathrm{[L^2]}$ is the permeability and $b\ \mathrm{[L]}$ is the aquifer thickness. From the Horner plot, we can infer a Horner straight line with a slope of $m=0.79$. Therefore the transmissivity of the aquifer can be calculated as $$\begin{aligned} \kappa b&=&0.1832\frac{Q\mu}{m}=0.1832\frac{((78000/3600)/78.68)\times1.295\times10^{4}}{7.09\times10^5}\&=&9.2\times10^{12}\ \mathrm{m^3}\end{aligned}$$ In addition, the straight line in the Horner plot can be extracted to a Horner time of 1, which corresponds to the infinite shutin time $(\Delta t)$. This leads to an extrapolated pressure $p_0$ of $67.5~\mathrm{kPa}$, which is the undisturbed reservoir pressure .
The benchmark project is defined in the input file buildup_test.prj
. It defines the process to
be solved as “LiquidFlow” and the primary variable is hence pressure.
The initial condition is set to $p_0=67.5\ \mathrm{bar}$ and the
undisturbed boundary is achieved by a large domain size
$(r=1000\ \mathrm{m})$. The timedependent source term is applied in this
benchmark. From the beginning until $t=424800$ sec, the pumping rate was
maintained at a constant rate. Afterwards, the well is shutin and pressure
starts to build up. The geometries used to specify the model domain, boundary
conditions, and source term can be found in line_1000_axi.gml
file.
The mesh is specified in line_1000_axi.vtu
, which is stored in the
VTK format and can be directly visualized in ParaView.
The pressure buildup test is comparable to a pumping recovery test as the extraction rate is first kept constant at $Q$, and then becomes zero at $t=t_p$. This benchmark then adopts the same assumptions as in the Theis' problem. The analytical solution of the pressure difference $\Delta p$ with respect to the initial pressure $p_0$ is the sum of two Theis curves: one starting at $t=0$ and another starting at $t=t_p$ but with an opposite extraction rate, i.e. for $t\leq t_p$, $$\Delta p=\rho g \frac{Q}{4\pi T}W\left(\frac{r^2S}{4Tt}\right)$$ and for $t>t_p$, $$\Delta p=\rho g \frac{Q}{4\pi T}W\left(\frac{r^2S}{4Tt}\right)+\rho g \frac{Q}{4\pi T}W\left(\frac{r^2S}{4T(tt_p)}\right)$$
The pressure evolution is simulated throughout the domain and the result is compared with the analytical solution at $r=10.287\ \mathrm{m}$. In Figure 2, it can be observed that the numerical model fits with the analytical solution very well. The absolute and relative error between the analytical and numerical solution is depicted in Figure 3.
Figure 2: OGS 6 result compared with analytical solution
Figure 3: Absolute and relative error
[1] RN Horne. Characterization, evaluation, and interpretation of well data. In: R DiPippo, editor,Geothermal Power Generation, chapter 6, pages 141–163.Elsevier, 2016.
\centering
$\Delta t$ (h)  $\Delta p$ (bar)  $\Delta t$ (h)  $\Delta p$ (bar) 

0.0024  0.174  0.1708  3.65 
0.0073  0.695  0.2442  4.00 
0.0098  1.13  0.3667  4.26 
0.0122  1.30  0.6111  5.13 
0.0171  1.57  0.8556  6.35 
0.022  1.74  1.2194  7.48 
0.0244  1.91  1.5861  8.17 
0.0292  2.00  1.8361  8.43 
0.0367  2.09  2.4417  9.22 
0.0414  2.17  3.4167  10.2 
0.0489  2.43  3.8611  10.4 
0.0586  2.61  6.3056  11.8 
0.0733  2.78  8.3056  12.6 
0.0856  2.96  10.9722  13.2 
0.0975  3.13 
: Pressure measurements during well buildup[]{label=“table:1”}
This article was written by Boyan Meng and Haibing Shao. If you are missing something or you find an error please let us know.
Generated with Hugo 0.101.0
in CI job 282302

Last revision: November 22, 2022
Commit: [PL/THM] Implement freezing for temperature eq. 68ebbec
 Edit this page on