Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Physics

formerly Central European Journal of Physics

Editor-in-Chief: Seidel, Sally

Managing Editor: Lesna-Szreter, Paulina

IMPACT FACTOR 2018: 1.005

CiteScore 2018: 1.01

SCImago Journal Rank (SJR) 2018: 0.237
Source Normalized Impact per Paper (SNIP) 2018: 0.541

ICV 2017: 162.45

Open Access
See all formats and pricing
More options …
Volume 15, Issue 1


Volume 13 (2015)

Research on carbonate reservoir interwell connectivity based on a modified diffusivity filter model

Wei Wang
  • Corresponding author
  • Research Centre of Multiphase Flow in Porous Media, China University of Petroleum (East China), Qingdao 266580, China
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Jun Yao
  • Research Centre of Multiphase Flow in Porous Media, China University of Petroleum (East China), Qingdao 266580, China
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Yang Li / Aimin Lv
  • Research Centre of Multiphase Flow in Porous Media, China University of Petroleum (East China), Qingdao 266580, China
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
Published Online: 2017-05-06 | DOI: https://doi.org/10.1515/phys-2017-0034


According to the solution of dual-porosity model, a diffusivity filter model of carbonate reservoir was established, which can effectively illustrate the injection signal attenuation and lag characteristic. The interwell dynamic connectivity inversion model combines a multivariate linear regression (MLR) analysis with a correction coefficient to eliminate the effect of fluctuating bottom-hole pressure (BHP). The modified MLR model was validated by synthetic field with fluctuating BHP. The method was applied to Tahe oilfield which showed that the inversion result was reliable. The interwell dynamic connectivity coefficients could reflect the real interwell connectivity of reservoir. The method is easy to use and proved to be effective in field applications.

Keywords: Carbonate reservoir; filter coefficient; interwell connectivity; bottom-hole pressure fluctuations

PACS: 47.10.-g; 47.56.+r

1 Introduction

Carbonate reservoir has various storage types that, are mainly cavities and fractures. It is difficult to describe the interwell connectivity due to its serious heterogeneity [13]. The evaluation of interwell connectivity is an important part of reservoir analysis, especially in carbonate reservoirs [4]. The interwell connectivity provides us a tool to divide fracture-vug units, identify location for new wells, and manage waterflood performance.

Both the injection and production rates are the most available data in any waterflooding reservoir. Numerous studies to investigate interwell connectivity using these data have been carried out. Heffer et al. [5] and Refunjol [6] used Spearman analysis to determine flow trends in a reservoir. Albertoni and Lake [7] used multivariate linear regression (MLR) analysis of the flow rates in a waterflooding reservoir to estimate interwell connectivity. They also introduced diffusivity filters to account for the delay where an injector signal reaches a producer. Yousef et al. [8] and Zhao et al. [9] introduced a capacitance model, which evaluated the connectivity by combining producer bottom-hole pressure (BHP) and flow rates.

There are many methods to analyze the interwell connectivity in a carbonate reservoir, such as pressure system analysis, interference test and similar interference test, fluid property contrast, water injection response, tracer agent, etc [10, 11]. Some researchers combined the dynamic and static data, including the palaeogeomorphology, formation structure, fracture, tracer and well interference, to estimate the interwell connectivity [12,13]. However, the required parameters in these methods are difficult to obtain and the corresponding methods are expensive. In addition, some methods need to shut-in wells, which will affect the normal production.

In this paper, a diffusivity filter model of carbonate reservoir is established based on the solution of an infinite reservoir with the double-porosity model proposed by Warren and Root [14]. In previous studies by Albertoni and Lake [7], MLR analysis was used to determine the interwell connectivity coefficients according to the injection and production data. We also propose a modified MLR model to account for the case of BHP fluctuating. The developed technique is applied to a synthetic field firstly and then applied to a real carbonate reservoir in Tahe oilfield.

2 Diffusivity filters of carbonate reservoir

2.1 Diffusivity filters development

According to the research results of reservoir studies, the fractured carbonate reservoir could be described by a double-porosity system: matrix system and fracture system (see Figure 1). The fracture system is the principle channel for fluid flow, and the matrix system is the major storage space. The fluid flow and pressure diffusion in such reservoirs are not like those in single porosity reservoirs [15]. In the double-porosity system, the fractures have higher permeability than the matrix. The fluids travel from the high permeability fractures system to the wellbore and the fracture pressure begins to decline. The hydrocarbons contained within the matrix begin to flow from matrix to nearby fractures.

The double-porosity model proposed by Warren and Root
Figure 1

The double-porosity model proposed by Warren and Root

Barenblatt et al. [16] suggested that the flow of a single phase slightly compressible fluid in a double-porosity medium could be expressed by the following equations: kfμ2pfr2+1rpfr+αμpmpf=ϕfCfpftαμpmpf+ϕmCmpmt=0(1) where p and µ are pressure and viscosity of fluid, respectively; k, ϕ, C are permeability, porosity and total compressibility, respectively; α is a geometric factor; and subscripts f and m denote the fracture system and matrix system, respectively.

Warren and Root [14] took the effects of porosity of fracture system and wellbore diameter into consideration and obtained an approximate solution for an infinite reservoir. According to the solution, the pressure change caused by a change of injection rate can be expressed as: Δp=pipfrw,t=qμ4πkfhlnηtrw2+Ei(at)Ei(aωt)+0.809=C1lnηtrw2+Ei(at)Ei(aωt)+0.809(2) Where C1=qμ4πkfh, η=kfμ(Cfϕf+Cmϕm), ω=CfϕfCfϕf+Cmϕm, a=λθω(1ω), λ=αrw2kf, θ=ηrw2; Δp is the pressure change; q is the well production rate; h is the formation thickness; Ei is the exponential integral function; rw is wellbore radius; η is the pressure conductivity coefficient.

Using the linear model of fluid productivity index and the superposition principle, the change in production rate caused by a unit injection impulse can be expressed as: Δq=C2lnηtd2+Ei(at)Ei(aωt)+0.809....t1C2lnηtd2+Ei(at)Ei(aωt)ln(t1)rw2Ei(a(t1))+Ei(aω(t1))...t>1(3) where Δq is production rate change; C2 is a new constant, C2 = C1J; J is productivity index; d is the distance between injection well and production well.

According to the definition of diffusivity filter proposed by Albertoni and Lake [7], the carbonate reservoir discrete diffusivity filter coefficients can be expressed as: α(n)=t=nt=n+1Δqdtt=0t=NΔqdtn=1,2,...,N(4) where α(n) is the diffusivity filter coefficient; N is the total number of filter coefficients. The normalization of the coefficients determines that the α(n) are less or equal to one, and the sum of all the coefficients is equal to one.

2.2 Validation of the diffusivity filters

Using two numerical simulation reservoir models to validate the diffusivity filters: one is a carbonate reservoir described by double-porosity model and the other is a sand reservoir described by single porosity model. Both are undersaturated oil reservoirs and the oil-water mobility ratio is equal to one. There is only one injector and one producer in the models. The water is injected only in the first month (one injection impulse) and the simulation extends for 12 months. Then the production rates are normalized to 0 and 1. If there was no dissipation in the reservoir, an impulse in the injector would be instantaneously produced in the producer. If there is some dissipation in the reservoir, the production response will remain for some months (see Figure 2).

The normalized production rate of numerical simulation model in carbonate reservoir and sand reservoir
Figure 2

The normalized production rate of numerical simulation model in carbonate reservoir and sand reservoir

From Figure 2, we know that the injection signal attenuation and time lag characteristics in carbonate reservoir are significantly different from those in the sand reservoir. There are lots of natural fractures in the carbonate reservoir and the permeability is larger than the sand reservoirs, so there is less dissipation in the carbonate reservoir. The production peak occurs in the first month and the impulse dissipation effect will remain for several months. There is a larger dissipation in the sand reservoir because of its lower permeability which makes the production peak much lower and occurs one month later than the impulse.

Figure 3 shows the diffusivity filter coefficients generated by Eq. 4 in the carbonate reservoir and sand reservoir. Comparing Figure 3 with Figure 2, the diffusivity filter coefficients are the same as the normalized production rate of the numerical simulation model. Therefore, the carbonate reservoir diffusivity filters can accurately describe the time lag and signal attenuation between injector and producer in the carbonate reservoir.

The diffusivity filter coefficients in carbonate reservoir and sand reservoir
Figure 3

The diffusivity filter coefficients in carbonate reservoir and sand reservoir

3 Build and solve the interwell dynamic connectivity model

3.1 Effective injection

Because of the dissipation of injection impulse, the production at time t is the sum of the effects of injection in the previous months. The effective injection rate at time t is expressed as: iijc(t)=n=0Nijαij(n)iij(tn)(5) where iijc is effective injection rate; iij is actual injection rate; and subscripts i and j denote injector index and producer index, respectively.

3.2 Interwell dynamic connectivity model

The most widely used interwell dynamic connectivity model is MLR model. In this model, the estimated production rate of a producer is given by: q^j(t)=β0j+i=1i=Iβij(iijc(t)i=1,2,3,,Ij=1,2,3,,J)(6) Where q^jt is the estimated production rate of producer j at time t; β0j is an additive constant term, when the waterflooding is balanced, β0j =0; I is the total number of injectors; J is the total number of producers; βij is the connectivity coefficient between injector i and producer j, the larger βij, the greater the connectivity.

We assume that the producer BHP is constant and the change in the production rate is exclusively caused by changes in the injection rates. Therefore, MLR model states that at any time the total production rate at well j is a linear combination of the rates of every injector plus a constant term.

The MLR model assumes that within the period of time selected for analysis, the producer BHP must be constant, no new wells and no new completions. In a real field, it is difficult to keep producer BHP constant. When the operation conditions of producer do not change, producer BHP is approximately constant.

Changing producer BHP would also affect production rates and then influence the injector-producer connectivity determination. By changing BHP of any producer, the production rates of the rest wells will be changed. When the producer BHP changes, we just shift the expected production rates by a constant rate [17]. Therefore, the modified MLR model can be rewritten as: q^j(t)=β0j+i=1i=Iβijiijc(t)+qBHPj(s)(7) Where qBHPj(s)=k=1k=KJpwfj(Ts1)pwfk(Ts)(8) is the correction coefficient; s is the time interval for BHP changes; Ts are the segmentation times; pwf is BHP of the producer; J is the liquid productivity index.

Albertoni and Lake [7] present the solution of a MLR problem. The minimization of the errors leads to the following set of I linear equations: σ112σ122σ1I2σ212σ222σ2I2σI12σI22σII2×β1jβ2jβIj=σ1j2σ2j2σIj2(9) which can be solved by standard means. The constant term is given by: β0j=p¯ji=1Iβijp¯i(10) where σii2 and σij2 are the injector-injector and injector-producer covariances, respectively.

3.3 Model validation

In order to verify the accuracy of the modified model, we use numerical simulation to generate a set of production rates from a given set of injection rates. There are five injectors and four producers in a homogeneous doubleporosity reservoir. The numerical simulation is run for 100 months under waterflooding and BHP of the producer P1 changes four times during the production period (Figure 4). Minimize the sum of squared error between actual and predicted production rates of all the producers which are calculated by the Eq. 7. Figure 5 shows a comparison between the production rate using MLR model, modified MLR model and the production rate observed in the numerical simulation. The coefficient of determination between MLR model production and numerical simulation production is R2 = 0.28, while the modified MLR model improves the coefficient of determination to R2 = 0.99 as shown in the Figure 5. Comparing the estimated β, the modified MLR model provides the most accurate estimation of the connectivity coefficient β (Figure 6). Therefore, the modified MLR model is able to predict the production rate accurately and it determines the true connectivity coefficients.

The qBHP changes for P1 well versus bottom-hole pressure changes
Figure 4

The qBHP changes for P1 well versus bottom-hole pressure changes

Production rate comparison between MLR model, modified MLR model and the numerical simulation
Figure 5

Production rate comparison between MLR model, modified MLR model and the numerical simulation

Correlation coefficient comparison between MLR and modified MLR model
Figure 6

Correlation coefficient comparison between MLR and modified MLR model

According to the actual production data of the oil field, the correction coefficient is often difficult to solve due to the lack of bottom-hole pressure data. However, we can use the time of producer production rate, which is changing dramatically as the time of BHP is fluctuating, which can be automatically identified by the computer algorithm, so an initial correction coefficient value can be given. Then the genetic algorithm is used to optimize the solution to minimize the error between the predicted production rate and actual production rate. After the time of BHP fluctuating is identified, the fluctuation of the production rate caused by the change of the production well’s own working system should also be excluded.

4 Field application

The modified technique was applied to the S48 unit in Tahe oilfield. The S48 unit is one of the main cave systems put into development firstly in Tahe oilfield. The storage space is composed of a large number of caves, fractures and pores which leads to a strong heterogeneity. The S48 unit started water waterflooding development since May 10, 2005 and the waterflood was balanced (the field injection rate is approximately equal to the total production rate). Due to the large number of fractures and vugs, the injection signal propagation is very fast in the reservoir, so the daily injection rate and daily production rate are selected as input signals. Meanwhile, the wellbore pressure fluctuation and shut-in situation happens frequently caused by water flooding method and many adjustment measures. Therefore, the modified MLR model which considers the fluctuation of bottom-hole pressure is applied.

The three-month production dynamic data was selected to study the interwell dynamic connectivity of the reservoir. The predicted production rate and the actual production rate curve of the block are shown in Figure 7, and the coefficient of determination between them is 88.7%. The average error of all production wells is 6.89%, which meets the requirement of field engineering. The S48 unit interwell dynamic connectivity diagram is plotted by using dynamic connectivity coefficient in Figure 8, where the reverse arrow points from the injector to the producer. Thus, the displacement direction is demonstrated intuitively in the reservoir. The injection rate is proportional to the dynamic connectivity coefficient of the well.

Production rate comparison between modified MLR simualtion result and the real field data
Figure 7

Production rate comparison between modified MLR simualtion result and the real field data

The interwell dynamic connectivity diagram of unit S48
Figure 8

The interwell dynamic connectivity diagram of unit S48

As shown in Figure 8, due to the strong heterogeneity of reservoir, dynamic connectivity between injectors and producers in the different directions are significantly different. The larger connectivity coefficients occur in the NW-SE trending, clearly showing the preferential permeability orientation. The main cave of unit S48, located in a NW-SE trending inherent karst highland, was formed by the first-order underground river [18]. Therefore, the injectors and producers in this direction are well-connected, which is consistent with the interwell connectivity distribution. Yang [19] has found TK410-TK425-S48-T401-TK412- T402 well groups are connected using the pressure transient analysis method before the unit had started water injection development. Many other researchers [20] also obtained the similar conclusions using other methods such as the pressure gradient method and the interwell energy interference analysis method.

After the waterflooding development, based on the injection and production performance data, interwell connectivity can be determined using our new method. The results show that TK412 is connected dynamically with well T402 and well T401, and TK425 is connected with well TK410 and well S48, which are consistent with the previous analysis. TK411 well has been injected with tracer during water injection. The tracer interpretation result demonstrates that the connectivity level of the wells S48, T401, TK467 and TK411 is the first order while the connectivity level of TK476, TK428CH and TK411 is the third order. The conclusion is consistent with the inversion result of this paper.

5 Conclusions

  1. The inject impulse dissipation in the carbonate reservoir is different from that in the sand reservoir. Due to fractures and vugs, the time lag and signal attenuation in the carbonate reservoir is smaller than that in the sand reservoir. Based on the typical pressure solution of doubleporosity reservoirs, diffusivity filter coefficient model has been built to model the propagation characteristics of injection signals in carbonate reservoirs, which can describe the attenuation and time lag of the injected signal accurately.

  2. The MLR model with fluctuating BHPs gives more reliable evaluations of interwell connectivity between injectors and producers, which greatly increases the amount of data available in practical application and improves the accuracy of inversion results.

  3. The inversion results of our method are consistent with the interwell connectivity of the reservoir. It can also characterize the connectivity between injectors and producers quantitatively as well as provide guidance for later injection. Our method is simple and convenient to apply. In addition, it doesn't influence the normal production of oil field.


This work was supported by the National Natural Science Foundation of China (No.51490654, No.51504276, No.51504277) and the Fundamental Research Funds for the Central Universities (No.16CX05018A).


  • [1]

    Yu J., A discussion on the distribution pattern of fissures, vugs and pores in Renqiu buried hill carbonate reservoir, Petroleum Exploration and Development, 1987, 14, 52-61. Google Scholar

  • [2]

    Zhou X., A primary discussion on the network-like oil and gas pools in carbonate rocks: taking the Lunnan Ordovician buried-hill pool in Tarim basin as an example, Petroleum Exploration and Development, 2000, 27, 5-8.Google Scholar

  • [3]

    Yang P., Sun Z., Liang X., et al., Seismic technology for predicting highly profitable wells in the fractured-vuggy carbonate reservoirs, Petroleum Exploration and Development, 2013, 40, 502-506. Google Scholar

  • [4]

    Li Y., The theory and method for development of carbonate fractured-cavity reservoirs in Tahe oilfield, Acta Petrolei Sinica,2013, 34, 115-121. Google Scholar

  • [5]

    Heffer K.J., Fox R.J., McGill C.A., Novel techniques show links between reservoir flow directionality, earth stress, fault structure and geomechanical changes in mature waterfloods, SPE 30711-PA, 1997. Google Scholar

  • [6]

    Refunjol B.T., Reservoir characterization of North Buck Draw field based on tracer response and production/injection analysis, MSc. Thesis, The University of Texas at Austin, Austin, USA, 1996. Google Scholar

  • [7]

    Albertoni A., Lake L.W., Inferring interwell connectivity from well-rate fluctuations in waterfloods, SPE 75225-MS, 2002. Google Scholar

  • [8]

    Yousef A.A., Lake L.W., Jensen J.L., Analysis and interpretation of interwell connectivity from production and injection rate fluctuations using a capacitance model, SPE 99998-MS, 2006. Google Scholar

  • [9]

    Zhao H., Yao J., Lv A., et al., Reservoir interwell dynamic connectivity inversion based on injection and production data, Journal of China University of Petroleum (Edition of Natural Science), 2010, 34, 91-94. Google Scholar

  • [10]

    Lv M., Yang Q., Chen K., Interwell communication in carbonate reservoir of Ordovician in Tahe oilfield, Tarim basin, Xinjiang Petroleum Geology, 2006, 27, 731-732. Google Scholar

  • [11]

    Wang X., Yan C., Yi X., et al., Well connectivity analysis of Ordovician carbonate reservoir in the fourth zone of Tahe oilfield, Journal of Chongqing University of Science and Technology (Natural Sciences Edition), 2010, 12, 52-54. Google Scholar

  • [12]

    Deng X., Gao P., Li S., et al., The method to distinguish the connectivity of the carbonate fractured-vuggy reservoir—taking the Lunguxi reservoir as example, Journal of Chongqing University of Science and Technology (Natural Science Edition.), 2012, 14, 71-74. Google Scholar

  • [13]

    Li Y., Quantitative characterization of Ordovician carbonate fracture cavity reservoir in Tahe oilfield, Journal of China University of Petroleum: Edition of Natural Science, 2012, 36, 1-8. Google Scholar

  • [14]

    Warren J.E., Root P.J., The behavior of naturally fractured reservoirs, SPE Journal, 1963, 245-255. Google Scholar

  • [15]

    Pirson S.J., Performance of fractured oil reservoirs, Aapg Bulletin American Association of Petroleum Geologists, 1953, 37, 232-244. Google Scholar

  • [16]

    Barenblatt G.I., Zheltov I.P., Kochina I.N., Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata], Journal of Applied Mathematics and Mechanics, 1960, 24, 1286-1303. CrossrefGoogle Scholar

  • [17]

    Kaviani D., Jensen J.L., Lake L.W., Estimation of interwell connectivity in the case of unmeasured fluctuating bottomhole pressures, Journal of Petroleum Science and Engineering, 2012, 90, 79-95. Google Scholar

  • [18]

    Lou Z., Zhu R., Yun L., et al., Distribution of the Ordovician fluid in the Tahe oilfield and dynamic response of cave system S48 to exploitation, Acta Geologica Sinica(English Edition), 2008, 82, 487-498. Google Scholar

  • [19]

    Yang M., Interwell communication in dissolved fracture-cavity type carbonate reservoir at block 4 in the Tahe oil fields, Xinjiang Geology, 2004, 22, 196-199. Google Scholar

  • [20]

    Zhou B., Cai Z., Li Q., A study on fluid connectivity of karst carbonate reservoir with static and dynamic data—an example of No. 4 district in Tahe oilfield, Tarim basin, Xinjiang Petroleum Geology, 2007, 28, 770-772. Google Scholar

About the article

Received: 2016-10-15

Accepted: 2017-01-05

Published Online: 2017-05-06

Citation Information: Open Physics, Volume 15, Issue 1, Pages 306–312, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2017-0034.

Export Citation

© 2017 Wei Wang et al.. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License. BY-NC-ND 3.0

Comments (0)

Please log in or register to comment.
Log in