^{1}

^{*}

^{1}

Oil field waterflooding is a complex man-controlled systematic behavior, and the related evaluation methods vary greatly. This paper put forward a fuzzy comprehensive method of evaluating controlled development level by analysis of the macroscopic evaluation to oil field waterflooding effect with combination of original reservoir geological state. This fuzzy evaluation technique bears unique advantages because there is little difference among evaluation indexes which represent the dynamic and static state of regional neighborhood of development units (blocks, Production Company , etc.). Not only the mathematical method for evaluating oil field waterflooding effect is set up, but also the method is applied in three blocks of D oil field. The calculated results show the effectiveness and practicability of the method.

Waterflooding as the dominant technology in oil field development is facing a universal problem of associated high water cut and how to improve the waterflooding efficiency during high water cut stage has become a highlight for all the researchers in the world. Many of the oil fields in China have entered extra-high water cut stage after several decades of development, which are now encountered by the conjuncture of intractable remaining oil exploitation, deficient reserve replacements and continuous production declination. Reasonable and proper development effect evaluation would help to determine the direction for remaining oil exploitation and adjustment and thus guiding more effective and scientific reservoir development, as well as providing reference to the composition of development plans and long-term decisions. The essence of evaluating the operational performance of any system is to judge the scientism, rationality and optimality of the system. Comprehensive evaluation of the operational performance of a system should set off from the perspective of systematology. Judging from the above point of view, during development effect evaluation, the first step is to find out evaluation indexes [

Many of the development effect evaluation methods usually lean on oil production mechanism [

On the other side, many single indexe [

V = { h z c } = { good normal poor } .

The universe of ambiguity factors could be defined as two aspects, first, factors that decide the original geologic property of the reservoir, and second, factors that dominate the development level of the oil field. To determine whether the water flooding effect is good or not mainly depends on the amount of oil production, which if explained from systematology is mainly controlled by two attributes, one is original reservoir geologic property and the other is controlled development level. Original reservoir geologic property refers to the natural oil productivity of the reservoir, while the controlled development level refers to treatments taken during water flooding process, including water injection pattern, well array, perforation, fracturing and acid treating, etc., some reservoirs with predominance natural oil productivity probably would not obtain satisfactory crude oil production after exploration, only by scientific and rational development control could the maximum crude oil recovery rate be reached. This is the purpose of comprehensive evaluation, to see if the reservoirs with different natural oil productivity are explored by reasonable development control and if different reservoirs are matched with corresponding development treatments. Water flooding effect evaluation would cover evaluations on both original reservoir geologic property and controlled development level and in addition the overall effects based on the former two evaluations.

Contexture of this paper follows the contents of regular comprehensive evaluation, first building an evaluation index system, and then calculating the evaluation index weight, successively, setting up unified standard function with fuzzy membership degree for each evaluation index, and finally the comprehensive evaluation.

Through system theory and control theory analysis, the structure of systematic behavior of oil field waterflooding can be found out. Water well system is defined as the group composed of all water wells, oil well system is defined as the group composed of all oil wells.

The index system of original reservoir property is the information indicator of reservoir system, while the index system of controlled development level is the information indicator of oil and water well system. Based on comprehensive evaluation of the oil field reality (D oil field), through oil field development engineering mechanism, the index system of the two oil field development attributes can be defined as follows by using quantitative and qualitative analysis (multiple correlation analysis).

All these indexes can be screen out as follows, the details are referenced to [

Factors for evaluating original reservoir geological property include:

1) channel sand proportion R h d ;

2) effective permeability K e ;

3) net thickness h e ;

4) coefficient of permeability variation V K ;

5) initial oil saturation S o ;

6) oil viscosity μ o ;

7) difference between reservoir pressure and saturation pressure P b ;

8) reserves abundance A q ;

9) displacement efficiency η D .

Factors for evaluating controlled development level include:

1) total pressure difference Δ P ;

2) oil and water well ratio R j ;

3) well spacing density A j ;

4) cumulative water storage rate γ w ;

5) cumulative injection production ratio R r c ;

6) waterflooding controlled degree C h ;

7) reserves recovery degree D h .

For comprehensive evaluation to any complex system, it is not easy to determine the index weight rationally after drawing up evaluation index system. For the systematic behavior of evaluation, if there are too many attributes index (for example, 16 indexes for waterflooding effect evaluation), only mathematic methods can be used to determine the weight. But this kind of quantitative method relies greatly on the samples need to be evaluated. There are three blocks (samples) in this D oil field. According to oil field waterflooding system attributes, the multiple correlation coefficient and variation coefficient are used to calculate index weight in this paper. Then, the final weight is got by geometric averaging the two weights.

Assume that the multiple correlation coefficients with m indexes are ρ 1 , ρ 2 , ⋯ , ρ m , then the number k coefficient reflects the ability of all the other indexes except k index to replace k index. Hence, the larger the k index, the smaller its function, so the reciprocal value of multiple correlation coefficient can be taken as weight. Then the weight formula for multiple correlation coefficients is developed,

v k 1 = ∏ i = 1 8 ρ i 2 ρ k 2 ∑ i = 1 8 ∏ j ≠ i , j = 1 8 ρ j 2 ( k = 1 , 2 , ⋯ , m )

Multiple correlation coefficient can be calculated from data in

The weight is calculated to be:

V 1 = ( v 1 1 , v 2 1 , ⋯ , v 9 1 ) T = ( 0.172 , 0.105 , 0.088 , 0.123 , 0.101 , 0.098 , 0.088 , 0.118 ) T (3.1)

In addition, from statistical theory we can know that, the larger the variance of each index, the stronger the corresponding attributes of each index for evaluating the targets, the more attention should be paid to this index, thus the larger of this index weight. This is the so called theory of determining index weight by variation coefficient.

As above mentioned, by using data in

V 2 = ( v 1 2 , v 2 2 , ⋯ , v 9 2 ) T = ( 0.102 , 0.081 , 0.186 , 0.053 , 0.026 , 0.067 , 0.221 , 0.216 , 0.048 ) T (3.2)

In order to avoid the unilateralism of different methods for calculating weights, comprehensive weight is obtained by geometric averaging the weights obtained from the two above methods.

ρ 1 2 | ρ 2 2 | ρ 3 2 | ρ 4 2 | ρ 5 2 | ρ 6 2 | ρ 7 2 | ρ 8 2 | ρ 9 2 |
---|---|---|---|---|---|---|---|---|

0.172 | 0.105 | 0.088 | 0.123 | 0.106 | 0.101 | 0.098 | 0.088 | 0.118 |

Formula for geometric average is:

V = V 1 ∘ V 2

i.e.

v k = v k 1 v k 2 , k = 1 , 2 , ⋯ , 9

Comprehensive weight is obtained from the above formula

V = ( v 1 , v 2 , ⋯ , v 9 ) T = ( 0.1428 , 0.0993 , 0.1377 , 0.0874 , 0.0567 , 0.0888 , 0.1583 , 0.1484 , 0.0806 ) (3.3)

In order to show the relationship between the indexes weight with the index itself, these data is tabularized in

Based on the same principle, the weight of control index can be obtained. The following formula (3.4)-(3.6) are weight from multiple correlation coefficient, weight from variation coefficient and comprehensive weight from geometric average respectively.

W 1 = ( w 1 1 , w 2 1 , ⋯ , w 7 1 ) T = ( 0.158 , 0.208 , 0.146 , 0.117 , 0.128 , 0.126 , 0.117 ) T (3.4)

W 2 = ( w 1 2 , w 2 2 , ⋯ , w 7 2 ) T = ( 0.315 , 0.102 , 0.248 , 0.148 , 0.073 , 0.039 , 0.039 ) T (3.5)

W = W 1 ∘ W 2 = ( w 1 , w 2 , ⋯ , w 7 ) T = ( 0.109 , 0.195 , 0.087 , 0.121 , 0.124 , 0.151 , 0.132 , 0.08 ) T (3.6)

These data is tabularized in

Channel sand proportion | Effective permeability | Net thickness | Permeability variation coefficient | Oil saturation |
---|---|---|---|---|

0.1428 | 0.0993 | 0.1377 | 0.0874 | 0.0567 |

Petroleum viscosity | Difference between reservoir pressure and saturation pressure | Reserves abybdabce | Displacement efficiengcy | |

0.0888 | 0.1583 | 0.1484 | 0.0806 |

Controll index | Total pressure difference | Oil and water well ratio | Well spacing density | Cumulative water storage rate |
---|---|---|---|---|

weight | 0.2516 | 0.1553 | 0.2028 | 0.1404 |

Controll index | Cumulative injectio production ratio | Waterflooding controlling degree | Reserve producing degree | |

weight | 0.103 | 0.0748 | 0.0722 |

A unified evaluation standard is needed in order to rationally evaluate each index; this is the function with fuzzy membership degree of each index. Based on characteristic analysis of 9 geological indexes and 7 controlling indexes, it is found that the comprehensive evaluation index can be classified into two types: one is tendency index, the other is moderate index.

For tendency index, there are two situations during evaluation: one is, the larger the index value, the larger the evaluation value (abbreviated as the larger the better); another one is, the larger the index value, the smaller the evaluation value (abbreviated as the smaller the better). Thus there are only two types of functions with fuzzy membership degree for tendency index. In this situation, assume the evaluation values corresponding to fuzzy remarks range

V = { h z c } = { good normal poor }

are x g , x n , x b respectively (subscript g stands for good, n for normal, b for poor).

For moderate index, it would be more complex when setting up corresponding function with fuzzy membership degree. First, the “moderate value x o ” (o stand for optimization) of moderate index should be calculated, or to be called “optimized as rational value”. Then set up function with fuzzy membership degree x * = x / x o . Apparently, now the more close of x * to 1, the larger the evaluation value, its standards values of “good, normal, poor” are all symmetrical to 1. In this situation, assume the evaluation values corresponding to fuzzy range

V = { h z c } = { good normal poor }

are x g + 1 , x g − 1 , x n + 1 , x n − 1 , x b + 1 , x b − 1 respectively, among them, superscript +1 stands for points larger than 1, −1 stands for points smaller than 1.

Based on above analysis, the broken line function with fuzzy membership degree of three types of index can be set up.

l Function with fuzzy membership degree of the larger the better index x

μ h ( x ) = { 0 , x ≤ x n x − x n x g − x n , x n < x < x g 1 , x ≥ x g (4.1)

μ z ( x ) = { 0 , x ≤ x b x − x b x n − x b , x b < x ≤ x n x g − x x g − x n , x n < x ≤ x g 0 , x ≥ x g (4.2)

μ c ( x ) = { 1 , x ≤ x b x n − x x n − x b , x b < x < x n 0 , x ≥ x n (4.3)

l Function with fuzzy membership degree of the smaller the better index x

μ h ( x ) = { 1 , x ≤ x b x n − x x n − x b , x b < x ≤ x n 0 , x > x n (4.4)

μ z ( x ) = { 0 , x ≤ x b x − x b x n − x b , x b < x ≤ x n x g − x x g − x n , x n < x ≤ x g 0 , x > x g (4.5)

μ c ( x ) = { 0 , x ≤ x n x − x n x g − x n , x n < x ≤ x g 1 , x > x g (4.6)

l Function with fuzzy membership degree of moderate index x^{*}

μ h ( x * ) = { 0 , x * ≤ x n − 1 x * − x n − 1 x g − 1 − x n − 1 , x n − 1 < x * < x g − 1 1 , x g − 1 ≤ x * ≤ x g + 1 x n + 1 − x * x n + 1 − x g + 1 , x g + 1 < x * < x n + 1 0 , x * ≥ x n + 1 (4.7)

μ z ( x * ) = { 0 , x * ≤ x b − 1 x * − x b − 1 x n − 1 − x b − 1 , x b − 1 < x * ≤ x n − 1 x g − 1 − x * x g − 1 − x n − 1 , x n − 1 ≤ x * < x g − 1 0 , x g − 1 ≤ x * ≤ x g + 1 x * − x g + 1 x n + 1 − x g + 1 , x g + 1 < x * ≤ x n + 1 x b + 1 − x * x b + 1 − x n + 1 , x n + 1 ≤ x * < x b + 1 0 , x * ≥ x b + 1 (4.8)

μ c ( x * ) = { 1 , x * ≤ x b − 1 x n − 1 − x * x n − 1 − x b − 1 , x b − 1 < x * < x n − 1 0 , x n − 1 ≤ x * ≤ x n + 1 x * − x n + 1 x b + 1 − x n + 1 , x n + 1 < x * < x b + 1 1 , x * ≥ x b + 1 (4.9)

There are two index system for oil field waterflooding effect evaluation: Original reservoir state evaluation; controlled development level evaluation. Comprehensive evaluation result can be got be combining the two results. Here formula (5.7) is used.

For oil field waterflooding effect evaluation, the 16 indexes can be classified into three types of functions with fuzzy membership degree, such as (4.1)-(4.3). The standing points for good, normal, poor are needed when we calculatethe fuzzy evaluation vectors V = { h z c } = { good normal poor } of three blocks in D oil field, i.e. evaluation standards x g , x n , x b or x g + 1 , x g − 1 , x n + 1 , x n − 1 , x b + 1 , x b − 1 . This evaluation index can be got at least by three methods: Oil displacement experimental calculation, mechanism calculation of production fluid flow through porous media, statistical calculation. There are three blocks need to be evaluated, all have certain representativeness. For example, according to

blocks | Index1 | Ind2 | Ind3 | Ind4 | Ind5 | Ind6 | Ind7 | Ind8 | Ind9 |
---|---|---|---|---|---|---|---|---|---|

block1 | 0.26 | 0.29 | 52.67 | 0.62 | 0.71 | 8.70 | 0.75 | 608.42 | 0.47 |

block2 | 0.34 | 0.40 | 66.91 | 0.64 | 0.71 | 8.70 | 1.78 | 753.87 | 0.47 |

block3 | 0.34 | 0.33 | 58.34 | 0.62 | 0.70 | 8.70 | 1.30 | 642.06 | 0.46 |

blocks | Index1 | Index2 | Index3 | Index4 | Index5 | Index6 | Index7 |
---|---|---|---|---|---|---|---|

block1 | −3.10 | 2.20 | 70.30 | 14.64 | 0.86 | 96.30 | 94.80 |

block2 | −3.48 | 2.06 | 79.10 | 33.68 | 1.07 | 93.80 | 92.40 |

block3 | −1.53 | 1.71 | 48.93 | 33.82 | 1.12 | 89.40 | 87.60 |

For 16 indexes for comprehensive evaluation, there are 10 the larger the better indexes: channel sand proportion R h d , net thickness h e , initial oil saturation S o , difference between reservoir pressure and saturation pressure P b , reserves abundance A q , displacement coefficient η D , total pressure difference Δ P , cumulative water storage rate γ w , waterflooding controlling degree C h , reserves recovery degree D h . Single index fuzzy evaluation vector can be calculated by formula (4.1)-(4.3).

There are 4 better indexes: effective permeability K e , coefficient of permeability variation V K , oil viscosity μ o , cumulative injection production ratio R r c . They can be obtained by formula (4.4)-(4.6).

There are two moderate indexes: oil water well ratio R j , well spacing density A j . For the two indexes, the rational oil and water well ratio and rational well spacing density of the three blocks should be calculated based on the technique, reservoir property, controlled development program, etc., of D oil field. Then use formula (4.7)-(4.9) to calculate.

The fuzzy evaluation vector of original reservoir geological property can be calculated by the above mentioned method, and it is signed as r i j k ( i = 1 , 2 , ⋯ , 9 ; j = 1 , 2 , 3 ; k = 1 , 2 , 3 ) , in it, subscript i = 1 , 2 , ⋯ , 9 stands for 9 reservoir geological property indexes, subscript j = 1 , 2 , 3 stands for good, normal, poor respectively, k = 1 , 2 , 3 stands for number of 3 blocks.

Similarly, the fuzzy evaluation vector of the controlled development level can be calculated by the above function with fuzzy membership degree

t i j k ( i = 1 , 2 , ⋯ , 7 ; j = 1 , 2 , 3 ; k = 1 , 2 , 3 ) .

The fuzzy evaluation result of each block is obtained by weighting the single index fuzzy vector. Synthetic operator M ( • , ⊕ ) [

D ˜ k = V ∘ R k = ( v 1 , v 2 , ⋯ , v 9 ) ∘ ( r i j k ) 9 × 3 (5.1)

where R k = ( r i j k ) 9 × 3 is the evaluation result of k block, synthetic fuzzy operator “ ∘ ” is

( v 1 , v 2 , ⋯ , v 9 ) ∘ ( r 11 k , r 21 k , ⋯ , r 91 k ) T = inf { ∑ i = 1 9 v i r i 1 k , 1 } (5.2)

when calculating the fuzzy evaluation result of reservoir geological property of three blocks using formula (5.1), since the weight V = ( v 1 , v 2 , ⋯ , v 9 ) T (formula (3.3)) and each factor in fuzzy matrix R k = ( r i j k ) 9 × 3 is smaller than 1, and the sum of R k corresponding column factors is 1, thus the use of synthetic fuzzy operator “ ∘ ” in formula (5.2) is the same with normal multiplication operator, then the fuzzy calculating formula (5.1) can be note as normal matrix multiplication, i.e.

D ˜ k = V T R k = ( v 1 , v 2 , ⋯ , v 9 ) ( r i j k ) 9 × 3

Or to be signed D ˜ = ( d 1 k , d 2 k , d 3 k ) (geological fuzzy evaluation vector), to k = 1 , 2 , 3

d i k = v 1 r 1 i k + v 2 r 2 i k + ⋯ + v 9 r 9 i k , i = 1 , 2 , 3 , k = 1 , 2 , 3 (5.3)

Similarly, the fuzzy evaluation vector for controlled development level can be got C ˜ = ( c 1 k , c 2 k , c 3 k ) T , among them

c i k = w 1 t 1 i k + w 2 t 2 i k + ⋯ + w 7 t 7 i k , i = 1 , 2 , 3 , k = 1 , 2 , 3 (5.4)

The fuzzy evaluation result of three blocks calculated by using formula (5.3) and (5.4) is listed in

For the convenience of application and comparison, we need assign a real number for these fuzzy vectors. Fuzzy levels are assigned as “good, normal, poor”, and the corresponding numerical values are “100, 50, 0” respectively for better and direct understanding, so the numerical value for fuzzy comprehensive evaluation vector ( r 1 , r 2 , r 3 ) T = ( good , mediate , poor ) T is

Y = 100 × r 1 + 50 × r 2 (5.5)

Fuzzy evaluation of vectors of three blocks is determined by using formula (5.5) and listed in

The scores for original reservoir geological property and controlled development level of three blocks in D oil field are calculated by using formula (5.3)-(5.5), the result is listed in attached C. The oil field waterflooding effect is mainly determined by controlled development level. Yet controlling development level depends greatly on original reservoir geological property. Therefore, the waterflooding effect can be defined as score for controlled development level divided by score for reservoir geological property, i.e.

score for comprehensive development effect = score for controlled development level score for reservoir geological property (5.6)

block | Reservoir geology evaluation | Evaluation of controlled development effect | Comprehensive evaluation | |||||||
---|---|---|---|---|---|---|---|---|---|---|

fuzzy | value | fuzzy | value | value | rank | |||||

good | normal | poor | good | normal | poor | |||||

1 | 0.38 | 0.068 | 0.547 | 47.8 | 0.385 | 0.068 | 0.547 | 41.9 | 0.766 | 2 |

2 | 0.17 | 0.218 | 0.609 | 75.1 | 0.174 | 0.218 | 0.609 | 28.3 | 0.626 | 3 |

3 | 0.11 | 0.589 | 0.298 | 64.8 | 0.114 | 0.589 | 0.298 | 40.8 | 0.791 | 1 |

In ideal condition, the formula bears “normal geological property (score 50) and normal control level (score 50)”. Hence, in order to rationally evaluate the comprehensive effect, the above formula should be adapted to the following to calculate comprehensive score:

Y x = Y c − Y ¯ c + 50 Y d − Y ¯ d + 50 (5.7)

Among them, Y x , stands for score for comprehensive development effect;

Y c , score for controlled development level;

Y ¯ c , score for controlled development level of three blocks in D oil field;

Y d , score for original reservoir geological property;

Y ¯ d , score for original reservoir geological property of three blocks in D oil field.

Thus, the comprehensive development effect of each block can be fully demonstrated by the comprehensive score. The comprehensive score of three blocks in D oil field can be obtained by using formula (5.7), which can be sequenced and be listed in

Based on score for original reservoir geological property obtained from formula (5.3) and (5.5), and score for controlled development level from formula (5.4) and (5.5), oil field waterflooding effect can be got after calculating the comprehensive effect score of each block area from formula (5.7). These scores can be comprehended as: development effect is mainly decided by two factors, one is original reservoir geological property, or the original oil production capacity of a block, this is the comprehensive score for reservoir property by fuzzy evaluation. The second one is if the controlled development is rational and scientific, this is the comprehensive score for controlled development level by fuzzy evaluation. Apparently, the better the original reservoir property, the more oil we can get, as long as the development is rationally and scientifically controlled. Even if the reservoir property is not very good, we can still obtain good results by rationally and scientifically making use of reservoir geology. Therefore, as defined in formula (5.7), it is reasonable to take the controlled development level score bears in reservoir property score as the “unified evaluation standard”.

From the result (

Research is partly supported by National Natural Science Foundation of China (Grant No. 12071218).

The authors declare no conflicts of interest regarding the publication of this paper.

Xiao, W. and Liang, Y.S. (2021) Fuzzy Comprehensive Evaluation of Oil Field Waterflooding Effect. International Journal of Geosciences, 12, 560-571. https://doi.org/10.4236/ijg.2021.126031