CN110410065B - Input parameter analysis method and device for dynamic connectivity model among multiple wells - Google Patents
Input parameter analysis method and device for dynamic connectivity model among multiple wells Download PDFInfo
- Publication number
- CN110410065B CN110410065B CN201910547438.4A CN201910547438A CN110410065B CN 110410065 B CN110410065 B CN 110410065B CN 201910547438 A CN201910547438 A CN 201910547438A CN 110410065 B CN110410065 B CN 110410065B
- Authority
- CN
- China
- Prior art keywords
- production
- wells
- coefficient
- multilayer
- inter
- 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.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 14
- 238000004519 manufacturing process Methods 0.000 claims abstract description 220
- 238000002347 injection Methods 0.000 claims abstract description 104
- 239000007924 injection Substances 0.000 claims abstract description 104
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 71
- 238000009826 distribution Methods 0.000 claims abstract description 13
- 239000010410 layer Substances 0.000 claims description 73
- 239000002356 single layer Substances 0.000 claims description 28
- 239000012530 fluid Substances 0.000 claims description 27
- 238000000034 method Methods 0.000 claims description 20
- 238000000605 extraction Methods 0.000 claims description 19
- 239000007788 liquid Substances 0.000 claims description 18
- 240000004282 Grewia occidentalis Species 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 abstract description 8
- 238000004364 calculation method Methods 0.000 description 7
- 101150004094 PRO2 gene Proteins 0.000 description 6
- 239000000243 solution Substances 0.000 description 6
- DPXHITFUCHFTKR-UHFFFAOYSA-L To-Pro-1 Chemical compound [I-].[I-].S1C2=CC=CC=C2[N+](C)=C1C=C1C2=CC=CC=C2N(CCC[N+](C)(C)C)C=C1 DPXHITFUCHFTKR-UHFFFAOYSA-L 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 101100084404 Mus musculus Prodh gene Proteins 0.000 description 3
- QHNORJFCVHUPNH-UHFFFAOYSA-L To-Pro-3 Chemical compound [I-].[I-].S1C2=CC=CC=C2[N+](C)=C1C=CC=C1C2=CC=CC=C2N(CCC[N+](C)(C)C)C=C1 QHNORJFCVHUPNH-UHFFFAOYSA-L 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000009472 formulation Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- IAZDPXIOMUYVGZ-UHFFFAOYSA-N Dimethylsulphoxide Chemical compound CS(C)=O IAZDPXIOMUYVGZ-UHFFFAOYSA-N 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides an input parameter analysis method of a multilayer interwell dynamic connectivity model, which comprises the following steps: determining the distribution mode of the production wells and the water injection wells in the current area, and calculating inter-well dynamic communication coefficients between the water injection wells and the production wells in the current area and multilayer asymmetric coefficients corresponding to the inter-well dynamic communication coefficients based on a multilayer inter-well dynamic communication model; determining an overdetermined coefficient corresponding to a dynamic connectivity model among the multiple layers of wells based on the production conditions and the production data of the production wells in the block; and selecting the optimal over-determined coefficient based on the multilayer asymmetric coefficients and the over-determined coefficient to determine the number of data points in the input parameters of the multilayer inter-well dynamic connectivity model, wherein the data points comprise production data of production wells in the block in a preset time interval. The invention can comprehensively and accurately reflect the formation information under the condition of limited data points, and can also truly reflect the communication coefficient between the current formation injection and production wells.
Description
Technical Field
The invention relates to the technical field of oil and gas field development, in particular to an input parameter analysis method and device of a multilayer inter-well dynamic connectivity model.
Background
The water flooding development is the most widely applied development mode of most oil fields, when the oil field is in the middle and later stages of development, the oil field enters a high water content and ultra-high water content stage, and because of reservoir heterogeneity, fluid fluidity difference, injection and production difference and the like, a water flow dominant channel is gradually formed in a stratum, so that a stratum seepage field is formed in a fixed position, the circulation of low-efficiency or ineffective water is serious, the water storage rate of the stratum is reduced, and the water flooding development effect is seriously influenced. The most important input parameters of the multi-layer interwell dynamic connectivity model are water injection amount of a water injection well and liquid production amount of a production well, the two parameters are basic data for inverting interwell dynamic connectivity, are original data of signals generated by water injection amount in reservoir propagation, and determine a final inversion result.
In applying the multi-layer inter-well dynamic connectivity model, a data point is taken at a predetermined time interval, and the data point contains all production data of the production well in the time period. The quantity of the data points directly determines the quality of an inversion result, if the quantity of the data points is too small, signals transmitted in the stratum cannot comprehensively and accurately reflect stratum information, the final solution result is influenced, and the obtained inter-well dynamic communication coefficient cannot reflect the real communication relation between injection wells and extraction wells. However, if too many data points are obtained, on one hand, a new requirement is provided for the storage life of the field production data, on the other hand, too many data points are obtained, which means that the field production history time is long, and during the whole production period, the inter-well dynamic communication relation between injection and production wells may change, and at this time, the inter-well dynamic communication coefficient obtained by inversion is influenced too much by the previous production data, and the communication coefficient between the injection and production wells in the current stratum cannot be truly reflected. Therefore, it is important to select a proper amount of data points.
Therefore, the invention provides an input parameter analysis method and device of a multilayer interwell dynamic connectivity model.
Disclosure of Invention
In order to solve the problems, the invention provides an input parameter analysis method of a multilayer interwell dynamic connectivity model, which comprises the following steps:
determining distribution modes of production wells and water injection wells in a current block, and calculating inter-well dynamic communication coefficients between the water injection wells and the production wells in the current block and multilayer asymmetric coefficients corresponding to the inter-well dynamic communication coefficients based on a multilayer inter-well dynamic communication model;
secondly, determining an overdetermined coefficient corresponding to the multi-layer inter-well dynamic connectivity model based on the production condition and the production data of the production wells in the block;
and thirdly, selecting an optimal overdetermined coefficient based on the multilayer asymmetric coefficients and the overdetermined coefficient to determine the number of data points in the input parameters of the multilayer inter-well dynamic connectivity model, wherein the data points comprise production data of production wells in the block in a preset time interval.
According to an embodiment of the present invention, when the distribution of the production wells and the water injection wells in the block is a multi-layer five-injection four-production model, the first step specifically comprises the following steps:
In a single oil reservoir layer, calculating to obtain a single-layer asymmetric coefficient based on the inter-well dynamic communication coefficient and the five-injection four-extraction model;
and in the multilayer oil reservoir layer, calculating the multilayer asymmetric coefficient by combining the five-injection four-extraction model based on the single-layer asymmetric coefficient.
According to one embodiment of the present invention, the single layer asymmetry coefficient is calculated by the following formula:
wherein A represents the single layer asymmetric coefficient, var1(λ) represents the variance of the inter-well dynamic communication coefficients between the central water injection well and the adjacent four production wells, var2(λ) represents the variance of the inter-well dynamic communication coefficients between the four-corner injection wells and the two adjacent production wells, var3And (lambda) represents the variance of the inter-well dynamic communication coefficients between the water injection wells at the four corners and the two production wells farther away.
According to one embodiment of the present invention, the multilayer asymmetry coefficient is calculated by the following formula:
wherein A' represents the multilayer asymmetric coefficient, L represents the number of reservoir layers, vari1(λ) represents the variance of the inter-well dynamic communication coefficients between the central water injection well of the ith zone and the adjacent four production wells, vari2(lambda) represents the variance of the inter-well dynamic communication coefficients between the water injection wells at the four corners of the ith layer and two adjacent production wells, var i3And (lambda) represents the variance of the inter-well dynamic communication coefficients between the water injection wells at the four corners of the ith layer and the two more distant production wells respectively.
According to an embodiment of the present invention, the second step specifically includes the following steps:
and judging whether the production well in the block has a liquid production profile or not based on the production condition of the production well in the block.
According to one embodiment of the invention, the over-determined factor is determined for a production well having a fluid production profile by the following formula:
wherein, Od pAnd representing the overdetermined coefficient corresponding to the dynamic connectivity model among the multilayer wells when the fluid production profile is known, wherein M represents the number of the data points, and I' represents the number of water injection wells communicated with the production wells.
According to one embodiment of the invention, the over-determined factor is determined for a production well without a fluid production profile by the following formula:
wherein, Od IAnd expressing the overdetermined coefficient corresponding to the dynamic connectivity model among the multilayer wells when the fluid production profile is unknown, and L expresses the number of the oil reservoir layers.
According to an embodiment of the present invention, the third step specifically includes the following steps:
and drawing a relation curve graph by taking the over-definite coefficient as a horizontal coordinate and the multilayer asymmetric coefficient as a vertical coordinate, and determining the optimal over-definite coefficient through the relation curve graph.
According to another aspect of the present invention, there is also provided an input parameter analysis apparatus for a multi-layer interwell dynamic connectivity model, the apparatus comprising:
the first module is used for determining the distribution mode of the production wells and the water injection wells in the current block, and calculating inter-well dynamic communication coefficients between the water injection wells and the production wells in the current block and multilayer asymmetric coefficients corresponding to the inter-well dynamic communication coefficients based on a multilayer inter-well dynamic communication model;
a second module for determining an over-determined coefficient corresponding to the multi-layer inter-well dynamic connectivity model based on production conditions and production data of production wells within the block;
and a third module for selecting an optimal overdetermined coefficient based on the multilayer asymmetric coefficients and the overdetermined coefficient to determine the number of data points in the input parameters of the multilayer inter-well dynamic connectivity model, wherein the data points contain production data of production wells in a block within a preset time interval.
According to one embodiment of the invention, when the production wells and the water injection wells in the block are distributed in a multi-layer five-injection four-production model, the first module comprises:
the single-layer asymmetric coefficient unit is used for calculating a single-layer asymmetric coefficient in a single oil reservoir layer based on the inter-well dynamic communication coefficient and the five-injection four-extraction model;
And the multilayer asymmetric coefficient unit is used for calculating the multilayer asymmetric coefficient by combining the five-injection four-extraction model based on the single-layer asymmetric coefficient in the multilayer oil reservoir layer.
The method and the device for analyzing the input parameters of the dynamic connectivity model among the multilayer wells determine the most important input parameters (water injection quantity of the water injection well and oil production quantity of the production well) of the dynamic connectivity model among the multilayer wells by comparing the asymmetric coefficient of the dynamic connectivity coefficient among the wells with the over-determined coefficient of the production well, thereby determining the final inversion result of the model and reasonably reflecting the formation information. The invention can comprehensively and accurately reflect the formation information under the condition of limited data points, and simultaneously, the proper amount of selected data points can also truly reflect the communication coefficient between the current formation injection and production wells, thereby fully considering the problem of actual parameter determination of the oil field site and having important reference value for the formulation of the later development scheme of the oil field in the high water cut period.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention and not to limit the invention. In the drawings:
FIG. 1 shows a flow diagram of a method for input parameter analysis of a multi-layer interwell dynamic connectivity model, according to one embodiment of the invention;
FIG. 2 is a plan view of an injection-production well for a five-injection four-production model;
FIG. 3 shows a flow diagram of a method for input parameter analysis of a multi-layer interwell dynamic connectivity model, according to another embodiment of the invention;
FIG. 4 shows a graph of simulated fluid production versus fitted fluid production for a Pro1 production well;
FIG. 5 shows a graph of simulated fluid production versus fitted fluid production for a Pro2 production well;
FIG. 6 shows a graph of simulated fluid production versus fitted fluid production for a Pro3 production well;
FIG. 7 shows a graph of simulated fluid production versus fitted fluid production for a Pro4 production well;
FIG. 8 shows a graph of the over-timed coefficient versus the multilayer asymmetry coefficient for a known production profile;
FIG. 9 shows a graph of the over-determined coefficient versus the multilayer asymmetry coefficient for unknown production profiles; and
FIG. 10 is a block diagram of an input parameter analysis apparatus for a multi-layer inter-well dynamic connectivity model according to an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention are described in further detail below with reference to the accompanying drawings.
In the prior art, characterization of fracture-cavity reservoir flow parameters is considered in a characterization method of inter-well communication channels and flow parameters of a fracture-cavity reservoir, and a detailed solving method is provided, but the method cannot optimize and analyze input parameters of an inter-well communication model, and has certain limitations.
FIG. 1 shows a flow chart of a method for analyzing input parameters of a multi-layer interwell dynamic connectivity model, according to one embodiment of the invention.
As shown in fig. 1, in step S101, the distribution modes of the production wells and the water injection wells in the current block are determined, and based on a multi-layer inter-well dynamic connectivity model, inter-well dynamic connectivity coefficients between the water injection wells and the production wells in the current block and multi-layer asymmetric coefficients corresponding to the inter-well dynamic connectivity coefficients are calculated.
Preferably, when the distribution of the production wells and the water injection wells in the block is a multi-layer five-injection four-production model (as shown in fig. 2), the step S101 specifically includes the following steps:
firstly, calculating a single-layer asymmetric coefficient of a single-layer oil reservoir based on an interwell dynamic communication coefficient and a five-injection four-extraction model.
Preferably, the single layer asymmetry factor is calculated by the following equation:
wherein A represents a single layer asymmetry coefficient, var1(λ) represents the variance of the inter-well dynamic communication coefficients between the central water injection well and the adjacent four production wells, var2(λ) represents the variance of the inter-well dynamic communication coefficients between the four-corner injection wells and the two adjacent production wells, var3And (lambda) represents the variance of the inter-well dynamic communication coefficients between the water injection wells at the four corners and the two production wells farther away.
And then, calculating the multilayer asymmetric coefficients of the multilayer oil reservoir by combining a five-injection four-extraction model based on the single-layer asymmetric coefficients.
Preferably, the multilayer asymmetry coefficient is calculated by the following formula:
wherein A' represents a multilayer asymmetric coefficient, LRepresents the number of reservoir layers, vari1(λ) represents the variance of the inter-well dynamic communication coefficients between the central water injection well of the ith zone and the adjacent four production wells, vari2(lambda) represents the variance of the inter-well dynamic communication coefficients between the water injection wells at the four corners of the ith layer and two adjacent production wells, vari3And (lambda) represents the variance of the inter-well dynamic communication coefficients between the water injection wells at the four corners of the ith layer and the two more distant production wells respectively.
Next, in step S102, the over-determined coefficients corresponding to the multi-layer inter-well dynamic connectivity model are determined based on the production conditions of the production wells within the block and the production data.
Preferably, step S102 specifically includes the following steps:
and judging whether the production well in the block has a liquid production profile or not based on the production condition of the production well in the block.
Preferably, for a production well having a fluid production profile, the over-determined factor is determined by the following equation:
wherein, Od pAnd representing the overdetermined coefficient corresponding to the dynamic connectivity model among the multilayer wells when the fluid production profile is known, wherein M represents the number of data points, and I' represents the number of water injection wells communicated with the production wells.
Preferably, for a production well without a fluid production profile, the over-determined factor is determined by the following formula:
wherein, Od IAnd (3) representing the overdetermined coefficient corresponding to the dynamic connectivity model among the multilayer wells when the fluid production profile is unknown, wherein L represents the number of the oil reservoir layers, and I' represents the number of the water injection wells communicated with the production wells.
Finally, in step S103, based on the multilayer asymmetric coefficients and the overdetermined coefficients, an optimal overdetermined coefficient is selected to determine the number of data points in the input parameters of the multilayer inter-well dynamic connectivity model, where the data points include production data of production wells in the block within a preset time interval.
Preferably, step S103 specifically includes the following steps:
and drawing a relation curve graph by taking the overdetermined coefficient as a horizontal coordinate and the multilayer asymmetric coefficient as a vertical coordinate, and determining the optimal overdetermined coefficient through the relation curve graph.
The method shown in figure 1 can fully utilize existing field data, invert the dynamic connectivity among the wells of the multi-layer oil reservoir as accurately as possible, and truly reflect the communication coefficient among injection wells and production wells in the stratum. The data change which can be provided by the oil reservoir site in the solving process is fully considered, the input parameters are analyzed and processed, and the method has important reference value for making the later development scheme of the oil field in the high water-cut period.
FIG. 2 shows a plan view of an injection-production well for a five-injection four-production model. Wherein Inj1, Inj2, Inj3, Inj4 and Inj5 are water injection wells, and Pro1, Pro2, Pro3 and Pro4 are production wells. According to the position relation and the space distance of the nine wells, the injection-production well pair can be divided into three types, namely:
the first type: correspondence between the intermediate water injection wells Inj1 and the adjacent four production wells (Pro1, Pro2, Pro3, Pro 4). At this time, the distances between the water injection well Inj1 and the adjacent four production wells are equal, and the four production wells are uniformly distributed, so the expected value of the communication coefficient of Inj1 to the adjacent four production wells is 0.25, that is, the injected water is uniformly distributed to the periphery. This category has a total of 4 injection and production well pairs.
The second type: and the corresponding relation between the water injection wells on the four corners and the two adjacent production wells. For example, Inj2 corresponds to Pro1 and Pro2, Inj3 corresponds to Pro1 and Pro3, Inj4 corresponds to Pro2 and Pro4, and Inj5 corresponds to Pro3 and Pro 4. This category has a total of 8 injection-production well pairs.
In the third category: the corresponding relation between the water injection wells on the four corners and the two distant production wells. For example, Inj2 corresponds to Pro3 and Pro4, Inj3 corresponds to Pro2 and Pro4, Inj4 corresponds to Pro1 and Pro3, and Inj5 corresponds to Pro1 and Pro 2. This category has a total of 8 injection-production well pairs.
As shown in FIG. 2, the production wells are symmetrically distributed and the distances between the injection and production wells in the above three types are the same, so that the inter-well dynamic communication coefficients of the injection and production well relation pairs in the above three types are the same.
FIG. 3 shows a flow chart of a method for analyzing input parameters of a multi-layer interwell dynamic connectivity model according to another embodiment of the invention.
Firstly, acquiring the data of oil deposit production wells in a current block, then determining the distribution mode of the production wells, and when the wells in the block adopt a five-injection four-extraction model as shown in figure 2: the asymmetry coefficients as well as the overdetermined coefficients need to be determined.
When the asymmetric coefficient is determined, in order to verify the accuracy of the solution result of the multi-layer inter-well dynamic connectivity model, the single-layer asymmetric coefficient A of the inter-well communication coefficient is determined, and the coefficient can accurately reflect the difference between the solution result and the ideal result.
Based on the above analysis, the single layer asymmetry factor A is defined as follows.
In the formula, var1(λ) -the variance of the inter-well dynamic connectivity coefficients between the central water injection well and the adjacent four production wells;
var2(lambda) -the variance of the inter-well dynamic communication coefficients between the four corner injection wells and the two adjacent production wells, respectively;
var3(lambda) -variance of inter-well dynamic communication coefficients between the four corner injection wells and the two more distant production wells, respectively.
According to the definition of the single-layer asymmetric coefficient A, the larger the value of the single-layer asymmetric coefficient A is, the larger the deviation of the calculation result from the expected value is, and the more inaccurate the calculation result is; the smaller the value of the calculation result is, the smaller the deviation of the calculation result from the expected value is, and the more accurate the calculation result is.
However, the parameter is only suitable for evaluating the solving precision of the dynamic connectivity model between single-layer oil reservoir wells and cannot meet the requirement of a multi-layer oil reservoir, so that a new parameter is required to be established for evaluating the solving precision of the dynamic connectivity model between the multi-layer oil reservoir wells.
On the basis of the single-layer asymmetric coefficient A, the multilayer asymmetric coefficient A' is formed by improving the multilayer oil deposit according to the characteristics of the multilayer oil deposit.
In the formula, L represents the number of oil reservoir layers;
vari1(λ) — the variance of the inter-well dynamic communication coefficients between the central water injection well of the ith zone and the adjacent four production wells;
vari2(lambda) -the variance of the inter-well dynamic communication coefficients between the water injection wells at the four corners of the ith layer and two adjacent production wells respectively;
vari3(lambda) -the variance of the inter-well dynamic communication coefficients between the water injection wells at the four corners of the ith layer and the two more distant production wells respectively.
The multilayer asymmetric coefficients eliminate the influence of the number of oil deposit layers on the result, so that the method is suitable for multilayer oil deposits, and the calculation results of the multilayer oil deposits can be compared in the same dimension.
When the over-determined coefficient is calculated, the ratio of the number of data points in the system to the number of unknown parameters in the model is defined as the over-determined coefficient, because the number of the unknown parameters corresponding to different models is different. The number of the unknown parameters of the dynamic connectivity model among the wells of the multilayer oil reservoir has great relation with the number of water injection well openings, the number of oil extraction well openings and the number of the oil reservoir layers of the whole oil reservoir, the number of the unknown parameters of the model is continuously increased along with the increase of the number of the production wells and the number of the oil reservoir layers, and the number of required data points is also continuously increased.
For a production well with a production profile, the overdetermined coefficients are defined as follows:
because the liquid production amount of each layer of the oil production well can be fitted in the process of solving the dynamic connectivity model among the multiple layers of wells, the over-determined coefficient corresponding to the model is not so sensitive to the number of the oil deposit layers.
In the formula, Od PAn overdetermined coefficient corresponding to a dynamic connectivity model among the multilayer wells when the fluid production profile is known;
M is the number of data points;
i' -the number of injection wells in communication with the production well.
For a production well without a fluid production profile, the overdetermined coefficients are defined as follows:
in the formula, Od IAn overdetermined coefficient corresponding to a dynamic connectivity model among the multilayer wells when the fluid production profile is unknown;
m is the number of data points;
l is the number of oil reservoir layers;
i' -the number of injection wells in communication with the production wells.
The overdetermined coefficient of the multi-layer interwell dynamic connectivity model has strong sensitivity to the layer number of the oil reservoir, and the difficulty of model solution is obviously increased due to the increase of the layer number.
The overdetermined coefficients under different conditions are obtained through data such as the multilayer oil deposit liquid production capacity, and the corresponding multilayer asymmetric coefficients are obtained by inverting the dynamic connectivity among the multilayer wells through the data points with different numbers. The overdetermined coefficient is used as an abscissa, the multilayer asymmetric coefficient is used as an ordinate, a relation curve of the overdetermined coefficient and the multilayer asymmetric coefficient is obtained, the value of the overdetermined coefficient is determined through an inflection point of the curve, and the minimum input quantity of data points needing to be input is further determined.
The method shown in FIG. 3 provides multilayer asymmetric coefficients, determines the minimum input quantity of model data points, and ensures the calculation accuracy of the multilayer inter-well dynamic connectivity model. Meanwhile, the operability is strong, the application range is wide, and the practical application of the oil field is facilitated.
In a fruitIn the embodiment, a three-layer five-injection four-extraction numerical simulation model is established, the distribution mode is shown as figure 2, the main parameters of the model are shown as table 1, and the horizontal permeability of the upper, middle and lower three-layer reservoirs is 1000 multiplied by 10 respectively-3μm2、500×10-3μm2、100×10-3μm2。
TABLE 1 three-layer five-injection four-extraction numerical simulation model principal parameters
The two models are used for continuously fitting the liquid production rate, because the multilayer model when the liquid production profile is unknown needs to fit the liquid production rate of each production well, and the multilayer model when the liquid production profile is known needs to fit the liquid production rate of each layer of each production well, the dynamic connectivity model between the multilayer wells when the liquid production profile is unknown is taken as an example for brief explanation. The fitting effect of each production well is shown in fig. 4-7.
Simulating different time lengths, namely different data points, by using a numerical model to obtain different over-determined coefficients, inverting the dynamic connectivity among the multiple layers of wells by using the data points with different numbers to obtain corresponding multilayer asymmetric coefficients A ', taking the over-determined coefficients as horizontal coordinates and the multilayer reservoir asymmetric coefficients A' as vertical coordinates to obtain a relation curve between the two, and explaining different dynamic connectivity models among the multiple layers of wells:
(1) for production wells with production profiles
As shown in fig. 8, for the multilayer interwell dynamic connectivity model when the fluid production profile is known, the asymmetry coefficient (multilayer asymmetry coefficient) is large when the over-determined coefficient is less than 7, and the asymmetry coefficient is small when the over-determined coefficient is greater than 7, and can be maintained at a low level. Therefore, when the model is applied, the accuracy of the model solution result can be ensured only by ensuring that the over-determined coefficient is not less than 7.
(2) For production wells without fluid production profile
As shown in fig. 9, for the multilayer interwell dynamic connectivity model when the fluid production profile is unknown, the asymmetry coefficient (multilayer asymmetry coefficient) is large when the over-determined coefficient is less than 6, and the asymmetry coefficient is small when the over-determined coefficient is greater than 6, and can be maintained at a low level. Therefore, when the model is applied, the accuracy of the model solving result can be ensured only by ensuring that the over-determined coefficient is not less than 6.
FIG. 10 is a block diagram of an input parameter analysis apparatus for a multi-layer inter-well dynamic connectivity model according to an embodiment of the invention. As shown in fig. 10, the input parameter analyzing apparatus 1000 includes a first module 1001, a second module 1002, and a third module 1003.
The first module 1001 is configured to determine distribution manners of the production wells and the water injection wells in the current block, and calculate inter-well dynamic communication coefficients between the water injection wells and the production wells in the current block and multi-layer asymmetric coefficients corresponding to the inter-well dynamic communication coefficients based on a multi-layer inter-well dynamic communication model.
Preferably, the first module 1001 comprises:
and the single-layer asymmetric coefficient unit is used for calculating to obtain a single-layer asymmetric coefficient in a single oil reservoir layer based on the inter-well dynamic communication coefficient and a five-injection four-extraction model.
And the multilayer asymmetric coefficient unit is used for calculating multilayer asymmetric coefficients in a multilayer oil reservoir layer by combining a five-injection four-extraction model based on single-layer asymmetric coefficients.
The second module 1002 is configured to determine an over-determined coefficient corresponding to a multi-layer inter-well dynamic connectivity model based on production conditions and production data of production wells within the block.
The third module 1003 is configured to select an optimal overdetermined coefficient based on the multi-layer asymmetric coefficients and the overdetermined coefficients to determine the number of data points in the input parameters of the multi-layer inter-well dynamic connectivity model, where the data points include production data of production wells in a block within a preset time interval.
It should be noted that the multi-layer inter-well dynamic connectivity model mentioned in the present invention can be as follows:
wherein q isjRepresents the predicted fluid production of the jth production well, qojRepresenting the injection-production imbalance coefficient, τjDenotes the time constant, NIIndicates the number of injection wells in the block, lambdaijRepresenting the inter-well dynamic communication coefficient, I 'between the ith water injection well and the jth production well' ijIndicates the water injection quantity from the i-th water injection well to the j-th production well, NpIndicates the number of production wells in the block, vljShowing the influence coefficient of the first production well on the liquid production capacity of the j production well,
represents the bottom hole flow pressure of the jth production well, IiRepresenting the water injection amount of the ith production well, n and m represent the time, n0Representing an initial time instant and deltan representing a time interval.
The dynamic connectivity model among the multiple layers of wells takes the minimum sum of squares of the difference between the predicted liquid production amount and the actual liquid production amount of the production wells in the block as an objective function, and the objective function is as follows:
wherein N istWhich represents the number of sampling time steps,
the actual liquid production of the jth production well is shown, and t is the liquid production time.
The constraints of the multi-layer interwell dynamic connectivity model are as follows:
when the production wells in the block are in a constant pressure production state, the dynamic connectivity model among the multilayer wells is as follows:
in conclusion, the input parameter analysis method and device for the multilayer inter-well dynamic connectivity model provided by the invention determine the most important input parameters (water injection amount of a water injection well and oil production amount of a production well) of the multilayer inter-well dynamic connectivity model by comparing the asymmetric coefficient of the inter-well dynamic connectivity coefficient with the over-determined coefficient of the production well, so as to determine the final inversion result of the model and reasonably reflect the formation information. The invention can comprehensively and accurately reflect the formation information under the condition of limited data points, and simultaneously, the proper amount of selected data points can also truly reflect the communication coefficient between the current formation injection and production wells, thereby fully considering the problem of actual parameter determination of the oil field site and having important reference value for the formulation of the later development scheme of the oil field in the high water cut period.
It is to be understood that the disclosed embodiments of the invention are not limited to the particular structures, process steps, or materials disclosed herein but are extended to equivalents thereof as would be understood by those ordinarily skilled in the relevant arts. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting.
Reference in the specification to "one embodiment" or "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment of the invention. Thus, the appearances of the phrase "one embodiment" or "an embodiment" in various places throughout this specification are not necessarily all referring to the same embodiment.
Although the embodiments of the present invention have been described above, the above description is only for the convenience of understanding the present invention, and is not intended to limit the present invention. It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (6)
1. A method for analyzing input parameters of a multi-layer interwell dynamic connectivity model, the method comprising the steps of:
determining distribution modes of production wells and water injection wells in a current block, and calculating inter-well dynamic communication coefficients between the water injection wells and the production wells in the current block and multilayer asymmetric coefficients corresponding to the inter-well dynamic communication coefficients based on a multilayer inter-well dynamic communication model;
secondly, determining an overdetermined coefficient corresponding to the multi-layer inter-well dynamic connectivity model based on the production condition and the production data of the production wells in the block;
selecting an optimal overdetermined coefficient based on the multilayer asymmetric coefficients and the overdetermined coefficient to determine the number of data points in input parameters of the multilayer inter-well dynamic connectivity model, wherein the data points comprise production data of production wells in a block in a preset time interval;
when the distribution mode of the production wells and the water injection wells in the block is a multilayer five-injection four-extraction model, in the multilayer oil reservoir layer, the multilayer asymmetric coefficients are calculated and obtained by combining the five-injection four-extraction model based on the single-layer asymmetric coefficients corresponding to the single oil reservoir layer through the following formula:
Wherein A' represents the multilayer asymmetric coefficient, L represents the number of reservoir layers, vari1(λ) represents the variance of the inter-well dynamic communication coefficients between the central water injection well in the ith layer and the adjacent four production wells, vari2(lambda) represents the variance of the inter-well dynamic communication coefficients between the water injection wells at the four corners of the ith layer and the two adjacent production wells, vari3(lambda) represents the variance of the inter-well dynamic communication coefficients between the water injection wells at the four corners of the ith layer and the two more distant production wells respectively;
judging whether the production well in the block has a liquid production profile or not based on the production condition of the production well in the block, and determining the over-determined coefficient by the following formula aiming at the production well with the liquid production profile:
wherein, Od pRepresenting the overdetermined coefficient corresponding to the dynamic connectivity model among the multilayer wells when the fluid production profile is known, wherein M represents the number of the data points, and I' represents the number of water injection wells communicated with the production wells;
for a production well that does not have a fluid production profile, the over-determined factor is determined by the following equation:
wherein, Od IAnd expressing the overdetermined coefficient corresponding to the dynamic connectivity model among the multilayer wells when the fluid production profile is unknown, and L expresses the number of the oil reservoir layers.
2. The method of claim 1, wherein when the production wells and the water injection wells within the block are distributed in a multi-zone five-injection four-production model, the first step comprises the following steps:
And in a single oil reservoir layer, calculating to obtain the single-layer asymmetric coefficient based on the inter-well dynamic communication coefficient and the five-injection four-extraction model.
3. The method of claim 2, wherein the single layer asymmetry factor is calculated by the formula:
wherein A represents the single layer asymmetry coefficient, var1(λ) represents the variance of the inter-well dynamic communication coefficients between the central water injection well and the adjacent four production wells, var2(λ) represents the variance of the inter-well dynamic communication coefficients between the four-corner injection wells and the two adjacent production wells, var3And (lambda) represents the variance of the inter-well dynamic communication coefficients between the water injection wells at the four corners and the two production wells farther away.
4. The method according to any one of claims 1 to 3, wherein step three comprises in particular the steps of:
and drawing a relation curve graph by taking the over-definite coefficient as a horizontal coordinate and the multilayer asymmetric coefficient as a vertical coordinate, and determining the optimal over-definite coefficient through the relation curve graph.
5. An input parameter analysis apparatus for a multi-layer inter-well dynamic connectivity model, wherein the method of any one of claims 1-4 is performed, the apparatus comprising:
The first module is used for determining the distribution mode of the production wells and the water injection wells in the current block, and calculating inter-well dynamic communication coefficients between the water injection wells and the production wells in the current block and multilayer asymmetric coefficients corresponding to the inter-well dynamic communication coefficients based on a multilayer inter-well dynamic communication model;
a second module for determining an over-determined coefficient corresponding to the multi-layer inter-well dynamic connectivity model based on production conditions of production wells within the block and production data;
a third module, configured to select an optimal overdetermined coefficient based on the multilayer asymmetric coefficients and the overdetermined coefficient to determine the number of data points in input parameters of the multilayer inter-well dynamic connectivity model, where the data points include production data of production wells in a block within a preset time interval;
when the distribution mode of the production wells and the water injection wells in the block is a multilayer five-injection four-extraction model, the first module comprises:
and the multilayer asymmetric coefficient unit is used for calculating the multilayer asymmetric coefficient by combining the five-injection four-extraction model based on the single-layer asymmetric coefficient corresponding to the single oil reservoir layer in the multilayer oil reservoir layer.
6. The apparatus of claim 5, wherein when the production wells and water injection wells within the block are distributed in a multi-zone five-injection four-production model, the first module comprises:
and the single-layer asymmetric coefficient unit is used for calculating and obtaining the single-layer asymmetric coefficient in a single oil reservoir layer based on the inter-well dynamic communication coefficient and the five-injection four-extraction model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910547438.4A CN110410065B (en) | 2019-06-24 | 2019-06-24 | Input parameter analysis method and device for dynamic connectivity model among multiple wells |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910547438.4A CN110410065B (en) | 2019-06-24 | 2019-06-24 | Input parameter analysis method and device for dynamic connectivity model among multiple wells |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110410065A CN110410065A (en) | 2019-11-05 |
| CN110410065B true CN110410065B (en) | 2022-05-03 |
Family
ID=68359553
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910547438.4A Active CN110410065B (en) | 2019-06-24 | 2019-06-24 | Input parameter analysis method and device for dynamic connectivity model among multiple wells |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110410065B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111506975B (en) * | 2020-04-10 | 2024-09-24 | 中国石油化工股份有限公司 | Method, device and computing equipment for judging communication relation between wells |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101725346A (en) * | 2009-12-15 | 2010-06-09 | 中国石油大学(华东) | Oil deposit inter-well dynamic connectivity inverting method |
| CN103670369A (en) * | 2013-12-12 | 2014-03-26 | 中国石油天然气股份有限公司 | Method and device for judging communication condition between injection wells and production wells |
| CN105089649A (en) * | 2014-04-17 | 2015-11-25 | 中国石油化工股份有限公司 | Method for establishing inter-well connectivity model |
| CN106837297A (en) * | 2016-12-22 | 2017-06-13 | 中国石油天然气股份有限公司 | Method for identifying connectivity among wells and predicting oil-water dynamic state |
| CN106845786A (en) * | 2016-12-26 | 2017-06-13 | 中国石油天然气股份有限公司 | Method and device for evaluating reservoir connectivity between injection well and production well |
| CN107291667A (en) * | 2016-04-01 | 2017-10-24 | 中国石油化工股份有限公司 | A kind of interwell communication degree determines method and system |
| CN108764760A (en) * | 2018-06-25 | 2018-11-06 | 成都北方石油勘探开发技术有限公司 | A kind of inter well connectivity analysis and method for early warning based on data mining |
| CN109447532A (en) * | 2018-12-28 | 2019-03-08 | 中国石油大学(华东) | A kind of oil reservoir inter well connectivity based on data-driven determines method |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8386226B2 (en) * | 2009-11-25 | 2013-02-26 | Halliburton Energy Services, Inc. | Probabilistic simulation of subterranean fracture propagation |
| GB2544922B (en) * | 2014-08-01 | 2022-01-12 | Geoquest Systems Bv | Methods and systems for simulating a hydrocarbon field using a multi-point well connection method |
| US10097404B2 (en) * | 2014-09-16 | 2018-10-09 | CloudGenix, Inc. | Methods and systems for time-based application domain classification and mapping |
| CN107038268B (en) * | 2016-10-13 | 2022-02-18 | 中国石油化工股份有限公司 | Method for determining water flooding wave sum coefficient of heterogeneous reservoir five-point well pattern |
-
2019
- 2019-06-24 CN CN201910547438.4A patent/CN110410065B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101725346A (en) * | 2009-12-15 | 2010-06-09 | 中国石油大学(华东) | Oil deposit inter-well dynamic connectivity inverting method |
| CN103670369A (en) * | 2013-12-12 | 2014-03-26 | 中国石油天然气股份有限公司 | Method and device for judging communication condition between injection wells and production wells |
| CN105089649A (en) * | 2014-04-17 | 2015-11-25 | 中国石油化工股份有限公司 | Method for establishing inter-well connectivity model |
| CN107291667A (en) * | 2016-04-01 | 2017-10-24 | 中国石油化工股份有限公司 | A kind of interwell communication degree determines method and system |
| CN106837297A (en) * | 2016-12-22 | 2017-06-13 | 中国石油天然气股份有限公司 | Method for identifying connectivity among wells and predicting oil-water dynamic state |
| CN106845786A (en) * | 2016-12-26 | 2017-06-13 | 中国石油天然气股份有限公司 | Method and device for evaluating reservoir connectivity between injection well and production well |
| CN108764760A (en) * | 2018-06-25 | 2018-11-06 | 成都北方石油勘探开发技术有限公司 | A kind of inter well connectivity analysis and method for early warning based on data mining |
| CN109447532A (en) * | 2018-12-28 | 2019-03-08 | 中国石油大学(华东) | A kind of oil reservoir inter well connectivity based on data-driven determines method |
Non-Patent Citations (2)
| Title |
|---|
| Inferring Interwell Connectivity from Well Bottomhole-Pressure Fluctuations in Waterfloods;Dihn Anh等;《SPE Reservoir Evaluation & Engineering》;20081231;第874-881页 * |
| Inferring interwell connectivity only from well-rate fluctuations in waterfloods;Albertoni,A.等;《SPE Reservoir Evaluation & Engineering》;20031231;第11-16页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110410065A (en) | 2019-11-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2549147C2 (en) | Systems, computer-implemented methods and computer-readable program products to compute approximate well drainage pressure for reservoir simulator | |
| CN105089649B (en) | A kind of inter well connectivity method for establishing model | |
| CN104060985B (en) | A kind of stratified deposit profile control and water plugging blocking agent penetration depth method of testing and system | |
| CN107291667B (en) | Method and system for determining communication degree between wells | |
| RU2709047C1 (en) | Method of adaptation of hydrodynamic model of productive formation of oil and gas condensate deposit taking into account uncertainty of geological structure | |
| CN110414048A (en) | Inter well connectivity analysis method and device | |
| CN110410065B (en) | Input parameter analysis method and device for dynamic connectivity model among multiple wells | |
| CN106522921B (en) | The stochastic modeling method and device of dynamic constrained | |
| Gao | Roughness and demand estimation in water distribution networks using head loss adjustment | |
| CN110566196B (en) | Reservoir connectivity analysis method | |
| CN112302594A (en) | Method, device and equipment for determining connection structure of heterogeneous water-drive reservoir | |
| CA3116482A1 (en) | Reservoir fluid property modeling using machine learning | |
| CA3095746C (en) | Optimized methodology for automatic history matching of a petroleum reservoir model with ensemble kalman filter | |
| CN113032996B (en) | Water channeling channel identification method for hypotonic fractured reservoir horizontal well | |
| CN112464552B (en) | Method and system suitable for predicting seepage field time sequence of underground water-sealed oil storage | |
| CN110991084B (en) | Reservoir permeability calculation method based on streamline numerical value well test | |
| CN112554854B (en) | System and method for determining water absorption profile | |
| CN114519260B (en) | Shale oil productivity prediction method, system, medium, equipment and terminal | |
| CN115062552B (en) | Method and system for predicting nitrogen huff and puff effect of fracture-cavity oil reservoir | |
| CN110135122B (en) | Low-permeability natural gas reservoir yield prediction method and system | |
| CN118296975B (en) | Inter-well connectivity prediction method based on knowledge interaction graph neural network | |
| CN114021297B (en) | Complex pipe network leakage positioning method based on echo state network | |
| US11808148B2 (en) | Systems and methods for back-allocation of oil produced by waterflooding | |
| CN113177319A (en) | Method and system for judging connectivity among wells based on neural network sensitivity analysis | |
| CN117468911A (en) | High-water saturation multi-layer production tight sandstone gas reservoir development potential evaluation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |