Abstract
The injection-production ratio of ultralow-permeability reservoirs is generally higher in the early stage of development because of the water absorption characteristics of transition layers and mudstone formations. In this paper, the water absorption characteristics of mudstone are experimentally studied, and the empirical function of the water absorption process is established. A new mathematical model of the whole lithology is established by applying the research results of mudstone water absorption characteristics. Combining the material balance method and finite difference method, the space terms in the basic differential equation are replaced by the material balance equation, and the finite difference in the time term is obtained. Then, the analytical solutions of the average pressures of the reservoir oil well area, reservoir water well area, transition layers, and mudstone formations are solved. Based on the static parameters of the reservoir in the Chaoyang Gou Oilfield of the ultralow-permeability reservoir in China, the new pressure prediction model is verified by the ideal model of numerical simulation and production data of the oil field. The experimental results show that the saturation water absorption rate of mudstone is 1.54-2.55%, and the water absorption process of mudstone cannot be described by the seepage equation of sandstone. The verification results of the numerical simulation show that the pressure of the transition layers and mudstone at the end of the water well gradually increases, while the pressure at the end of the oil well basically remains unchanged, which is consistent with the assumptions of the model. The verification results of the oilfield production data show that the water well static pressure and oil well static pressure calculated by the new model are highly consistent with the actual values, which well explains the phenomenon of the low reservoir pressure level under the condition of a high injection production ratio in an ultralow-permeability reservoir.
1. Introduction
Due to the water absorption characteristics of transition layers and overlying mudstone formations, the injection-production ratio of ultralow-permeability reservoirs is usually greater than 1 at the initial stage of development, and the change law of reservoir pressure is quite different from that of conventional reservoirs [1–4]. The material balance method and numerical simulation method are commonly used research methods to solve these problems. Liu et al. and Fan et al. studied the variation of reservoir pressure around water wells in ultralow-permeability reservoirs at different injection-production ratio conditions using numerical simulation [5, 6]. Johnson and Rodgerson theoretically deduced the formula to calculate the pressure field around the crack in infinite elasticity [7]. Yin et al. used the finite element method to simulate the borehole wallrock stress field and pressure field evolution in the relief process [8]. Zou et al. deduced the injection-production ratio model of low-permeability reservoirs for balanced and unbalanced systems based on the extreme value principle [9]. Wu qualitatively analysed that the abnormally high injection-production ratio in low-permeability reservoirs was caused by ineffective water injection and presented the ineffective water absorption volume that accounted for 28.3% of the total water volume according to experience [10]. However, these scholars did not quantitatively introduce the water absorption characteristics of mudstone into the calculation model.
Scholars at home and abroad have conducted many experiments on water absorption characteristics of mudstone [2, 11–17]. Erguler and Ulusay selected different types of mudstone cores from different parts of Turkey to conduct mudstone water absorption experiments. The water absorption process of mudstone was divided into the isokinetic water absorption stage, accelerated water absorption stage, and isokinetic water absorption stage, but the function of the water absorption process was not provided [18]. Jiang et al. verified the water absorption characteristics of mudstone and analysed the corresponding influencing factors using computerized tomographic scanning [19]. Li et al. selected mudstone samples from Tibet for water absorption experiments and scanning electron microscopy (SEM) tests. They presented an inflection point on the water absorption curve [20]. Guo et al. studied the water absorption capacity of mudstone using the self-developed intelligent testing system for water absorption of deep soft rock [21]. He et al. also presented that the mudstone with high water invasion pressure (the situation of high injection-production ratio at the early stage of development of ultralow-permeability reservoirs) had stronger water absorption capacity than nonpressure water absorption [22].
In this paper, the water absorption characteristics of mudstone are analysed by taking mudstone water absorption experiments, and the relation between the mudstone water absorption rate and water invasion pressure is quantified by establishing an empirical function of the water absorption process. Furthermore, a mathematical model of whole lithologic pressure diffusion considering the transition reservoir and mudstone formation is established. The model can reflect the effect of water absorption characteristics of the transition layers and overlying mudstone formations on pressure in ultralow-permeability reservoirs. Moreover, it can more accurately predict the average formation pressure to determine a reasonable injection-production ratio.
2. Water Absorption Characteristics of Mudstone
The water absorption characteristics of mudstone at different water invasion pressures are studied by water absorption experiments of mudstone. The empirical function for the water absorption process of mudstone is established based on the experimental results.
2.1. Water Absorption Experiment of Pure Mudstone
The water absorption characteristics of mudstone at different water invasion pressures are studied by taking the mudstone cores of inspection wells in Chaoyang Gou Oilfield of an ultralow-permeability reservoir as the experimental object.
2.1.1. Experimental Materials and Processes
The core samples in the experiment are taken from the inspection wells in Block S of ultralow-permeability reservoirs in China. The X-ray diffraction analysis results of the inspection well core samples are statistically analysed; the basic parameters and compositions of the core samples are shown in Table 1.
The flow chart of the mudstone water absorption experiment is shown in Figure 1. The specific flow chart is as follows.
S01: core samples and determination of basic parameters. The drilled core samples were placed in a constant temperature box at 60°C for two weeks until the difference between weighings was less than 1%. Then, the permeability was measured with nitrogen.
S02: the core sample is put into the core holder, and the core holder is placed on the electronic scale, which is convenient for measuring the water absorption flux at different times. After connecting the pipeline, formation water is recharged from the outlet to the core, then the back-pressure valve is closed, and the pressure gauge is connected at the outlet.
S03: the flow pump is set in constant pressure displacement mode. The same confining pressure is loaded by the confining pressure controller, and the displacement pressure is set to the pump according to the experimental scheme. The experiment is started, and the accumulated water absorption flux is recorded at different times until the accumulated flow flux is basically unchanged, and the experiment is finished.
S04: in the experimental results, the cumulative water absorption flux is calculated and recorded according to the weight of the electronic scales at different times. At the initial stage of mudstone water absorption, data are recorded every 2 hours; after the mudstone water absorption is stable, data are recorded every 6 hours. After the experiment, the maximum water absorption flux of mudstone with constant water invasion pressure is recorded, and the saturation water absorption is calculated.
2.1.2. Experimental Result
The saturation water absorption flux (maximum water absorption flux with constant water invasion pressure) and saturation water absorption rate (ratio of saturation water absorption flux to core sample volume with constant water invasion pressure) are obtained by mudstone water absorption experiments. To provide an experimental data basis for establishing an empirical function of the mudstone water absorption process, 7 types of water invasion pressures are set, and 3 cores are made for each type of water invasion pressure. The experimental results are shown in Table 2. Taking one core as an example for each type of water invasion pressure, the curves of water absorption characteristics are shown in Figure 2.
According to the experimental results, the water absorption characteristics of mudstone are as follows. (i) When the water invasion pressure is constant, the water absorption rate increases with the increase in water absorption time, but the increase rate gradually decreases and finally remains basically unchanged. (ii) A higher water invasion pressure corresponds to higher saturation water absorption and more time required for water absorption to reach saturation. (iii) During the experiment, the pressure at the end of the core sample remains at the initial pressure of 0.1 MPa, which indicates that the water absorption velocity through the cross-section of the mudstone core is very low, and no water flows through the end of the core until the water absorption reaches saturation. Therefore, the water absorption process of mudstone is different from that of sandstone seepage, and the seepage equation cannot be used to describe the water absorption characteristics of mudstone.
2.2. Empirical Function of Water Absorption Process
Based on the experimental results of mudstone water absorption and the references [23, 24], the water absorption process function is determined. According to Figure 2, the water absorption characteristic curve of mudstone shows a negative exponential complex function relation. where is the water absorption rate, which varies with the water invasion pressure and time, %; is the saturation water absorption rate, %; is the water absorption attenuation coefficient, d-1; and is the water absorption time, d.
2.2.1. Relation between Saturation Water Absorption Rate and Water Invasion Pressure
The logarithmic relation between saturation water absorption rate and water invasion pressure is determined by fitting the statistical data of the saturation water absorption rate at different water invasion pressures. The results are shown in Figure 3. where is the water invasion pressure, MPa; is the initial pressure, MPa.
2.2.2. Relation between Water Absorption Attenuation Coefficient and Water Invasion Pressure
The power function relation between the water absorption attenuation coefficient and water invasion pressure is determined by fitting the statistical data of the water absorption attenuation coefficient at different water invasion pressures. The results are shown in Figure 4.
Therefore, the empirical formulas of cumulative water influx () and instantaneous water influx () are obtained.
Then, the velocity of water invasion per unit cross-sectional area is where is the cross-sectional area of mudstone core, m2; is the distance perpendicular to the cross-section, m. According to Equation (6), the velocity of water invasion per unit cross-sectional area of mudstone is related to the pressure of water invasion and the thickness of the mudstone layer, which provides a basis for solving the whole lithologic pressure diffusion mathematical model.
3. Mathematical Model of Pressure Diffusion for Whole Lithology
Taking an injection-production unit as the research object, i.e., a group of wells consisting of a central injection well and surrounding oil-water wells, we provide an example of the inverse nine-point well pattern.
Vertically, considering the water absorption characteristics of transition layers and mudstone formations, the geological model is divided into three parts: sandstone reservoir, transition layers, and mudstone formations. The cylindrical coordinate system is established with the central well as its pole, as shown in Figure 5.
The mathematical model established in this paper makes the following assumptions: (1) For reservoirs, there are a horizontal radial flow to oil wells and a vertical flow to the transition layers in the control area of water wells, while there is only a horizontal radial flow to oil wells in the control area of oil wells. (2) For transition layers and mudstone formations, there is only a vertical flow of injected water. (3) Both transition lithology and sandy mudstone have permeability, and the flow flux is calculated according to Darcy’s formula, while the water absorption flux is calculated according to the established empirical formula. (4) In the early stage of ultralow-permeability reservoir development, a clay stabilizer is injected at the same time of the water injection to prevent the clay from swelling and hydration migration from blocking the throat, so it is assumed that the permeability does not change during water injection.
Taking the reservoir as the research object, the continuity equation is established as follows considering the point source of the water well, point sink of the oil well, and transfer between reservoir and transition layers.
Taking the transition lithology in transition layers as the research object, the fluid continuity equation is
Taking sandy mudstone in transition layers as the research object, the fluid continuity equation is as follows:
The seepage equation considering the start-up pressure gradient is as follows:
The pressure diffusion mathematical model considering the water absorption characteristics of the transition layers and overlying mudstone is obtained by substituting Equations (10a) and (10b) into Equations (7), (8), and (9).
For reservoirs,
For transition lithology,
For sandy mudstone,
The initial condition is
The boundary condition is
4. Solution of the Mathematical Model
The pressure diffusion mathematical model is a set of partial differential Equation (5), which cannot directly derive analytical solutions. To solve this type of equations, there are two main methods: numerical difference method and simplifying the treatment to obtain analytical formulas. The former must compile a complex simulation software, while the current commercial numerical simulation software cannot equivalently handle mudstone water absorption characteristics. In this paper, the second method is used to simplify the geological model into two zones: the seepage zone in the control area of the water well and the seepage zone in the control area oil well. The schematic diagram is shown in Figure 6.
After oil wells are put into production, the reservoir pressure in the area controlled by oil wells decreases and the increases, which will cause the underground fluid transfer from the water well to the reservoir controlled by oil wells and transition lithologic. Combined with the material balance method, the partial differential Equations (11a), (11b) and (11c) can be rewritten into ordinary differential equations.
For the controlled area of water well in the reservoir,
For the controlled area of oil well in the reservoir,
For transition lithology in the transition layers,
For sandy mudstones in the transition layers
The volume flow velocity () from the reservoir to the transition lithology in the water well controlled area is
The volume flow velocity (qz2) from the reservoir to the transition lithology in the water well controlled area is
According to the empirical Equation (6) of mudstone water absorption, the volume flow velocity from sandy mudstone to mudstone in the water well controlled area is where the conductivity () can be expressed as
The intermediate transfer cross-section area in an oil well controlled area can be expressed as
The average permeability can be obtained by weighing the reservoir thickness.
The average threshold pressure gradient can be obtained by weighing the reservoir thickness.
The average distance of longitudinal intermediate transfer can be expressed as
Let the parameter group be
After substituting Equations (15a), (15b), (15c), (16), (17), (18a), (18b), (19a), (19b), (20a), (20b) and (20c) and the parameter groups (Equations (21), (22) and (23)) into the simplified pressure diffusion model (Equations (14a), (14b), (14c) and (14d)), the equations can be simplified.
The differential term at the left end of Equations (24a), (24b), (24c) and (24d) is replaced by the difference term, and the pressure at the right end of the equation is replaced by the variable of time .
The equation group about , , , and is obtained by substituting Equationa (25a) and (25b) into Equations (24a), (24b), (24c) and (24d), which has exactly four unknown variables. The average pressure in different parts can be obtained by solving the equation group.
Let the constant terms be
The above equation can be rewritten as follows:
The average reservoir pressure in the water well controlled area is calculated as follows:
The average reservoir pressure in the oil well controlled area is calculated as follows:
The average pressure of transition lithology in the water well controlled area is calculated as follows:
The average pressure of mudstone formations in the water well controlled area is calculated as follows:
The coefficients are , when calculating , , , and . Therefore, the coefficients of the above formula have definite formulas. The average pressure everywhere can be calculated from time , from which the remainder (, , ,…) can be deduced.
5. Model Validation and Application
5.1. Numerical Simulation Verification of the New Model
Although the commercial software Eclipse cannot simulate the water absorption characteristics of mudstone formations, the accurate pressure value predicted by Eclipse can be determined to be correct by setting an analytical water body to approximate the overflow and combining the historical fitting technology. Therefore, the numerical simulation results can be used as the evaluation criteria of the new model.
In this paper, the pressure prediction formula Equations (29a), (29b), (29c) and (29d) for an ultralow-permeability reservoir is validated by establishing an ideal model of the whole lithology numerical simulation. An ideal model containing 9 layers is established using the Eclipse numerical simulation software. The grid parameters are shown in Table 3, and the static reservoir parameters are shown in Table 4.
The ideal model is set to a constant flow production mode with a flow rate of 4.8 L/h, and the initial pressure of the reservoir is set at 30 atm. The comparison between the initial pressure profile of the ideal model and the pressure profile after stabilization is shown in Figure 7. The results of the Eclipse numerical simulation software and those of the new model are compared in Figure 8.
(a)
(b)
The verification results of the ideal model show that (1) under the condition of constant flow rate production, the pressure profile shows that the pressure of transition layers and mudstone at the end of the water well gradually increases, while the pressure at the end of the oil well basically remains unchanged, which is consistent with the hypothesis of the mathematical model established in this paper. (2) A higher injection-production ratio in the early stage of development is required to increase the pressure of the transition layers and mudstone formations at the water well end and gradually spread to the oil well end. (3) The pressure calculation results show that the new pressure prediction model established in this paper is consistent with the numerical simulation results.
5.2. Production Data Validation of the New Model
We take the actual production data of Chaoyang Gou Oilfield in China as the verification standard and use the average value of reservoir physical property parameters as the calculation basis of the theoretical model, as shown in Table 5. Under the conditions of different injection production ratios, the predicted average pressure in different areas of the whole lithologic reservoir is shown in Figure 9.
(a) Reservoir oil well area
(b) Reservoir injection well area
(c) Transition layers
(d) Mudstone formations
The results show that the average pressure of the oil well area and average pressure of the water well area calculated by the theoretical model are very consistent with the actual monitoring values; the average error is only 9.3%. The minimum oil well static pressure of the oil well area in the reservoir is 5.7 MPa, the maximum is 6.01 MPa, and the range of pressure change is small. However, the minimum injection well static pressure is 20.1 MPa, the maximum is 26.8 MPa, the overall range is maintained at 22.7 MPa, and the range of pressure change is large. The phenomena are mainly due to the extremely poor physical properties of the ultralow-permeability reservoir and poor connectivity between injection and oil wells. Even if the injection production ratio changes, the static pressure in the oil well area will not significantly. Simultaneously, due to the serious pressure build-up in the water well area, the initial stage of water injection will maintain a high pressure, and a large amount of injected water under the effect of high pressure in the reservoir will successively enter the transition layers and mudstone formations.
6. Summary and Conclusions
(1)The water absorption characteristics of mudstone formations in ultralow-permeability reservoirs are defined through water absorption experiments of mudstone. The water absorption rate increases with the increase in water absorption time with the condition of constant water invasion pressure, but the increase gradually decreases. The saturation water absorption of mudstone increases with the increase in water invasion pressure, but the increase gradually decreases. The water absorption process of mudstone is quite different from the sandstone seepage process, so it cannot be described by the seepage equation(2)An empirical function of negative exponential for the mudstone water absorption process is established based on the experimental results. The relation between saturation water absorption rate and water invasion pressure is logarithmic in the empirical function, and the relation between water absorption attenuation coefficient and water invasion pressure is a power function(3)A new mathematical model of the whole lithologic is established. Combining the finite difference method and material balance method, the space terms in the basic differential equation are replaced by the material balance equation, and the finite difference of the time term is made. Then, the analytical solutions of the average pressure of different regions are solved by applying the research results of mudstone water absorption characteristics(4)An ideal numerical simulation model is established, and the pressure profile shows that under the condition of constant flow rate production, the pressure of the transition layers and mudstone at the end of the water well gradually increases, while the pressure at the end of the oil well basically remains unchanged, which is consistent with the hypothesis of the mathematical model established in this paper. The verification results of oilfield production data show that the water well static pressure and oil well static pressure calculated by the new model are highly consistent with the actual value, which well explains the low reservoir pressure level at high injection production ratios in ultralow-permeability reservoirs
Nomenclature
| Rws: | Water absorption rate varying with water invasion pressure and time (%) |
| Rwsm: | Saturation water absorption (%) |
| δ: | Water absorption attenuation coefficient (d-1) |
| t: | Water absorption time (d) |
| P: | Water invasion pressure (MPa) |
| Pi: | Initial pressure (MPa) |
| Q: | Cumulative water influx (m3·s-1) |
| q: | Water invasion velocity per unit section area (m·s-1) |
| S: | Section area (m2) |
| d: | Water invasion distance perpendicular to cross section (m) |
| φ1: | Porosity of reservoir |
| p1, p2, p3: | Fluid pressure of reservoir, transition lithology, and sandy mudstone (MPa) |
| Dp: | Threshold pressure gradient (MPa·m-1) |
| , : | Average formation pressure of reservoir in the area controlled by water and oil wells (MPa) |
| , : | Average pressure of transition layers and mudstone formations in the well control area (MPa) |
| Vw1, Vo1: | Volume of reservoir within the control range of water and oil wells (m3) |
| Vw2, Vw3: | Volume of transitional lithology and sandy mudstone in the controlled range of the water well (m3) |
| qw, qL: | Injection rate of the water well and fluid production of the oil well (m3·d-1) |
| qwi, qlj: | Water injection intensity of water well and fluid production intensity of oil well (m3·(d·m)-1) |
| qk1, qk2, qk3: | Water absorption intensity of transitional lithology, sandy mudstone, and mudstone (m3·(d·m2)-1) |
| qk4: | Transfer intensity of unclosed spillover of sandy mudstone (m3·(d·m2)-1) |
| λ1, λz1, λz2, λz3: | Conductivity in reservoir, conductivity from reservoir to transition lithology, conductivity from transition lithology to sandy mudstone, and conductivity from sandy mudstone to mudstone |
| , : | Average permeability from the middle of the reservoir to the middle of the transition lithology and average permeability from the middle transition lithology to the middle of the sandy mudstone (10-3 μm2) |
| , : | Average threshold pressure gradient from the middle reservoir to the middle of the transition lithology and average threshold pressure gradient from the middle of the transition lithology to the middle of the sandy mudstone (MPa/m) |
| Ce1, Ce2, Ce3: | Total compressibility of the reservoir, transition lithology and sandy mudstone (MPa-1) |
| A: | Cross-section area of flow (m2) |
| L: | Distance between the center of the water well controlled area and the center of the oil well controlled area (m) |
| s1, s2, s3: | Cross-section area of transition lithology, sandy mudstone, and mudstone (m2) |
| h1, h2, h3: | Thickness of the reservoir, transition lithology, and sandy mudstone (m) |
| d1, d2, d3: | Distance from the middle of the reservoir to the middle of the transition lithology, from the middle of the transition lithology to the middle of the sandy mudstone, and from the top of the sandy mudstone to the middle of the mudstone (m) |
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was funded by the Ministry of Science and Technology of the People's Republic of China, and the name of the project is “Demonstration Project for Development of Dense Oil in Jiyang Depression, Bohai Bay Basin (2017ZX05072)” of the National Science and Technology Major Special Project.