# Optimization of Logging Interpretation Parameters Based on Committee Machine

Authors: Ziyi Na; Lixia Dang; Mingyang Chang; Rui Deng
Journal: Mathematical Problems in Engineering
 (2022)

---

## Abstract

Aiming at the problems of low precision and low applicability of the selection method of macroscopic capture cross-section parameters of rock skeleton, formation water, oil and gas, and mud in the volume model of remaining oil saturation logging, this study proposes a method for optimizing the capture section parameters based on the committee machine regression model. To begin, we select well sections from well logging data in a reasonable manner, allocate well sections with different component parameters in the same proportion, and construct sample datasets. The basic experts are then selected as multiple regression models, particle swarm optimization, and robust regression methods to train and learn the input parameters. By combining multiple experts, the regression committee machine improves the overall performance of the intelligent model. Finally, the genetic algorithm is used as a combiner to determine the contribution of each basic expert network in the final output, and the optimized parameters are obtained, which are then fed into the volume model to calculate the remaining oil saturation of the newly developed production wells, guiding the perforation and development work. The model is used to evaluate the remaining oil in theX oilfield, and the calculated water saturation matches the oil test results, proving the model’s accuracy and availability. The use of real-world data demonstrates that this method can effectively characterize the four parameter values in the volume model and provide reliable geophysical technical support for the evaluation of remaining oil.


---

## Body

## 1. Introduction

Most of China’s oilfields have been developed by water injection and have gradually entered the later stage of development, the evaluation of residual oil saturation is becoming more and more difficult. Residual oil logging evaluation is currently the main means of macroscopic evaluation of residual oil. In the process of oil and gas reservoir development, the geophysical characteristics will change accordingly, and the study of changes in logging response during the development process can guide the residual oil tapping. The methods of identifying and evaluating the water flooded layer can be mainly divided into three types: one is the laboratory geological analysis technology; the second is the reservoir geochemical technology; the third is the geophysical technology [1], among which the neutron lifetime logging is an important means of identifying and evaluating water flooded formations in geophysical logging. Due to the high sensitivity of neutron lifetime logging technology in water flooded layers with high salinity, this technology is usually used in oil fields to qualitatively identify and quantitatively evaluate water flooded layers, mainly by measuring the macro capture section of the formation and calculating the remaining oil saturation of the reservoir [2].Based on the petrophysical volume model, the relationship between the capture section and water saturation can be established [3]. The rock volume model refers to dividing the rock volume into different parts according to the composition of the reservoir and its physical properties, and regards the logging results as the sum of the contributions of each part. The formation porosity and shale content required to obtain water saturation by using the volume model can be obtained from the open hole logging data [4], and the reliability is high. The porosity and shale content of the formation required for water saturation can be obtained with high reliability from open hole logging data; however, the selection of macro capture cross sections, such as rock skeleton, formation water, oil and gas, and mud, is based on empirical values, which do not accurately reflect formation characteristics and affect residual oil saturation calculation accuracy.Artificial neural networks and geophysical logging have a growing number of interdisciplinary applications as big data and machine learning technology become more advanced. Predecessors have conducted a great deal of exploratory research and made some progress on the topic of calculating reservoir physical parameters. Liu et al. used the hilchie index when calculating the formation mud content [5]. Kang et al. used the principle of normalized nonlinear least square method to optimize the parameters of regional formation components [6]. An et al. established LSTM cyclic neural network method to predict mud content and porosity [7]. Zhai and Dong used the density regression method to calculate the density value of the skeleton [8].The majority of the currently used machine learning applications for parameter optimization of volume models only take into account one model, and frequently only one of the four parameters is used as the research subject without explaining the entire set. This study applies the committee machine model in light of this issue. By combining multiple experts, the regression committee machine enhances the performance of the intelligent model overall. It selects the multiple regression model, particle swarm optimization algorithm, and robust regression model as the expert network to optimize the four parameters of rock skeleton, formation water, oil and gas, and mud macrocapture cross section in the volume model, and determines the weight of each expert through genetic algorithm. Finally, the optimized four capture cross-section values are output and brought into the water saturation calculation formula to calculate the remaining oil saturation of the newly developed production wells, so as to improve the accuracy of interpretation of water flooded zones in this area.


## 2. The Principle of Parameter Optimization

According to the petrophysical volume model, the reservoir can be regarded as a simple structure composed of mud, framework, and pores. The framework often includes different lithologic components, and the pores contain oil, gas, water, and other fluids [9], as shown in Figure 1.Figure 1
Schematic diagram of the volume model.Taking the capture cross-section value measured by neutron logging instrument as the comprehensive response of macrocapture cross sections of skeleton, mud, water, and hydrocarbon, the calculation formula of water saturation calculated by volume model method is as follows:(1)Sw=Σlog−Σma−ϕΣh−ΣmaϕΣw−Σh−VshΣsh−ΣmaϕΣw−Σh,where Σlog is the macroscopic capture section of the formation, Σw is the macroscopic capture section of water, Σma is the macroscopic capture section of all dry solids (framework, silt, and dry clay colloid), Σsh is the macroscopic capture section of hydrocarbons, Σh is the macroscopic capture section of the formation mud, Vsh is the shale content of the formation, ϕ is the total porosity of the formation (including shale irreducible water), and Sw is the water saturation.Apart from the precise macrocapture cross section, the formula requires the determination of six parameters, includingVsh, Σma, Σsh, Σh, Σw, and ϕ. Formation porosity and shale content can be obtained from open hole well logging data with high reliability. In previous studies, the combination of empirical values and theoretical values is usually used to determine the four parameters: Σma, Σsh, Σh, and Σw, but the thermal neutron capture of each component in the formation is not constant [10], and it varies depending on salinity and hydrogen level in the formation. The parameters chosen based on empirical values are entered into the model, and the resulting water saturation does not correctly reflect the formation’s genuine state [11]. In this instance, the parameters of each component skeleton can be inversed according to the actual measured values utilizing the parameter optimization approach, allowing the selected parameters to reflect the formation’s actual conditions.The optimization method of the capture cross-section parameters of rock components is to take the parameters ofΣma, Σsh, Σh, and Σw as unknown quantities, and use logging data to obtain porosity, shale content, and the total capture cross-section values in the formation as known quantities to construct a linear system of equations::(2)∑z=Vsh∑sh+Vma∑ma+Vw∑w+Vo∑o,s·t·Vsh+Vma+Vw+Vo=1,where ∑z is the total trapping section of the formation obtained by measurement, ∑sh, ∑ma, ∑w, and ∑o are the parameters to be determined, and the constraint condition is that the sum of the volume content of each component is 1.


## 3. Committee Machine Model

Nilsson proposed the concept of committee machine in 1965, which is actually a two-layer network, but only needs to adjust the weights of one-layer network during training. The committee machine is usually divided into a neural network ensemble, which is made up of several expert machine with independent prediction or classification functions. The local results obtained by each expert machine are combined and decision-making is made, and then the global decision-making model of the entire committee machine system is obtained through some synthetic mechanism.Different networks in the committee machine are different experts. The committee machine requires the experts to have necessary differences. The performance of the committee machine can be improved by integrating different experts through the committee machine. Expert differences can be obtained by varying the initial training settings and training sets of experts, or by modifying the network differences [12].Haykin divides committee machine into two categories [13], one is static structure and the other is dynamic structure. The integrated decision-making method of committee machine output to each expert does not involve the committee machine’s input in the static structure; instead, the committee machine’s input is used to determine how to combine the partial solutions of each expert’s output. The committee machine model presented in this research is classified as a static structure, and the model flow chart is shown in Figure 2. The multiple linear regression model, the particle swarm optimization algorithm, and the robust regression model are the three different networks selected in the field of experts. The three methods of network training are different, which meet the different conditions of the experts, and finally output the results through the combiner.Figure 2
Schematic diagram of the committee machine model.


### 3.1. Expert Selection

#### 3.1.1. Linear Regression Model

Linear regression is one of the most well-known modeling techniques and one of the first-choice techniques people use when learning predictive models. In this technique, the dependent variable is continuous, the independent variable can be either continuous or discrete [14], and the regression line is linear in nature.The condition that the multiple linear regression model needs to meet is that there is a linear relationship between the explanatory variableY and multiple explanatory variables x1,x2,…,xk [15], which is a multiple linear function of the explanatory variables, namely,(3)Y=β0+β1X1+β2X2+⋯+βkXk+μ,where β0,β1⋯βk are k + 1 unknown parameters, β0 is known as regression constants, in which β0,β1⋯βk are regression coefficients, y is called as explained variables, and x1,x2,…,xk are k explanatory variables, in which μ known as random errors.Corresponding to the model in parametric inversion,(4)∑z=Vsh∑sh+Vma∑ma+Vw∑w+Vh∑h.In this case,k = 4, and the constant term β0 is 0, β1 for Vsh, β2 for Vma, β3 for Vw, β4 for Vh, X1 for ∑sh, X2 for ∑ma, X3 for ∑w, X4 for ∑h, and Y for ∑Z, to satisfy the conditions for a multiple regression model.The use of multiple linear regression models requires a goodness of fit test [16]; that is, for formula (3), the regression equation is obtained as follows:(5)Y^=β^0+β^1X1+β^2X2+⋯+β^KXk.For the evaluation of the fitting effect of the regression model, it is usually selectedR¯2 as the test index [17], and its value is between 0 and 1. The larger the value, the better the fitting effect of the equation. The calculation formula is as follows:(6)R¯2=1−SSE/n−p−1SST/n−1=1−∑Yi−Y^2/n−p−1∑Yi−Y¯2/n−1,where n − p − 1 and N − 1 are the degrees of freedom of SSE and SST.


#### 3.1.2. Particle Swarm Optimization Algorithm Model

The particle swarm optimization algorithm model is derived from simulations of bird migration and aggregation in the foraging process [18], and the birds achieve the optimal goal through a cooperative group, by following the optimal value of the current global optimization search, with fewer parameters and faster convergence [19].The mathematical description of the particle swarm optimization algorithm model is that, for a population, given an N-dimensional search space, a total number ofM particles are generated in the process of its movement, and the population is set as X=x1,x2,⋯,xm, the corresponding position vector of each particle is partxi=xi1,xi2,⋯,xiN, velocity vector partvk=xi1,xi2,⋯,xim, i=1,2,⋯M, position vector represents the potential solution of the objective function, and the velocity vector can control the size of the step size when searching for the solution. The optimal position of each particle in the pursuit of its own history is set as poi=pk1,pk2,⋯,pkm; the optimal solution found by the pursuit group is set as goi=g1,g1,⋯,gn. After n + 1 iterations, the change of the Kth dimension velocity and position can be expressed by the following formula:(7)viki+1=W×Viki+c1×rand1×pbiki−Piki+c2×rand2×gbki−Piki,Piki+1=Piki+Viki+1,where i represents the ith particle, k represents the dimension of the current solution, t is the number of iterations, c1 and c2 are the learning factors, rand1 and rand2 represent random numbers, and W represents the weight.Given different learning factors, the particles can have different learning abilities, and the learning factors can follow a certain change rule, which plays an important role in regulating the learning method of the particles [20]. For example, making the parameters decrease linearly from the maximum value indicates that the particles are less dependent on the population in the later stage of the iteration. By making the two learning factors change asynchronously, particles can take different optimization forms in the early and late optimization process [21].The weight plays an important role in particle solution finding [22]. By changing the size of the weight, you can decide whether to let it solve in the original direction. The larger the weight value, the faster the particle moves in the original direction, and the smaller the correction amount of the speed, so that it has a smaller convergence speed. When the weight value decreases, it means to increase the amount of correction to the particle direction, so that it can find a solution in the new direction.


#### 3.1.3. Robust Regression

At present, most of the regression models for regression analysis of logging data are implemented by the least squares method. However, during the logging process, due to the influence of the logging instrument and accuracy, there will be abnormal points in the logging data, and these abnormal points will affect the regression model. There is a great influence, and the regression model obtained using the least squares method is quite different from the actual situation.Robust regression is a method used for regression when the least squares method has abnormal points in the data samples [23], which can detect abnormal points or sample values that have the greatest impact on the model. The most commonly used in robust regression is the maximum likelihood M estimation [24], and the general model of M estimation is as follows:(8)Yi=∑j=1pxijβj+ei,i=1,⋯,n,where β1,⋯,βp is the regression coefficient, e1,⋯,ei is independent and identically distributed, and the variance is 0. The purpose of the least square method is to obtain the value of the regression coefficient that minimizes the value of the objective function. The idea of using the M estimation of robust regression is to iteratively weight the coefficients of the least squares method [25], and then combine the value of the regression residuals to determine the weight of each point, so as to make the function robust. The optimized objective function is as follows:(9)∑j=1nwiYi−∑j=1pxijβj2=min.M estimation reduces the effect of outliers by controlling the weight of outliers. By giving smaller weights to points with large residual values and larger weights to points with small residuals, the purpose of screening data is achieved and the minimum value is established. The optimal solution is obtained by quadratic estimation, and the weight coefficients are improved by iterative iterations. There are many ways to construct weights. In this experiment, the Biweight method [26] provided by DPS is adopted:(10)wi=0,uic>1,1−uicc22,uic≤1,where ui=ei/s is the residual index and s is the residual scale.


### 3.2. Combiners

The simple average method or the weighted average method is often chosen as the combination strategy in the regression committee model. Based on the weighted average method, the final output of the committee machine is obtained by weighted sum of the output of each expert. The weight can be selected by experience method, artificial method, or by optimization algorithm. In this paper, genetic algorithm is used as a weight calculator to search the global optimal solution randomly.Under the condition of producing a certain number of individuals first [27], the genetic algorithm produces more individuals through breeding, crossing and variation, and screening these individuals according to certain conditions, and in the screening process, keep the number of individuals at a constant level and repeat the selection process to find the optimal solution [28]. The specific operation steps encode the actual problem as a chromosome, randomly generate a certain number of individuals, that is, a population, calculate the fitness value of each individual, and judge whether the initial solution is the optimal solution through a given termination condition, and then choose whether to output. If the result is not its optimal solution, a new population is generated through the genetic operator, and the above operations are repeated. Real-coded genetic algorithms (GAs) are used by Roy et al. to train ANNs, and they produce effective results. [29].The optimization objective function designed by genetic algorithm is as follows:(11)E=∑j=1N−In∑i=13wiexp−12Yhatj−Y0j2.The restrictions are as follows:(12)∑i=13wi=1,0≤wi≤1,where w1, w2, and w3 are the weights of multiple linear regression, robust regression, and particle swarm optimization, respectively, Yhat is the prediction value, Y0 is the target value, and N is the quantity of training data.


## 3.1. Expert Selection

### 3.1.1. Linear Regression Model

Linear regression is one of the most well-known modeling techniques and one of the first-choice techniques people use when learning predictive models. In this technique, the dependent variable is continuous, the independent variable can be either continuous or discrete [14], and the regression line is linear in nature.The condition that the multiple linear regression model needs to meet is that there is a linear relationship between the explanatory variableY and multiple explanatory variables x1,x2,…,xk [15], which is a multiple linear function of the explanatory variables, namely,(3)Y=β0+β1X1+β2X2+⋯+βkXk+μ,where β0,β1⋯βk are k + 1 unknown parameters, β0 is known as regression constants, in which β0,β1⋯βk are regression coefficients, y is called as explained variables, and x1,x2,…,xk are k explanatory variables, in which μ known as random errors.Corresponding to the model in parametric inversion,(4)∑z=Vsh∑sh+Vma∑ma+Vw∑w+Vh∑h.In this case,k = 4, and the constant term β0 is 0, β1 for Vsh, β2 for Vma, β3 for Vw, β4 for Vh, X1 for ∑sh, X2 for ∑ma, X3 for ∑w, X4 for ∑h, and Y for ∑Z, to satisfy the conditions for a multiple regression model.The use of multiple linear regression models requires a goodness of fit test [16]; that is, for formula (3), the regression equation is obtained as follows:(5)Y^=β^0+β^1X1+β^2X2+⋯+β^KXk.For the evaluation of the fitting effect of the regression model, it is usually selectedR¯2 as the test index [17], and its value is between 0 and 1. The larger the value, the better the fitting effect of the equation. The calculation formula is as follows:(6)R¯2=1−SSE/n−p−1SST/n−1=1−∑Yi−Y^2/n−p−1∑Yi−Y¯2/n−1,where n − p − 1 and N − 1 are the degrees of freedom of SSE and SST.


### 3.1.2. Particle Swarm Optimization Algorithm Model

The particle swarm optimization algorithm model is derived from simulations of bird migration and aggregation in the foraging process [18], and the birds achieve the optimal goal through a cooperative group, by following the optimal value of the current global optimization search, with fewer parameters and faster convergence [19].The mathematical description of the particle swarm optimization algorithm model is that, for a population, given an N-dimensional search space, a total number ofM particles are generated in the process of its movement, and the population is set as X=x1,x2,⋯,xm, the corresponding position vector of each particle is partxi=xi1,xi2,⋯,xiN, velocity vector partvk=xi1,xi2,⋯,xim, i=1,2,⋯M, position vector represents the potential solution of the objective function, and the velocity vector can control the size of the step size when searching for the solution. The optimal position of each particle in the pursuit of its own history is set as poi=pk1,pk2,⋯,pkm; the optimal solution found by the pursuit group is set as goi=g1,g1,⋯,gn. After n + 1 iterations, the change of the Kth dimension velocity and position can be expressed by the following formula:(7)viki+1=W×Viki+c1×rand1×pbiki−Piki+c2×rand2×gbki−Piki,Piki+1=Piki+Viki+1,where i represents the ith particle, k represents the dimension of the current solution, t is the number of iterations, c1 and c2 are the learning factors, rand1 and rand2 represent random numbers, and W represents the weight.Given different learning factors, the particles can have different learning abilities, and the learning factors can follow a certain change rule, which plays an important role in regulating the learning method of the particles [20]. For example, making the parameters decrease linearly from the maximum value indicates that the particles are less dependent on the population in the later stage of the iteration. By making the two learning factors change asynchronously, particles can take different optimization forms in the early and late optimization process [21].The weight plays an important role in particle solution finding [22]. By changing the size of the weight, you can decide whether to let it solve in the original direction. The larger the weight value, the faster the particle moves in the original direction, and the smaller the correction amount of the speed, so that it has a smaller convergence speed. When the weight value decreases, it means to increase the amount of correction to the particle direction, so that it can find a solution in the new direction.


### 3.1.3. Robust Regression

At present, most of the regression models for regression analysis of logging data are implemented by the least squares method. However, during the logging process, due to the influence of the logging instrument and accuracy, there will be abnormal points in the logging data, and these abnormal points will affect the regression model. There is a great influence, and the regression model obtained using the least squares method is quite different from the actual situation.Robust regression is a method used for regression when the least squares method has abnormal points in the data samples [23], which can detect abnormal points or sample values that have the greatest impact on the model. The most commonly used in robust regression is the maximum likelihood M estimation [24], and the general model of M estimation is as follows:(8)Yi=∑j=1pxijβj+ei,i=1,⋯,n,where β1,⋯,βp is the regression coefficient, e1,⋯,ei is independent and identically distributed, and the variance is 0. The purpose of the least square method is to obtain the value of the regression coefficient that minimizes the value of the objective function. The idea of using the M estimation of robust regression is to iteratively weight the coefficients of the least squares method [25], and then combine the value of the regression residuals to determine the weight of each point, so as to make the function robust. The optimized objective function is as follows:(9)∑j=1nwiYi−∑j=1pxijβj2=min.M estimation reduces the effect of outliers by controlling the weight of outliers. By giving smaller weights to points with large residual values and larger weights to points with small residuals, the purpose of screening data is achieved and the minimum value is established. The optimal solution is obtained by quadratic estimation, and the weight coefficients are improved by iterative iterations. There are many ways to construct weights. In this experiment, the Biweight method [26] provided by DPS is adopted:(10)wi=0,uic>1,1−uicc22,uic≤1,where ui=ei/s is the residual index and s is the residual scale.


## 3.1.1. Linear Regression Model

Linear regression is one of the most well-known modeling techniques and one of the first-choice techniques people use when learning predictive models. In this technique, the dependent variable is continuous, the independent variable can be either continuous or discrete [14], and the regression line is linear in nature.The condition that the multiple linear regression model needs to meet is that there is a linear relationship between the explanatory variableY and multiple explanatory variables x1,x2,…,xk [15], which is a multiple linear function of the explanatory variables, namely,(3)Y=β0+β1X1+β2X2+⋯+βkXk+μ,where β0,β1⋯βk are k + 1 unknown parameters, β0 is known as regression constants, in which β0,β1⋯βk are regression coefficients, y is called as explained variables, and x1,x2,…,xk are k explanatory variables, in which μ known as random errors.Corresponding to the model in parametric inversion,(4)∑z=Vsh∑sh+Vma∑ma+Vw∑w+Vh∑h.In this case,k = 4, and the constant term β0 is 0, β1 for Vsh, β2 for Vma, β3 for Vw, β4 for Vh, X1 for ∑sh, X2 for ∑ma, X3 for ∑w, X4 for ∑h, and Y for ∑Z, to satisfy the conditions for a multiple regression model.The use of multiple linear regression models requires a goodness of fit test [16]; that is, for formula (3), the regression equation is obtained as follows:(5)Y^=β^0+β^1X1+β^2X2+⋯+β^KXk.For the evaluation of the fitting effect of the regression model, it is usually selectedR¯2 as the test index [17], and its value is between 0 and 1. The larger the value, the better the fitting effect of the equation. The calculation formula is as follows:(6)R¯2=1−SSE/n−p−1SST/n−1=1−∑Yi−Y^2/n−p−1∑Yi−Y¯2/n−1,where n − p − 1 and N − 1 are the degrees of freedom of SSE and SST.


## 3.1.2. Particle Swarm Optimization Algorithm Model

The particle swarm optimization algorithm model is derived from simulations of bird migration and aggregation in the foraging process [18], and the birds achieve the optimal goal through a cooperative group, by following the optimal value of the current global optimization search, with fewer parameters and faster convergence [19].The mathematical description of the particle swarm optimization algorithm model is that, for a population, given an N-dimensional search space, a total number ofM particles are generated in the process of its movement, and the population is set as X=x1,x2,⋯,xm, the corresponding position vector of each particle is partxi=xi1,xi2,⋯,xiN, velocity vector partvk=xi1,xi2,⋯,xim, i=1,2,⋯M, position vector represents the potential solution of the objective function, and the velocity vector can control the size of the step size when searching for the solution. The optimal position of each particle in the pursuit of its own history is set as poi=pk1,pk2,⋯,pkm; the optimal solution found by the pursuit group is set as goi=g1,g1,⋯,gn. After n + 1 iterations, the change of the Kth dimension velocity and position can be expressed by the following formula:(7)viki+1=W×Viki+c1×rand1×pbiki−Piki+c2×rand2×gbki−Piki,Piki+1=Piki+Viki+1,where i represents the ith particle, k represents the dimension of the current solution, t is the number of iterations, c1 and c2 are the learning factors, rand1 and rand2 represent random numbers, and W represents the weight.Given different learning factors, the particles can have different learning abilities, and the learning factors can follow a certain change rule, which plays an important role in regulating the learning method of the particles [20]. For example, making the parameters decrease linearly from the maximum value indicates that the particles are less dependent on the population in the later stage of the iteration. By making the two learning factors change asynchronously, particles can take different optimization forms in the early and late optimization process [21].The weight plays an important role in particle solution finding [22]. By changing the size of the weight, you can decide whether to let it solve in the original direction. The larger the weight value, the faster the particle moves in the original direction, and the smaller the correction amount of the speed, so that it has a smaller convergence speed. When the weight value decreases, it means to increase the amount of correction to the particle direction, so that it can find a solution in the new direction.


## 3.1.3. Robust Regression

At present, most of the regression models for regression analysis of logging data are implemented by the least squares method. However, during the logging process, due to the influence of the logging instrument and accuracy, there will be abnormal points in the logging data, and these abnormal points will affect the regression model. There is a great influence, and the regression model obtained using the least squares method is quite different from the actual situation.Robust regression is a method used for regression when the least squares method has abnormal points in the data samples [23], which can detect abnormal points or sample values that have the greatest impact on the model. The most commonly used in robust regression is the maximum likelihood M estimation [24], and the general model of M estimation is as follows:(8)Yi=∑j=1pxijβj+ei,i=1,⋯,n,where β1,⋯,βp is the regression coefficient, e1,⋯,ei is independent and identically distributed, and the variance is 0. The purpose of the least square method is to obtain the value of the regression coefficient that minimizes the value of the objective function. The idea of using the M estimation of robust regression is to iteratively weight the coefficients of the least squares method [25], and then combine the value of the regression residuals to determine the weight of each point, so as to make the function robust. The optimized objective function is as follows:(9)∑j=1nwiYi−∑j=1pxijβj2=min.M estimation reduces the effect of outliers by controlling the weight of outliers. By giving smaller weights to points with large residual values and larger weights to points with small residuals, the purpose of screening data is achieved and the minimum value is established. The optimal solution is obtained by quadratic estimation, and the weight coefficients are improved by iterative iterations. There are many ways to construct weights. In this experiment, the Biweight method [26] provided by DPS is adopted:(10)wi=0,uic>1,1−uicc22,uic≤1,where ui=ei/s is the residual index and s is the residual scale.


## 3.2. Combiners

The simple average method or the weighted average method is often chosen as the combination strategy in the regression committee model. Based on the weighted average method, the final output of the committee machine is obtained by weighted sum of the output of each expert. The weight can be selected by experience method, artificial method, or by optimization algorithm. In this paper, genetic algorithm is used as a weight calculator to search the global optimal solution randomly.Under the condition of producing a certain number of individuals first [27], the genetic algorithm produces more individuals through breeding, crossing and variation, and screening these individuals according to certain conditions, and in the screening process, keep the number of individuals at a constant level and repeat the selection process to find the optimal solution [28]. The specific operation steps encode the actual problem as a chromosome, randomly generate a certain number of individuals, that is, a population, calculate the fitness value of each individual, and judge whether the initial solution is the optimal solution through a given termination condition, and then choose whether to output. If the result is not its optimal solution, a new population is generated through the genetic operator, and the above operations are repeated. Real-coded genetic algorithms (GAs) are used by Roy et al. to train ANNs, and they produce effective results. [29].The optimization objective function designed by genetic algorithm is as follows:(11)E=∑j=1N−In∑i=13wiexp−12Yhatj−Y0j2.The restrictions are as follows:(12)∑i=13wi=1,0≤wi≤1,where w1, w2, and w3 are the weights of multiple linear regression, robust regression, and particle swarm optimization, respectively, Yhat is the prediction value, Y0 is the target value, and N is the quantity of training data.


## 4. Interpretation of Residual Oil Saturation of the Committee Machine Model

X oilfield has gradually entered the aquatic production period of medium and high water content. The mineralization of distinct formations in the vertical direction is complex and variable due to the uneven advancement and alternating mixed injection of clean water and sewage induced by plane water injection. The saturation logging interpretation settings are chosen using a traditional statistical method that does not adjust to the present scenario. Early open hole interpretation and logging can be used to classify wells with comparable lithology and water injection distribution as one type of interpretation parameter, avoiding unpredictability in interpretation. In study area A of X oilfield, we select the wells with the same lithology and formation, as well as the horizon of interpretation results and perforation verification. Finally, 998 logging data groups from 30 horizons were obtained. The committee machine model’s training target horizon types include oil layer, oil-water layer, water layer, and dry layer. The interpretation parameters with popularization and application value in the block are obtained through the unified selection of the component parameters of different reservoir types. The committee machine is used to select the interpretation parameters of block A of the X oilfield, output the parameter values with high accuracy, and then feed them into the volume model to calculate the remaining oil saturation of the newly developed production wells, thereby guiding perforation and development work.


### 4.1. Data Processing

First, we perform correlation analysis onVsh, Vma, Vw, and Vh to avoid the problem of multicollinearity and calculate the correlation coefficient between the four variables. The correlation coefficient is a statistical indicator of the closeness of the relationship between the response variables. The value range of the correlation coefficient is between −1 and 1. 1 means that the two variables are perfectly linearly correlated, -1 means that the two variables are completely negatively correlated, and 0 means that the two variables are not correlated. The closer the data is to 0, the weaker the correlation is.The calculation results are shown in Figure3. The absolute value of the correlation coefficient between the four variables does not exceed 0.7. It can be considered that the correlation between the variables is not strong and multicollinearity processing is not required.Figure 3
Correlation coefficient analysis diagram.


### 4.2. Parameter Optimization

Based on the above regression, the committee machine built the parameter optimization model of the rock component capture section in well XII in theX oilfield, input the test set’s logging data into the committee machine model, first output three groups of parameter values through the calculation of three experts, and then optimize the combination weight of the three groups of parameter values via the genetic algorithm combiner to obtain the committee machine out. Finally, the results are entered into the verification set to validate the model’s effect.The performance evaluation indicators of the committee’s machine model are shown in Table1. The average absolute error (MAE) and the adjusted goodness of fit (R¯2) are selected as the two indicators of model performance evaluation. The MAE is a measure of the error between pairs of observations expressing the same phenomenon. Outliers have little influence when using the MAE loss function, and the fitting straight line can better represent the distribution of normal data. The adjusted goodness of fit is used to assess the overall fitness of the regression equation and to express the overall relationship between the dependent and independent variables. In the training set, the mean absolute error (MAE) between the formation capture cross-section value calculated by the parameters output by the multiple linear regression expert and the actual formation capture cross-section value is 1.0171, the error is the largest, the MAE of the committee machine is 1.0100, and the error is the smallest. The MAE of the robust regression expert and the particle swarm optimization expert are 1.0118 and 1.0161, respectively. We verify that the MAE value of the centralized committee machine is 1.0529, with the minimum error. The adjusted goodness of fit is the highest value of the committee machine in both the training and validation sets, indicating that it performs better than a single expert network.Table 1
Committee machine model performance evaluation index.
Training setMAER¯2Validation setMAER¯2Multiple regression1.01710.6820Multiple regression1.05980.6328Robust regression1.01180.6632Robust regression1.05400.6262Particle swarm optimization1.01610.6123Particle swarm optimization1.05550.5119Committee machine1.01000.6825Committee machine1.05290.6377It can be seen from Figure4 that the data point distribution of the committee’s machine model and the actual data is the most uniform and the convergence degree is the best with respect to the 45 degree line distribution.Figure 4
Correlation coefficient analysis diagram.Using the committee machine model constructed above, the parameters of each component capture section of wellX in oilfield are optimized. The parameter values of ∑sh, ∑ma, ∑w, and ∑o obtained by the committee machine are 33.1747c.u, 11.4581c.u, 50.3596c.u, and 10.3007c.u, respectively. As shown in Table 2, the capture cross-section values calculated by the committee machine model itself and the three expert systems included in it are within the theoretical range.Table 2
Committee machine model output parameter results.
∑shc.u.∑mac.u.∑wc.u.∑oc.u.Particle swarm optimization10.932630.873352.6328.9641Multiple regression11.492934.109748.80588.9307Robust regression11.426232.684051.024210.9495Committee machine11.458133.174750.359610.3007


### 4.3. Residual Oil Saturation Interpretation

The four capture section parameters output by the committee machine are applied to equation (1) to obtain the water saturation value after parameter optimization, and the following interpretation results of well XII can be obtained. In Figure 5, SWPNN is the water saturation curve calculated by using the empirical value of the previous application area, and SWCAL is the water saturation curve calculated by using the interpretation parameters output by the above committee machine. It is obvious that the water saturation curve calculated after the optimized capture section parameters is obviously smooth and consistent with the shape corresponding to other curves. The correlation shows that the accuracy of this interpretation has been improved. It can be seen from the interpretation results that the remaining oil saturation of layers A1 and A2 is about 60%, which has development and production potential. The remaining saturation of layers A3 and A4 is lower than 40%, and the porosity is smaller than that of layers A1 and A2. The degree of water flooded layer is high. It is basically an oil-bearing water layer with low exploitation potential.Figure 5
Interpretation results of well XII.After the parameter optimization interpretation, the study area conducted perforation analysis on the potential A1 and A2 layers of the well according to the interpretation results and verified the interpretation profile of the oil production rate of the interpreted layer. The interpretation conclusions are shown in Table3.Table 3
Calculation of oil production rate.
HorizonSWPNN (%)SWCAL (%)Perforation water production rate (%)Interpretation conclusionA161.642.848.5Low water flooded layerA258.741.246.9Low water flooded layerThe results show that the accuracy of the committee machine output parameters in the calculation of residual oil saturation has improved, and in terms of water production rate, the value calculated after parameter optimization is significantly closer to the value of water production rate, indicating that the regression method of the capture section parameter value in the block through the committee machine model is reliable and can be popularized and applied.


## 4.1. Data Processing

First, we perform correlation analysis onVsh, Vma, Vw, and Vh to avoid the problem of multicollinearity and calculate the correlation coefficient between the four variables. The correlation coefficient is a statistical indicator of the closeness of the relationship between the response variables. The value range of the correlation coefficient is between −1 and 1. 1 means that the two variables are perfectly linearly correlated, -1 means that the two variables are completely negatively correlated, and 0 means that the two variables are not correlated. The closer the data is to 0, the weaker the correlation is.The calculation results are shown in Figure3. The absolute value of the correlation coefficient between the four variables does not exceed 0.7. It can be considered that the correlation between the variables is not strong and multicollinearity processing is not required.Figure 3
Correlation coefficient analysis diagram.


## 4.2. Parameter Optimization

Based on the above regression, the committee machine built the parameter optimization model of the rock component capture section in well XII in theX oilfield, input the test set’s logging data into the committee machine model, first output three groups of parameter values through the calculation of three experts, and then optimize the combination weight of the three groups of parameter values via the genetic algorithm combiner to obtain the committee machine out. Finally, the results are entered into the verification set to validate the model’s effect.The performance evaluation indicators of the committee’s machine model are shown in Table1. The average absolute error (MAE) and the adjusted goodness of fit (R¯2) are selected as the two indicators of model performance evaluation. The MAE is a measure of the error between pairs of observations expressing the same phenomenon. Outliers have little influence when using the MAE loss function, and the fitting straight line can better represent the distribution of normal data. The adjusted goodness of fit is used to assess the overall fitness of the regression equation and to express the overall relationship between the dependent and independent variables. In the training set, the mean absolute error (MAE) between the formation capture cross-section value calculated by the parameters output by the multiple linear regression expert and the actual formation capture cross-section value is 1.0171, the error is the largest, the MAE of the committee machine is 1.0100, and the error is the smallest. The MAE of the robust regression expert and the particle swarm optimization expert are 1.0118 and 1.0161, respectively. We verify that the MAE value of the centralized committee machine is 1.0529, with the minimum error. The adjusted goodness of fit is the highest value of the committee machine in both the training and validation sets, indicating that it performs better than a single expert network.Table 1
Committee machine model performance evaluation index.
Training setMAER¯2Validation setMAER¯2Multiple regression1.01710.6820Multiple regression1.05980.6328Robust regression1.01180.6632Robust regression1.05400.6262Particle swarm optimization1.01610.6123Particle swarm optimization1.05550.5119Committee machine1.01000.6825Committee machine1.05290.6377It can be seen from Figure4 that the data point distribution of the committee’s machine model and the actual data is the most uniform and the convergence degree is the best with respect to the 45 degree line distribution.Figure 4
Correlation coefficient analysis diagram.Using the committee machine model constructed above, the parameters of each component capture section of wellX in oilfield are optimized. The parameter values of ∑sh, ∑ma, ∑w, and ∑o obtained by the committee machine are 33.1747c.u, 11.4581c.u, 50.3596c.u, and 10.3007c.u, respectively. As shown in Table 2, the capture cross-section values calculated by the committee machine model itself and the three expert systems included in it are within the theoretical range.Table 2
Committee machine model output parameter results.
∑shc.u.∑mac.u.∑wc.u.∑oc.u.Particle swarm optimization10.932630.873352.6328.9641Multiple regression11.492934.109748.80588.9307Robust regression11.426232.684051.024210.9495Committee machine11.458133.174750.359610.3007


## 4.3. Residual Oil Saturation Interpretation

The four capture section parameters output by the committee machine are applied to equation (1) to obtain the water saturation value after parameter optimization, and the following interpretation results of well XII can be obtained. In Figure 5, SWPNN is the water saturation curve calculated by using the empirical value of the previous application area, and SWCAL is the water saturation curve calculated by using the interpretation parameters output by the above committee machine. It is obvious that the water saturation curve calculated after the optimized capture section parameters is obviously smooth and consistent with the shape corresponding to other curves. The correlation shows that the accuracy of this interpretation has been improved. It can be seen from the interpretation results that the remaining oil saturation of layers A1 and A2 is about 60%, which has development and production potential. The remaining saturation of layers A3 and A4 is lower than 40%, and the porosity is smaller than that of layers A1 and A2. The degree of water flooded layer is high. It is basically an oil-bearing water layer with low exploitation potential.Figure 5
Interpretation results of well XII.After the parameter optimization interpretation, the study area conducted perforation analysis on the potential A1 and A2 layers of the well according to the interpretation results and verified the interpretation profile of the oil production rate of the interpreted layer. The interpretation conclusions are shown in Table3.Table 3
Calculation of oil production rate.
HorizonSWPNN (%)SWCAL (%)Perforation water production rate (%)Interpretation conclusionA161.642.848.5Low water flooded layerA258.741.246.9Low water flooded layerThe results show that the accuracy of the committee machine output parameters in the calculation of residual oil saturation has improved, and in terms of water production rate, the value calculated after parameter optimization is significantly closer to the value of water production rate, indicating that the regression method of the capture section parameter value in the block through the committee machine model is reliable and can be popularized and applied.


## 5. Conclusion

Aiming at the problem that the capture cross-section value of each component of the formation cannot be accurately determined in the process of obtaining water saturation from the volume model, a method for parameter inversion based on the committee machine method is proposed. A static committee machine is constructed by combining multiple experts using genetic algorithm. This method is used to explain and analyze the water flooding situation of some wells inX oilfield and to verify the results of actual measures. It is confirmed that the method has high accuracy in practical application.(1)
The three experts selected by the committee machine are different from each other. The combination of multiple experts can avoid the problem of unreliable prediction results of a single algorithm. The goodness of fit of the results output by the committee machine reaches 0.6825, which is close to the actual situation of the formation. Using a genetic algorithm as a combiner can better express the differences between experts than a simple average method, resulting in a better final output effect.(2)
The macrocapture cross section after parameter optimization through the committee machine model is more scientific and accurate than that selected according to the theoretical empirical value. The water saturation obtained in the volume model is closer to the perforation water production rate, which provides a new idea for the quantitative interpretation of residual oil saturation logging and has guiding significance for the selection of the next key production target layer in the oilfield.(3)
In residual oil saturation logging, the optimization of macrocapture cross-section parameters of rock skeleton, formation water, oil and gas, and mud in the volume model is essentially a regression problem. The model’s application effect is directly related to data quality. The fitting effect of the model output can be improved by increasing the number of actual logging data and enriching the sample types. Furthermore, how to choose the basic expert network with the best performance to combine and improve the model’s generalization ability requires more research.



---

*Source: 6440600-2022-11-17.xml*