Log‐based rock compressibility estimation for Asmari carbonate formation

Rock compressibility has many applications in the upstream petroleum industry, for example, reservoir material balance calculations and geomechanics, as related to possible formation compaction and subsidence. When measurements on core are not available, empirical correlations may need to be considered, derived for specific fields. In the study presented, the suitability of the log-based method for evaluating rock compressibility is investigated, with the aim of deriving a suitable correlation for Asmari carbonate, Southwest Iran. A stepwise log-based method of rock compressibility determination is presented, including the construction of a geomechanical earth model, MEM. In comparing a log approach with core-derived results, it was found that the log-based method showed a rather reliable estimation of pore compressibility for carbonates studied, except for intervals with extremely large wellbore washout. Although the selected procedure is similar to that used for sandstone by other researchers, involving in some cases complex procedures, a relatively good matching result could be obtained with a relatively simple procedure. To validate the results from this study, a comparison of compressibility values is made with several industry correlations. An application of formation compressibility for estimating possible formation compaction and associated subsidence is also presented. However, a minimal effect is indicated in this case, as the carbonates studied are rather consolidated. In summary, this study presents a new correlation for rock compressibility for Southwestern Iranian carbonate oilfields (Asmari carbonate), validating the applicability of a log-based technique.


Introduction
Rock compressibility is a significant parameter in reservoir studies, for example, material balance as related to estimating (remaining) reserves. Another application involves the prediction of possible formation compaction and associated surface subsidence, a complex geomechanical problem. Rock compressibility has also application in well completions and associated sand control operations (Ghalambor et al. 1994). Core-based compressibility measurements, especially for consolidated rock, are always the first choice. However, when such data are not available, alternative techniques need to be considered.
This study presents an investigation of a cost-effective log-based method of compressibility estimation for Asmari carbonate formation. The approach is validated by comparison with core-derived measurements. Compressibility results were also used to evaluate the possibility of compaction-related subsidence.
It is common practice by the industry to use published, empirical correlations of rock compressibility, presented at net over-burden pressure, NOBP versus atmospheric porosity, in cases when core measurements are unavailable. Correlations are derived for specific NOBP conditions (reservoirs, fields or regions) and therefore may not be directly applicable for another region. Ideally, correlations should be derived for specific fields, to minimize possible errors. This study fills the gap in considering Iranian, Asmari carbonate formation, as demonstrated by comparison with existing correlations reported in the literature. For each existing industry correlation, deviations for the study dataset were determined, resulting in the identification of the most applicable correlation. The identified correlation can then be deployed for other fields in the region by simply knowing atmospheric porosity.

3
A concern is that core compressibility measurements may not always be reliable, specifically as derived from hydrostatic tests. There is also the possibility of invalid, failure-related results from uniaxial testing. The less expensive, log-based technique presented here may prove to be advantageous, especially when the formation being considered is relatively uniform, as is the case for the carbonate formation under consideration. The methodology selected is similar to more standard techniques used for sandstone formations; more complex methods such as those proposed by Saxena (2011) were not considered. Log-based estimates are compared with core-based evaluation, imposing a cutoff to identify unreliable results.
Finally, using the compressibility values derived for the reservoir under consideration, following a rigorous stepwise mathematical description, the compaction-related subsidence for the reservoir is investigated, following decades of production.

Review of rock compressibility definitions and applications
Rock stress and related strain, and their relations have been in use for considerable time. In solid mechanics, for elastic materials, the following relationship holds between strain, and the effective stress, σ (Hakiki and Shidqi 2018): where E is Young's modulus of elasticity.
In fluid mechanics, the following relationship holds between shear stress, τ and shear rate (strain), γ: where μ is fluid viscosity.
For almost all reservoir rocks which are elastically deformable, the relationship between volume and loading is nonlinear (van der Knaap 1959). The coefficient of compressibility has been defined as the strain divided by stress, which is comparable to the reciprocal of Young's modulus in solid mechanics and to the reciprocal of viscosity in fluid mechanics. The coefficient of isothermal compressibility is generally defined as the relative volume change of a material (i.e., strain) per unit change in effective stress under conditions of constant temperature. The compressibility coefficients (which are always positive values) are defined separately for the bulk, rock matrix and pore space as the pore pressure varies (after Zimmerman et al. 1986): is the coefficient of bulk compressibility as the pore pressure drops [psi −1 ]; the first subscript denotes the type of compressibility and the second subscript denotes the changing pressure, C mp (or simply C m ) is the matrix compressibility coefficient [psi −1 ], C pp (or simply C p ) is the pore volume compressibility coefficient [psi −1 ], V b , V m , and V p are the bulk, matrix and pore space volumes [cm 3 ], respectively, The subscript "c" indicates that confining pressure is constant; "p" indicates that the pore pressure is constant.
NOBP is the net over-burden pressure or effective stress [psi].
The net over-burden pressure, NOBP is the effective stress, with the definition first given by Terzaghi (1936), and later modified by Biot (1941) to: σ V is the vertical or overburden stress (psi −1 ), which can be found from bulk density logs (Ramiah et al. 2019) or from the estimated overburden stress gradient in a region, α is Biot's coefficient [dimensionless], and P p is pore pressure [psi −1 ], which is typically measured by repeat formation tester (RFT), modular dynamics testing (MDT), well test analysis, logging while drilling sensor (Ashena et al. 2020), or can be estimated using Dc-exponent or kick data from daily drilling reports (Ashena et al. 2020;Farsimadan et al. 2020).
Comparing different compressibility terms, the following relationship holds between the mentioned compressibility coefficients (Domenico 1977): As C m is an order of magnitude smaller than C p , it may be ignored (Laurent et al. 1993): Combining Eqs. 7 and 8, C b can be related to C p by: Alternatively, it can be expressed as: Pore volume compressibility, C p , which is also named rock compressibility C r , has two main applications in the petroleum industry: • Reservoir simulation: In any (dynamic) reservoir simulation software, rock compressibility is traditionally required as the only geomechanical input parameter which can have significant effect on the results (Falcao et al. 2015). The effect is included in material balance calculations and the evaluation of original-oil-in-place, especially in undersaturated volumetric reservoirs (i.e., oil reservoirs without any gas caps) and also when the limits of the field are unknown or undefined. In simple simulation models, a single, constant compressibility value is entered for the entire reservoir simulation model. However, using a relationship for compressibility would better represent the formation stress-strain behavior for a particular reservoir situation (Falcao et al. 2015). In any (dynamic) reservoir simulation software, rock compressibility is traditionally required as the only geomechanical input parameter which can have significant effect on the results (Falcao et al. 2015). The effect is included in material balance calculations and the evaluation of original-oil-in-place, especially in under-saturated volumetric reservoirs (i.e., oil reservoirs without any gas caps) and also when the limits of the field are unknown or undefined. In simple simulation models, a single, constant compressibility value is entered for the entire reservoir simulation model. However, using a relationship for compressibility would better represent the formation stress-strain behavior for a particular reservoir situation (Falcao et al. 2015). • Compaction and subsidence: Having determined pore compressibility, possible formation compaction and associated subsidence at surface can be estimated. Using core and log-based estimates of rock compressibility, this study shows how the magnitude of compaction is estimated for a specific case involving Asmari carbonate. Having determined pore compressibility, possible formation compaction and associated subsidence at surface can be estimated. Using core and log-based estimates of rock compressibility, this study shows how the magnitude of compaction is estimated for a specific case involving Asmari carbonate.

Methods for estimating rock compressibility
There are two methods discussed in this paper to estimate rock compressibility:

Selection of best-matching industry correlation
Core compressibility measurements may be unavailable for a region under study, or the core compressibility test may be technically difficult to execute, produce questionable results or is very costly, e.g., unconsolidated formations. In such cases, published correlations of rock (or pore) compressibility C r at initial net over-burden-pressure, NOBP i vs atmospheric porosity, i.e., rock porosity at the surface, is typically used by the industry (Hall 1953;Newman 1973;Khatchikian 1995). Table 1 shows several published correlations for C r for carbonates as a function of atmospheric porosity. Most of these empirical correlations have been developed only for well-consolidated rock, except for Horne (1990) and modified Horne correlations (Jalalh 2006), derived for both, consolidated and unconsolidated sandstone, and also carbonates. For very shallow sandstone formations, containing soil or clay, soil compressibility models should be considered (Santamarina et al. 2019). For consolidated formations, the lower the porosity, the greater the pore compressibility (van der Knaap 1959). However, this trend may be the opposite for friable and unconsolidated formations. Deviation from expected trends, e.g., the modified Horne correlation may be related to other factors such as in situ depth and effective stress, as well as specific rock type; if such factors are ignored compressibility estimates are often in error (Newman 1973).
Generally, to estimate compressibility for a field or region in a cost-effective manner, it is recommended to first take core samples and conduct compressibility laboratory tests in some wells of the field. Next, using the available laboratory data, it is recommended to find the best-matching correlation for the region, among the available correlations. This study follows such recommendation, comparing estimates of compressibility derived from measurement with various suitable correlations on a comparative basis. Error statistics have been used to gauge the goodness of fit, including Mean Square Error, MSE, and Average Absolute Percent Error, AAPE. The best-matching correlation is selected as the one with the lowest error indicator, particularly MSE.

Log-based method
The second approach used in this study to estimate compressibility is the so-called dynamic measurement approach, the use of petrophysical well log data. As well logs, including sonic are acquired under in situ reservoir conditions of stress, essentially uniaxial, log-derived compressibility should be equivalent to uniaxial strain compressibility, which is an advantage (Khatchikian 1995). The log-based method can also offer an option of early pore compressibility estimation for field development without extensive SCAL at high cost (Saxena 2011). Finally, the log-based method has the added benefit of providing a continuous profile of the computed pore compressibility with depth which allows for convenient compaction calculations (Wolfe et al. 2005). Despite all these advantages, dynamic pore compressibility tends to underestimate compressibility due to the existence of microcracks (Geertsma 1957), unless the dynamic effect is considered.
Various studies for the log-based approach have been reported in the literature, but most of these are for sandstones, for example, Ong et al (2001) and Wolfe et al (2005), using a mathematical model comprising computation of rock mechanical properties. There are a few literature references for carbonates, e.g., Saxena (2011), using a differential effective medium (DEM) model, applicable for complex carbonates. However, the methodology applied in the study reported here was not the complex methodology of Saxena but rather similar to that of Ong et al and Wolfe et al., assuming that matrix compressibility can be assumed to be negligible for carbonates. The procedure is outlined below.
After gathering petrophysical well logs and related information: gamma ray, dipole sonic slowness, neutron density, interpreted lithology and porosity, a quality check was executed first. Comparing caliper log readings with hole size, wellbore washouts due to sloughing of the fractured carbonate rocks were identified, with extreme washouts resulting in log readings deemed unreliable. Next, the available laboratory-measured, geomechanical data on cores from this well were collected.
In order to make a valid comparison of compressibility, log-derived with that from core, NOBP values should be evaluated using Eq. 6, which requires vertical, overburden stress, pore pressure and Biot's coefficient. For the logged interval of 8383-9187 ft (2555-2800 m), the vertical stress, σ V is found using the bulk density log. For the non-logged section, above 8383 ft (2555 m), the overburden pressure gradient of 1.1 psi/ft was chosen, commonly in use for the region. The pore pressure, P p was obtained from wireline logs. Biot's coefficient, α was assumed equal to 0.9, to be compatible with the value considered in the core-based method. Having gathered all of the required data, the following procedure was followed to find the estimated log-based pore compressibility:

Log-based procedure
To estimate C p vs NOBP, a geomechanical earth model (MEM) was essentially constructed, using as input the petrophysical log data and the pore pressure profile. Assuming zero matrix compressibility and the elastic rock, the following five steps were taken to find the pore compressibility for each NOBP:

Estimation of K b,d
The bulk modulus represents the effective pressure, NOBP change with respect to the relative bulk volume change, Using the following log parameters: acoustic compressional wave velocity, V p , shear slowness, V s , and the density, ρ b , the dynamic bulk modulus, K b,d was estimated (Gatens et al. 1990;Yu and Smith 2011;Munir et al. 2011): Equation 12 may be alternatively written as shown in Eq. 13, using the reciprocal relation of wave velocity and slowness ( V p = 1 Δt p and V s = 1 Δt s ):

Conversion of K b,d to K b,s
The next step involves the conversion of the dynamic bulk modulus (K b,d ) to the static one (K b,s ). Based on Fjaer et al (2008), the bulk modulus is related to Young's modulus, E and Poisson's ratio, ν as follows: Writing Eq. 14 both for the static and dynamic cases and taking the ratio gives: Assuming ν s = ν d in Eq. 15, appropriate for an elastic medium (Larsen et al. 2000): Alternatively: Some log-based compressibility estimation approaches in the literature have attempted converting the dynamic to static moduli (Khatchikian 1995) by considering only drained conditions, i.e., converting undrained conditions to the drained dry situations using Gassman's equation (Gassmann 1951;Domenico 1977). However, these references are less comprehensive as they do not consider all the reasons for arising discrepancies, as differences between dynamic and static moduli may be attributable to several causes. Discrepancies arise due to different drainage conditions, different strain amplitudes and rates, possible rock heterogeneity and anisotropy (Fjaer 2019). Firstly, comparing drainage conditions gives an indication that deformations induced by elastic dynamic waves in a fully saturated rock are always undrained, while "static" deformations are often drained (Gassmann 1951;Biot 1955). Secondly, different frequency or strain amplitudes need to be considered (Plona and Cook 1995), where strain amplitudes are typically of the order of 10 −7 for dynamic tests and 10 −2 to 10 −3 for static tests. A closed crack or an uncemented grain contact may remain immobilized during the oscillating, low amplitude deformation associated with a passing elastic wave, whereas it may become mobilized as a result of static loading (Walsh 1965). This phenomenon was confirmed by other researchers (Cheng and Johnston 1981;Jizba and Nur 1990;Khatchikian 1995;Fjaer 1999), who found that the difference decreases at high confining stresses, inferring that microcracks close under confining pressures (Geertsma 1957). As an example, Cheng and Johnston (1981) found that the ratio of static to dynamic modulus for some sandstone samples varied from 0.5 at atmospheric conditions to 0.9 at 2800 psi. Based on Eissa and Kazi (1988), the ratio of E s /E d usually ranges from 0.5 to 1. As To consider the above issues and possible discrepancies, it is recommended to use core laboratory data to relate the dynamic (log-based) E d to the static Young's modulus from core, E s . In this approach, the dynamic and static bulk moduli are also correlated. There are several correlations between dynamic and static Young's moduli (Elkatatny et al. 2018). The approach followed here was that used by Wang (2000) and Ashena et al. (2020): Therefore:

Conversion of K b,s to C b,s
Generally, the bulk compressibility coefficient, C b is defined as the reciprocal of the bulk modulus, K b . The static bulk compressibility, C b,s can then be expressed as: Similarly,

Conversion of C b,s to C p
Static bulk compressibility data (C b,s ) is essentially converted to static pore compressibility (C p,s ), simply named pore compressibility (C p ). Considering Eq. 10, C p (or C r ) is found for the log depth interval:

Final calibration
During the process of determining C p values using petrophysical well logs and static Young's moduli from core data were considered, to find a calibrated correlation for conversion of the dynamic moduli to static ones, Eq. 18. In comparing log and core-derived values for C p , and in case there is a mismatch between the log-estimated value and the measured value on core, the correlation converting the dynamic to static moduli should be updated until the best match is obtained. Based on the results, if the log-based C p values are found near the core measurements, this stage may be ignored, as was the case with this case study. It should be noted that some discrepancy between the laboratory results and the log-based compressibility estimates is considered natural due to: • a possible error in the conversion of the dynamic to static moduli using empirical correlations in the log-based method, • differences between the hydrostatic nature of experimentally derived compressibility and real anisotropic in situ stress conditions, and • differences in the fluid type in pores.
Different fluid contained in the core samples (gas, oil or salt water) influence measured pore compressibility. This variation is often ignored by petrophysicists and core lab technicians. Depending on the type of measuring apparatus, the fluid type in the core samples is different: • In the Core Compressibility Measurement System (CMS), dry (air-bearing) core samples are used, resulting in measured Cr values that are not fully representative of oil-bearing core samples. This situation causes a difference between the core laboratory compressibility measurement and the real one. In this study, laboratory compressibility data were obtained using the CMS apparatus. • In the Rock Compressibility System (RCS), brine-containing core samples are used. The fluid error is much lower using RCS measurements compared to CMS measurements.

Best-matching correlation
To develop rock compressibility correlations, 45 laboratorymeasured compressibility data values were utilized from an Asmari carbonate formation, representative of Southwest Iran. These measurements were plotted on a graph with various industry correlations in the background, Fig. 1. These correlations are: Hall (1953); Newman (1973) and modified Horne (Jalalh 2006). Comparing measurements with each of the correlations and using the error metrics described earlier, Table 2 shows that the modified Horne correlation gives the best overall match to the data, with the lowest MSE error of 3.125, the Absolute Root MSE (ARMSE) of 1.767, and an AAPE of 24.44%. Conversely, Newman's equation was identified as the poorest matching correlation. It is, therefore, recommended that the modified Horne correlation should be used to estimate rock compressibility for other Asmari carbonate formations in the region. Equations for error statistics used are listed in "Appendix 1". Taking a closer look at the distribution of data points in Fig. 1, it can be observed that all correlations converge at higher porosity values where variations in lithology are less sensitive. Conversely, lithology effects tend to be amplified for lower porosity values, particularly in the presence of micro-fractures/ cracks and vugs due to leached marine shell fragments, commonly encountered in Asmari carbonate formations. More elevated compressibility values observed ( Fig. 1) at lower porosity values are probably related to the lithology variations mentioned. It should also be mentioned that Newman's correlation was derived for less consolidated carbonates, compared to Hall's correlation, where it can be seen that Newman's correlation is a good representation of the more elevated compressibility values observed for samples from this study.  Log-based method Figure 2 shows petrophysical well logs, consisting of gamma ray, dipole sonic slowness, neutron, density, interpreted lithology and porosity for the well under study. Comparing the caliper log readings with the hole size, one can observe the occurrence of wellbore washout due to sloughing of the fractured carbonate rocks. The adverse effect of washout has been compensated in most intervals, except at some intervals (indicated in the same figure) with extremely high washout. For these intervals, the log readings are deemed to be unreliable. Figure 3 shows the constructed MEM showing the vertical stress, and current and initial pore pressures (PP and PP i , Fig. 2 Petrophysical well logs for the well under study, including interpreted lithology and porosity required for estimation of rock moduli, strength and pore compressibility of the formation for the interval 8383-9187 ft (2555-2800 m). Due to some zones with excessive washout of the wellbore, which can be observed by comparing the caliper and the hole size, log measurements for these zones could not be corrected and are unreliable. Legend: TVD, true verti-cal depth; Az, azimuth; Inc, inclination angle; DEG, degree; CGR, compensated gamma ray; GAPI, gamma ray unit; Cal, caliper; BS, bit size; IN, inch; DTCO, acoustic compressional slowness; DTSM, acoustic shear slowness; US/FT, micro-second per foot; NPHI, neutron porosity; V/V, (pore) volume to (bulk) volume ratio; RHOB, bulk density, G/CC, gram per cubic centimeter respectively), the current and initial NOBP, rock strength (UCS), the dynamic and static Young's moduli (E d and E s from core data), the dynamic and static bulk moduli (K b,d and K b,s , respectively), static bulk compressibility (C b,s from core data), porosity and pore compressibility (C p ). Table 3 shows the available geomechanical laboratorymeasured data of core samples for this well under study, which consist of two points for Young's modulus and three points for core-based pore compressibility. As shown in Table 4, the log-based pore compressibility values show some discrepancy from the core laboratory data: the greatest difference being 58% (absolute percent error, AAPE) at the measured depth, MD of 8,384.64 m, with the lowest difference of 19.7% at the MD of 8,461.15 m. As noted above, Legend: Sv, Vertical Stress; PP_i, initial formation pore pressure; PP, current pore pressure; NOBP, net over-burden pressure; NOBP_i, ini-tial net over-burden pressure; UCS, uniaxial compressive strength; Ed, dynamic Young's modulus; Es, static Young's modulus; Kbd, dynamic bulk modulus; Kbs, static bulk modulus; Cbs, static bulk compressibility; Cp, pore compressibility; Cp (core), pore compressibility (from measured core samples) this discrepancy can be attributed to a possible combination of effects: log values due to washouts, the hydrostatic nature of the experimental pore compressibility measurements, the discrepancy in the fluid type, and the possible error in the conversion of the dynamic to static moduli. In layers with unreliable well log measurements, particularly the ones indicated in Figs. 2 or 3, extremely high pore compressibility values are indicated, even greater than 50 × 10 −6 /psi, which may be invalid. For layers with extremely low porosity values, in the order of 10 −4 , typically occurring in thin and tight carbonate layers, extremely high compressibility values can be observed. Considering the possibility of low porosity values in carbonates, a cutoff for pore compressibility may be selected, say 15-20 × 10 −6 /psi may be more appropriate, based on this study.

Application to compaction and subsidence
A geomechanical application of compressibility data is the estimation of possible formation compaction and associated subsidence. There are two extreme end-members when considering surface subsidence due to reservoir compaction: (1) plastic flow, the result of competent zones between the reservoir and the surface, and (2) extreme elastic behavior, where there are competent zones present. In the first mentioned case, associated with shallow reservoirs, the surface subsidence bowl is situated directly above the reservoir and its outline mirrors that of the field boundaries below, where the subsidence volume would equate to reservoir voidage using material balance. In the second case, where competent members, e.g., substantial shales are present, the subsidence bowl is aerially significantly larger than the size of the field due to the shales taking up some of the energy through deformation, resulting in a subsidence bowl volume which is significantly less than reservoir voidage. During hydrocarbon production, the depletion of the reservoir pore pressure leads to significant decrease in the total horizontal stresses (or increase in the effective horizontal stresses), whereas the vertical, overburden stress remains unchanged (Wolfe et al. 2005). Based on the amount of depletion-induced reservoir compaction, adverse effects, such as subsidence, casing deformation and seismicity may occur. Nagel (2001) gives an overview of compaction and subsidence in the petroleum industry with several field examples.
Many compaction and subsidence studies (e.g. Raghavan et al. 1972;Jelmert and Toverud 2018;Jelmert 2019) are based on certain assumptions to estimate compaction-related subsidence, including constant grain volume, the validity of Darcy's law, instantaneous response to fluid pressure changes, constant overburden stress, and that flow in the vertical direction due to compaction may be neglected. An interesting subsidence study was undertaken for the large Groningen gas field, Netherlands by Hettema et al (2000), utilizing a relatively simple approach to estimating subsidence was applied, relating reduction in reservoir thickness over several decades of production to the subsidence bowl at surface. Considering this approach, an equation can be derived in terms of (absolute) change in reservoir thickness, Eq. 22. For derivation of this equation, see "Appendix 2".
where φ 2 is the current porosity at the time of the study; ΔP p = P p,2 − ΔP p,1 (subtraction of the initial pore pressure from the current pore pressure, resulting in negative values).
The above equation for subsidence estimation is comparable with the developed equation by Fjaer et al (2008) assuming linear poroelasticity for a homogenous, isotropic reservoir rock: where E and ΔP p should be both in the same units, e.g., psi. The negative sign has been used behind ΔP p to make it positive, and thus, the change in reservoir thickness or ground level is in absolute terms, |Δh|.
To estimate the change in reservoir thickness with both methods, a stepwise procedure is outlined. Starting with reservoir thickness, h 1 is the difference in true vertical-depth from the mean-sea-level (MSL), TVD MSL (7,471-8,078ft or 2,277-2,462 m), average values for φ 2 and ΔP p were considered for use in Eq. 22; average values for E s , ν and α were considered in Eq. 23. The average value for φ 2 was found by taking the average of all log-based porosity values (Fig. 3); the average ΔP p in the interval was found by subtracting the initial reservoir pressure from the current pore pressure (a negative value); the C p values were used from the available core laboratory measurements (Tables 3 or 4); the average values for E and ν were found by averaging from Fig. 3, ignoring values for unreliable layers; the average value for α was assumed to be 0.9.  Table 5 summarizes all parameters and also gives the determined subsidence, utilizing both formulations. A surface subsidence of less than 0.4 ft (0.12 m) was estimated to have occurred following 70 years of production history. The reason for this relatively small amount of subsidence, despite an excessive reservoir depletion of 2330 psi, is that the carbonate rock is rather consolidated and the reservoir thickness is just ~ 600 ft. Recently, using gas injection, the pore pressure has been raised to prevent the continuity of this subsidence phenomenon.

Conclusions
Two methods for obtaining formation compressibility were presented, considering an example for Asmari carbonate formation, Southwest Iran. A comparison was made between laboratory-derived core measurements and values obtained from open-hole logs. These estimates where subsequently used to derive estimates of compaction and related subsidence, based on two alternative formulations.
• To validate laboratory compressibility measurements, it is good practice to compare these against industry correlations. For the Asmari carbonate under study, the modified Horne correlation was found to be the best match to specific laboratory measurements. • Matching of log-based compressibility with specific corederived measurements, in this case, Asmari carbonate, is a suitable approach for determining the best overall compressibility for a reservoir. The log-based method utilized was similar to that used for sandstones, following construction of a geomechanical earth model (MEM), resulting in a log of compressibility versus depth. Comparison of the logbased and core data gave a similar result, omitting intervals with large wellbore washout. Some discrepancy was also attributed to the hydrostatic nature of core compressibility measurements. Another reason for some elevated carbonate values observed is related to the diagenetic nature of this, typically micro-fractures and vugs.
• The magnitude of possible compaction related subsidence, assuming plastic flow was estimated for a compaction related reservoir thickness change, assuming reservoir depletion resulting in a pressure decrease of 2330 psi, following 70 years of production. An evaluation was performed utilizing two alternative formulations. In this case, due to the competent nature of the carbonate studied, possible subsidence is rather minimal, of the order of 0.4 ft (0.12 m).
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.

Appendix 2: Derivation of subsidence equation
The derivation is given in three steps, as follows:

Pore volume change
During hydrocarbon production and assuming reservoir pressure depletion, NOBP increases proportionally, where dNOBP = −dP p . Using this equality and integrating both sides of Eq. 5 for C p , gives the following for V p : The above equation can be rewritten as: where ΔP p = P 2 − P 1 , which is negative.

Porosity change
To find the change in porosity, Eq. 5 is considered in terms of C p , where V p /V b is replaced by φ, and dNOBP is replaced by −dP p . resulting in: Simplifying: Rearranging and integrating Eq. 30 gives the following result: where ΔP p = P p,2 − P p,1 , which is negative.
Alternatively, Eq. 31 can be obtained by dividing both sides of Eq. 29 by the bulk volume, V b , assuming the volume remains constant during compaction. However, ignoring the bulk volume change resultants in an unacceptable, final porosity, φ 2 . To find a more accurate equation for φ with respect to pore pressure, a bulk volume change with pressure is considered. Assuming the matrix volume change to be zero during a pressure change: Dividing both sides of Eq. 29 by the two sides of Eq. 32 gives: Multiplying and dividing both sides of the above equation by the corresponding bulk volumes (V b,2 and V b,1 ), gives: Considering the definition of porosity gives: Rearranging gives the final result:

Subsidence
To find the equation for subsidence or ground level change, and assuming plastic flow, the ratio of heights (h 2 /h 1 ) is proportional to the ratio of bulk volumes as follows: Assuming the matrix volume remains unchanged during compaction, and multiplying and dividing the right hand side of Eq. 35 by the equivalent matrix volumes in Eq. 32 gives: The right-hand side of the above equation can be rearranged as follows: Simplifying gives: Combining Eqs. 36 and 33 gives the following: Replacing 1 2 in Eq. 37 with its equivalent in Eq. 34 results in the following: Using Eq. 37, the change in reservoir thickness, Δh can be expressed as: And in absolute terms: Equation 38 can then be rewritten as: If the original porosity before any substantial production is not available, φ 1 in Eq. 41 is replaced by its equivalent in terms of φ 2 , taken from Eq. 34, resulting in the alternative expression: where φ 2 is the porosity after compaction.