The definition of productive reservoirs by the method of fuzzy inference
S. V. Denisov^{1} , V. E. Lyalin^{2} , K. S. Chernyshev^{3} , M. V. Danilov^{4} , V. P. Tsepelev^{5}
^{1}Ufa State Petroleum Technological University, Ufa, Russia
^{2, 3, 4, 5}Kalashnikov Izhevsk State Technical University, Izhevsk, Russia
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 19, 2018, p. 241245.
https://doi.org/10.21595/vp.2018.20060
Received 3 July 2018; accepted 16 September 2018; published 24 September 2018
The lithological classification and allocation of reservoirs are based on the difference of physical and geophysical parameters of rocks. Finding the values of physical and geophysical parameters in some ranges makes it possible to predict the lithology of the formation. The class of information systems based on fuzzy logic provides an effective use of existing knowledge about a certain object. This allows to apply fuzzy logic in intelligent systems of interpretation of well geophysical research materials.
Keywords: well logging, productive collector, fuzzy logic, interpretation of materials.
1. Introduction
The lithological classification and allocation of reservoirs are based on the difference of physical and geophysical parameters of rocks. Finding the values of physical and geophysical parameters in some ranges makes it possible to predict the lithology of the formation. Since the ranges overlap for different rocks, it is necessary to identify the lithology of rocks by a set of reservoir, physical and geophysical parameters.
In automated systems, interpretation of geophysical research wells used algorithm lithological division with a scan to dam [1, 2]. The algorithm consists in compiling a complex code according to the data of geophysical surveys of wells and comparing it with predictive diagnostic codes representing the geophysical characteristics of rocks in the section. Approximate interval estimates are used instead of numerical values. For example, small values are assigned the code 00, medium – 01, large – 10, very large – 11. The set of different geophysical measuring methods, recorded in the accepted order, forms a complex code of lithological type. The comparison of the complex formation code with the diagnostic one makes it possible to determine the type of lithological rock in this formation. We note two main drawbacks of this method. The first is that in case of mismatch of the actual code with any diagnostic, lithological type of formation remains uncertain. The second disadvantage is the need to evaluate the boundary values of the geophysical methods in the preparation of complex codes. From the method of division of changing parameters on the intervals strongly depends on the form of a comprehensive code.
The further development of this method is the algorithm of lithological dissection of the section with probability estimation [8]. For this formation, the most probable lithological type to which it should be assigned is determined by the complex of indications of geophysical methods [1]. This algorithm eliminates the first drawback of the method of lithological dismemberment by diagnostic codes, but is not exempt from the need to specify intervals. In addition, the assumption is made about the normal distribution of the indications of the methods of geophysical studies of wells.
These algorithms are not representatives of intelligent systems, as they are not selfadjusting and selflearning. Class of information systems based on fuzzy logic [9] provides a more effective use of existing knowledge about a certain object. This allows to apply fuzzy logic in intelligent systems of interpretation of well geophysical research materials.
The advantages of fuzzy expert systems are: the possibility of parallel implementation of the existing rules and forecasting new States of the system; the plurality of interpretations of the values of variables, as well as the fact that the description of the problem and the rules is conducted in natural language using linguistic variables [6].
2. The model of the system
The method of fuzzy inference is based on the introduction of some rulesstatements that make the result dependent on the conditions [3, 10]. Conditions and results are described by linguistic variables. In the same way as in the algorithm lithological division at the diagnostic codes, for the variables of the system (indications of geophysical methods) is selected resolution form and entered the encoding of the values of linguistic variables in the form of numbers of terms of the set. As a variable, depending on the indications of geophysical methods, the lithological type acts. The rules necessary for fuzzy selection can be formulated using the table of diagnostic complex codes [1]. For example, the complex codes for the terrigenous section given in [1] are replaced by the rules in Table 1.
Table 1. Rules for the replacement of complex codes for the terrigenous section
${X}_{1}$

${X}_{2}$

${X}_{3}$

${X}_{4}$

${X}_{5}$

${X}_{6}$

$Y$

L

L

L

M

L

H

0

M

L

L

M

L

H

0

L

L

H

M

L

H

1

M

L

H

M

L

H

1

L

L

M

M

L

H

2

M

L

M

M

L

H

2

L

M

H

M

M

M

3

M

M

H

M

M

M

3

L

M

M

M

M

M

3

M

H

H

H

L

H

4

M

M

H

H

L

H

4

M

H

H

H

L

H

4

H

H

L

L

H

L

5

H

M

L

L

H

L

5

M

H

L

L

H

L

5

M

H

L

M

H

L

5

$Y\in \{0;1;2;3;4;5\}$ – lithological types (waterbearing Sandstone, oilbearing Sandstone, Sandstone with undefined saturation, clay Sandstone, dense rock, clay, respectively);
${X}_{i}$, $i=\overline{\mathrm{1,6}}$ – readings of geophysical methods caliper logging [4, 5], selfpotential, resistivity, microsonde, gamma radiation, neutron radiation, respectively. Linguistic variables have basic values: “low” (L), “medium” (M), “high” (H)

To rules (Table 1) it is necessary to add additional rules received on the basis of the decoded data on interpretation of materials of geophysical researches of wells (GRW).
3. Description of the algorithm
Let’s consider the algorithm of recognition of lithology of layers on the basis of fuzzy logic with the use of the known base of rules that provides more effective use of the available knowledge about some object.
The preliminary stage for this algorithm is the procedure of splitting the well into layers, which determines the composition of the layers on the basis of preliminary data on the characteristics of the wells of this field.
In the layers, which can be obtained by means of neural network separation [7, 11], are the average values of the parameters: $\mathbf{p}\left(h\right)=\underset{{h}_{BL}\left(h\right)}{\overset{{h}_{BR}\left(h\right)}{\int}}\mathbf{P}\left(\xi \right)d\xi /{h}_{BR}\left(h\right){h}_{BL}\left(h\right)$, where $h$ – depth; ${h}_{BL}\left(h\right)$, ${h}_{BR}\left(h\right)$ – left and right boundaries of the formation, lying at a depth of $h$, respectively; $\mathbf{P}\left(h\right)$ – the values of the parameters measured by any method.
In order to train the system for the subsequent forecasting of reservoir productivity, it is necessary to involve the results of processing and interpretation of several wells. For training and forecasting, consider the fuzzy inference method. We consider a system corresponding to a certain set of selected layers $G$. the Properties characterizing this class of layers are denoted by ${v}_{i}$, $i=\overline{1,m}$. the properties become variables of the system under consideration, which is a model of the class $G$. one of the properties is the lithological type. Other properties are the average integral values of the parameters $p\left(h\right)$. The selected properties are characterized by a set of linguistic variables ${x}_{i}$, $i=\overline{1,m}$. Linguistic variables can have the main values: “low” (L), “medium” (M), “high” (H), and, if necessary, some intermediate {L, LM, M, MN, H}. We write $n$ rules linking the state variables of the system, in the form of:
where ${Z}_{ij}$ is some value from the term set {L, LM, M, MN, H}.
The rule is read: variable ${x}_{k}$ takes the corresponding value from the set $\left\{L,...,H\right\}$, if other variables ${x}_{i}$, $i=1,...,m$, $i\ne k$ have some values from this set.
Rules (1) may be established by experts. If there is a data system corresponding to the objects $G$, the rules can be obtained by the classification learning algorithm.
We introduce the encoding of the values of linguistic variables in the form of the numbers of terms of the set ${x}_{i}=k$, $k=\overline{0,{K}_{i}}$, $i=\overline{1,m}$. Each term is characterized by its membership function ${\mu}_{ik}\left({v}_{i}\right)$. Let the objects ${g}_{j}\in G$, $j=\overline{1,n}$ correspond to the data system. The values of the system variables ${x}_{i}$, $i=\overline{1,m}$ form the objectproperty matrix $Q=\Vert {q}_{ij}\Vert $, $i=\overline{1,m}$, $j=\overline{1,n}$. The columns of the matrix $Q$ is formed by the values of the variables ${x}_{i}=k$, $k=\overline{0,{K}_{i}}$, $i=\overline{1,m}$. Thus, the elements of the matrix $Q$ are the values of ${Z}_{ij}$ from Eq. (1). For certainty, we believe that the lithological type corresponds to the first variable ${v}_{1}$. States $c\in C$ correspond to combinations of values of variables ${v}_{i}$, $i=\overline{1,m}$. Each state will be considered as a rule. If the sample is small, the expert can add a number of additional rules to the obtained rules, obtained on the basis of decoded data on the interpretation of GRW materials.
To apply the technique of fuzzy logic required is the assignment of membership functions for each term. If a data system is available, the optimal parameters of the fuzzy system can be obtained as a result of training. You can use a point estimation algorithm to do this. There is a data system in which the training sample is allocated. After selecting the resolution form, the rules are determined (1). Then, the intersection of fuzzy sets operation calculates the degree of execution of each rule. For the right part of the $j$th rule is the degree of ownership ${\mu}_{j}\left(Y\right)=\underset{i}{\mathrm{m}\mathrm{i}\mathrm{n}}\left[{\mu}_{i}\left({X}_{i}\right)\right]$, where ${X}_{i}$ – the base variable for the linguistic variable $Y$, included in the left part of the rule ($i=1,...,m$, $i\ne k$), and ${X}_{i}$ – the base variable for the linguistic variable ${x}_{k}$, included in the right part of the rule. The membership function corresponding to the variable $Y$ is determined by the ratio:
where ${\mu}_{k}\left(Y\right)={b}_{k}^{{L}_{k}}\left(Y\right)\text{,}$${b}_{k}^{l}\left(Y\right)={b}_{k}^{l1}\left(Y\right)+{\mu}_{k}^{l}\left(Y\right){b}_{k}^{l1}{\mu}_{k}^{l}\left(Y\right)\text{,}$$l=\overline{1,{L}_{k}}\text{;}$${b}_{k}^{0}=0$ using a probabilistic $Y$conorm, where ${L}_{k}$ is the number of rules with the same right part.
The function Eq. (2) calculates the point estimate of the variable $Y$:
The functions of the trapezoidal membership can be described using the coordinates of the corner points of the trapezoids for each term. We denote these parameters ${x}_{pi}^{t}$, $i=\overline{1,m}$, $t=\overline{1,{T}_{i}}$, $p=\overline{\mathrm{1,4}}$, where ${T}_{i}$ is the number of terms for the variable ${x}_{i}$. the form of the rule also determines the left and right boundaries of terms $B{L}_{i}^{t}$, $B{R}_{i}^{t}$, $i=\overline{1,m}$, $t=\overline{1,{T}_{i}}$, which generally do not coincide with ${x}_{1i}^{t}$, ${x}_{4i}^{t}$.
The result of the point estimation algorithm is the result ${Y}_{0}$, depending on the rules (1) and parameters ${x}_{pi}^{t}$, $B{L}_{i}^{t}$, $B{R}_{i}^{t}$, $i=\overline{1,m}$, $t=\overline{1,{T}_{i}}$, $p=\overline{\mathrm{1,4}}$. Therefore, we can assume that the function is defined: $Y=F(\mathbf{X},\mathbf{Z})$, where $\mathbf{X}=\left[{X}_{i}\right]$ is the vector of input variables; $\mathbf{Z}$ is a vector consisting of variables (parameters) ${x}_{pi}^{t}$, $B{L}_{i}^{t}$, $B{R}_{i}^{t}$, $i=\overline{1,m}$, $t=\overline{1,{T}_{i}}$, $p=\overline{\mathrm{1,4}}$.
For training, we define the target function:
where $\left\{\u27e8{\mathbf{X}}^{s},{Y}^{s}\u27e9\right\}$, $\left(s=\overline{1,h}\right)$ is the data system used for training (training sample); $h$ is the amount of training sample.
4. Results and discussions
The minimization of the objective function Eq. (3) gives the optimal values of the parameters of the fuzzy system $Z=\{{x}_{pi}^{t},B{L}_{i}^{t},B{R}_{i}^{t},i=\overline{1,m};t=\overline{1,{T}_{i}};p=\overline{\mathrm{1,4}}\}$. the Effective solution of the problem (4) is again carried out using a hybrid genetic algorithm BGAVM. The optimal membership functions obtained as a result of training on the rules (1) and corresponding to the considered geophysical methods are used further for the prediction and identification of lithological types of formations in new wells.
Fig. 1. The results of forecasting of productivity of the formation
a)
b)
In Fig. 1, a and b for two wells shows the forecast of reservoir productivity in comparison with the data decoding. The forecast corresponds to the dark bars, decoding – light. It is seen that for both wells, the productive formations corresponding to oil are recognized by 100 %. Water layers are recognized at 95 and 97 %, respectively. Layers of type oilwater are recognized worse: the system took 5 layers of oilwater to oil and 1 layer to the water.
5. Conclusions
A method has been developed to determine the informativeness of logging methods when recognizing the lithologic structure of a borehole section by the NL model, depending on the number of logging methods. It is important that the trained fuzzy system does not pass the productive formations of the oil type, and therefore the method can be used for rapid analysis of new wells.
References
 Dyakonova T. F. Primenenie EVM pri interpretatsii dannyih geofizicheskih issledovaniy skvazhin. Nedra, Moscow, 1991, p. 220, (in Russian). [Search CrossRef]
 Lyalin V. E. Intellektualnyie informatsionnoizmeritelnyie tehnologii i programmnoapparatnyie kompleksyi dlya avtomatizatsii geofizicheskih issledovaniy skvazhin. Monograph, Izdvo IE UrO RAN, Ekaterinburg, Izhevsk, 2001, (in Russian). [Search CrossRef]
 Casillas J., Cordon O., Herrera F. Learning fuzzy rulebased systems using ant colony optimization algorithms. 2nd International Workshop on Ant Algorithms, Brussels, Belgium, 2000. [Search CrossRef]
 Klir G. J., Yuan Bo Fuzzy Sets and Fuzzy Logic: Theory and Applications. Upper Saddle River, Prentice Hall, 1997. [Search CrossRef]
 Chang Hsien Cheng, Chen Hui Chuan, Fang Jen Ho Lithology determination from well logs with fuzzy associative memory neural network. IEEE Transactions on Geoscience and Remote Sensing, Vol. 35, 1997, p. 773780. [Publisher]
 Batycky R. P. A ThreeDimensional TwoPhase Field Scale Streamline Simulator. Ph.D. Thesis, Stanford University, 1997. [Search CrossRef]
 Saggaf M., Nebrija Ed L. A fuzzy logic approach for the estimation of facies from wireline logs. AAPG Bulletin, Vol. 87, 2003, p. 12231240. [Publisher]
 Hsieh Bieng Zih, Lewis Charles, Lin Zsay Shing Lithology identification of aquifers from geophysical well logs and fuzzy logic analysis: ShuiLin Area, Taiwan. Computers and Geosciences, Vol. 31, Issue 3, 2005, p. 263275. [Publisher]
 Saggaf M., Nebrija Ed L. Estimation of missing logs by regularized neural networks. AAPG Bulletin, Vol. 87, Issue 8, 2003, p. 13771389. [Publisher]
 Danilov V. A., Senilov M. A., Lyalin V. E. Interpretation of data of geophysical well logging with application of rank algorithm of an adaptive choice of subclasses. 5th International Congress of Mathematical Modeling, Dubna, 2002. [Search CrossRef]
 Bosch David, Ledo J., Queralt Pilar Fuzzy logic determination of lithologies from well log data: application to the KTB project data set (Germany). Surveys in Geophysics, Vol. 34, Issue 4, 2013, p. 413439. [Publisher]