CN103670369A - Method and device for judging communication condition between injection wells and production wells - Google Patents

Method and device for judging communication condition between injection wells and production wells Download PDF

Info

Publication number
CN103670369A
CN103670369A CN201310681831.5A CN201310681831A CN103670369A CN 103670369 A CN103670369 A CN 103670369A CN 201310681831 A CN201310681831 A CN 201310681831A CN 103670369 A CN103670369 A CN 103670369A
Authority
CN
China
Prior art keywords
mrow
msub
production
injection
wells
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201310681831.5A
Other languages
Chinese (zh)
Other versions
CN103670369B (en
Inventor
杨超
许晓明
李嘉
李彦兰
李华
何辉
孙景民
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201310681831.5A priority Critical patent/CN103670369B/en
Publication of CN103670369A publication Critical patent/CN103670369A/en
Application granted granted Critical
Publication of CN103670369B publication Critical patent/CN103670369B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Injection Moulding Of Plastics Or The Like (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Meat, Egg Or Seafood Products (AREA)

Abstract

The invention discloses a method and a device for judging the communication condition between injection wells and production wells. The method for judging the communication condition between injection wells and production wells comprises the following steps: s1: data points of the oilfield block to be determined are optimized unreasonably; s2: according to the optimized data of S1, calculating an initial value of the injection-production interwell communication coefficient and an initial value of a time lag constant of the oil field block to be determined; s3: substituting the optimized data of S1 and the initial value of the communication coefficient between injection wells and the production well calculated in S2 into an objective function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the objective function value is smaller than the fitting error; s4: and judging the communication condition between the injection wells and the production wells by utilizing the communication coefficient between the injection wells and the production wells, the time lag constant and the interference constant between the production wells which are calculated by the S3. The method and the device are beneficial to improving the rationality of judging the communication condition between the injection wells and the production wells of the high-water-content oil field.

Description

Method and device for judging communication condition between injection wells and production wells
Technical Field
The invention belongs to the field of oil field injection and production, and particularly relates to a method and a device for judging the communication condition between injection wells and production wells.
Background
The process of reservoir exploitation is often accompanied by the continuous reduction of reservoir pressure, and the exploitation work of the reservoir is influenced when the reservoir pressure is reduced to a certain degree. In order to reduce the pressure drop of the reservoir, the reservoir pressure is supplemented by a water flooding development method. When water injection development is implemented, a plurality of water injection wells and production wells are often required to form an injection-production well pattern. Because the geological condition of the oil reservoir is very complicated, the injected water of one water injection well has different driving degrees on different production wells, namely the water driving degrees are different. The water flooding degree has great influence on the yield of the production well, and the yield of the production well can be continuously stable or gradually increased if the water flooding degree is good; if the water flooding level is poor, the production well production will gradually decrease. The communication condition between injection and production wells is an important factor influencing the water flooding degree, so that the judgment of the communication condition between the injection and production wells is very important for the exploitation of the whole oil reservoir.
The existing method for judging the communication condition between injection wells and production wells generally adopts a capacitance resistance method.
The capacitance resistance method is characterized in that the communication condition between injection wells and production wells is represented by an injection-production inter-well communication coefficient (the injection-production inter-well communication coefficient reflects the strength of the communication condition between an injection well and a surrounding production well) and a time lag constant (the time lag constant represents the time required by injected water to reach the production well) which are calculated by establishing an objective function. The established objective function is represented by the sum of the squares of the differences between the calculated fluid production and the actual fluid production for the production well. The objective function considers the medium between an injection well, a production well and an injection-production well in the oil deposit as a complete system. In the system, the injected water is influenced by the communication condition between the injection wells and the production wells when the injected water is spread between the injection wells and the production wells, and the influence is directly reflected in the actual liquid production amount of the production wells.
The capacitance resistance method takes each production well as an independent individual, the calculated liquid production capacity of the production well in the established objective function at a certain time step is represented as the influence of the calculated liquid production capacity of the previous time step, the communication coefficient among injection and production wells and the time lag constant, and the objective function is expressed by a function formula formed by the initial calculated liquid production capacity of the production well, the communication coefficient among injection and production wells and the time lag constant through iterative calculation of the calculated liquid production capacity of the previous time step. The initial calculated fluid production is the same as the initial actual fluid production of the production well. By substituting the production data of the oil field, the communication coefficient and the time lag constant between the injection wells and the production wells can be calculated by a capacitance resistance method, and the communication condition between the injection wells and the production wells can be reflected by the communication coefficient and the time lag constant between the injection wells and the production wells.
In the prior art, the objective function established by the capacitance-resistance method does not take the interferences of other production wells into consideration. The high water content oil field is generally developed by adopting a dense well pattern, the interference among production wells has large influence on the communication condition among injection and production wells and cannot be ignored, so that the capacitance resistance method has certain irrationality in judging the communication condition among the injection and production wells of the high water content oil field.
Disclosure of Invention
The invention aims to provide a method for judging the communication condition between injection wells and production wells, so as to judge the communication condition between injection wells and production wells of a high water-cut oil field more reasonably.
The invention provides a method for judging the communication condition between injection wells and production wells, which is characterized by comprising the following steps of:
s1: data points for which data optimization of the oilfield block to be determined is not reasonable, the data of the oilfield block to be determined comprising: the number of injection wells, the number of production wells, the water injection amount of each injection well in each time step, the liquid production amount of each production well in each time step, the average permeability among injection and production wells and the injection and production well spacing;
s2: according to the optimized data of S1, calculating an initial value of the injection-production interwell communication coefficient and an initial value of a time lag constant of the oil field block to be determined;
s3: substituting the optimized data of S1 and the initial value of the communication coefficient between injection wells and the production well calculated in S2 into an objective function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the objective function value is smaller than the fitting error; the objective function is a function for calculating the minimum value of the superposition value of the sum of squares of the difference values of the calculated liquid production rate and the actual liquid production rate of each production well in the oil field block to be determined in all time steps;
s4: judging the strength of the communication condition between the injection wells and the production wells in the oil field block to be determined by utilizing the communication coefficient between the injection wells and the production wells calculated in the S3;
judging the shortest time for the injected water in the oil field block to reach each production well to be determined by utilizing the time lag constant calculated in the S3;
and judging that each production well in the oil field block to be determined is influenced by the interference of other production wells by using the interference constant among the production wells calculated by the S3.
The method for judging the communication condition between the injection wells and the production wells further adopts the following scheme:
the data of the oilfield block to be determined in S1 further includes: the included angle between the connecting line between the injection wells and the source direction and the average permeability between the injection wells and the production wells in the X, Y direction;
the step S4 is implemented by the following steps:
a. judging the initial value of the communication coefficient between the injection wells and the production wells and the initial value of the time lag constant calculated in the S2 through constraint conditions, and rejecting unreasonable judging data; the constraint conditions comprise an injection-production inter-well communication coefficient constraint condition, a time lag constant constraint condition and a production inter-well interference constant constraint condition and are used for constraining an injection-production inter-well communication coefficient, a time lag constant and a production inter-well interference constant;
b. substituting the optimized data of S1 and the initial value of the communication coefficient between injection wells and the initial value of the time lag constant after the judgment of a into an objective function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the objective function value is smaller than a preset value;
c. and c, judging the communication coefficient among the injection wells, the time lag constant and the interference constant among the production wells obtained by calculation in the step b through the constraint condition, eliminating unreasonable judgment data, and obtaining the communication coefficient among the injection wells, the time lag constant and the interference constant among the production wells which are finally used for judging the communication condition among the injection wells.
The data points for which the S1 optimization is not reasonable include:
setting the communication coefficient and the time lag constant between injection wells and production wells of the invalid well to be 0;
and setting the communication coefficient of the injection-production conversion well to be 0.
The step of calculating the initial value of the communication coefficient between the injection wells and the production wells of the oil field block to be determined in the step of S2 includes:
calculating to obtain an initial value of the communication coefficient between the injection wells and the production wells by using the following formula:
<math> <mrow> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mover> <mi>k</mi> <mo>&OverBar;</mo> </mover> <mi>ij</mi> </msub> <mo>/</mo> <msub> <mi>L</mi> <mi>ij</mi> </msub> </mrow> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </msubsup> <msub> <mover> <mi>k</mi> <mo>&OverBar;</mo> </mover> <mi>ij</mi> </msub> <mo>/</mo> <msub> <mi>L</mi> <mi>ij</mi> </msub> </mrow> </mfrac> </mrow> </math>
wherein f isijRepresenting the communication coefficient between the injection well of the ith hole and the production well of the jth hole;
Figure BDA0000436033460000032
the average permeability between the injection well of the ith and the production well of the jth is expressed by 10-3Square micron; l isijThe well spacing of the ith injection well and the jth production well is expressed in meters; n is a radical ofproRepresenting the number of production wells in ports; wherein N ispro
Figure BDA0000436033460000033
LijIs a known value.
The initial value of the time lag constant is the average time lag constant of the oilfield block to be determined;
the calculation step of the average time lag constant of the oilfield block to be determined comprises the following steps:
substituting the optimized data of S1 into an average time lag constant formula and calculating the average time lag constant when the formula value of the average time lag constant is smaller than a preset value;
the average time lag constant is expressed as:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>[</mo> <mi>q</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>q</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k</mi> <mi>&tau;</mi> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mo>[</mo> <mrow> <mo>(</mo> <msub> <mi>e</mi> <mi>w</mi> </msub> <mo>+</mo> <mi>fI</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mi>n</mi> <mo>-</mo> <mi>k</mi> </mrow> <mi>&tau;</mi> </mfrac> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mi>&tau;</mi> </mfrac> </mrow> </msup> <mo>)</mo> </mrow> <mo>]</mo> <mo>]</mo> </mrow> <mn>2</mn> </msup> </mrow> </math>
wherein n represents the total calculation time step in months; i (k) represents the water injection quantity of the block at the k time step, and the unit is cubic meter per day; q (k) represents the oil production of the block at the kth time step in cubic meters per day; q (0) represents the overall initial production of the block in cubic meters per day; f represents the block overall connectivity coefficient; τ represents the mean time lag constant in days; e.g. of the typewRepresents external energy replenishment in joules; wherein I (k), q (k) are known amounts.
The fitting error is 0.1.
The objective function in S3 is:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </munderover> <msup> <mrow> <mo>{</mo> <msub> <mi>q</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mo>{</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>l</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>[</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> </mrow> <mn>2</mn> </msup> </mrow> </math>
wherein q isj(k) Representing the actual liquid production amount of the jth production well in k time step, wherein the unit is cubic meter per day;
Figure BDA0000436033460000042
representing the initial yield of the jth production well in cubic meters per day; Δ t represents a time step in months; f. ofijRepresenting the communication coefficient of the ith injection well and the jth production well; alpha is alphajRepresenting the interference constant of other production forming wells to the jth production well; i isi(p) represents the water injection quantity of the injection well at the ith hole in p time steps, and the unit is cubic meter per day; tau isjRepresents the time lag constant for the jth production well in days; n is a radical ofiIndicating the number of injection wells in units of openings; n is a radical ofproRepresenting the number of production wells in ports; n represents the total calculation time step in months; wherein q isj(k)、
Figure BDA0000436033460000043
Δt、Ii(p)、Ni、NproAnd n is a known amount.
The constraint conditions are as follows:
<math> <mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> <mo>&le;</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <msub> <mi>a</mi> <mi>ij</mi> </msub> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mo>&le;</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>=</mo> <mn>1,2</mn> <mo>,</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>,</mo> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0.145</mn> <mi>&Delta;t</mi> <mo>&le;</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>&le;</mo> <mn>1000</mn> <mi>&Delta;t</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>&GreaterEqual;</mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
wherein,
Figure BDA0000436033460000045
j=1,2,...,Nprois a constraint condition of communication coefficient between injection wells and production wells, aijIn order to influence the coefficients of the effects,
<math> <mrow> <msub> <mi>a</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>></mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>&le;</mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow> </math> said rijBy passing <math> <mrow> <msub> <mi>r</mi> <mi>ij</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>a</mi> <mrow> <mn>1</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>b</mi> <mrow> <mn>1</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>a</mi> <mrow> <mn>2</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>b</mi> <mrow> <mn>2</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>,</mo> </mrow> </math>
a 1 ij b 1 ij = a 2 ij b 2 ij = k xij k yij Calculating to obtain;
0.145Δt≤τjthe time lag constant constraint condition is not more than 1000 delta t;
αjthe constraint condition of the interference constant between production wells is more than or equal to 0;
wherein d isijThe distance between the ith injection well and the jth production well is expressed in meters; r isijThe limit influence distance between the ith injection well and the jth production well is expressed in meters; a is1ij、a2ijRespectively representing the major semi-axes of ellipses formed by taking the injection well as the production well as the center of the well pair formed by the ith injection well and the jth production well, and the major semi-axes are meters; b1ij、b2ijRespectively representing the minor semi-axes of ellipses formed by taking the injection well as the production well as the center of the well pair formed by the ith injection well and the jth production well, and the unit is meter; thetaijThe included angle between the direction of the source and the connecting line of the ith injection well and the jth production well is expressed in degrees; Δ t represents a time step in months; f. ofijRepresenting the communication coefficient of the ith injection well and the jth production well; alpha is alphajRepresenting the interference constant of other production forming wells to the jth production well; tau isj(ii) a Represents the time lag constant for the jth production well in days; n is a radical ofproRepresenting the total number of producing wells in units of ports; wherein, b1ij、b2ijSimulating one half of the comprehensively determined reasonable well spacing by adopting an oil reservoir engineering numerical value; thetaijAnd Δ t is a known quantity.
The invention also provides a device for judging the communication condition between injection wells and production wells, which is characterized by comprising
The data preprocessing module is used for optimizing unreasonable data points of the oilfield block to be determined, and the data of the oilfield block to be determined comprises: the number of injection wells, the number of production wells, the water injection amount of each injection well in each time step, the liquid production amount of each production well in each time step, the average permeability among injection and production wells and the injection and production well spacing;
the initial value calculation module is used for calculating an initial value of the communication coefficient between the injection wells and the production wells of the oil field block to be determined and an initial value of a time lag constant according to the data optimized by the data preprocessing module;
the result calculation module is used for substituting the data optimized by the data preprocessing module and the initial value of the communication coefficient between injection wells and the initial value of the time lag constant calculated by the initial value calculation module into an objective function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the objective function value is smaller than the fitting error; the objective function is a function for calculating the minimum value of the superposition value of the sum of squares of the difference values of the calculated liquid production rate and the actual liquid production rate of each production well in the oil field block to be determined in all time steps;
and the communication condition judging module is used for judging the strength of the communication condition between the injection wells and the production wells in the oil field block to be determined by utilizing the communication coefficient between the injection wells and the production wells calculated in the S3, judging the shortest time for the injection water to reach each production well by utilizing the time lag constant calculated in the S3 and judging the interference influence of each production well on other production wells by utilizing the interference constant between the production wells calculated in the S3.
The device for judging the communication condition between the injection and production wells further adopts the following scheme:
the data of the oilfield block to be determined in the data preprocessing module further comprises an included angle between a connecting line between injection wells and a source direction and an average permeability between the injection wells and the production wells in the X, Y direction;
the result calculation module further comprises the following modules:
the constraint module is used for judging the initial value of the communication coefficient between the injection wells and the production wells and the initial value of the time lag constant calculated in the introduced initial value calculation module through constraint conditions and rejecting unreasonable judgment data; the constraint conditions comprise an injection-production inter-well communication coefficient constraint condition, a time lag constant constraint condition and a production inter-well interference constant constraint condition and are used for constraining an injection-production inter-well communication coefficient, a time lag constant and a production inter-well interference constant;
the calculation module is used for substituting the data optimized by the data preprocessing module and the initial value of the communication coefficient between injection wells and the initial value of the time lag constant judged by the constraint module into a target function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the target function value is smaller than the fitting error;
and the judging module is used for judging the communication coefficient among the injection wells, the time lag constant and the interference constant among the production wells which are calculated in the calculating module through the constraint condition, eliminating unreasonable judging data and obtaining the communication coefficient among the injection wells, the time lag constant and the interference constant among the production wells which are finally used for judging the communication condition among the injection wells.
The invention has the beneficial effects that:
1. the method considers the influence of the production well-to-well interference on the communication condition between the injection wells and the production wells, and finally expresses the communication condition between the injection wells and the production wells through the communication coefficient between the injection wells and the production well, the time lag constant and the interference constant between the production wells, so that the rationality of judging the communication condition between the injection wells and the production wells of the high-water-content oil field is enhanced.
2. According to the method, constraint conditions are introduced, the calculated communication coefficient between injection wells and production wells, the time lag constant and the interference constant between production wells are judged, unreasonable judgment data are eliminated, the judgment data are made to be more in line with geological understanding of an oil field block to be determined, the rationality of a target function is further improved, and finally the communication condition between injection wells and production wells of the high-water-content oil field can be clearly expressed through the optimal solution of the calculated target function so as to guide the fine injection and production of the oil field.
Drawings
FIG. 1 is a schematic view of a one-injection-one-mining limit influence distance;
FIG. 2 is a flowchart of a method for determining a communication status between injection wells and production wells according to a first embodiment of the present invention;
FIG. 3 is a flowchart of step S4 of a method for determining the communication status between injection wells and production wells according to a second embodiment of the present invention;
FIG. 4 is a schematic view of a first embodiment of a communication status determination device for injection wells and production wells according to the present invention;
FIG. 5 is a schematic diagram of a result calculation module according to a second embodiment of the apparatus for determining the communication status between injection wells and production wells of the present invention;
FIG. 6 is a diagram of the distribution of the initial injection-production inter-well communication coefficients;
FIG. 7 is a distribution diagram of the communication coefficient between injection wells and production wells.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. It is to be understood that the embodiments described are only a few embodiments of the present application and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 2, a first embodiment of the method for determining the communication status between injection wells and production wells in the present invention includes the following steps:
s1: data points for which data optimization of the oilfield block to be determined is not reasonable, the data of the oilfield block to be determined comprising: the number of injection wells, the number of production wells, the water injection amount of each injection well in each time step, the liquid production amount of each production well in each time step, the average permeability among injection and production wells and the injection and production well spacing;
the data points for which the S1 optimization is not reasonable include:
setting the communication coefficient and the time lag constant between injection wells and production wells of the invalid well to be 0; the invalid well is a production stopping well, and the production data is less than the average single well production data of the block;
setting the communication coefficient of the injection-production conversion well to be 0;
s2: according to the optimized data of S1, calculating an initial value of the injection-production interwell communication coefficient and an initial value of a time lag constant of the oil field block to be determined;
calculating to obtain an initial value of the communication coefficient between the injection wells and the production wells by using the following formula: :
<math> <mrow> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mover> <msub> <mi>k</mi> <mi>ij</mi> </msub> <mo>&OverBar;</mo> </mover> <mo>/</mo> <msub> <mi>L</mi> <mi>ij</mi> </msub> </mrow> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </msubsup> <mover> <msub> <mi>k</mi> <mi>ij</mi> </msub> <mo>&OverBar;</mo> </mover> <mo>/</mo> <msub> <mi>L</mi> <mi>ij</mi> </msub> </mrow> </mfrac> </mrow> </math>
wherein f isijEighthly, representing the communication coefficient between the injection well at the ith hole and the production well at the jth hole;
Figure BDA0000436033460000072
the average permeability between the injection well of the ith and the production well of the jth is expressed by 10-3Square micron; l isijThe well spacing of the ith injection well and the jth production well is expressed in meters; n is a radical ofproRepresenting the number of production wells in ports; wherein,
Figure BDA0000436033460000073
Lijare known;
the initial value of the time lag constant adopts the average time lag constant of the oilfield block to be determined;
the calculation step of the average time lag constant of the oilfield block to be determined comprises the following steps:
substituting the optimized data of S1 into an average time lag constant formula and calculating the average time lag constant when the formula value of the average time lag constant is smaller than a preset value;
the average time lag constant is expressed as:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>[</mo> <mi>q</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>q</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k</mi> <mi>&tau;</mi> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mo>[</mo> <mrow> <mo>(</mo> <msub> <mi>e</mi> <mi>w</mi> </msub> <mo>+</mo> <mi>fI</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mi>n</mi> <mo>-</mo> <mi>k</mi> </mrow> <mi>&tau;</mi> </mfrac> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mi>&tau;</mi> </mfrac> </mrow> </msup> <mo>)</mo> </mrow> <mo>]</mo> <mo>]</mo> </mrow> <mn>2</mn> </msup> </mrow> </math>
wherein n represents the total calculation time step in months; i (k) represents the water injection quantity of the block at the k time step, and the unit is cubic meter per day; q (k) represents the oil production of the block at the kth time step in cubic meters per day; q (0) represents the overall initial production of the block in cubic meters per day; f represents the block overall connectivity coefficient; τ represents the mean time lag constant in days; e.g. of the typewRepresents external energy replenishment in joules; wherein I (k), q (k) are known amounts
S3: substituting the optimized data of S1 and the initial value of the communication coefficient between injection wells and the production well calculated in S2 into an objective function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the objective function value is smaller than the fitting error; the objective function is a function for calculating the minimum value of the superposition value of the sum of squares of the difference values of the calculated liquid production rate and the actual liquid production rate of each production well in the oil field block to be determined in all time steps;
the objective function is:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </munderover> <mo>{</mo> <msub> <mi>q</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mo>{</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>[</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> </mrow> </math>
wherein q isj(k) Represents the actual fluid production at k time step in cubic meters per day;representing the initial production of the production well in cubic meters per day; Δ t represents the time step, month; f. ofijRepresenting the communication coefficient of the ith injection well and the jth production well; alpha is alphajRepresenting the interference constant of other production forming wells to the jth production well; i isi(p) represents the water injection quantity of the injection well at the ith hole in p time steps, and the unit is cubic meter per day; tau isjRepresents the time lag constant for the jth production well in days; n is a radical ofiIndicating the number of injection wells in units of openings; n is a radical ofproRepresenting the number of production wells in ports; n represents the total calculation time step in months; wherein q isj(k)、
Figure BDA0000436033460000083
Δt、Ii(p)、Ni、NproAnd n is a known amount.
S4: judging the strength of the communication condition between the injection wells and the production wells in the oil field block to be determined by utilizing the communication coefficient between the injection wells and the production wells calculated in the S3;
judging the shortest time for the injected water in the oil field block to reach each production well to be determined by utilizing the time lag constant calculated in the S3;
and judging that each production well in the oil field block to be determined is influenced by the interference of other production wells by using the interference constant among the production wells calculated by the S3.
In the embodiment, the influence of the interference between the production wells on the communication condition between the injection wells and the production wells is considered, the interference constant between the production wells is introduced into the establishment of the objective function, and the communication condition between the injection wells and the production wells is judged by combining the communication coefficient between the injection wells and the production wells, the time lag constant and the interference constant between the production wells, so that the rationality of judging the communication condition between the injection wells and the production wells of the high-water-content oil field is enhanced, and the production work of the oil field is guided more favorably.
Steps S2 and S4 in the second embodiment of the method for determining the communication status between injection wells and production wells according to the present invention are the same as those in the first embodiment.
Step S1 in the second embodiment of the present invention is: data points for which data optimization of the oilfield block to be determined is not reasonable, the data of the oilfield block to be determined comprising: the number of injection wells, the number of production wells, the water injection amount of each injection well at each time step, the liquid production amount of each production well at each time step, the average permeability among injection and production wells, the injection and production well spacing, the included angle between the connecting line among the injection and production wells and the material source direction and the average permeability among the injection and production wells in the X, Y direction;
fig. 3 shows the specific steps included in S3 in the second embodiment. As shown in fig. 3, the S3 is implemented by the following three steps:
a. judging the initial value of the communication coefficient between the injection wells and the production wells and the initial value of the time lag constant calculated in the S2 through constraint conditions, and rejecting unreasonable judging data; the constraint conditions comprise an injection-production inter-well communication coefficient constraint condition, a time lag constant constraint condition and a production inter-well interference constant constraint condition and are used for constraining an injection-production inter-well communication coefficient, a time lag constant and a production inter-well interference constant;
the constraint conditions are as follows:
<math> <mrow> <mi>s</mi> <mo></mo> <mo>.</mo> <mi>t</mi> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> <mo>&le;</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <msub> <mi>a</mi> <mi>ij</mi> </msub> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mo>&le;</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>=</mo> <mn>1,2</mn> <mo>,</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>,</mo> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0.145</mn> <mi>&Delta;t</mi> <mo>&le;</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>&le;</mo> <mn>1000</mn> <mi>&Delta;t</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>&GreaterEqual;</mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
wherein,
Figure BDA0000436033460000092
j=1,2,...,Nprois a constraint condition of communication coefficient between injection wells and production wells, aijIn order to influence the coefficients of the effects,
<math> <mrow> <msub> <mi>a</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>></mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>&le;</mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow> </math> rijby passing <math> <mrow> <msub> <mi>r</mi> <mi>ij</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>a</mi> <mrow> <mn>1</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>b</mi> <mrow> <mn>1</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>a</mi> <mrow> <mn>2</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>b</mi> <mrow> <mn>2</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>,</mo> </mrow> </math> a 1 ij b 1 ij = a 2 ij b 2 ij = k xij k yij Calculating to obtain;
0.145Δt≤τjthe time lag constant constraint condition is not more than 1000 delta t;
αjthe constraint condition of the interference constant between production wells is more than or equal to 0;
wherein d isijThe distance between the ith injection well and the jth production well is expressed in meters; r isijThe limit influence distance between the ith injection well and the jth production well is expressed in meters; a is1ij、a2ijRespectively representing the major semi-axes of ellipses formed by taking the injection well as the production well as the center of the well pair formed by the ith injection well and the jth production well, and the major semi-axes are meters; b1ij、b2ijRespectively representing the minor semi-axes of ellipses formed by taking the injection well as the production well as the center of the well pair formed by the ith injection well and the jth production well, and the unit is meter; thetaijThe included angle between the direction of the source and the connecting line of the ith injection well and the jth production well is expressed in degrees; Δ t represents a time step in months; f. ofijRepresenting the communication coefficient of the ith injection well and the jth production well; alpha is alphajRepresenting the interference constant of other production forming wells to the jth production well; tau isj(ii) a Represents the time lag constant for the jth production well in days; n is a radical ofproRepresenting the total number of producing wells in units of ports; wherein, b1ij、b2ijSimulating one half of the comprehensively determined reasonable well spacing by adopting an oil reservoir engineering numerical value; thetaijΔ t is a known quantity;
b. substituting the optimized data of S1 and the initial value of the communication coefficient between injection wells and the initial value of the time lag constant after the judgment of a into an objective function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the objective function value is smaller than the fitting error;
the objective function is:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </munderover> <mo>{</mo> <msub> <mi>q</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mo>{</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <msub> <mi>a</mi> <mi>j</mi> </msub> <mo>[</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> </mrow> </math>
wherein q isj(k) Represents the actual fluid production at k time step in cubic meters per day;
Figure BDA0000436033460000101
representing the initial production of the production well in cubic meters per day; Δ t represents the time step, month; f. ofijRepresenting the communication coefficient of the ith injection well and the jth production well; alpha is alphajRepresenting the interference constant of other production forming wells to the jth production well; i isi(p) represents the water injection quantity of the injection well at the ith hole in p time steps, and the unit is cubic meter per day; tau isjRepresents the time lag constant for the jth production well in days; n is a radical ofiIndicating the number of injection wells in units of openings; n is a radical ofproRepresenting the number of production wells in ports; n represents the total calculation time step in months; wherein q isj(k)、
Figure BDA0000436033460000102
Δt、Ii(p)、Ni、NproAnd n is a known amount.
And the fitting error adopts a communication coefficient among injection wells and production wells, a time lag constant and an interference constant among production wells as an optimal solution of a target function when the minf (x) is less than or equal to 0.1.
And when the calculated objective function value is smaller than the fitting error, the injection-production inter-well communication coefficient, the time lag constant and the production inter-well interference constant adopt a sequential quadratic programming method.
c. And c, judging the communication coefficient among the injection wells, the time lag constant and the interference constant among the production wells obtained by calculation in the step b through the constraint condition, eliminating unreasonable judgment data, and obtaining the communication coefficient among the injection wells, the time lag constant and the interference constant among the production wells which are finally used for judging the communication condition among the injection wells.
In the second embodiment of the invention, the influence of the interference between the production wells on the communication condition between the injection wells is considered, the interference constant between the production wells is introduced into the establishment of the objective function, and the communication condition between the injection wells and the production wells is judged by combining the communication coefficient between the injection wells and the production wells, the time lag constant and the interference constant between the production wells, so that the rationality for judging the communication condition between the injection wells and the production wells of the high-water-content oil field is greatly enhanced, and the production work of the high-water-content oil field is guided more favorably.
Considering that the injection-production communication relationship is greatly influenced by geological factors, and the injection-production communication coefficient, the time lag constant and the production inter-well interference constant are consistent with geological knowledge, for example, a relatively large communication coefficient should not exist between injection-production wells with relatively long distances, the production inter-well interference constant should not be a negative value and the like, the second embodiment of the invention introduces a constraint condition into the established objective function, discriminates the calculated injection-production inter-well communication coefficient, the time lag constant and the production inter-well interference constant, rejects unreasonable discrimination data, further improves the rationality of the calculated discrimination data, finally achieves the purpose of secondarily improving the rationality of discrimination of the injection-production inter-well communication condition of the high water-cut oil field, and is more favorable for guiding the fine injection-production of the high water-cut oil field.
The second embodiment of the method for determining the communication status between injection wells and production wells according to the present invention is further described below by taking the application of the present invention to a certain block of a high water content oil field as an example.
At present, a certain block of a certain high-water-cut oil field is developed for 20 years, the water content reaches more than 80%, after the injection and production well pattern of the oil field is adjusted for many times, the water drive control degree (only 55%) and the water drive use degree (20% -35%) are still low, the flow direction of injected water is unclear, the injection and production relation is incomplete, fine water injection regulation and control are urgently needed to be implemented, and the communication condition between injection and production wells is judged to be the premise of implementing the fine water injection regulation and control.
The total number of the wells in a certain block of the high water-cut oil field is 28, wherein: 17 production wells, 5 injection wells and 6 transfer wells, the production starting time is 1 month in 1992, the data acquisition cutoff calculation time is 3 months in 2012, and the time step is 1 month. And calculating and counting the water injection amount of each injection well, the liquid production amount of each production well, the included angle between the material source direction (the direction of sand body deposition in the oil reservoir) and the connecting line between the injection and production wells, the distance between the injection and production wells, the average permeability in the X, Y direction between the injection and production wells, the average permeability between the injection and production wells and the like according to the production data and the geological static data of each well.
Aiming at the communication condition among injection and production wells of a certain block of the high water-cut oil field, the invention comprises the following steps:
s1: data points of the oilfield block to be determined are optimized unreasonably;
before judging the communication condition between injection wells and production wells, firstly preprocessing statistical data and the like, optimizing unreasonable data points and improving the fitting effect between the actual liquid production amount and the calculated liquid production amount.
Optimizing unreasonable data points includes:
setting the communication coefficient and the time lag constant between injection wells and production wells of the invalid well to be 0;
setting the communication coefficient of the injection-production conversion well to be 0;
and obtaining the required data through the steps. Due to the large amount of data required, only a part of them, i.e., the data shown in tables 1 to 3, will be listed below. The statistical data are as follows:
table 1 injection well parameter data table
Figure BDA0000436033460000111
TABLE 2 production well parameter data sheet
Figure BDA0000436033460000121
TABLE 3 parameter data sheet between injection wells and production wells
Figure BDA0000436033460000122
S2: according to the optimized data of S1, calculating an initial value of the injection-production interwell communication coefficient and an initial value of a time lag constant of the oil field block to be determined;
firstly, calculating the communication coefficient f between injection wells and production wellsijIs started. According to the equilibrium state multivariable linear regression model, assuming that the flow reaches a steady state in the acquired time step, and obtaining the communication coefficient between the injection well at the ith port and the production well at the jth port as follows:
<math> <mrow> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mi>ij</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mi>wfi</mi> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>wfj</mi> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </munderover> <msub> <mi>T</mi> <mi>ij</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mi>wfi</mi> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>wfj</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein, TijThe conductivity coefficient between the injection well of the ith hole and the production well of the jth hole is expressed by 10-3Square micron meter/(millipascal second); p is a radical ofwfi、pwfjRespectively represent the i-th port injectionThe bottom hole flow pressure of the well and the j production well is measured in millipascals.
And (3) obtaining the conductivity coefficient between the injection well of the ith opening and the production well of the jth opening according to a steady-state linear flow equation:
<math> <mrow> <msub> <mi>T</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mover> <mi>k</mi> <mo>&OverBar;</mo> </mover> <mi>ij</mi> </msub> <msub> <mi>A</mi> <mi>ij</mi> </msub> </mrow> <mrow> <mi>&mu;</mi> <msub> <mi>L</mi> <mi>ij</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein,
Figure BDA0000436033460000132
the average permeability between the injection well of the ith and the production well of the jth is expressed by 10-3Square micron; a. theijRepresents the cross-sectional area of flow in square meters; l isijThe well spacing of the ith injection well and the jth production well is expressed in meters.
Substituting the formula (2) into the formula (1), assuming that the bottom flowing pressure of the production well is the same, and assuming that the flowing sectional area of each injection well to the surrounding production well is the same, obtaining a communication coefficient initial value calculation formula:
<math> <mrow> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mover> <mi>k</mi> <mo>&OverBar;</mo> </mover> <mi>ij</mi> </msub> <mo>/</mo> <msub> <mi>L</mi> <mi>ij</mi> </msub> </mrow> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </munderover> <msub> <mover> <mi>k</mi> <mo>&OverBar;</mo> </mover> <mi>ij</mi> </msub> <mo>/</mo> <msub> <mi>L</mi> <mi>ij</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </math>
it can be seen from the formula (3) that the initial value of the communication coefficient can be obtained by substituting the average permeability between injection wells and the distance between injection wells and production wells into the formula (3).
The initial value of the time lag constant adopts the average time lag constant of the oilfield block to be determined.
The average time lag constant of the oilfield block to be determined is obtained by substituting the optimized data of S1 into an average time lag constant formula and calculating the average time lag constant when the value of the average time lag constant formula is smaller than a preset value;
the average time lag constant is expressed as:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>[</mo> <mi>q</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>q</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k</mi> <mi>&tau;</mi> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mo>[</mo> <mrow> <mo>(</mo> <msub> <mi>e</mi> <mi>w</mi> </msub> <mo>+</mo> <mi>fI</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>n</mi> <mo>-</mo> <mi>k</mi> </mrow> <mi>&tau;</mi> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mi>&tau;</mi> </mfrac> </mrow> </msup> <mo>)</mo> </mrow> <mo>]</mo> <mo>]</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein n represents the total calculation time step in months; i (k) represents the water injection quantity of the block at the k time step, and the unit is cubic meter per day; q (k) represents the oil production of the block at the kth time step in cubic meters per day; q (0) represents the overall initial production of the block in cubic meters per day; f represents the block overall connectivity coefficient; τ represents the mean time lag constant in days; e.g. of the typewRepresents external energy replenishment in joules; wherein I (k), q (k) are known amounts.
And (4) substituting the optimized data of S1 into the formula (4) to calculate the average time lag constant tau of the oilfield block to be determined to be 322 days by a sequential quadratic programming method.
S3: substituting the optimized data of S1 and the initial value of the communication coefficient between injection wells and the production well calculated in S2 into an objective function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the objective function value is smaller than the fitting error; the objective function is a function for calculating the minimum value of the superposition value of the sum of squares of the difference values of the calculated liquid production rate and the actual liquid production rate of each production well in the oil field block to be determined in all time steps;
the S3 is realized by the following three steps:
a. judging the initial value of the communication coefficient between the injection wells and the production wells and the initial value of the time lag constant calculated in the S2 through constraint conditions, and rejecting unreasonable judging data; the constraint conditions comprise an injection-production inter-well communication coefficient constraint condition, a time lag constant constraint condition and a production inter-well interference constant constraint condition and are used for constraining an injection-production inter-well communication coefficient, a time lag constant and a production inter-well interference constant;
the communication condition of injection and production communication is greatly influenced by geological factors, and the communication coefficient, the time lag constant and the interference constant between production wells conform to geological knowledge, for example, the communication coefficient between injection and production wells which are far away should not be larger, the interference constant between production wells should not be negative, and the like. In order to make the obtained data conform to the geological knowledge of the oilfield block to be determined and further improve the rationality of the judgment data, constraint conditions need to be established to judge the rationality of the obtained data.
The established constraint conditions comprise an injection-production inter-well communication coefficient constraint condition, a time lag constant constraint condition and a production inter-well interference constant constraint condition.
Firstly, establishing a communication coefficient constraint condition between injection wells and production wells;
the method for establishing the communication coefficient constraint condition between the injection wells and the production wells comprises the following steps:
A. acquiring a limit influence distance formula between injection wells and production wells;
and the limit influence distance between the injection wells and the production wells is the corresponding well distance when the elliptic isobars formed in the limit wave area are tangent. In the limit wave area of the ellipse formed by one well, the direction with larger radius of the ellipse is along the direction of the object source, and the direction with smaller radius is vertical to the direction of the object source. If the well spacing exceeds the limit influence distance, the communication relation between the two wells is considered to be weak, and the communication coefficient between the injection wells and the production wells is 0.
Fig. 1 is a schematic diagram of the injection-production limit influence distance, and the injection-production limit influence distance shown in fig. 1 is calculated by taking an injection-production in an oil field block to be determined as an example. The upper left ellipse represents the limiting swept area of the injection well, with the major semi-axis a1The minor semi-axis is b1(ii) a The lower right ellipse represents the limiting swept area of the producing well, with the semi-major axis a2The minor semi-axis is b2(ii) a Theta is an included angle between the direction of the source and a connecting line between the injection well and the production well; r isijThe method is characterized in that the method is a connecting line between two points in the graph, namely the limit influence distance between injection wells and production wells, and the method is calculated according to an elliptic formula:
<math> <mrow> <msub> <mi>r</mi> <mi>ij</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mi>coa</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>a</mi> <mrow> <mn>1</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>b</mi> <mrow> <mn>1</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>a</mi> <mrow> <mn>2</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>b</mi> <mrow> <mn>2</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein r isijThe limit influence distance between the ith injection well and the jth production well is represented by meter; a is1ij、a2ijRespectively showing the major semi-axis and the minor semi-axis of an ellipse formed by taking the injection well as the production well as the center of the ellipse in the well pair formed by the ith injection well and the jth production well; b1ij、b2ijThe minor semi-axis and the meter of an ellipse formed by taking the injection well as the production well as the center of the well pair formed by the ith injection well and the jth production well respectively are shown; thetaijThe included angle between the direction of the source and the connecting line of the injection well of the ith opening and the production well of the jth opening is expressed as degree; wherein, b1ij、b2ijOne half of the reasonable well spacing, theta, which can be comprehensively determined by numerical simulation of oil reservoir engineeringijIs a known amount;
B. calculating a in step A1ij、a2ij
According to the actual condition of an oil reservoir, the stratum permeability is a main factor influencing the limiting radius between injection wells and production wells, and is deduced according to an anisotropic stratum elliptic flow equation, and the following relations exist:
a 1 ij b 1 ij = a 2 ij b 2 ij = k xij k yij - - - ( 6 )
wherein k isxijAnd kyijAverage permeability in the direction of X, Y for the ith injection well and the jth production well, respectively, 10-3Square micron; a is1ij、a2ijThe major semi-axis and the minor semi-axis of an ellipse formed by taking the injection well as the center and taking the production well as the center in a well pair formed by the ith injection well and the jth production well respectively are meters; b1ij、b2ijA minor semi-axis, m, of an ellipse formed by centering the injection well and centering the production well in a well pair formed by the injection well and the production well of the ith port and the jth production well respectively; wherein k isxij、kyij、b1ij、b2ijIn known amounts.
C. Establishing a limit influence distance matrix:
<math> <mrow> <mi>R</mi> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>r</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>r</mi> <mn>12</mn> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>r</mi> <mrow> <mn>1</mn> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>r</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>r</mi> <mn>22</mn> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>r</mi> <mrow> <mn>2</mn> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>r</mi> <mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>r</mi> <mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>r</mi> <mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> <msub> <mo>)</mo> <mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> <mo>&times;</mo> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mrow> </msub> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein r isijThe limit influence distance between the ith injection well and the jth production well is represented by meter; n is a radical ofiIndicating the number of injection wells, the mouth; n is a radical ofproIndicates the number of production wells, ports; wherein N isi、NproIs a known amount, rijCalculated by the step A, B;
D. establishing an injection-production well spacing matrix:
<math> <mrow> <mi>D</mi> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>d</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>d</mi> <mn>12</mn> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>d</mi> <mrow> <mn>1</mn> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>d</mi> <mn>12</mn> </msub> </mtd> <mtd> <msub> <mi>d</mi> <mn>22</mn> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>d</mi> <mrow> <mn>2</mn> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>d</mi> <mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>d</mi> <mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>d</mi> <mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>)</mo> </mrow> <mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein d isijRepresenting the distance between the ith injection well and the jth production well in meters; n is a radical ofiIndicating the number of injection wells, the mouth; n is a radical ofproRepresenting production wellsNumber, mouth; wherein d isij、Ni、NproIn known amounts.
E. Establishing a calculation influence coefficient matrix:
<math> <mrow> <mi>A</mi> <mo>=</mo> <msub> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mi>ij</mi> </msub> <mo>)</mo> </mrow> <mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> <mo>&times;</mo> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein, <math> <mrow> <msub> <mi>a</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>></mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>&le;</mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> </math> dij、rijNiindicating the number of injection wells, the mouth; n is a radical ofproIndicating the number of production wells, ports.
F. Establishing a communication coefficient constraint condition between injection wells and production wells;
by a in formula (9)ijRestricting the communication coefficient between injection wells and production wells, and establishing the restriction condition of the communication coefficient between injection wells and production wells, namely:
<math> <mrow> <mn>0</mn> <mo>&le;</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <msub> <mi>a</mi> <mi>ij</mi> </msub> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mo>&le;</mo> <mn>1</mn> </mrow> </math> j=1,2,...,Npro <math> <mrow> <msub> <mi>a</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>></mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>&le;</mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </math>
secondly, establishing a time lag constant constraint condition;
to ensure computation of weight terms in an objective function
Figure BDA0000436033460000166
And
Figure BDA0000436033460000167
in the (0, 1) range, the invention derives a reasonable range of time lag constants: 0.145 Deltat is less than or equal to tauj1000 Δ t, the range of time lag constants being related to the time step Δ t.
Finally, establishing constraint conditions of interference constants among production wells;
since the interference between production wells cannot be negative, the interference constant of other production forming wells on the jth production well cannot be less than 0, i.e., αj≥0。
The constraint condition is finally obtained as follows:
<math> <mrow> <mi>s</mi> <mo>.</mo> <mo></mo> <mi>t</mi> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> <mo>&le;</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <msub> <mi>a</mi> <mi>ij</mi> </msub> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mo>&le;</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>=</mo> <mn>1,2</mn> <mo>,</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>,</mo> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0.145</mn> <mi>&Delta;t</mi> <mo>&le;</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>&le;</mo> <mn>1000</mn> <mi>&Delta;t</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>&GreaterEqual;</mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow> </math>
substituting the optimized data of S1 into a formula (3) to obtain a result, and then judging the result according to a formula (11) of a communication coefficient constraint condition to obtain an initial value matrix of the communication coefficient between injection wells, wherein the initial value matrix is as follows:
TABLE 4 initial value matrix of communication coefficient between injection wells and production wells
Figure BDA0000436033460000165
Figure BDA0000436033460000171
The initial value of the time lag constant is equal to or more than 0.145 and equal to or more than 322 and equal to or more than 1000, and the initial value of the time lag constant is 322 days because the constraint condition is met.
As shown in fig. 6, the injection and production wells are connected on the injection and production well distribution map by lines with different thickness degrees according to the data in table 4, and the larger the initial connectivity coefficient is, the thicker the lines are, and the initial injection and production inter-well connectivity coefficient distribution map is further drawn.
b. Substituting the optimized data of S1 and the initial value of the communication coefficient between injection wells and the initial value of the time lag constant after the judgment of a into an objective function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the objective function value is smaller than the fitting error;
when the calculated fluid production capacity and the actual fluid production capacity are well fitted, the sum of squares of the difference values of the calculated fluid production capacity and the actual fluid production capacity approaches to 0, so that the objective function adopts a function for calculating the minimum value of the superposition value of the sum of squares of the difference values of the calculated fluid production capacity and the actual fluid production capacity of each production well in the oil field block to be determined in all time steps. Because the calculated liquid production amount and the actual liquid production amount cannot be completely fitted, a certain fitting error exists, and the fitting error is 0.1 in the process of distinguishing the high-water-content oil field block.
The objective function is the following equation:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </munderover> <msup> <mrow> <mo>[</mo> <msub> <mi>q</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>]</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein q isj(k) Representing the actual liquid production amount of the jth production well at the K time step, wherein the cubic meter is per day;representing the calculated liquid production amount of the jth production well at the K time step, wherein the cubic meter is per day; n represents the total calculation time step, month; n is a radical ofproRepresenting the total number of producing wells, ports.
The high water content oil field has the characteristics of dense injection and production well pattern and serious interference between production wells, so that the interference constant alpha between production wells is introduced when calculating the liquid production amountjAnd represents the interference of other production forming wells to the jth production well.
Introduction of interference constant alpha between production wellsjThe post-calculated fluid production is expressed as:
<math> <mrow> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein,representing the calculated liquid production capacity of the jth production well on the k time step, cubic meters per day; Δ t represents the time step, month; f. ofijRepresenting the communication coefficient of the ith injection well and the jth production well; alpha is alphajIndicating the interference of other production forming wells on the jth production well; tau isjRepresents the time lag constant, days, for the jth producer well; i isi(k) And the water injection quantity of the injection well at the ith hole at the k time step is shown, and the cubic meter is taken every day.
Thus, the fluid production at the kth time point of the jth production well is expressed as the combined effect of the fluid production at the previous time point, the water injection at the surrounding injection well, and the interference at the surrounding production well.
The recursive operation is performed on equation (6) as follows:
<math> <mrow> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> </mrow> </math>
<math> <mrow> <mo>=</mo> <mo>[</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> <mo>]</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> </mrow> </math>
<math> <mrow> <mo>=</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> </mrow> </math>
<math> <mrow> <mo>=</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>3</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>3</mn> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <msup> <mrow> <mo>-</mo> <mi>e</mi> </mrow> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> <mo>+</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> </mrow> </math>
<math> <mrow> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> </mrow> </math>
......
<math> <mrow> <mo>=</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <mrow> <mi></mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> <mo>+</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>n</mi> <mo>+</mo> <mn>2</mn> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> </mrow> </math>
<math> <mrow> <mo>+</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> </mrow> </math>
<math> <mrow> <mo>=</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>]</mo> </mrow> </math>
equation (6) becomes, through the above operation:
<math> <mrow> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mo>{</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </mrow> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>[</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </mrow> </msup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein,
Figure BDA00004360334600001810
the calculated liquid production amount of the jth production well on the k time step is expressed, and the unit is cubic meter per day;
Figure BDA00004360334600001811
representing the initial liquid production of a jth production well, and the unit is cubic meter per day; Δ t represents a time step in months; f. ofijRepresenting the communication coefficient of the ith injection well and the jth production well; alpha is alphajIndicating interference of wells formed by other production with the jth production well; i isi(p) represents the water injection quantity of the injection well at the ith hole in p time steps, and the unit is cubic meter per day; tau isjTime lag constant for the jth production well in days; n is a radical ofiIndicating the number of injection wells in units of openings.
Substituting the formula (14) into the formula (12) to obtain a judging objective function of the communication condition between the injection wells and the production wells:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </munderover> <msup> <mrow> <mo>{</mo> <msub> <mi>q</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mo>{</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>[</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow> </math>
and (3) substituting the data of the number of the injection wells and the production wells optimized by the S1, the total calculation time step, the water injection quantity of the injection wells, the liquid production quantity of the production wells, the injection-production well spacing and the like, and the initial values of the communication coefficient between the injection wells and the production wells and the initial values of the time lag constant in the table 4 into a formula (15), wherein the fitting error adopts 0.1, and then solving the communication coefficient between the injection wells and the production wells, the time lag constant and the interference constant between the production wells by using a sequential quadratic programming method.
c. And c, judging the communication coefficient among the injection wells, the time lag constant and the interference constant among the production wells obtained by calculation in the step b through the constraint condition, eliminating unreasonable judgment data, and obtaining the communication coefficient among the injection wells, the time lag constant and the interference constant among the production wells which are finally used for judging the communication condition among the injection wells.
And (c) judging the communication coefficient among injection wells, the time lag constant and the interference constant among production wells obtained by calculation in the step b by a constraint condition formula (11), and finally obtaining an injection-production inter-well communication condition matrix consisting of the communication coefficient among injection wells, the time lag constant and the interference constant among production wells, wherein the table is as follows:
TABLE 5 matrix of communication conditions between injection wells and production wells
Figure BDA0000436033460000191
S4: judging the strength of the communication condition between the injection wells and the production wells in the oil field block to be determined by utilizing the communication coefficient between the injection wells and the production wells calculated in the S3;
as shown in fig. 7, the data in table 5 are plotted on the injection and production well distribution map by connecting the injection and production wells by lines having different thickness degrees, and the larger the communication coefficient between the injection and production wells is, the thicker the lines are, so as to obtain the injection and production well communication coefficient distribution map shown in fig. 8.
The line thickness degree in FIG. 8 shows that there are communication conditions between G2 and the surrounding production wells G5-14, G28-24, G5-22 and G5-32, wherein the connectivity with G5-22 is strongest. The water well G28-28 has communication conditions with the surrounding water wells G28-26, G5-36, G5-34, G5-52, G32-28, G5-54, G32-30 and G5-58, wherein the communication conditions with a plurality of production wells G5-32, G5-34 and G32-28 are stronger. The G30-30 wells are communicated with the G5-54 wells, the G32-30 wells, the G5-58 wells and the G5-64 wells, wherein the G5-58 wells are communicated with the G5-58 wells. G5-1 is communicated with a plurality of wells such as G5-52, G32-28, G5-54, G32-30 and G5-64, wherein the communication with G5-52 is strong. The well G5 has communication conditions with G5-52, G32-30, G5-64, G34-32 and G36-32, wherein the communication conditions with G5-64 and G5-52 are stronger.
Judging the shortest time for the injected water to reach each production well by using the time lag constant calculated in the step S3;
the time for the injected water to reach each injection well fastest can be clearly analyzed through the time lag constant of each production well listed in table 5;
and judging that each production well is influenced by the interference of other production wells by using the interference constant among the production wells calculated by the S3.
The interference effect of each production well from other production wells can be clearly analyzed by the interference constants between production wells listed in table 5.
Fig. 4 is a schematic diagram of a first embodiment of a device for determining a communication status between injection wells and production wells according to the present invention. As shown in fig. 4, a first embodiment of the apparatus for determining the communication status between injection wells and production wells in the present invention includes the following modules:
a data preprocessing module P1 for data points for which data optimization of the oilfield block to be determined is not reasonable, the data of the oilfield block to be determined including: the number of injection wells, the number of production wells, the water injection amount of each injection well in each time step, the liquid production amount of each production well in each time step, the average permeability among injection and production wells and the injection and production well spacing; (ii) a
The data preprocessing module P1 optimizes unreasonable data points including:
setting the communication coefficient and the time lag constant between injection wells and production wells of the invalid well to be 0;
and setting the communication coefficient of the injection-production conversion well to be 0.
An initial value calculation module P2, which is used for calculating the initial value of the injection-production inter-well communication coefficient and the initial value of the time lag constant of the oil field block to be determined according to the optimized data of S1;
the initial value calculation module P2 calculates the initial value of the communication coefficient between the injection wells and the production wells by using the following formula:
<math> <mrow> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mover> <msub> <mi>k</mi> <mi>ij</mi> </msub> <mo>&OverBar;</mo> </mover> <mo>/</mo> <msub> <mi>L</mi> <mi>ij</mi> </msub> </mrow> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </msubsup> <mover> <msub> <mi>k</mi> <mi>ij</mi> </msub> <mo>&OverBar;</mo> </mover> <mo>/</mo> <msub> <mi>L</mi> <mi>ij</mi> </msub> </mrow> </mfrac> </mrow> </math>
wherein f isijRepresenting the communication coefficient between the injection well of the ith hole and the production well of the jth hole;
Figure BDA0000436033460000212
the average permeability between the injection well of the ith and the production well of the jth is expressed by 10-3Square micron; l isijThe well spacing of the ith injection well and the jth production well is expressed in meters; n is a radical ofproRepresenting the number of production wells in ports; wherein,
Figure BDA0000436033460000213
Lijare known;
the initial value of the time lag constant adopts the average time lag constant of the oilfield block to be determined;
the calculation step of the average time lag constant of the oilfield block to be determined comprises the following steps:
substituting the optimized data of the data preprocessing module P1 into an average time lag constant formula and calculating the average time lag constant when the formula value of the average time lag constant is smaller than a preset value;
the average time lag constant is expressed as:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>[</mo> <mi>q</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>q</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k</mi> <mi>&tau;</mi> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mo>[</mo> <mrow> <mo>(</mo> <msub> <mi>e</mi> <mi>w</mi> </msub> <mo>+</mo> <mi>fI</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mi>n</mi> <mo>-</mo> <mi>k</mi> </mrow> <mi>&tau;</mi> </mfrac> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mi>&tau;</mi> </mfrac> </mrow> </msup> <mo>)</mo> </mrow> <mo>]</mo> <mo>]</mo> </mrow> <mn>2</mn> </msup> </mrow> </math>
wherein n represents the total calculation time step in months; i (k) represents the water injection quantity of the block at the k time step, and the unit is cubic meter per day; q (k) represents the oil production of the block at the kth time step in cubic meters per day; q (o) represents the overall initial production of the block in cubic meters per day; f represents the block overall connectivity coefficient; τ represents the mean time lag constant in days; e.g. of the typewRepresents external energy replenishment in joules; it is composed ofWherein I (k), q (k) are known amounts
And a result calculation module P3, configured to substitute the optimized data of the data preprocessing module P1 and the initial value of the communication coefficient between injection and production wells and the initial value of the time lag constant calculated in the initial value calculation module P2 into an objective function, and calculate the communication coefficient between injection and production wells, the time lag constant, and the interference constant between production wells when the objective function value is smaller than the fitting error. The objective function is a function for calculating the minimum value of the superposition value of the sum of squares of the difference values of the calculated liquid production rate and the actual liquid production rate of each production well in the oil field block to be determined in all time steps;
the objective function is:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </munderover> <mo>{</mo> <msub> <mi>q</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mo>{</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>[</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> </mrow> </math>
wherein q isj(k) Represents the actual fluid production at k time step in cubic meters per day;
Figure BDA0000436033460000216
representing the initial production of the production well in cubic meters per day; Δ t represents the time step, month; f. ofijRepresenting the communication coefficient of the ith injection well and the jth production well; alpha is alphajRepresenting the interference constant of other production forming wells to the jth production well; i isi(p) represents the water injection quantity of the injection well at the ith hole in p time steps, and the unit is cubic meter per day; tau isjRepresents the time lag constant for the jth production well in days; n is a radical ofiIndicating the number of injection wells in units of openings; n is a radical ofproRepresenting the number of production wells in ports; n represents the total calculation time step in months; wherein q isj(k)、
Figure BDA0000436033460000221
Δt、Ii(p)、Ni、NproAnd n is a known amount.
And the injection-production inter-well communication condition judgment module P4 is used for judging the strength of the communication condition between the injection wells and the production wells in the oil field block to be determined by using the communication coefficient between the injection wells and the production wells calculated by the result calculation module P3, judging the shortest time for the injection water in the oil field block to be determined to reach each production well by using the time lag constant calculated by the result calculation module P3, and judging the influence of the interference of other production wells on each production well in the oil field block to be determined by using the interference constant between the production wells calculated by the result calculation module P3.
In the embodiment, the influence of the interference among the production wells on the communication condition among the injection and production wells is considered, the interference constant among the production wells is introduced into the objective function, and the communication condition among the injection and production wells is judged by combining the communication coefficient among the injection and production wells, the time lag constant and the interference constant among the production wells, so that the rationality of judging the communication condition among the injection and production wells of the high-water-content oil field is enhanced, and the production work of the oil field is guided more favorably.
The initial value calculation module P2 and the injection-production well communication condition determination module P4 in the second embodiment of the injection-production well communication condition determination device of the present invention are the same as those in the first embodiment.
The data preprocessing module P1 in the second embodiment of the present invention is used to optimize unreasonable data points for the data of the oilfield block to be determined. The data of the oilfield block to be determined comprises: the number of injection wells, the number of production wells, the water injection amount of each injection well at each time step, the liquid production amount of each production well at each time step, the average permeability among injection and production wells, the injection and production well spacing, the included angle between the connecting line among the injection and production wells and the material source direction and the average permeability among the injection and production wells in the X, Y direction;
fig. 5 is a diagram illustrating the result calculation module P3 according to the second embodiment. As shown in fig. 5, the result calculation module P3 includes the following modules:
and the constraint module P3.1 is used for judging the initial value of the communication coefficient between the injection wells and the production wells and the initial value of the time lag constant calculated by the initial value calculation module P2 through constraint conditions and rejecting unreasonable judgment data.
The constraint conditions comprise an injection-production inter-well communication coefficient constraint condition, a time lag constant constraint condition and a production inter-well interference constant constraint condition and are used for constraining an injection-production inter-well communication coefficient, a time lag constant and a production inter-well interference constant;
the constraint conditions are as follows:
<math> <mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> <mo>&le;</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <msub> <mi>a</mi> <mi>ij</mi> </msub> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mo>&le;</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>=</mo> <mn>1,2</mn> <mo>,</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>,</mo> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0.145</mn> <mi>&Delta;t</mi> <mo>&le;</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>&le;</mo> <mn>1000</mn> <mi>&Delta;t</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>&GreaterEqual;</mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
wherein,
Figure BDA0000436033460000231
j=1,2,...,NproIs a constraint condition of communication coefficient between injection wells and production wells, aijIn order to influence the coefficients of the effects, <math> <mrow> <msub> <mi>a</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>></mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>&le;</mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow> </math> rijby passing <math> <mrow> <msub> <mi>r</mi> <mi>ij</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>a</mi> <mrow> <mn>1</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>b</mi> <mrow> <mn>1</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>a</mi> <mrow> <mn>2</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>b</mi> <mrow> <mn>2</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>,</mo> </mrow> </math> a 1 ij b 1 ij = a 2 ij b 2 ij = k xij k yij Calculating to obtain;
0.145Δt≤τjthe time lag constant constraint condition is not more than 1000 delta t;
αjthe constraint condition of the interference constant between production wells is more than or equal to 0;
wherein d isijThe distance between the ith injection well and the jth production well is expressed in meters; r isijThe limit influence distance between the ith injection well and the jth production well is expressed in meters; a is1ij、a2ijRespectively representing the major semi-axes of ellipses formed by taking the injection well as the production well as the center of the well pair formed by the ith injection well and the jth production well, and the major semi-axes are meters; b1ij、b2ijRespectively representing the minor semi-axes of ellipses formed by taking the injection well as the production well as the center of the well pair formed by the ith injection well and the jth production well, and the unit is meter; thetaijThe included angle between the direction of the source and the connecting line of the ith injection well and the jth production well is expressed in degrees; Δ t represents a time step in months; f. ofijRepresenting the communication coefficient of the ith injection well and the jth production well; alpha is alphajRepresenting the interference constant of other production forming wells to the jth production well; tau isjRepresents the time lag constant for the jth production well in days; n is a radical ofproRepresenting the total number of producing wells in units of ports; wherein, b1ij、b2ijSimulating one half of the comprehensively determined reasonable well spacing by adopting an oil reservoir engineering numerical value; thetaijAnd Δ t is a known quantity.
And a calculation module P3.2 for substituting the initial value of the communication coefficient between injection wells and production wells and the initial value of the time lag constant after the data are optimized by the data preprocessing module P1 and the initial value of the communication coefficient between injection wells and the production wells are judged by the initial value calculation module P2 into an objective function, and calculating the communication coefficient between injection wells and the production wells, the time lag constant and the interference constant between production wells when the objective function value is smaller than the fitting error.
The objective function is:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </munderover> <mo>{</mo> <msub> <mi>q</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mo>{</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>[</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> </mrow> </math>
wherein q isj(k) Represents the actual fluid production at k time step in cubic meters per day;
Figure BDA0000436033460000236
representing the initial production of the production well in cubic meters per day; Δ t represents the time step, month; f. ofijRepresenting the communication coefficient of the ith injection well and the jth production well; alpha is alphajRepresenting the interference constant of other production forming wells to the jth production well; i isi(p) represents the water injection quantity of the injection well at the ith hole in p time steps, and the unit is cubic meter per day; tau isjRepresents the time lag constant for the jth production well in days; n is a radical ofiIndicating the number of injection wells in units of openings; n is a radical ofproRepresenting the number of production wells in ports; n represents the total calculation time step in months; wherein q isj(k)、Δt、Ii(p)、Ni、NproAnd n is a known amount.
And the fitting error adopts a communication coefficient among injection wells and production wells, a time lag constant and an interference constant among production wells as an optimal solution of a target function when the minf (x) is less than or equal to 0.1.
And when the calculated objective function value is smaller than the fitting error, the injection-production inter-well communication coefficient, the time lag constant and the production inter-well interference constant adopt a sequential quadratic programming method.
And a judging module P3.3 for judging the communication coefficient, the time lag constant and the interference constant among the injection wells and the production wells calculated in the calculating module according to the constraint condition, eliminating unreasonable judging data and obtaining the communication coefficient, the time lag constant and the interference constant among the injection wells and the production wells which are finally used for judging the communication condition among the injection wells and the production wells.
In the second embodiment of the invention, the influence of the interference between the production wells on the communication condition between the injection wells and the production wells is considered, the interference constant between the production wells is introduced into the objective function, and the communication condition between the injection wells and the production wells is judged by combining the communication coefficient between the injection wells and the production wells, the time lag constant and the interference constant between the production wells, so that the rationality for judging the communication condition between the injection wells and the production wells of the high-water-content oil field is greatly enhanced, and the production work of the high-water-content oil field is guided more favorably.
Considering that the injection-production communication relationship is greatly influenced by geological factors, and the injection-production communication coefficient, the time lag constant and the production inter-well interference constant are consistent with geological knowledge, for example, a relatively large communication coefficient should not exist between injection-production wells with relatively long distances, the production inter-well interference constant should not be a negative value and the like, the second embodiment of the invention introduces a constraint condition into the established objective function, discriminates the calculated injection-production inter-well communication coefficient, the time lag constant and the production inter-well interference constant, eliminates unreasonable discrimination data, further improves the rationality of the calculated result, finally secondarily improves the rationality of discrimination of the injection-production inter-well communication condition of the high water-cut oil field, and is more favorable for guiding the fine injection-production of the high water-cut oil field.
The foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be appreciated by those skilled in the art that the present invention is not limited by the embodiments described above, which are described in the specification only to illustrate the principles of the invention. The present invention is subject to various changes and modifications without departing from the spirit and scope of the invention, and such changes and modifications are intended to be included within the scope of the appended claims. The scope of the invention is defined by the claims in this specification and their equivalents.

Claims (10)

1. A method for judging the communication condition between injection wells and production wells is characterized by comprising the following steps:
s1: data points for which data optimization of the oilfield block to be determined is not reasonable, the data of the oilfield block to be determined comprising: the number of injection wells, the number of production wells, the water injection amount of each injection well in each time step, the liquid production amount of each production well in each time step, the average permeability among injection and production wells and the injection and production well spacing;
s2: according to the optimized data of S1, calculating an initial value of the injection-production interwell communication coefficient and an initial value of a time lag constant of the oil field block to be determined;
s3: substituting the optimized data of S1 and the initial value of the communication coefficient between injection wells and the production well calculated in S2 into an objective function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the objective function value is smaller than the fitting error; the objective function is a function for calculating the minimum value of the superposition value of the sum of squares of the difference values of the calculated liquid production rate and the actual liquid production rate of each production well in the oil field block to be determined in all time steps;
s4: judging the strength of the communication condition between the injection wells and the production wells in the oil field block to be determined by utilizing the communication coefficient between the injection wells and the production wells calculated in the S3;
judging the shortest time for the injected water in the oil field block to reach each production well to be determined by utilizing the time lag constant calculated in the S3;
and judging that each production well in the oil field block to be determined is influenced by the interference of other production wells by using the interference constant among the production wells calculated by the S3.
2. The method for determining the communication status between injection wells and production wells according to claim 1,
the data of the oilfield block to be determined in S1 further includes: the included angle between the connecting line between the injection wells and the source direction and the average permeability between the injection wells and the production wells in the X, Y direction;
the step S4 is implemented by the following steps:
a. judging the initial value of the communication coefficient between the injection wells and the production wells and the initial value of the time lag constant calculated in the S2 through constraint conditions, and rejecting unreasonable judging data; the constraint conditions comprise an injection-production inter-well communication coefficient constraint condition, a time lag constant constraint condition and a production inter-well interference constant constraint condition and are used for constraining an injection-production inter-well communication coefficient, a time lag constant and a production inter-well interference constant;
b. substituting the optimized data of S1 and the initial value of the communication coefficient between injection wells and the initial value of the time lag constant after the judgment of a into an objective function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the objective function value is smaller than a preset value;
c. and c, judging the communication coefficient among the injection wells, the time lag constant and the interference constant among the production wells obtained by calculation in the step b through the constraint condition, eliminating unreasonable judgment data, and obtaining the communication coefficient among the injection wells, the time lag constant and the interference constant among the production wells which are finally used for judging the communication condition among the injection wells.
3. A method for determining the communication status between injection and production wells according to claim 1 or 2, wherein the step S1 of optimizing unreasonable data points includes:
setting the communication coefficient and the time lag constant between injection wells and production wells of the invalid well to be 0;
and setting the communication coefficient of the injection-production conversion well to be 0.
4. The method for determining the communication status between injection wells and production wells according to claim 1 or 2, wherein the step of calculating the initial value of the communication coefficient between injection wells and production wells of the oilfield block to be determined in S2 includes:
calculating to obtain an initial value of the communication coefficient between the injection wells and the production wells by using the following formula:
<math> <mrow> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mover> <msub> <mi>k</mi> <mi>ij</mi> </msub> <mo>&OverBar;</mo> </mover> <mo>/</mo> <msub> <mi>L</mi> <mi>ij</mi> </msub> </mrow> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </msubsup> <mover> <msub> <mi>k</mi> <mi>ij</mi> </msub> <mo>&OverBar;</mo> </mover> <mo>/</mo> <msub> <mi>L</mi> <mi>ij</mi> </msub> </mrow> </mfrac> </mrow> </math>
wherein f isijEighthly, representing the communication coefficient between the injection well at the ith hole and the production well at the jth hole;
Figure FDA0000436033450000022
the average permeability between the injection well of the ith and the production well of the jth is expressed by 10-3Square micron; l isijThe well spacing of the ith injection well and the jth production well is expressed in meters; n is a radical ofproRepresenting the number of production wells in ports; wherein N ispro
Figure FDA0000436033450000023
LijIs a known value.
5. The method for judging the communication condition between the injection and production wells according to claim 1 or 2, wherein the initial value of the time lag constant is the average time lag constant of the oilfield block to be determined;
the calculation step of the average time lag constant of the oilfield block to be determined comprises the following steps:
substituting the optimized data of S1 into an average time lag constant formula and calculating the average time lag constant when the formula value of the average time lag constant is smaller than a preset value;
the average time lag constant is expressed as:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>[</mo> <mi>q</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>q</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>k</mi> <mi>&tau;</mi> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mo>[</mo> <mrow> <mo>(</mo> <msub> <mi>e</mi> <mi>w</mi> </msub> <mo>+</mo> <mi>fI</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>n</mi> <mo>-</mo> <mi>k</mi> </mrow> <mi>&tau;</mi> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mi>&tau;</mi> </mfrac> </mrow> </msup> <mo>)</mo> </mrow> <mo>]</mo> <mo>]</mo> </mrow> <mn>2</mn> </msup> </mrow> </math>
wherein n represents the total calculation time step in months; i (k) indicates that the block is at the k-th time stepThe unit of water injection is cubic meter per day; q (k) represents the oil production of the block at the kth time step in cubic meters per day; q (0) represents the overall initial production of the block in cubic meters per day; f represents the block overall connectivity coefficient; τ represents the mean time lag constant in days; e.g. of the typewRepresents external energy replenishment in joules; wherein I (k), q (k) are known amounts.
6. The method for judging the communication condition between injection and production wells according to claim 1 or 2, wherein the fitting error is 0.1.
7. A method for determining the communication status between injection and production wells according to claim 1 or 2, wherein the objective function in S3 is:
<math> <mrow> <mi>min</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>pro</mi> </msub> </munderover> <mo>{</mo> <msub> <mi>q</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mover> <mi>q</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>k</mi> <mi>&Delta;t</mi> </mrow> <mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mi></mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mo>{</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>[</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>ij</mi> </msub> <msub> <mi>I</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>p</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <msup> <mrow> <mo>[</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> <mi>&Delta;t</mi> </mrow> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </mrow> </msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mfrac> <mi>&Delta;t</mi> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> </mfrac> </msup> <mo>)</mo> </mrow> <mo>]</mo> <mo>}</mo> </mrow> <mn>2</mn> </msup> </mrow> </math>
wherein q isj(k) Representing the actual liquid production amount of the jth production well in k time step, wherein the unit is cubic meter per day;
Figure FDA0000436033450000032
representing the initial yield of the jth production well in cubic meters per day; Δ t represents a time step in months; f. ofijRepresenting the communication coefficient of the ith injection well and the jth production well; alpha is alphajRepresenting the interference constant of other production forming wells to the jth production well; i isi(p) represents the water injection quantity of the injection well at the ith hole in p time steps, and the unit is cubic meter per day; tau isjRepresents the time lag constant for the jth production well in days; n is a radical ofiIndicating the number of injection wells in units of openings; n is a radical ofproRepresenting the number of production wells in ports; n represents the total calculation time step in months; wherein q isj(k)、
Figure FDA0000436033450000033
Δt、Ii(p)、Ni、NproAnd n is a known amount.
8. The method for judging the communication condition between the injection and production wells according to claim 2, wherein the constraint condition is as follows:
<math> <mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> <mo>&le;</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>i</mi> </msub> </munderover> <msub> <mi>&alpha;</mi> <mi>ij</mi> </msub> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mo>&le;</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>=</mo> <mn>1,2</mn> <mo>,</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>,</mo> <msub> <mi>N</mi> <mi>pro</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0.145</mn> <mi>&Delta;t</mi> <mo>&le;</mo> <msub> <mi>&tau;</mi> <mi>j</mi> </msub> <mo>&le;</mo> <mn>1000</mn> <mi>&Delta;t</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&alpha;</mi> <mi>j</mi> </msub> <mo>&GreaterEqual;</mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
wherein,j=1,2,...,Nprois a constraint condition of communication coefficient between injection wells and production wells, aijIn order to influence the coefficients of the effects,
<math> <mrow> <msub> <mi>a</mi> <mi>ij</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>></mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>d</mi> <mi>ij</mi> </msub> <mo>&le;</mo> <msub> <mi>r</mi> <mi>ij</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow> </math> said rijBy passing <math> <mrow> <msub> <mi>r</mi> <mi>ij</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>a</mi> <mrow> <mn>1</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>b</mi> <mrow> <mn>1</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>a</mi> <mrow> <mn>2</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mrow> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mi>ij</mi> </msub> </mrow> <msup> <msub> <mi>b</mi> <mrow> <mn>2</mn> <mi>ij</mi> </mrow> </msub> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msup> <mo>,</mo> </mrow> </math>
a 1 ij b 1 ij = a 2 ij b 2 ij = k xij k yij Calculating to obtain;
0.145Δt≤τjthe time lag constant constraint condition is not more than 1000 delta t;
αjthe constraint condition of the interference constant between production wells is more than or equal to 0;
wherein d isijThe distance between the ith injection well and the jth production well is expressed in meters; r isijThe limit influence distance between the ith injection well and the jth production well is expressed in meters; a is1ij、a2ijRespectively representing the major semi-axes of ellipses formed by taking the injection well as the production well as the center of the well pair formed by the ith injection well and the jth production well, and the major semi-axes are meters; b1ij、b2ijRespectively representing the minor semi-axes of ellipses formed by taking the injection well as the production well as the center of the well pair formed by the ith injection well and the jth production well, and the unit is meter; thetaijThe included angle between the direction of the source and the connecting line of the ith injection well and the jth production well is expressed in degrees; Δ t represents a time step in months; f. ofijRepresenting the communication coefficient of the ith injection well and the jth production well; alpha is alphajRepresenting the interference constant of other production forming wells to the jth production well; tau isj(ii) a Represents the time lag constant for the jth production well in days; n is a radical ofproRepresenting the total number of producing wells in units of ports; wherein, b1ij、b2ijSimulating one half of the comprehensively determined reasonable well spacing by adopting an oil reservoir engineering numerical value; thetaijAnd Δ t is a known quantity.
9. A device for judging the communication condition between injection wells and production wells is characterized by comprising
The data preprocessing module is used for optimizing unreasonable data points of the oilfield block to be determined, and the data of the oilfield block to be determined comprises: the number of injection wells, the number of production wells, the water injection amount of each injection well in each time step, the liquid production amount of each production well in each time step, the average permeability among injection and production wells and the injection and production well spacing;
the initial value calculation module is used for calculating an initial value of the communication coefficient between the injection wells and the production wells of the oil field block to be determined and an initial value of a time lag constant according to the data optimized by the data preprocessing module;
the result calculation module is used for substituting the data optimized by the data preprocessing module and the initial value of the communication coefficient between injection wells and the initial value of the time lag constant calculated by the initial value calculation module into an objective function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the objective function value is smaller than the fitting error; the objective function is a function for calculating the minimum value of the superposition value of the sum of squares of the difference values of the calculated liquid production rate and the actual liquid production rate of each production well in the oil field block to be determined in all time steps;
and the communication condition judging module is used for judging the strength of the communication condition between the injection wells and the production wells in the oil field block to be determined by utilizing the communication coefficient between the injection wells and the production wells calculated in the S3, judging the shortest time for the injection water to reach each production well by utilizing the time lag constant calculated in the S3 and judging the interference influence of each production well on other production wells by utilizing the interference constant between the production wells calculated in the S3.
10. The apparatus according to claim 9, wherein the data of the field block to be determined in the data preprocessing module further includes an angle between a connecting line between injection wells and a source direction and an average permeability between injection wells in a direction of X, Y;
the result calculation module further comprises the following modules:
the constraint module is used for judging the initial value of the communication coefficient between the injection wells and the production wells and the initial value of the time lag constant calculated in the introduced initial value calculation module through constraint conditions and rejecting unreasonable judgment data; the constraint conditions comprise an injection-production inter-well communication coefficient constraint condition, a time lag constant constraint condition and a production inter-well interference constant constraint condition and are used for constraining an injection-production inter-well communication coefficient, a time lag constant and a production inter-well interference constant;
the calculation module is used for substituting the data optimized by the data preprocessing module and the initial value of the communication coefficient between injection wells and the initial value of the time lag constant judged by the constraint module into a target function, and calculating the communication coefficient between injection wells, the time lag constant and the interference constant between production wells when the target function value is smaller than the fitting error;
and the judging module is used for judging the communication coefficient among the injection wells, the time lag constant and the interference constant among the production wells which are calculated in the calculating module through the constraint condition, eliminating unreasonable judging data and obtaining the communication coefficient among the injection wells, the time lag constant and the interference constant among the production wells which are finally used for judging the communication condition among the injection wells.
CN201310681831.5A 2013-12-12 2013-12-12 Method and device for judging communication condition between injection wells and production wells Active CN103670369B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310681831.5A CN103670369B (en) 2013-12-12 2013-12-12 Method and device for judging communication condition between injection wells and production wells

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310681831.5A CN103670369B (en) 2013-12-12 2013-12-12 Method and device for judging communication condition between injection wells and production wells

Publications (2)

Publication Number Publication Date
CN103670369A true CN103670369A (en) 2014-03-26
CN103670369B CN103670369B (en) 2016-08-17

Family

ID=50309059

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310681831.5A Active CN103670369B (en) 2013-12-12 2013-12-12 Method and device for judging communication condition between injection wells and production wells

Country Status (1)

Country Link
CN (1) CN103670369B (en)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105089649A (en) * 2014-04-17 2015-11-25 中国石油化工股份有限公司 Method for establishing inter-well connectivity model
CN105404726A (en) * 2015-11-02 2016-03-16 中国地质大学(武汉) Method and device for reverting inter-well connectivity through Gaussian distribution-based capacitance model
CN107291667A (en) * 2016-04-01 2017-10-24 中国石油化工股份有限公司 A kind of interwell communication degree determines method and system
CN105971577B (en) * 2016-07-08 2018-09-04 中国石油天然气股份有限公司 Method and device for improving connectivity of fireflooding injection well and production well
CN110410065A (en) * 2019-06-24 2019-11-05 中国石油化工股份有限公司 The input parameters analysis method and device of multilayer inter-well dynamic connectivity model
CN111119808A (en) * 2018-10-31 2020-05-08 北京国双科技有限公司 Construction method of oil field injection-production relation model and identification method of oil field injection-production relation
CN111126650A (en) * 2018-10-31 2020-05-08 北京国双科技有限公司 Model construction method and device for oil field simulation system
CN111364953A (en) * 2020-03-09 2020-07-03 中国石油大学(华东) Method and system for identifying dynamic connectivity among wells based on graph neural network
CN111364970A (en) * 2020-03-30 2020-07-03 海默潘多拉数据科技(深圳)有限公司 A Quantification Method of Interwell Connectivity Coefficient
CN111997575A (en) * 2020-08-21 2020-11-27 中国海洋石油集团有限公司 Irregular well pattern optimization water injection method based on planar displacement component
CN112302594A (en) * 2020-12-01 2021-02-02 中国石油大学(北京) Method, Apparatus and Equipment for Determining Connecting Structure of Heterogeneous Water Drive Oil Reservoir
CN116305963A (en) * 2023-03-23 2023-06-23 西南石油大学 Inter-well connectivity evaluation method based on time lag correlation analysis
CN116341406A (en) * 2023-03-03 2023-06-27 常州大学 A Calculation Method for Determining the Timing of Injection-production Interference in Inclined Reservoirs
CN117684929A (en) * 2022-12-14 2024-03-12 中国科学院沈阳自动化研究所 A global optimization control method for energy consumption in oil and water well systems based on inter-well connectivity

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107130951B (en) * 2017-05-17 2019-09-10 中国石油天然气股份有限公司 Method and system for monitoring communication condition between steam flooding wells

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5679893A (en) * 1996-02-12 1997-10-21 Chevron, U.S.A., Inc. Determining the oil decline for an oil-producing reservoir
RU2278250C1 (en) * 2005-03-09 2006-06-20 Государственное образовательное учреждение высшего профессионального образования "Санкт-Петербургский государственный горный институт им. Г.В. Плеханова (технический университет)" Oil field development
WO2011037574A1 (en) * 2009-09-25 2011-03-31 Hewlett-Packard Development Company, L.P. Silicon-germanium, quantum-well, light-emitting diode
CN102146788A (en) * 2011-03-30 2011-08-10 中国石油化工股份有限公司 Method for enhancing recovery ratio through three-dimensional streamline regulation of water drive reservoir
CN102419787A (en) * 2011-10-19 2012-04-18 中国石油天然气股份有限公司 Method and system for identifying water flow dominant channel of high-water-content oil field
CN102704924A (en) * 2012-06-05 2012-10-03 中国石油天然气股份有限公司 Method and device for determining effective dry layer
CN103422849A (en) * 2013-07-18 2013-12-04 中国石油天然气股份有限公司 Method and device for reconstructing well pattern of old oil field developed by water injection

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5679893A (en) * 1996-02-12 1997-10-21 Chevron, U.S.A., Inc. Determining the oil decline for an oil-producing reservoir
RU2278250C1 (en) * 2005-03-09 2006-06-20 Государственное образовательное учреждение высшего профессионального образования "Санкт-Петербургский государственный горный институт им. Г.В. Плеханова (технический университет)" Oil field development
WO2011037574A1 (en) * 2009-09-25 2011-03-31 Hewlett-Packard Development Company, L.P. Silicon-germanium, quantum-well, light-emitting diode
CN102146788A (en) * 2011-03-30 2011-08-10 中国石油化工股份有限公司 Method for enhancing recovery ratio through three-dimensional streamline regulation of water drive reservoir
CN102419787A (en) * 2011-10-19 2012-04-18 中国石油天然气股份有限公司 Method and system for identifying water flow dominant channel of high-water-content oil field
CN102704924A (en) * 2012-06-05 2012-10-03 中国石油天然气股份有限公司 Method and device for determining effective dry layer
CN103422849A (en) * 2013-07-18 2013-12-04 中国石油天然气股份有限公司 Method and device for reconstructing well pattern of old oil field developed by water injection

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
冯其红等: "低渗透油藏井间动态连通性研究方法", 《特种油气藏》 *
尚福华等: "基于决策树的注采连通关系判别研究", 《计算机应用研究》 *
朱国金等: "注采井动态连通关系计算方法研究及应用", 《中国海上油气》 *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105089649B (en) * 2014-04-17 2018-03-13 中国石油化工股份有限公司 A kind of inter well connectivity method for establishing model
CN105089649A (en) * 2014-04-17 2015-11-25 中国石油化工股份有限公司 Method for establishing inter-well connectivity model
CN105404726A (en) * 2015-11-02 2016-03-16 中国地质大学(武汉) Method and device for reverting inter-well connectivity through Gaussian distribution-based capacitance model
CN105404726B (en) * 2015-11-02 2018-02-13 中国地质大学(武汉) A kind of capacitor model inverting inter well connectivity method and device based on Gaussian Profile
CN107291667B (en) * 2016-04-01 2020-11-13 中国石油化工股份有限公司 Method and system for determining communication degree between wells
CN107291667A (en) * 2016-04-01 2017-10-24 中国石油化工股份有限公司 A kind of interwell communication degree determines method and system
CN105971577B (en) * 2016-07-08 2018-09-04 中国石油天然气股份有限公司 Method and device for improving connectivity of fireflooding injection well and production well
CN111119808A (en) * 2018-10-31 2020-05-08 北京国双科技有限公司 Construction method of oil field injection-production relation model and identification method of oil field injection-production relation
CN111126650A (en) * 2018-10-31 2020-05-08 北京国双科技有限公司 Model construction method and device for oil field simulation system
CN111126650B (en) * 2018-10-31 2022-05-27 北京国双科技有限公司 Model construction method and device for oil field simulation system
CN111119808B (en) * 2018-10-31 2022-04-12 北京国双科技有限公司 Construction method of oil field injection-production relation model and identification method of oil field injection-production relation
CN110410065A (en) * 2019-06-24 2019-11-05 中国石油化工股份有限公司 The input parameters analysis method and device of multilayer inter-well dynamic connectivity model
CN110410065B (en) * 2019-06-24 2022-05-03 中国石油化工股份有限公司 Input parameter analysis method and device for dynamic connectivity model among multiple wells
CN111364953A (en) * 2020-03-09 2020-07-03 中国石油大学(华东) Method and system for identifying dynamic connectivity among wells based on graph neural network
CN111364953B (en) * 2020-03-09 2022-02-22 中国石油大学(华东) Method and system for identifying dynamic connectivity among wells based on graph neural network
CN111364970B (en) * 2020-03-30 2022-03-01 海默潘多拉数据科技(深圳)有限公司 Method for quantizing inter-well communication coefficient
CN111364970A (en) * 2020-03-30 2020-07-03 海默潘多拉数据科技(深圳)有限公司 A Quantification Method of Interwell Connectivity Coefficient
CN111997575A (en) * 2020-08-21 2020-11-27 中国海洋石油集团有限公司 Irregular well pattern optimization water injection method based on planar displacement component
CN111997575B (en) * 2020-08-21 2022-07-29 中国海洋石油集团有限公司 Irregular well pattern optimization water injection method based on planar displacement component
CN112302594A (en) * 2020-12-01 2021-02-02 中国石油大学(北京) Method, Apparatus and Equipment for Determining Connecting Structure of Heterogeneous Water Drive Oil Reservoir
CN117684929A (en) * 2022-12-14 2024-03-12 中国科学院沈阳自动化研究所 A global optimization control method for energy consumption in oil and water well systems based on inter-well connectivity
CN116341406A (en) * 2023-03-03 2023-06-27 常州大学 A Calculation Method for Determining the Timing of Injection-production Interference in Inclined Reservoirs
CN116341406B (en) * 2023-03-03 2024-04-23 常州大学 Calculation method for determining injection and production interference time of inclined oil reservoir
CN116305963A (en) * 2023-03-23 2023-06-23 西南石油大学 Inter-well connectivity evaluation method based on time lag correlation analysis
CN116305963B (en) * 2023-03-23 2025-03-11 西南石油大学 Inter-well connectivity evaluation method based on time lag correlation analysis

Also Published As

Publication number Publication date
CN103670369B (en) 2016-08-17

Similar Documents

Publication Publication Date Title
CN103670369B (en) Method and device for judging communication condition between injection wells and production wells
CN105117511B (en) A kind of characterizing method of fracture hole oil reservoir interwell communication passage and flow parameter
CN103939066B (en) An a kind of note adopts the method that well group determines water injection rate determination oilwell produced fluid amount more
CN109948272A (en) Tune based on inter well connectivity blocks up dynamic prediction method and system
CN108984886B (en) Method for inverting dynamic connectivity among wells of multi-layer oil reservoir
CN104989341B (en) Method for determining effective displacement injection-production well spacing of low-permeability oil reservoir
CN110147561B (en) Method for predicting volume fracture network of tight oil and gas reservoir containing natural fracture
CN106599374B (en) A kind of river facies relation and its derivation method suitable for runoff tidal estuary section
CN109577929B (en) Quantitative evaluation method for establishing effective displacement of ultra-low permeability tight reservoir horizontal well
CN101725346A (en) Oil deposit inter-well dynamic connectivity inverting method
CN104060985B (en) A kind of stratified deposit profile control and water plugging blocking agent penetration depth method of testing and system
CN104879104B (en) Oil reservoir water injection method
CN105631529A (en) Edge water gas reservoir water breakthrough time prediction method
CN113836695B (en) Oil reservoir numerical simulation method based on gridless connecting element
CN114357592B (en) Analytical method for water inflow of suspended water-stop curtain in circular foundation pit of confined water layer
CN104915530A (en) Method for establishing communication relationship between oil deposit wells
CN112836438B (en) Method and device for predicting water breakthrough time of gas well of heterogeneous edge water-gas reservoir
CN104407397A (en) Method for determining lower limit of physical property of compact oil and gas reservoir
CN111927410A (en) Well-reservoir collaborative separate injection test regulation design method
CN117235963A (en) Injection and production unit water drive volume sweep coefficient calculation method
CN106682333A (en) Method for determining rational well spacing of low-permeability reservoir
CN105631763A (en) Underground water type drinking water source protection area division method
CN113343434B (en) A development method using large pores in oil fields
CN116703072A (en) Water Quantity-Water Quality Comprehensive Dispatching Method in Areas Connected by Rivers and Reservoirs
Lian et al. Evaluation and applicability study on prediction methods of water inflow in mines

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant