## Abstract

Aperture, which refers to the opening size of a fracture, is a critical parameter controlling rock mass permeability. Moreover, distribution of permeability within the reservoir is commonly affected by natural fracture occurrences. In a water-based mud environment, borehole-imaging tools are able to identify both location and aperture size of the intersected fractures, whereas in oil-based environment, due to invasion of resistive mud into the fractures, this technique is impractical.

#### Authors

Hesam Zarehparvar Ghoochaninejad^{1}, Mohammad Reza Asef^{1}, Seyed Ali Moallemi^{2}

^{1}Faculty of Earth Sciences Kharazmi UniversityTehran Iran. ^{2}National Iranian Oil Company, Exploration DirectorateTehran Iran

Received: 19 April 2017 / Accepted: 25 September 2017 / Published online: 4 October 2017

Recently, some artificial intelligence techniques facilitated reliable estimations of reservoir parameters. In this paper, a teaching–learning-based optimization algorithm (TLBO) trained an initial fuzzy inference system to estimate hydraulic aperture of detected fractures using well logs responses. Comparing the results with real measurements revealed that the model can provide reliable estimations in both conductive and resistive mud environments, wherever the aperture size is unknown. TLBO, besides of its easier application, outperformed earlier optimization algorithms, which were used to evaluate the method effectiveness.

## Introduction

Natural fractures can affect essential reservoir characteristics tremendously. As a fluid conduit, natural open fractures can control the reservoir permeability (Boro et al. 2014). Conventional techniques of oil recovery may drain fractured reservoirs inefficiently, leading to early breakthrough in the secondary recovery (Alvarado and Manrique 2010; Speight 2016). In these cases, where production substantially takes place through fracture networks, aperture size is one of the basic parameters controlling reservoir permeability (Nelson 2001; Tiab and Donaldson 2015).

Mechanically, aperture is defined as the physical separation of the fracture walls (Fossen 2010). Moreover, the equivalent open part of the aperture is an indication of roughness characteristics, through which fluid can flow, and it is assumed as the hydraulic aperture as illustrated in Fig. 1 (Klimczak et al. 2010; Dardashti and Ajalloeian 2015).

**Fig. 1** A schematic illustration of fracture hydraulic and mechanical apertures with flow (Q) passing through them

Witherspoon et al. (1980) and Klimczak et al. (2010) indicated that the cubic law accurately describes the relationship between fluid flow and aperture size, considering fracture walls as smooth parallel plates. Accordingly, fracture permeability is related to the square of aperture size. Nevertheless, asperity and aperture size variation are commonly observed (Klimczak et al. 2010). Thus, the equivalent hydraulic aperture can express the relationship more accurately as follows (Sarkar et al. 2004; Watanabe et al. 2008; Ishibashi et al. 2012):

where K is the fracture permeability (millidarcy) and w_{h }is the hydraulic aperture size (microns).

During fracture characterization, the aperture size is commonly characterized with the image logs. In conductive mud environments, imaging tools can determine fracture aperture on the order of microns (Tiab and Donaldson 2015), whereas in oil-based mud environments, due to the invasion of resistive mud into the fracture planes, the imaging tools such as Oil-Base Micro Imager (OBMI) display the same appearance for both the open and healed fractures (Bishop 2015). In addition, the vast majority of natural fractures have an aperture size less than the OBMI tool pixel width (Cheung et al. 2001). Statistics have shown that the most frequent range of aperture size is between 10 and 40 microns (Golf-Racht 1982). Therefore, the presence of a quick and cost-effective method to estimate aperture size of detected fractures can be useful, wherever the aperture size is unknown.

Nakashima and Kikuchi (2007) performed a quantitative estimation of fracture apertures, of a size higher than 0.2 mm, using NMR logging data. They showed that there is a linear relationship between fracture aperture size and bulk-fluid porosity. Wu (2013) proposed a workflow to determine the relationship between aperture size and a fracture indicator value by using measurements from multi-axial electromagnetic induction well logging instruments. Ramandi et al. (2016) used X-ray micro-computed tomography to measure fracture apertures in a fractured coal sample from digital images. Shalaby and Islam (2017) have shown that fractured zones can be detected by using conventional well logs.

Recently, artificial intelligence techniques were utilized successfully for prediction of different reservoir parameters (Ahmadi 2011; Wang et al. 2013; Anifowose et al. 2017). El Ouahed et al. (2005) developed a 2-D fracture network map and fracture intensity map in Hassi Messaoud oil field by using fuzzy logic and artificial neural network. Ja’fari et al. (2012) suggested a model to estimate fracture intensity from conventional well logs by employing the adaptive neuro-fuzzy inference system. Accordingly, a fuzzy system maps input space to an output space using a set of fuzzy rules.

In this paper, a new method is introduced to estimate aperture size from well log data. The proposed method is an integration of fuzzy inference system and teaching–learning-based optimization algorithm (TLBO), in which the TLBO is utilized to train an initial fuzzy structure. TLBO is one of the most recent algorithms (Rao et al. 2011), which stands out as a powerful technique for optimization problems (Satapathy and Naik 2015; Kothavade and Deshpande 2016).

In order to evaluate the method effectiveness, various optimization algorithms available in the literature were employed to train the same fuzzy structure. Comparing the results indicated that the hybrid TLBO-Fuzzy Inference System (TLBOFIS) provides a higher level of accuracy. Finally, by estimating aperture size of detected fractures in both conductive and resistive mud environments, validation of the model was verified.

## Methodology

### Fuzzy inference system

Fuzzy set theory was first introduced by Zadeh (1965) as a method for handling data uncertainty. Unlike the Boolean logic, where variables take only a true or false value (1 or 0), fuzzy logic deals with the degree of truth and variables may have any value between 0 and 1. Capability of fuzzy sets to express gradual transitions from membership to non-membership provides a powerful representation of measurement uncertainties (Klir and Yuan 1995).

In fuzzy logic, the fuzzy inference system performs a formulation between the given input–output data, using a series of fuzzy if–then rules. As illustrated in Fig. 2, fuzzy inference system comprises five functional blocks: rule base, database (or dictionary), fuzzifier, inference engine, and defuzzifier. The rule base contains a selection of fuzzy rules, and the database defines their membership functions. Rule base and database together are usually known as the knowledge base (Jang 1993). By using membership functions in the fuzzifier block, the crisp (non-fuzzy) input data are transformed into the fuzzy sets; the fuzzy systems require a fuzzy input to be applied. The membership function defines how each point in the input space is mapped to a degree of membership between 0 and 1 (Matlab Users Guide 2013). Consequently, a set of if–then rules will be constructed, which relates the input and output datasets. The inference engine performs inference operation on the fuzzy rules. Finally, if a crisp output is required, defuzzifier transforms the fuzzy conclusions into the crisp values.

**Fig. 2** Basic structure of fuzzy inference system

There are three types of fuzzy inference systems: Mamdani (Mamdani and Assilian 1975), Tsukamoto (Tsukamoto1979), and Takagi–Sugeno–Kang (Sugeno 1985). Sugeno system is a more compact and computationally efficient model, in which its output membership functions are either constant or linear (Matlab Users Guide 2013). Considering x and y as the input values and z as the crisp output, a typical rule in the Sugeno model has the form of:

**If x is A and y is B, then z = f(x,y)**

where A and B are fuzzy sets and f(x,y) is a function in consequent. Once the fuzzy structure initialized, the membership function parameters can be trained by employing optimization algorithms.

### Teaching–learning-based optimization algorithm (TLBO)

TLBO is a metaheuristic algorithm, recently developed by Rao et al. (2011) as a method for constrained optimization problems. The TLBO inspired by the process of teaching and learning within a classroom works on the effect of influence of a teacher on learners and learners on each other. This algorithm comprises two parts: “teacher phase” and “learner phase.” During the teacher phase, a teacher tries to increase the mean result of the class from the initial level to his own level (Rao 2015). The algorithm considers the most knowledgeable learner as the teacher. In this phase, existing solution updates according to the following expression:

_{j,k,i}is the updated value of X

_{j,k,i }in subject j, learners k, and iteration i. D

_{j,k,i }is the difference between existing mean result and corresponding result of teacher in each subject.

Local search of the TLBO occurs during the learner phase. This phase uses accepted function values at the end of teacher phase as its own input. In this phase, learners enhance their knowledge by interacting randomly with themselves. A learner enhances knowledge of other learners, when he or she has more knowledge than they do.

The TLBO is a population-based evolutionary algorithm, which uses a population of solutions to proceed to the global solution (Ganesh and Reddy 2014). It is interesting to note that the TLBO does not require any algorithm-specific parameter. Common controlling parameters, such as the population size and number of iterations, are sufficient for its performance (Rao 2015). This ability of TLBO makes it an efficient tool, over the other algorithms, especially when dealing with multi-objective optimization problems (Zuperl and Cus 2016; Rai 2017).

### Hybrid TLBO-Fuzzy Inference System (TLBOFIS)

Training of the initial fuzzy inference system may be assumed as an optimization problem to be solved using optimization techniques. Firstly, by applying subtractive clustering, an initial Sugeno fuzzy system is generated. Then, the TLBO is employed to train the initialized fuzzy model. Figure 3 illustrates the flowchart of the proposed hybrid TLBO-Fuzzy Inference System approach (TLBOFIS). The methodology may be divided into three stages: (1) data preparation; (2) initializing fuzzy inference system; and (3) optimizing the fuzzy system using TLBO. The following subsections explain each stage, respectively.

**Fig. 3** Flowchart of the proposed hybrid TLBOFIS approach

### Data preparation

Here, Electrical Micro Imaging (EMI) logs along with the conventional well logs from the same intervals were the used data for initializing fuzzy model. These data were acquired from five wells drilled into two different hydrocarbon fields. The study area lies in the Zagros Basin in the Fars Province of Iran. The first field contains three wells with recorded well logs ranging in depths between 3050 and 3500 m. In the second field, which includes the other two wells, studied depths range between 2500 and 2850 m. Lithology of the logged intervals consists of anhydrite, bedded dolomite, limestone, and shale, assigned to the Permo-Triassic Dehram Group.

EMI is a water-based formation-imaging tool, which performs a qualitative micro-resistivity reading around the borehole. Six imaging pads of the EMI, in contrast to four-arm tools, provide more accurate hole volumes and diameter computations (Changxiong 2004).

Available image logs were processed using CIFLog-GeoMatrix software, in order to obtain static and dynamic wellbore images. Using these images, several features such as the natural fractures, bed boundaries, and faults can be discriminated. Figure 4 displays EMI static and dynamic images in cross section of a studied well. Commonly, higher conductivities are illustrated by darker colors; therefore, the open fractures appear as dark-colored sinusoidal features, due to the mud conductivity (Fig. 4).

**Fig. 4** A cross section of the studied wells, showing the EMI static and dynamic images. Additionally, a conductive fracture is specified

By processing available image logs, a total number of 132 conductive fractures were detected and their hydraulic aperture sizes were calculated. In the first studied hydrocarbon field, fractures show a NEE-SWW dominant strike, southeast azimuth, and dip inclination ranging between 75° and 85°. Fracture orientations in the second hydrocarbon field are disordered. Nevertheless, NE–SW strike, mean azimuth of 110° and mean dip of 78°, is dominant.

Luthi and Souhaité (1990) and Faivre (1993) showed that in a conductive mud environment, imaging tools measure the fracture aperture size as follows:

where W refers to fracture width (i.e., aperture); R_{m} and R_{xo }indicate the resistivity of mud and shallow rock, respectively; c is tool constant; and b is a value slightly smaller than 1. Excess current A is caused by the presence of fracture.

In the literature, various empirical models were proposed in order to relate the mechanical and hydraulic apertures (Renshaw 1995; Zimmerman and Bodvarsson 1996; Li and Jiang 2013). The CIFLog-GeoMatrix computes the average hydraulic aperture as the cubic mean value of fracture trace aperture.

In the next step, image logs were correlated against petrophysical well log and responses at the exact depth of each fracture were specified. In Fig. 5, measured aperture values are compared against acoustic travel time (DT), bulk density (RHOB), deep laterolog (LLD), gamma ray (GR), neutron porosity (NPHI), and photoelectric factor (PEF) logs. Among the available logs, LLD, NPHI, DT, and RHOB showed the highest correlations with aperture size and therefore, they were chosen as the input of the fuzzy model:

**Fig. 5** Cross-plots of well logs responses versus the measured hydraulic aperture (H. Aperture) for detected fractures in the studied wells. R^{2} indicates the coefficient of determination

-LLD measures formation deep resistivity and, among other well logs, has the highest correlation with aperture size. Generally, wider apertures show lower resistivities (Fig. 5).

-NPHI tracks hydrogen concentration in the formation and shows a linear relationship with measured aperture sizes, as illustrated in Fig. 5.

-DT displays travel time of elastic wave through the formation rock. Generally, void spaces reduce the speed of the sound (Serra 2008). Bigger apertures accompany with slower velocity records.

-RHOB is a record of formation bulk density and is related to the rock porosity. As seen in Fig. 5, lower densities imply bigger aperture sizes.

Final data sets, consisting of measured hydraulic aperture and well log responses, were merged and divided into two groups: a group including 98 sets for training the model and another group with 34 sets for testing.

In order to remove the effect value dimensions and achieve a better estimation, all input data were normalized between 0 and 1. Normalization was performed by the following expression:

where x_{i} is the data point to be normalized and min(x_{i}) and max(x_{i}) are the minimum and maximum values.

## Initializing fuzzy inference system

In order to generate the initial fuzzy inference system and define membership functions, a fuzzy clustering technique may be used. Subtractive clustering and Fuzzy C-Means are two basic methods of fuzzy clustering. Here, a subtractive clustering method was employed to determine the fuzzy rules. In a set of data, subtractive clustering defines the number of clusters and cluster centers and returns a fuzzy structure that contains fuzzy if–then rules. This method eliminates the need of specifying grid resolution, in which the agreement between the accuracy and computational complexity must be regarded (Bataineh et al. 2011).

Membership functions define each crisp data with a degree of membership. The membership functions can have different forms, such as the triangular, trapezoidal, and Gaussian. Here, the Gaussian membership function shows the best fit to the input well logs (Fig. 6). Gaussian membership function is represented by center and width values:

where *c*_{i} and *σ*_{i} are the mean (center) and standard deviation (width) of the *i*_{th} fuzzy set, respectively.

**Fig. 6** Histograms showing the fitted Gaussian membership function to the RHOB, NPHI, DT, and LLD well logs