A New Quantitative Inversion Method for Reservoir Interwell Dynamic Connectivity at Extra High Water Cut Stage

Research Article

Austin J Hydrol. 2015; 2(1): 1015.

A New Quantitative Inversion Method for Reservoir Interwell Dynamic Connectivity at Extra High Water Cut Stage

Wen H¹*, Sun N², Liu Y¹, Lu J³ and Li B³

¹Petroleum Engineering Institute, Northeast Petroleum University, China

²Gas Recovery Plant of Songyuan, Jinlin Oilfield Company, China

35th Oil Recovery Plant, Daqing Oilfield Limited Company, China

*Corresponding author: Hua Wen, Petroleum Engineering Institute, Northeast Petroleum University, Daqing 163318, China

Received: July 25, 2015; Accepted: October 05, 2015; Published: October 14, 2015


The reservoir interwell water flooding status and oil-water movement direction at extra high water cut stage are complex, and inefficient and ineffective circling are becoming increasingly serious. The strengthening inference of reservoir interwell dynamic connectivity has been an important part of reservoir evaluation and dynamic analysis. According to the system analysis method, a weighted grey correlation model based on optimized entropy is proposed, and the inversion decision model for characterizing interwell dynamic connectivity of injection-production system is established, and reservoir interwell dynamic connectivity is inversed quantitatively with injection-production dynamic data. The interwell dynamic connectivity in the G20-39 district of Daqing LaSaXing oil field is inversed by using the established model; the interwell dynamic connectivity map is drawn. By adjusting the measure schemes for oil and water wells, a better development effect can be obtained, which verifies the validity of the presented method. The practice shows that this method has strong adaptability and accurate inversion result, the interwell dynamic connectivity inversed by this method is accord with the actual situation of reservoir, and can provide the foundation for the implementation of development adjustment measure of LaSaXing oilfield at extra high water cut stage.

Keywords: Grey system; Grey correlation degree; Grey entropy; Interwell dynamic connectivity; Reservoir dynamic; Inversion


Strengthening the estimation of d reservoir interwell dynamic connectivity and identifying the Injected water flow direction have been an important part of reservoir evaluation and dynamic analysis at extra high water cut stage, and it is important significance for implementing the scheme and enhancing oil recovery [1-2]. At present, the frequently-used methods of reservoir interwell dynamic connectivity mainly include geochemical method, tracer test, well testing, numerical simulation, and soon [3-8]. The implementations of these methods are more complex, expensive, and can affect the normal production of oilfield. In recent years, it has been mainly popular that the interwell dynamic connectivity is inversed through using production data analyzed by statistical methods at home and abroad, which has been an important method for research the reservoir connectivity, the method models mainly include correlation analysis model [9-10], multiple linear regression model [11], capacitance model [12-14], system analysis model [15-16], and soon. The correlation analysis method is one of the reservoir engineering methods for reservoir dynamic analysis. The grey correlation analysis in grey theory is the most widely used and applied to the reservoir dynamic analysis because of its simple calculation, clear mechanism, advantages of clear required less data and strong adaptability. Wang [17] and Liu [18] separately calculated the correlation of interwell dynamic data by using the rate correlation degree and gray correlation analysis method, then the strong or weak dynamic connectivity between oil and water well is quantitatively determined through the correlation value. But both the rate correlation degree and the grey correlation analysis method take the simple arithmetic average value of the grey correlation coefficient at each time point as the grey correlation degree, which may cause the loss of information and local point correlation value controls the whole grey correlation trend. Multiple linear regression model, capacitance model and system analysis model have disadvantages of complex calculation, difficult application to the inversion for interwell dynamic connectivity, and they cannot inverse the interwell connectivity between production wells without considering interwell fluid interference between production wells. Aiming at above problems, this paper views the reservoir as a complicated nonlinear grey system, grey correlation coefficient is unequally weighted average to calculate grey correlation degree through combining with the grey sequence entropy analysis method, the weighted grey correlation model based on optimized entropy is established, and the inversion decision model for characterizing interwell dynamic connectivity of injection-production system is established, and reservoir interwell dynamic connectivity is inversed quantitatively with injection-production dynamic data.

Interwell Dynamic Connectivity Model

The principle

Reservoir is a complicated nonlinear dynamic balance system, according to the systematology idea, injection well, production well and interwell pore throat passage of reservoir are viewed as a whole grey system. Practice has proved that the production well performance of reservoir at extra high water cut stage is mainly affected by water flooding situation of corresponding injection well and the fluid disturbance around the production well. After water is injected into the reservoir, the crude oil and part of the inefficient, ineffective water are displaced to oil well through the good connectivity reservoir, during displacement, the change of water injection can make dynamic index of corresponding production well fluctuate, or the enhanced liquid, controlled liquid measures of a production may make dynamic index of around production well fluctuate, which is a feature reflection of reservoir connectivity between wells. The interwell grey correlation degree reflects the strength of the interaction between wells. The interwell connectivity and oil-water movement direction can be determined through calculating the grey correlation degree between the input signal (water injection) of injection well and the output signal of corresponding production well (liquid production, oil production, water cut, bottom hole flowing pressure) in a certain period of time, or through calculating the correlation degree between the mass production of a production well and the output signal of the surrounding production well. Therefore, injection-production dynamic data are used to inverse the reservoir interwell dynamic connectivity, which not only can obtain the more accurate and effective result, also is a relatively simpler and cheaper method. If the grey correlation between the water injection of a injection well and the dynamic index of a surrounding production well is high, it shows that the water flooding degree of the injection well to the production well is strong, the reservoir interwell connectivity is good, the water flow direction mainly points to the production well which is mainly affected by the injection well. If the grey correlation degree between the liquid production of a production well and the output signal of a surrounding production well is high, it shows that the connectivity between the production well is good, the interwell interference degree is strong, water flow direction occurs steering, the surrounding production well is mainly controlled by the production well, clarifying these problems can formulate the feasible development adjustment measures, that is the directions between injection and production wells in which is the increasing water injection is to increase oil production, which is the increasing injection is to increase water production (It indicates that the preponderant flow channel has formed ), which is the increasing of oil production and water production, and the directions between production wells in which should be to enhance liquid, and which should be to control liquid, thus it provides reliable basis for adjusting the injection-production relationship and liquid production structure and changing the liquid flow direction.

Establishment of inversion model

According to the above principle, Supposing there are n injection wells w1,w2…, wn and m production wells O1,O2,…,Om around a production well in system, eight correlations are established between the oil production, liquid production, water cut, Bottom Hole Flowing Pressure (BHP) of a production well and the water injection of corresponding injection well and the liquid production of surrounding production well, which are water-oil, waterliquid, water-water cut, water-BHP, liquid-oil, liquid-liquid, liquidwatercut, liquid-BHP, theirs correlation degrees are respectively calculated to construct an attribute matrix E:

E=    oil    liq        f w    p wf w 1 w n o 1 o m ( γ 11 γ 12 γ 13 γ 14 γ n1 γ n2 γ n3 γ n4 γ n+11 γ n+12 γ n+13 γ n+14 γ n+m1 γ n+m2 γ n+m3 γ n+m4 )                 (1) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyraiabg2da9uaabeqaceaaaeaafaqabeqaeaaaaeaacaqGGaGaaeiiaiaabccacaWGVbGaamyAaiaadYgaaeaacaqGGaGaaeiiaiaabccacaWGSbGaamyAaiaadghaaeaacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaWGMbWaaSbaaSqaaiaadEhaaeqaaaGcbaGaaeiiaiaabccacaWGWbWaaSbaaSqaaiaadEhacaWGMbaabeaaaaaakeaafaqabeGbbaaaaeaacaqG3bWaaSbaaSqaaiaaigdaaeqaaaGcbaGaeSO7I0eabaGaam4DamaaBaaaleaacaWGUbaabeaaaOqaaiaad+gadaWgaaWcbaGaaGymaaqabaaakeaacqWIUlstaeaacaWGVbWaaSbaaSqaaiaad2gaaeqaaaaakmaabmaabaqbaeqabyabaaaaaeaacqaHZoWzdaWgaaWcbaGaaGymaiaaigdaaeqaaaGcbaGaeq4SdC2aaSbaaSqaaiaaigdacaaIYaaabeaaaOqaaiabeo7aNnaaBaaaleaacaaIXaGaaG4maaqabaaakeaacqaHZoWzdaWgaaWcbaGaaGymaiaaisdaaeqaaaGcbaGaeSO7I0eabaGaeSO7I0eabaGaeSO7I0eabaGaeSO7I0eabaGaeq4SdC2aaSbaaSqaaiaad6gacaaIXaaabeaaaOqaaiabeo7aNnaaBaaaleaacaWGUbGaaGOmaaqabaaakeaacqaHZoWzdaWgaaWcbaGaamOBaiaaiodaaeqaaaGcbaGaeq4SdC2aaSbaaSqaaiaad6gacaaI0aaabeaaaOqaaiabeo7aNnaaBaaaleaacaWGUbGaae4kaiaaigdacaaIXaaabeaaaOqaaiabeo7aNnaaBaaaleaacaWGUbGaae4kaiaaigdacaaIYaaabeaaaOqaaiabeo7aNnaaBaaaleaacaWGUbGaae4kaiaaigdacaaIZaaabeaaaOqaaiabeo7aNnaaBaaaleaacaWGUbGaae4kaiaaigdacaaI0aaabeaaaOqaaiabl6Uinbqaaiabl6Uinbqaaiabl6Uinbqaaiabl6Uinbqaaiabeo7aNnaaBaaaleaacaWGUbGaae4kaiaad2gacaaIXaaabeaaaOqaaiabeo7aNnaaBaaaleaacaWGUbGaae4kaiaad2gacaaIYaaabeaaaOqaaiabeo7aNnaaBaaaleaacaWGUbGaae4kaiaad2gacaaIZaaabeaaaOqaaiabeo7aNnaaBaaaleaacaWGUbGaae4kaiaad2gacaaI0aaabeaaaaaakiaawIcacaGLPaaaaaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabIcacaqGXaGaaeykaaaa@B90D@

A comprehensive correlation matrix is calculated and obtained through comprehensively combining with the weight matrix Q of oil production, liquid production, water cut and BHP, that is the inversion model of the interwell dynamic connectivity coefficient T between the surrounding production well, injection well and center production well:

T=E Q T = w 1 w n o 1 o m ( T 1 T n T n+1 T n+m )                             (2) MathType@MTEF@5@5@+=feaaguart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamivaiabg2da9iaadweacqGHxkcXcaWGrbWaaWbaaSqabeaacaWGubaaaOGaeyypa0tbaeqabyqaaaaabaGaam4DamaaBaaaleaacaaIXaaabeaaaOqaaiabl6UinbqaaiaadEhadaWgaaWcbaGaamOBaaqabaaakeaacaWGVbWaaSbaaSqaaiaaigdaaeqaaaGcbaGaeSO7I0eabaGaam4BamaaBaaaleaacaWGTbaabeaaaaGcdaqadaqaauaabeqageaaaaqaaiaadsfadaWgaaWcbaGaaGymaaqabaaakeaacqWIUlstaeaacaWGubWaaSbaaSqaaiaad6gaaeqaaaGcbaGaamivamaaBaaaleaacaWGUbGaae4kaiaaigdaaeqaaaGcbaGaeSO7I0eabaGaamivamaaBaaaleaacaWGUbGaae4kaiaad2gaaeqaaaaaaOGaayjkaiaawMcaaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeikaiaabkdacaqGPaaaaa@6E6B@

In the formula, the weight matrix Q = (0.3419, 0.2503, 0.3253, 0.0826).

The correlation degree γ in the matrix E is solved by using the weighted grey correlation model based on the grey entropy optimization. Supposing the oil production or liquid production, water injection and so on as the system time sequence:

3rd pending

Where, X0 is the main behavior sequence, Xi is the comparison sequence in system. For ρ€(0,1), there is

γ( x 0 (k), x i (k))= min i min k | x 0 (k) x i (k) |+ρ max i max k | x 0 (k) x i (k) | | x 0 (k) x i (k) |+ρ max i max k | x 0 (k) x i (k) |                (4) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4SdCMaaiikaiaadIhadaWgaaWcbaGaaGimaaqabaGccaGGOaGaam4AaiaacMcacaGGSaGaamiEamaaBaaaleaacaWGPbaabeaakiaacIcacaWGRbGaaiykaiaacMcacqGH9aqpdaWcaaqaamaaxababaGaciyBaiaacMgacaGGUbaaleaacaWGPbaabeaakmaaxababaGaciyBaiaacMgacaGGUbaaleaacaWGRbaabeaakmaaemaabaGaamiEamaaBaaaleaacaaIWaaabeaakiaacIcacaWGRbGaaiykaiabgkHiTiaadIhadaWgaaWcbaGaamyAaaqabaGccaGGOaGaam4AaiaacMcaaiaawEa7caGLiWoacqGHRaWkcqaHbpGCdaWfqaqaaiGac2gacaGGHbGaaiiEaaWcbaGaamyAaaqabaGcdaWfqaqaaiGac2gacaGGHbGaaiiEaaWcbaGaam4AaaqabaGcdaabdaqaaiaadIhadaWgaaWcbaGaaGimaaqabaGccaGGOaGaam4AaiaacMcacqGHsislcaWG4bWaaSbaaSqaaiaadMgaaeqaaOGaaiikaiaadUgacaGGPaaacaGLhWUaayjcSdaabaWaaqWaaeaacaWG4bWaaSbaaSqaaiaaicdaaeqaaOGaaiikaiaadUgacaGGPaGaeyOeI0IaamiEamaaBaaaleaacaWGPbaabeaakiaacIcacaWGRbGaaiykaaGaay5bSlaawIa7aiabgUcaRiabeg8aYnaaxababaGaciyBaiaacggacaGG4baaleaacaWGPbaabeaakmaaxababaGaciyBaiaacggacaGG4baaleaacaWGRbaabeaakmaaemaabaGaamiEamaaBaaaleaacaaIWaaabeaakiaacIcacaWGRbGaaiykaiabgkHiTiaadIhadaWgaaWcbaGaamyAaaqabaGccaGGOaGaam4AaiaacMcaaiaawEa7caGLiWoaaaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeikaiaabsdacaqGPaaaaa@9F13@

γ( X 0 , X i )= k=1 n ω k γ( x 0 (k), x i (k))                    (5) MathType@MTEF@5@5@+=feaaguart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4SdCMaaiikaiaadIfadaWgaaWcbaGaaGimaaqabaGccaGGSaGaamiwamaaBaaaleaacaWGPbaabeaakiaacMcacqGH9aqpdaaeWbqaaiabeM8a3naaBaaaleaacaWGRbaabeaakiabeo7aNjaacIcacaWG4bWaaSbaaSqaaiaaicdaaeqaaOGaaiikaiaadUgacaGGPaGaaiilaiaadIhadaWgaaWcbaGaamyAaaqabaGccaGGOaGaam4AaiaacMcacaGGPaaaleaacaWGRbGaeyypa0JaaGymaaqaaiaad6gaa0GaeyyeIuoakiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGOaGaaeynaiaabMcaaaa@61D2@

γ(x0(k),xi(k)) is called as grey correlation coefficient, γ(X0,Xi) is the weighted grey correlation degree between X0 and Xi. Where, ρ is the resolution coefficient is equal to 0.5; ωk is the weight of point?, k=1 n ω k =1(0< ω k <1) MathType@MTEF@5@5@+=feaaguart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaabCaeaacqaHjpWDdaWgaaWcbaGaam4AaaqabaGccqGH9aqpcaaIXaGaaiikaiaaicdacqGH8aapcqaHjpWDdaWgaaWcbaGaam4AaaqabaGccqGH8aapcaaIXaGaaiykaaWcbaGaam4Aaiabg2da9iaaigdaaeaacaWGUbaaniabggHiLdaaaa@4856@ ωkIn order to solve the model parameter ?k, the grey entropy analysis method is used to determine the weight of correlation coefficient [19].

is the model parameter to be determined.

2 Solution of model parameters

In order to solve the model parameter ωk, the grey entropy analysis method is used to determine the weight of correlation coefficient [19].

The specific process is as follows.

Map : RjPj

Definition 1, calling the mapping:

p k = ω k γ( x 0 (k), x i (k)) k=1 n ω k γ( x 0 (k), x i (k))              (6) MathType@MTEF@5@5@+=feaaguart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiCamaaBaaaleaacaWGRbaabeaakiabg2da9maalaaabaGaeqyYdC3aaSbaaSqaaiaadUgaaeqaaOGaeq4SdCMaaiikaiaadIhadaWgaaWcbaGaaGimaaqabaGccaGGOaGaam4AaiaacMcacaGGSaGaamiEamaaBaaaleaacaWGPbaabeaakiaacIcacaWGRbGaaiykaiaacMcaaeaadaaeWbqaaiabeM8a3naaBaaaleaacaWGRbaabeaakiabeo7aNjaacIcacaWG4bWaaSbaaSqaaiaaicdaaeqaaOGaaiikaiaadUgacaGGPaGaaiilaiaadIhadaWgaaWcbaGaamyAaaqabaGccaGGOaGaam4AaiaacMcacaGGPaaaleaacaWGRbGaeyypa0JaaGymaaqaaiaad6gaa0GaeyyeIuoaaaGccaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeikaiaabAdacaqGPaaaaa@67F1@

Above formula is the distribution mapping of the weighted grey correlation coefficient, the mapping ?? is the distribution density of the weighted grey correlation coefficient, and meets the k=1 n p k =1 MathType@MTEF@5@5@+=feaaguart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaabCaeaacaWGWbWaaSbaaSqaaiaadUgaaeqaaOGaeyypa0JaaGymaaWcbaGaam4Aaiabg2da9iaaigdaaeaacaWGUbaaniabggHiLdaaaa@3FB5@

Definition 2, calling function:

H ( R i )= k=1 n p k ln p k                   (7) MathType@MTEF@5@5@+=feaaguart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamisamaaBaaaleaacqGHxkcXaeqaaOGaaiikaiaadkfadaWgaaWcbaGaamyAaaqabaGccaGGPaGaeyypa0JaeyOeI0YaaabCaeaacaWGWbWaaSbaaSqaaiaadUgaaeqaaOGaciiBaiaac6gacaWGWbWaaSbaaSqaaiaadUgaaeqaaaqaaiaadUgacqGH9aqpcaaIXaaabaGaamOBaaqdcqGHris5aOGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeikaiaabEdacaqGPaaaaa@57C2@

Above formula is the weighted grey correlation entropy of Xi, any drives ?? to equally change will make the correlation entropy H (Ri) increase.

Considering the influence of the comparison sequence on the main behavior sequence should be remain stable, that is the weighted correlation coefficient should maintain balanced, ?? should maintain balanced, so there is:

max H ( R i )= k=1 n p k ln p k                            (8) s.t. p k 0 s.t. k=1 n p k =1 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaaciGGTbGaaiyyaiaacIhacaWGibWaaSbaaSqaaiabgEPiedqabaGccaGGOaGaamOuamaaBaaaleaacaWGPbaabeaakiaacMcacqGH9aqpcqGHsisldaaeWbqaaiaadchadaWgaaWcbaGaam4AaaqabaGcciGGSbGaaiOBaiaadchadaWgaaWcbaGaam4AaaqabaaabaGaam4Aaiabg2da9iaaigdaaeaacaWGUbaaniabggHiLdGccaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGOaGaaeioaiaabMcaaeaacaWGZbGaaiOlaiaadshacaGGUaGaaWiiaiaadchadaWgaaWcbaGaam4AaaqabaGccqGHLjYScaaIWaaabaGaam4Caiaac6cacaWG0bGaaiOlaiaamccadaaeWbqaaiaadchadaWgaaWcbaGaam4AaaqabaGccqGH9aqpcaaIXaaaleaacaWGRbGaeyypa0JaaGymaaqaaiaad6gaa0GaeyyeIuoaaaaa@76BC@

For the extremum problem of above function, the weight vector is obtained through constructing the Lagrange function and operating matrix under the constraint condition of grey correlation entropy maximization:

Ω = Г-1 b (9)


Γ=[ γ 1 γ 2 γ 2 γ 3 γ n1 γ n 1 1 1 1 1 ] MathType@MTEF@5@5@+=feaaguart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeu4KdCKaeyypa0ZaamWaaeaafaqabeqbfaaaaaqaaiabeo7aNnaaBaaaleaacaaIXaaabeaaaOqaaiabgkHiTiabeo7aNnaaBaaaleaacaaIYaaabeaaaOqaaaqaaaqaaaqaaaqaaiabeo7aNnaaBeaaleaacaaIYaaabeaaaOqaaiabgkHiTiabeo7aNnaaBaaaleaacaaIZaaabeaaaOqaaaqaaaqaaaqaaaqaaiablgVipbqaaaqaaaqaaaqaaaqaaaqaaiabeo7aNnaaBaaaleaacaWGUbGaeyOeI0IaaGymaaqabaaakeaacqaHZoWzdaWgaaWcbaGaamOBaaqabaaakeaacaaIXaaabaGaaGymaaqaaiaaigdaaeaacaaIXaaabaGaaGymaaaaaiaawUfacaGLDbaaaaa@539F@ Ω=[ ω 1 ω 2 ω n1 ω n ] MathType@MTEF@5@5@+=feaaguart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuyQdCLaeyypa0ZaamWaaeaafaqabeqbbaaaaeaacqaHjpWDdaWgaaWcbaGaaGymaaqabaaakeaacqaHjpWDdaWgaaWcbaGaaGOmaaqabaaakeaacqWIUlstaeaacqaHjpWDdaWgaaWcbaGaamOBaiabgkHiTiaaigdaaeqaaaGcbaGaeqyYdC3aaSbaaSqaaiaad6gaaeqaaaaaaOGaay5waiaaw2faaaaa@498B@ b=[ 0 0 0 1 ]                     (10) MathType@MTEF@5@5@+=feaaguart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOyaiabg2da9maadmaabaqbaeqabuqaaaaabaGaaGimaaqaaiaaicdaaeaacqWIUlstaeaacaaIWaaabaGaaGymaaaaaiaawUfacaGLDbaacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGOaGaaeymaiaabcdacaqGPaaaaa@4ED9@

The model with higher resolution can well control the influence of fluctuation of the partial correlation coefficient on overall correlation. In view of the advantages of the model can accurately calculate the weight ωk, and then inverse the reservoir interwell dynamic connectivity.

Example Analysis

The G20-39 block in Daqing LaSaXing oilfield was selected to apply the above model. This block mainly develops PI reservoir group, and there are 13 sedimentary units. The main reservoir sand body distribution is striped and sheet like, and non-essential reservoir sand body distribution is narrow stripe, scattered and poorly developed. The block was put into operation in 1983, it applied the four-spot water flooding pattern, and well spacing was 350 meters. Since 1999 the block has entered into extra high water cut stage, the production decreasing magnitude is large, the water cut is rising rapidly, it is urgent need to adjust the current water injection scheme, and select the applicable measure for injection well and production well, thus, it is very important to carry out the interwell dynamic connectivity research. The interwell dynamic connectivity in the block was inversed according to its computation steps.

Understanding of interwell dynamic connectivity

The first thing is to understand the interwell connectivity relationship, the connectivity degree can be determined through the correlation degree value between injection well and production well in the each injection-production system. According to the evaluation standard definition of reservoir dynamic grey correlation degree in LaSaXing oilfield, if the correlation degree T is greater than or equal to 0.8, the connectivity is good; if the correlation degree T is between 0.6 and 0.8, the connectivity is little good; if the correlation degree T is less than or equal to 0.6, the connectivity is poor. The model was used to analyze the interwell correlation degree respectively between the monthly water injection of injection well and monthly oil production, monthly liquid production, monthly water cut, BHP of surrounding production well, respectively between the monthly liquid production and monthly oil production, monthly liquid production, monthly water cut, bottom hole flowing pressure of surrounding production well during 2013 in the G20-39 block, by which could determine the connectivity between injection well and production well, and between production well. Seen from the (Table 1) and (Figure 1), the results show that the connectivity degree between injection well and the surrounding production well are not uniform. The correlation degree analysis results are in good agreement with the actual connectivity, which prove that the inversion analysis is effective and accurate. Seen from the results, the interwell connectivities in the G20-39 block belong to the little good class.