CN113738319A - Water drive invalid cycle identification method, device and medium based on Lorentz curve - Google Patents

Water drive invalid cycle identification method, device and medium based on Lorentz curve Download PDF

Info

Publication number
CN113738319A
CN113738319A CN202111040985.7A CN202111040985A CN113738319A CN 113738319 A CN113738319 A CN 113738319A CN 202111040985 A CN202111040985 A CN 202111040985A CN 113738319 A CN113738319 A CN 113738319A
Authority
CN
China
Prior art keywords
specified
water injection
well
coefficient
layer
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
CN202111040985.7A
Other languages
Chinese (zh)
Other versions
CN113738319B (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.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN202111040985.7A priority Critical patent/CN113738319B/en
Publication of CN113738319A publication Critical patent/CN113738319A/en
Application granted granted Critical
Publication of CN113738319B publication Critical patent/CN113738319B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
    • E21B43/20—Displacing by water
    • 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
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/14—Obtaining from a multiple-zone well
    • 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
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/30—Specific pattern of wells, e.g. optimising the spacing of wells
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Fluid Mechanics (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The embodiment of the specification discloses a method, a device and a medium for identifying water drive ineffective circulation based on a Lorentz curve, wherein the method comprises the steps of fully mastering the communication rule among wells of an oil reservoir, and calculating the water absorption capacity of each small layer of an injection and production unit to be researched according to the physical property parameters and development dynamics of the oil reservoir; and then, determining a water absorption coefficient of variation based on the Lorentz curve, judging the water flooding effect according to the water absorption coefficient of variation, and judging the invalid circulation horizon according to a critical value of the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient. And adjusting according to the water-drive effect and the identified ineffective circulation horizon, and then evaluating after adjusting by utilizing a Lorentz curve and a water-drive ineffective circulation horizon determination method, thereby laying a foundation for better guiding the water-drive ineffective circulation identification and measure adjustment work of the oil field.

Description

Water drive invalid cycle identification method, device and medium based on Lorentz curve
Technical Field
The specification relates to the technical field of oil reservoir development, in particular to a method, a device and a medium for identifying water drive invalid circulation based on a Lorentz curve.
Background
The sandstone oil reservoir water injection development is a main technology of oil field development, and can make most oil fields enter a high water-cut period, even an ultrahigh water-cut period, and because the sandstone oil reservoir has reservoir heterogeneity, the water injection amount of the water injection well at different positions and different directions has obvious difference, the single-layer water injection amount heterogeneity degree is high, the single-layer outburst is serious, the injected water is low-efficiency and ineffective in circulation, and the oil stabilization and water control become the focus problem of the high water-cut oil field development. Therefore, the water drive quantitative characterization and the water drive low-efficiency ineffective cycle identification of the sandstone reservoir are the key points for improving the current situation of high water content of the reservoir to realize water control and yield increase.
At present, the identification methods for water injection low-efficiency or ineffective circulation at home and abroad mainly comprise a direct observation method, an interwell tracer method, a well testing data method, a well logging data interpretation method, a production data method, a method combining mathematics and oil reservoir engineering and the like, but each identification method not only needs to consume a large amount of time and material cost, but also can only realize the qualitative analysis of the water drive effect of the sandstone oil reservoir, realize the quantitative representation of the water drive development effect, and efficiently and accurately identify a water injection low-efficiency or ineffective circulation channel, which is still a great technical problem in the field of oil and gas field development.
Disclosure of Invention
The embodiment of the specification aims to provide a method, a device and a medium for identifying water drive invalid circulation based on a Lorentz curve, which can more accurately determine the development effect and the invalid circulation horizon and lay a foundation for better guiding the identification of the water drive invalid circulation of an oil field and the adjustment of measures.
The specification provides a method, a device and a medium for identifying water drive invalid circulation based on a Lorentz curve, which are realized by the following modes:
the specification provides a water drive invalid circulation identification method based on a Lorentz curve, which comprises the following steps: calculating single-layer unidirectional water injection quantity distributed to the direction of the specified oil well in the specified small layer by the specified water injection well according to the single-well water injection quantity, the plane distribution coefficient and the longitudinal distribution coefficient of the specified water injection well in the target oil reservoir; wherein the plane distribution coefficient is determined according to the seepage resistance coefficient of the specified small layer between the specified water injection well and the specified oil well and the pressure of the specified water injection well in the specified small layer; the longitudinal distribution coefficient is determined at least according to the seepage resistance coefficient, the average permeability and the reservoir transformation coefficient of the specified small layer between the specified water injection well and the specified oil well; calculating the ratio of the single-layer unidirectional water injection quantity to the total water injection quantity to be used as a dimensionless water injection coefficient of the specified water injection well in the direction from a specified small layer to a specified oil well; the total water injection amount is the accumulated value of the water injection amount of each water injection well; calculating the ratio of the average horizon thickness of the specified small layer between the specified water injection well and the specified oil well to the total horizon thickness as a dimensionless effective thickness coefficient of the specified small layer between the specified water injection well and the specified oil well; the total layer thickness comprises an accumulated value of the layer average thickness corresponding to each appointed small layer; and determining whether a seepage channel of the specified water injection well in the specified small layer to the direction of the specified oil well is an invalid circulation channel or not according to the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient.
Optionally, the longitudinal distribution coefficient is determined in a coefficient manner:
Figure BDA0003249189700000021
wherein, Ci,kRepresents the longitudinal distribution coefficient, H, of a given water injection well i to a given sub-zone ki,j,kRepresenting the average thickness of the horizon of the specified small layer K between the specified water injection well i and the specified oil well j, Ki,j,kRepresents the average permeability, M, of the given fraction k between the given water injection well i and the given oil well ji,j,kRepresenting the reservoir transformation coefficient, R, of the specified small layer k between the specified water injection well i and the specified oil well ji,j,kRepresenting the seepage resistance coefficient of the specified small layer k between the specified water injection well i and the specified oil well j;
Figure BDA0003249189700000022
to specify the number of wells in the small layer k, S represents the number of small layers obtained after the target reservoir formation is partitioned.
Optionally, the seepage resistance coefficient of the specified water injection well i in the direction from the specified small layer k to the specified oil well j is determined by using the following calculation model:
Figure BDA0003249189700000023
Figure BDA0003249189700000024
wherein λ isi,j,k、Li,j,kRespectively the apparent viscosity and the effective well spacing between a designated water injection well i and a designated oil well j of a designated small layer Kro、Krw、μo、μwThe relative permeability of the oil phase, the relative permeability of the water phase, the viscosity of the crude oil and the viscosity of water are respectively.
Optionally, the method further includes: determining the water absorption capacity variation coefficient of the target oil reservoir by using the following calculation model:
Figure BDA0003249189700000031
wherein V is the coefficient of variation of water absorption, WI,J,K、ZI,J,KRespectively representing a dimensionless water injection coefficient and a dimensionless effective thickness coefficient, I, K, J representing numbers corresponding to a water injection well, an oil well and a small layer, and I, K, J values are configured in a circular nesting mode; s represents the number of small layers obtained after the target oil reservoir stratum is divided, N represents the number of oil wells in the target oil reservoir, and M represents the number of water injection wells in the target oil reservoir; n, M and S represent values of I, J, K in the sequential accumulation process, M is more than or equal to 1 and less than or equal to M, N is more than or equal to 1 and less than or equal to N, and S is more than or equal to 1 and less than or equal to S; wherein, the loop nesting mode comprises the following steps: taking I as an outermost layer, J as a secondary outer layer and K as an innermost layer, wherein the cycle corresponding to each layer is that the value of the corresponding element of each layer is sequentially increased from 1 to the maximum value of the corresponding element according to the sequence of the dimensionless water injection coefficient from small to large, and the step length of each increase is 1;
and analyzing the development effect of the target oil reservoir by using the water absorption variation coefficient.
Optionally, under the condition that the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient is greater than an invalid circulation threshold, determining that a seepage channel of the specified water injection well in the specified small layer towards the specified oil well is an invalid circulation channel; and determining the invalid circulation threshold according to a data set consisting of the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient before the target oil reservoir is plugged and reformed.
Optionally, the method further includes: plugging or fracturing the target oil reservoir according to the ineffective circulation channel; and updating the reservoir transformation coefficients of the specified small layers of the target oil reservoir between the specified water injection well and the specified oil well after plugging or fracturing treatment so as to identify the invalid circulation channel of the target oil reservoir after the current plugging by using the updated reservoir transformation coefficients.
In another aspect, the present specification further provides a water drive invalid circulation identification apparatus based on a lorentz curve, including: the water injection rate determining module is used for calculating single-layer unidirectional water injection rate distributed to the direction of the specified oil well by the specified water injection well in the specified small layer according to the single-well water injection rate, the plane distribution coefficient and the longitudinal distribution coefficient of the specified water injection well in the target oil reservoir; wherein the plane distribution coefficient is determined according to the seepage resistance coefficient of the specified small layer between the specified water injection well and the specified oil well and the pressure of the specified water injection well in the specified small layer; the longitudinal distribution coefficient is determined at least according to the seepage resistance coefficient, the average permeability and the reservoir transformation coefficient of the specified small layer between the specified water injection well and the specified oil well; the water injection coefficient calculation module is used for calculating the ratio of the single-layer unidirectional water injection quantity to the total water injection quantity to be used as a dimensionless water injection coefficient of the specified water injection well in the direction of the specified oil well in the specified small layer; the total water injection amount is the accumulated value of the water injection amount of each water injection well; the thickness coefficient calculation module is used for calculating the ratio of the average thickness of the layer position of the specified small layer between the specified water injection well and the specified oil well to the total layer position thickness as a dimensionless effective thickness coefficient of the specified small layer between the specified water injection well and the specified oil well; the total layer thickness comprises an accumulated value of the layer average thickness corresponding to each appointed small layer; and the invalid circulation channel identification module is used for identifying whether the seepage channel of the specified water injection well in the specified small layer to the specified oil well direction is an invalid circulation channel according to the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient.
Optionally, the apparatus further includes a vertical distribution coefficient calculation module, where the vertical distribution coefficient calculation module is configured to calculate the vertical distribution coefficient by using the following method:
Figure BDA0003249189700000041
wherein, Ci,kRepresents the longitudinal distribution coefficient, H, of a given water injection well i vertically to a given small layer ki,j,kRepresenting the average thickness of the horizon of the specified small layer K between the specified water injection well i and the specified oil well j, Ki,j,kRepresents the average permeability, M, of the given fraction k between the given water injection well i and the given oil well ji,j,kRepresenting the reservoir transformation coefficient, R, of the specified small layer k between the specified water injection well i and the specified oil well ji,j,kRepresenting the seepage resistance coefficient of the specified small layer k between the specified water injection well i and the specified oil well j;
Figure BDA0003249189700000042
to specify the number of wells in the small layer k, S represents the number of small layers obtained after the target reservoir formation is partitioned.
Optionally, the apparatus further includes a seepage resistance coefficient calculation module, where the seepage resistance coefficient calculation module is configured to calculate the seepage resistance coefficient by using the following method:
Figure BDA0003249189700000043
Figure BDA0003249189700000044
wherein λ isi,j,k、Li,j,kRespectively the apparent viscosity and the effective well spacing between a designated water injection well i and a designated oil well j of a designated small layer Kro、Krw、μo、μwThe relative permeability of the oil phase, the relative permeability of the water phase, the viscosity of the crude oil and the viscosity of water are respectively.
In another aspect, the present specification also provides a computer readable storage medium having stored thereon computer instructions which, when executed, implement the steps of the method of any one or more of the above embodiments.
According to the method, the device and the medium for identifying the water drive ineffective cycle based on the Lorentz curve, which are provided by one or more embodiments of the specification, the water absorption capacity of each small layer of the researched injection and production unit is calculated by fully mastering the communication rule among wells of an oil reservoir and according to the physical property parameters and development dynamics of the actual oil reservoir; and then, determining a water absorption coefficient of variation based on a Lorentz curve, judging a water drive effect according to the water absorption coefficient of variation, and judging an invalid circulation layer position according to a ratio of the dimensionless water injection coefficient and the dimensionless effective thickness coefficient. Meanwhile, a reservoir transformation coefficient is further determined by combining the reservoir transformation degree, and the communication rule of the transformed reservoir is indirectly represented by combining the reservoir transformation coefficient, so that the reservoir development effect and whether the seepage channel is an ineffective circulation channel or not are evaluated, the determination and evaluation of the reservoir ineffective circulation channel can be more in line with the modified reservoir communication rule, the understanding of the reservoir water injection development effect is further improved, and the design of the reservoir water injection development scheme is better guided.
Drawings
In order to more clearly illustrate the embodiments of the present specification or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments described in the present specification, and for those skilled in the art, other drawings can be obtained according to the drawings without any creative effort. In the drawings:
FIG. 1 is a schematic flow chart diagram illustrating an embodiment of a method for identifying a water flooding invalid cycle provided herein;
FIG. 2 is a schematic diagram of the distribution and communication characteristics of single sand bodies and the numbering of each communication layer provided in the present specification;
FIG. 3 is a schematic diagram of the injection and production well unit distribution provided in the present specification;
FIG. 4 is a graphical representation of pressure data for each of the oil and water injection wells provided herein;
FIG. 5 is a schematic Lorentzian curve provided herein before adjustment;
FIG. 6 is a graphical representation of pressure data for each of the oil and water injection wells provided herein;
FIG. 7 is a graph illustrating an adjusted Lorentzian curve provided herein;
fig. 8 is a schematic block structure diagram of an embodiment of a water flooding invalid cycle identifying apparatus provided in this specification.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in the present specification, the technical solutions in one or more embodiments of the present specification will be clearly and completely described below with reference to the drawings in one or more embodiments of the present specification, and it is obvious that the described embodiments are only a part of the embodiments of the specification, and not all embodiments. All other embodiments obtained by a person skilled in the art based on one or more embodiments of the present specification without making any creative effort shall fall within the protection scope of the embodiments of the present specification.
Fig. 1 is a schematic flow chart of an embodiment of a water drive invalid cycle identification method based on a lorentz curve provided in this specification. Although the present specification provides the method steps or apparatus structures as shown in the following examples or figures, more or less steps or modules may be included in the method or apparatus structures based on conventional or non-inventive efforts. In the case of steps or structures which do not logically have the necessary cause and effect relationship, the execution order of the steps or the block structure of the apparatus is not limited to the execution order or the block structure shown in the embodiments or the drawings of the present specification. When the described method or module structure is applied to a device, a server or an end product in practice, the method or module structure according to the embodiment or the figures may be executed sequentially or in parallel (for example, in a parallel processor or multi-thread processing environment, or even in an implementation environment including distributed processing and server clustering). In a specific embodiment, as shown in fig. 1, in an embodiment of the method for identifying a water flooding invalid cycle based on a lorentz curve provided in this specification, the method may include:
s20: calculating single-layer unidirectional water injection quantity distributed to the direction of the specified oil well in the specified small layer by the specified water injection well according to the single-well water injection quantity, the plane distribution coefficient and the longitudinal distribution coefficient of the specified water injection well in the target oil reservoir; wherein the plane distribution coefficient is determined according to the seepage resistance coefficient of the specified small layer between the specified water injection well and the specified oil well and the pressure of the specified water injection well in the specified small layer; the longitudinal distribution coefficient is determined at least according to the seepage resistance coefficient, the average permeability and the reservoir transformation coefficient of the designated small layer between the designated water injection well and the designated oil well.
Multiple water injection wells and oil wells can be distributed in the target oil reservoir, and the water injection wells can inject water into the stratum so as to recover oil in a water-drive mode. In some embodiments, the reservoir may be divided into a plurality of small layers by analyzing the interbed heterogeneity characteristics of the target reservoir, and each small layer may be analyzed separately. The designated water injection well is any water injection well in the target oil reservoir. The designated oil well is any oil well communicated with the designated water injection well in the target oil reservoir. The designated small layer can be any small layer after the reservoir is divided.
The water injection quantity of each water injection well in each small layer along the direction of the oil well communicated with the water injection well can be calculated according to the physical property parameters and the development dynamics of the actual oil deposit by analyzing the communication rule among all wells of the target oil deposit. In some embodiments, the injection-production pressure difference can be considered, the water injection plane distribution coefficient and the water injection longitudinal distribution coefficient are introduced, and the water injection quantity of each layer in the direction of each affected oil well communicated with the single well water injection quantity distribution can be obtained by combining the production dynamic data.
The single-layer unidirectional water injection quantity distributed to the direction of the specified oil well in the specified small layer of the specified water injection well can be calculated according to the single-well water injection quantity, the plane distribution coefficient and the longitudinal distribution coefficient of the specified water injection well in the target oil reservoir. The planar distribution coefficient may be determined from a seepage resistance coefficient of the designated fraction between the designated water injection well and the designated oil well and a pressure of the designated water injection well at the designated fraction. The longitudinal distribution coefficient may be determined based at least on a seepage resistance coefficient, an average permeability, and a reservoir reformation coefficient of the designated interval between the designated water injection well and the designated oil well.
The vertical distribution coefficient may be determined in the following manner:
Figure BDA0003249189700000071
wherein, Ci,kRepresents the longitudinal distribution coefficient, H, of a given water injection well i vertically to a given small layer ki,j,kRepresenting the average thickness of the horizon of the specified small layer K between the specified water injection well i and the specified oil well j, Ki,j,kRepresents the average permeability, M, of the given fraction k between the given water injection well i and the given oil well ji,j,kRepresenting the reservoir transformation coefficient, R, of the specified small layer k between the specified water injection well i and the specified oil well ji,j,kRepresenting the seepage resistance coefficient of the specified small layer k between the specified water injection well i and the specified oil well j; n is the number of wells in the designated small layer k, and S represents the number of small layers obtained after the target reservoir stratum is divided.
The planar distribution coefficient may be determined in the following manner:
Figure BDA0003249189700000072
αi,j,kdenotes the planar distribution coefficient, Pi,kIndicates the pressure, P, of a given water injection well i at a given interval kj,kIndicating the pressure of the designated injection well j in the designated small layer k, and N indicating the number of wells of the target reservoir.
Wherein, the seepage resistance coefficient of the specified water injection well i in the direction from the specified small layer k to the specified oil well j can be determined by utilizing the following calculation model:
Figure BDA0003249189700000073
Figure BDA0003249189700000074
wherein λ isi,j,k、Li,j,kRespectively the apparent viscosity and the effective well spacing between a designated water injection well i and a designated oil well j of a designated small layer Kro、Krw、μo、μwThe relative permeability of the oil phase, the relative permeability of the water phase, the viscosity of the crude oil,Viscosity of water.
The single-layer unidirectional water injection quantity can be as follows:
Qi,j,k＝αi,j,k×Ci,k×Qwi
wherein Q isi,j,kRepresents the single-layer unidirectional water injection quantity Q distributed from the designated water injection well i to the designated oil well j in the designated small layer kwiIndicating the single well water injection amount of the designated water injection well i.
In the above embodiment, Hi,j,k、Ki,j,k、λi,j,k、Li,j,k、Kro、Krw、μo、μwThe corresponding data may be obtained from the target reservoir test prior to the target reservoir modification, or may be calculated based on data obtained from the target reservoir test. Mi,j,kIs used to characterize the extent and manner of modification to the target reservoir. For example, before reforming a target reservoir, Mi,j,kThe value may be 1. Assuming that the area of the designated small layer k between the designated water injection well i and the designated oil well j is one injection-production unit (i, j, k), M is carried out after the transformation of one injection-production unit (i, j, k) of the target oil reservoiri,j,kThe value of (c) can be determined according to the modification degree and modification mode of the injection-production unit (i, j, k). For example, when the injection-production unit (i, j, k) is subjected to plugging treatment, the seepage resistance of the injection-production unit (i, j, k) should be larger, and the permeability should be smaller; accordingly, the reservoir transformation coefficient M of the injection-production unit (i, j, k) can be seti,j,kIs a value less than 1, and M can be determined according to the plugging ratio of the injection-production cell (i, j, k), the material used, or the likei,j,kThe specific size of (a); the greater the plugging degree, the greater Mi,j,kThe larger the adjustment, i.e. Mi,j,kThe smaller the value of (c). If the injection-production unit (i, j, k) is subjected to fracturing treatment, the seepage resistance and the permeability of the injection-production unit (i, j, k) should be smaller and larger; accordingly, the reservoir transformation coefficient M of the injection-production unit (i, j, k) can be seti,j,kIs a value greater than 1, and M can be determined according to the fracturing ratio, fracturing manner, etc. of the injection-production unit (i, j, k)i,j,kThe specific size of (a); the greater the degree of fracturing, Mi,j,kThe greater the amount of adjustment of (c) is,namely Mi,j,kThe larger the value of (c).
In the simulation analysis process of the overall mining influence of the reservoir transformation mode of the target oil reservoir on the target oil reservoir, parameters such as permeability and the like can only be measured before simulation test, and after the simulation reservoir transformation is carried out on the target oil reservoir, the data of the target oil reservoir after the simulation transformation cannot be obtained through the test, so that the actual influence of the reservoir transformation on the target oil reservoir cannot be accurately given only on the basis of the initially measured parameters such as permeability and the like. In the embodiment of the specification, the influence of reservoir transformation on parameters such as reservoir permeability, reservoir seepage resistance and the like is indirectly represented by using the reservoir transformation coefficient, and on the basis, the invalid circulation channel and the overall development effect after the reservoir transformation are determined, so that the reservoir transformation simulation analysis can better accord with the actual target reservoir development, and the design of a reservoir water injection development scheme can be better guided.
Of course, M is as defined abovei,j,kThe values are preferably given as examples, but do not constitute a direct limitation to the embodiments of the present disclosure, and those skilled in the art can make appropriate adjustments in the above configurations and adjustment modes, but they should be covered within the protection scope of the present disclosure as long as the functions and effects achieved by the present disclosure are the same or similar. For example, Mi,j,kThe initial value of (1) may be other than 1, and other modification methods may be adopted for reservoir modification.
In other embodiments, the communication condition of each small layer of each injection-production well group can be judged in advance according to the geological data of the well group, and the data processing efficiency is further improved when only the connected injection-production units are subjected to invalid circulation channel analysis. The left side of figure 2 shows the single sand body distribution and the intercommunication characteristic in experimental zone, and the right side of figure 2 shows that each intercommunication layer in experimental zone is unified to be numbered. In FIG. 2, I1、I2Indicating a water injection well, P1To P6The production wells are shown, and the (1) to (8) respectively show the corresponding injection and production units between one water injection well and one injection and production well. And (3) judging whether the (1) to (8) are disconnected or connected by referring to the left image of the figure 2 and the related geological materials, wherein the connected value is 1, and the disconnected value is 0. When the invalid circulation channel is analyzed, the method can be only carried out on the injection and production unit with the value 1And performing calculation and analysis to further improve the data processing efficiency. Accordingly, when only the connected injection and production units are subjected to the invalid circulation channel analysis, the water injection wells and the oil wells related to the connected injection and production units in the calculation formulas of the embodiments of the present specification are all the water injection wells and the oil wells related to the connected injection and production units, and the number of the water injection wells and the number of the oil wells in the target oil reservoir are also all the number of the water injection wells and the number of the oil wells related to the connected injection and production units.
Of course, those skilled in the art may also have other modifications to the calculation method of the above-mentioned planar distribution coefficient, longitudinal distribution coefficient and single-layer unidirectional water injection amount, such as configuring weights for each factor in the above calculation process according to the communication rule between wells, etc., according to the technical spirit of the embodiments of the present disclosure. But all that is required is to cover the scope of the present specification as long as the functions and effects achieved by the present specification are the same or similar.
S22: calculating the ratio of the single-layer unidirectional water injection quantity to the total water injection quantity to be used as a dimensionless water injection coefficient of the specified water injection well in the direction from a specified small layer to a specified oil well; the total water injection amount is the accumulated value of the water injection amount of each water injection well.
S24: calculating the ratio of the average horizon thickness of the specified small layer between the specified water injection well and the specified oil well to the total horizon thickness as a dimensionless effective thickness coefficient of the specified small layer between the specified water injection well and the specified oil well; the total layer thickness comprises an accumulated value of the layer average thickness corresponding to each designated small layer.
S26: and determining whether a seepage channel of the specified water injection well in the specified small layer to the direction of the specified oil well is an invalid circulation channel or not according to the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient.
The single-layer unidirectional water injection quantity Q can be calculatedi,j,kRelative total water injection
Figure BDA0003249189700000091
Ratio of
Figure BDA0003249189700000092
As a dimensionless water injection coefficient W of a specified water injection well i in the direction from a specified small layer k to a specified oil well ji,j,k. And calculating the average thickness H of the layer of the designated small layer k between the designated water injection well i and the designated oil well ji,j,kRelative total layer thickness
Figure BDA0003249189700000101
Ratio of
Figure BDA0003249189700000102
Dimensionless effective thickness coefficient Z between designated water injection well i and designated oil well j as designated zone ki,j,k. Wherein S represents the number of small layers obtained after the target oil reservoir stratum is divided, N represents the number of oil wells in the target oil reservoir, and M represents the number of water injection wells in the target oil reservoir. Thus, a data pair consisting of a dimensionless water injection coefficient and a dimensionless effective thickness coefficient can be obtained: (W)i,j,k，Zi,j,k)。
Then, the ratio of each dimensionless water injection coefficient to the corresponding dimensionless effective thickness coefficient can be calculated. And determining whether the seepage channel of the specified water injection well in the specified small layer towards the specified oil well is an invalid circulation channel according to the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient. If the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient is larger than the invalid circulation threshold, determining the seepage channel of the designated water injection well in the designated stratum towards the designated oil well as the invalid circulation channel. The invalid cycle threshold value can be determined according to the communication rule of the target oil reservoir. Preferably, in some embodiments, the ineffective circulation threshold value may be determined according to a data set composed of a ratio of a non-dimensional water injection coefficient to a non-dimensional effective thickness coefficient before the target reservoir is plugged and transformed, so that the ineffective circulation threshold value is more consistent with geological features of the target reservoir.
Fig. 3 shows a schematic diagram of injection and production well distribution, Pro1, Pro2, Pro3 in fig. 3 respectively show water injection wells numbered 1, 2, 3, and Inj1 shows an oil well numbered 1. The mode can be used for judging whether the seepage channel of any water injection well in any small layer to any oil well direction is an invalid circulation channel. Supposing that the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient is greater than the water injection well Pro3, the oil well Inj1 and the 2 nd small layer of the invalid circulation threshold, it can be shown that the water injection well Pro3 is most obvious to the oil well Inj1 in the preferential seepage passage of the corresponding small layer 2, the water absorption capacity is integrally larger, and the main position for forming the preferential seepage passage is the main position, namely the position of the small layer 2 between the water injection well Pro3 and the oil well Inj1 is more likely to be the invalid water circulation layer. Therefore, the method provided by the embodiment can more directly reflect the relative water absorption of the small layers among the injection wells and more directly reflect the relative difference and relative dominance of the small layers, so that the scheme provided by the embodiment can be used for more accurately and quantitatively determining the ineffective circulation channel, and a foundation is laid for better guiding the adjustment work of oilfield measures.
The above data pair (W) may also be usedi,j,k，Zi,j,k) According to the dimensionless water injection coefficient Wi,j,kIs arranged based on the order from small to large. Then, a data pair (X) can be constructed based on the sorted data pairsm,n,s,Ym,n,s)，Xm,n,sTo accumulate the effective thickness factor, Xm,n,sIs the cumulative water injection coefficient.
May be based on the ordered (W)i,j,k，Zi,j,k) And (4) data pairs, namely, renumbering the water injection well, the oil well and the small layer. For the purpose of distinguishing and representing, the dimensionless water injection coefficient and the dimensionless effective thickness coefficient after renumbering can be represented as WI,J,K、ZI,J,KI, J, K has a value range of [1, M]、[1,N]、[1,S]. The new numbers of the water injection well, the oil well and the small layer can be associated with the original numbers of the water injection well, the oil well and the small layer so as to facilitate the tracing.
In some embodiments, the value of I, J, K may be configured in a loop nesting manner. Specifically, I may be used as the outermost layer, J may be used as the second outermost layer, and K may be used as the innermost layer, where the corresponding cycle of each layer is: the value of the corresponding element of each layer is as follows (W)i,j,k，Zi,j,k) The sequence of the data pairs is increased from 1 to the maximum value of the corresponding element in turn, and each time the value is increasedThe stride is 1. Correspondingly, after the innermost layer completes one cycle each time, adding 1 to the secondary outer layer; when the secondary outer layer also completes one cycle, adding 1 to the outermost layer; and then, executing the innermost layer cycle, and adding 1 to the next outer layer after the innermost layer completes one cycle, and adding 1 to the outermost layer after the completion. And the rest is repeated until the outermost layer also completes the circulation.
For example, starting from I, J, K all taking the value of 1, the value of I, J is kept unchanged, and the value of K is increased; after K is increased to S, adding 1 to the value of J, keeping the value of I, J unchanged, and increasing the value of K; after K is increased to S, adding 1 to the value of J, keeping the value of I, J unchanged, and increasing the value of K; until J is increased to N and K is also increased to S, the value of I is added by 1; and repeating the steps until the value of I is increased to M, J is increased to N, and K is increased to S, and then the renumbering of the injection well, the oil well and the small layer in the target oil deposit is completed.
Then, a data pair (X) can be constructed in the following mannerm,n,s,Ym,n,s)：
Figure BDA0003249189700000111
Figure BDA0003249189700000112
Wherein, WI,J,K、ZI,J,KRespectively representing a dimensionless water injection coefficient and a dimensionless effective thickness coefficient, I, K, J representing numbers corresponding to a water injection well, an oil well and a small layer, and I, K, J values are configured in a circular nesting mode; n, M and S represent values of I, J, K in the sequential accumulation process, M is more than or equal to 1 and less than or equal to M, N is more than or equal to 1 and less than or equal to N, and S is more than or equal to 1 and less than or equal to S; wherein, the loop nesting mode comprises the following steps: taking I as an outermost layer, J as a secondary outer layer and K as an innermost layer, wherein the cycle corresponding to each layer is that the value of the corresponding element of each layer is sequentially increased from 1 to the maximum value of the corresponding element according to the sequence from small to large of the dimensionless water injection coefficient, and the step of each increase is 1.
Then, can be(Xm,n,s,Ym,n,s) And projecting the water injection quantity to a rectangular Cartesian coordinate system to obtain a Lorentz curve for evaluating the heterogeneity of the water injection quantity of the target oil reservoir. The water absorption coefficient of variation of the lorentz curve may be calculated as the water absorption coefficient of variation of the target reservoir.
In some embodiments of the present disclosure, the coefficient of variation of water absorption can be quantitatively calculated by:
Figure BDA0003249189700000121
wherein V is the water absorption variation coefficient of the target oil reservoir.
The development effect of the target reservoir can be evaluated according to the Lorentz curve and/or the water absorption coefficient of variation determined in the above embodiment. The Lorentz curve can be used for evaluating the water absorption nonuniformity among injection and production wells of each small layer, and the water absorption variation coefficient can be used for evaluating the water drive ineffective circulation degree of a target oil reservoir. The smaller the area enclosed by the Lorentz curve and the coordinate axis is, the smaller the water absorption coefficient of variation is, and the better the development effect is. If the development effect does not meet the actual production requirement, the target oil reservoir can be subjected to plugging or fracturing treatment in a targeted manner according to the determined ineffective circulation channel in the embodiment.
For example, a certain ineffective circulation channel can be blocked. And the reservoir transformation coefficient of each designated small layer in the target oil reservoir between the designated water injection well and the designated oil well can be updated based on the plugging degree, the seepage resistance coefficient is calculated by utilizing the updated reservoir transformation coefficient, and the water injection quantity of the single well is quantitatively determined and distributed to each layer of water injection quantity in the direction of each effective oil well communicated with the single well. And determining whether the seepage channel of the specified water injection well in the specified small layer towards the direction of the specified oil well is an invalid circulation channel or not based on the water injection quantity distribution determined quantitatively. And evaluating the development effect of the target oil reservoir after plugging treatment, and if the development effect does not meet the requirement, further plugging treatment is carried out based on the determined ineffective circulation channel until the development effect meets the requirement, so that the modification and adjustment measures of the target oil reservoir are more targeted, and the development effect of the target oil reservoir is improved. Meanwhile, a reservoir transformation coefficient is further determined by combining the transformation degree, and the communication rule of the transformed reservoir is indirectly represented by combining the reservoir transformation coefficient, so that the reservoir development effect and whether the seepage channel is an ineffective circulation channel or not are evaluated, the determination and evaluation of the reservoir ineffective circulation channel can be more in line with the reservoir communication rule after transformation, the knowledge of the reservoir water injection development effect is further improved, and the design of a reservoir water injection development scheme is better guided.
Based on the scheme provided by the above embodiment, the embodiment of the present specification further provides an example of identifying and evaluating the flooding ineffective cycle of the injection and production unit as follows.
(1) The experimental area of the example comprises 4 oil production wells, 2 water injection wells and 3 total layers of oil reservoirs.
(2) According to the hydroelectric similarity principle, the physical parameters of each small layer, the injection-production well distance and the transformation condition are utilized, the seepage resistance coefficient is introduced to evaluate the single-layer seepage resistance among wells in the injection-production units, the fluid viscosity, the oil-water two-phase seepage relation, the average effective thickness of the reservoir, the average permeability, the effective distance among wells and the reservoir transformation condition are fully considered, and the seepage resistance coefficient of each injection-production unit is calculated by utilizing the scheme provided by the embodiment.
(3) And (3) taking the injection-production pressure difference into consideration, introducing a water injection amount plane distribution coefficient and a water injection amount longitudinal distribution coefficient, and calculating to obtain the water injection amount of each layer of the single-well in the direction of each affected well communicated with the single-well by combining the physical property parameters of the injection-production well and the production dynamic data.
(4) Determining an invalid circulation horizon by using a ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient; and the heterogeneity of the water absorption of each small layer of the injection and production well can be evaluated by utilizing a Lorentz curve, the water drive invalid circulation degree can be evaluated by utilizing the variation coefficient of the water absorption, and the quantitative research on the development effect of heterogeneous sandstone reservoirs with different characteristics can be realized.
Supposing that the thickness of each layer in the experimental area is 5m, the permeability from the first layer to the third layer is 50mD, 100mD and 150mD in sequence, the reservoir stratum of the oil reservoir is not reformed in the primary development, and the coefficient of the reformation of the reservoir stratum isMi,j,kThe value is 1, the viscosity of crude oil is 15 mPas, the viscosity of water is 1 mPas, and the water injection well 1 co-injects 1000m of water3And the water injection well 2 injects 2000m of water3The pressure measurements at each of the production and injection wells (i.e., injection wells) are shown in fig. 4. Pressure data for each of the production and injection wells prior to the adjustment measure is shown in FIG. 4. Then, the adjusted lorentz curve can be calculated by the method provided in the above embodiment, as shown in fig. 5. The coefficient of variation of water absorption before adjustment was 0.297469. The ratio of the dimensionless water injection coefficient and the dimensionless effective thickness coefficient of each layer of the injection-production well before adjustment is shown in table 1.
TABLE 1
Figure BDA0003249189700000131
Figure BDA0003249189700000141
According to the ratio of the dimensionless water injection coefficient and the dimensionless effective thickness coefficient of each layer of the injection and production well before adjustment, the critical value is determined to be 2, if the ratio is higher than the critical value, the seepage channel of the injection and production well at the position is regarded as a water drive ineffective circulation channel, according to the initial development effect of the oil reservoir, the water drive ineffective circulation channel marked in the table 1 can be blocked in a targeted manner, the transformation coefficient of the reservoir after blocking is shown in the table 2, and the pressure detection data of each layer of the production well and the water injection well after adjustment is shown in the figure 6.
TABLE 2
i,j, k 1,1,1 1,2,1 1,3,1 1,4,1 1,1,2 1,2,2 1,3,2 1,4,2
Coefficient of reservoir reformation (M) 1 1 1 1 1 1 1 1
i,j, k 1,1,3 1,2,3 1,3,3 1,4,3 2,1,1 2,2,1 2,3,1 2,4,1
Coefficient of reservoir reformation (M) 0.8 1 0.8 1 1 1 1 1
i,j, k 2,1,2 2,2,2 2,3,2 2,4,2 2,1,3 2,2,3 2,3,3 2,4,3
Coefficient of reservoir reformation (M) 1 1 1 1 0.7 0.6 0.7 1
The adjusted lorentz curve calculated by the method provided in the above embodiment is shown in fig. 7. The coefficient of variation of the water absorption capacity after adjustment is 0.253069. The ratio of the dimensionless water injection coefficient and the dimensionless effective thickness coefficient of each layer of the injection and production well after adjustment is shown in table 3.
TABLE 3
Figure BDA0003249189700000142
Figure BDA0003249189700000151
The non-dimensional water injection coefficient and the non-dimensional effective thickness coefficient ratio of each layer of the injection well before adjustment can be obtained according to the table 1, and the water injection well 1-the production well 1-the third layer, the water injection well 1-the production well 3-the third layer, the water injection well 2-the production well 2-the third layer and the water injection well 2-the production well 3-the third layer are water drive ineffective circulation channels before adjustment. Comparing fig. 5 and fig. 7, it can be seen that the enclosed area of the lorentz curve and the coordinate axis after plugging is smaller. The coefficient of variation of the water absorption before and after the scheme is adjusted can be obtained, and the coefficient of variation of the water absorption after the scheme is adjusted is smaller. As can be seen from the ratio of the dimensionless water injection coefficient and the dimensionless effective thickness coefficient of each layer of the injection and production well after adjustment in the table 3, the number of the reservoir water drive ineffective circulation channels after plugging is reduced, and the adjusted water injection well 2-the production well 3-the third layer and the water injection well 2-the production well 4-the third layer are water drive ineffective circulation channels which are key plugging parts in the next step.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. For details, reference may be made to the description of the related embodiments of the related processing, and details are not repeated herein.
The foregoing description has been directed to specific embodiments of this disclosure. Other embodiments are within the scope of the following claims. In some cases, the actions or steps recited in the claims may be performed in a different order than in the embodiments and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing may also be possible or may be advantageous.
Based on the method for identifying the water drive ineffective circulation based on the lorentz curve, one or more embodiments of the specification further provide a device for identifying the water drive ineffective circulation based on the lorentz curve. The apparatus may include systems, software (applications), modules, components, servers, etc. that utilize the methods described in the embodiments of the present specification in conjunction with hardware implementations as necessary. Based on the same innovative conception, embodiments of the present specification provide an apparatus as described in the following embodiments. Since the implementation scheme of the apparatus for solving the problem is similar to that of the method, the specific implementation of the apparatus in the embodiment of the present specification may refer to the implementation of the foregoing method, and repeated details are not repeated. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated. Specifically, fig. 8 is a schematic block diagram of an embodiment of a water drive invalid circulation identification apparatus based on a lorentz curve provided in the specification, and as shown in fig. 8, the apparatus may include:
the water injection rate determining module 80 is used for calculating single-layer unidirectional water injection rate distributed to the direction of the specified oil well by the specified water injection well in the specified small layer according to the single-well water injection rate, the plane distribution coefficient and the longitudinal distribution coefficient of the specified water injection well in the target oil reservoir; wherein the plane distribution coefficient is determined according to the seepage resistance coefficient of the specified small layer between the specified water injection well and the specified oil well and the pressure of the specified water injection well in the specified small layer; the longitudinal distribution coefficient is determined at least according to the seepage resistance coefficient, the average permeability and the reservoir transformation coefficient of the specified small layer between the specified water injection well and the specified oil well;
a water injection coefficient calculation module 82, configured to calculate a ratio of the single-layer unidirectional water injection amount to the total water injection amount, where the ratio is used as a dimensionless water injection coefficient of the specified water injection well in a direction from a specified small layer to a specified oil well; the total water injection amount is the accumulated value of the water injection amount of each water injection well;
a thickness coefficient calculation module 84, configured to calculate a ratio of the average horizon thickness of the specified small layer between the specified water injection well and the specified oil well to the total horizon thickness as a dimensionless effective thickness coefficient of the specified small layer between the specified water injection well and the specified oil well; the total layer thickness comprises an accumulated value of the layer average thickness corresponding to each appointed small layer;
and the invalid circulation channel identification module 86 is used for identifying whether the seepage channel of the specified water injection well in the direction of the specified oil well in the specified small layer is an invalid circulation channel according to the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient.
In other embodiments of the present description, the apparatus further includes a vertical distribution coefficient calculation module, where the vertical distribution coefficient calculation module is configured to calculate the vertical distribution coefficient by:
Figure BDA0003249189700000171
wherein, Ci,kRepresents the longitudinal distribution coefficient, H, of a given water injection well i vertically to a given small layer ki,j,kRepresenting the average thickness of the horizon of the specified small layer K between the specified water injection well i and the specified oil well j, Ki,j,kRepresents the average permeability, M, of the given fraction k between the given water injection well i and the given oil well ji,j,kRepresenting the reservoir transformation coefficient, R, of the specified small layer k between the specified water injection well i and the specified oil well ji,j,kRepresenting the seepage resistance coefficient of the specified small layer k between the specified water injection well i and the specified oil well j;
Figure BDA0003249189700000172
to specify the number of wells in the small layer k, S represents the number of small layers obtained after the target reservoir formation is partitioned.
In other embodiments of the present description, the apparatus further comprises a seepage resistance coefficient calculation module, and the seepage resistance coefficient calculation module is configured to calculate the seepage resistance coefficient by:
Figure BDA0003249189700000173
Figure BDA0003249189700000174
wherein λ isi,j,k、Li,j,kRespectively the apparent viscosity and the effective well spacing between a designated water injection well i and a designated oil well j of a designated small layer Kro、Krw、μo、μwThe relative permeability of the oil phase, the relative permeability of the water phase, the viscosity of the crude oil and the viscosity of water are respectively.
It should be noted that the above-described apparatus may also include other embodiments according to the description of the method embodiment. The specific implementation manner may refer to the description of the related method embodiment, and is not described in detail herein.
Based on the method for identifying the water drive invalid cycle based on the lorentz curve, the specification further provides a computer-readable storage medium, on which computer instructions are stored, and the computer instructions, when executed, implement the steps of the method according to any one or more of the embodiments. The storage medium may include a physical device for storing information, and typically, the information is digitized and then stored using an electrical, magnetic, or optical media. The storage medium may include: devices that store information using electrical energy, such as various types of memory, e.g., RAM, ROM, etc.; devices that store information using magnetic energy, such as hard disks, floppy disks, tapes, core memories, bubble memories, and usb disks; devices that store information optically, such as CDs or DVDs. Of course, there are other ways of storing media that can be read, such as quantum memory, graphene memory, and so forth.
It should be noted that the embodiments of the present disclosure are not limited to the cases where the data model/template is necessarily compliant with the standard data model/template or the description of the embodiments of the present disclosure. Certain industry standards, or implementations modified slightly from those described using custom modes or examples, may also achieve the same, equivalent, or similar, or other, contemplated implementations of the above-described examples. The embodiments using these modified or transformed data acquisition, storage, judgment, processing, etc. may still fall within the scope of the alternative embodiments of the present description.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment. In the description of the specification, reference to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the specification. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.
The above description is only an example of the present specification, and is not intended to limit the present specification. Various modifications and alterations to this description will become apparent to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present specification should be included in the scope of the claims of the present specification.

Claims (10)

1. A water drive invalid cycle identification method based on a Lorentz curve is characterized by comprising the following steps:
calculating single-layer unidirectional water injection quantity distributed to the direction of the specified oil well in the specified small layer by the specified water injection well according to the single-well water injection quantity, the plane distribution coefficient and the longitudinal distribution coefficient of the specified water injection well in the target oil reservoir; wherein the plane distribution coefficient is determined according to the seepage resistance coefficient of the specified small layer between the specified water injection well and the specified oil well and the pressure of the specified water injection well in the specified small layer; the longitudinal distribution coefficient is determined at least according to the seepage resistance coefficient, the average permeability and the reservoir transformation coefficient of the specified small layer between the specified water injection well and the specified oil well;
calculating the ratio of the single-layer unidirectional water injection quantity to the total water injection quantity to be used as a dimensionless water injection coefficient of the specified water injection well in the direction from a specified small layer to a specified oil well; the total water injection amount is the accumulated value of the water injection amount of each water injection well;
calculating the ratio of the average horizon thickness of the specified small layer between the specified water injection well and the specified oil well to the total horizon thickness as a dimensionless effective thickness coefficient of the specified small layer between the specified water injection well and the specified oil well; the total layer thickness comprises an accumulated value of the layer average thickness corresponding to each appointed small layer;
and determining whether a seepage channel of the specified water injection well in the specified small layer to the direction of the specified oil well is an invalid circulation channel or not according to the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient.
2. The method of claim 1, wherein the vertical distribution coefficient is determined in a coefficient manner:
Figure FDA0003249189690000011
wherein, Ci,kRepresents the longitudinal distribution coefficient, H, of a given water injection well i to a given sub-zone ki,j,kRepresenting the average thickness of the horizon of the specified small layer K between the specified water injection well i and the specified oil well j, Ki,j,kRepresents the average permeability, M, of the given fraction k between the given water injection well i and the given oil well ji,j,kRepresenting the reservoir transformation coefficient, R, of the specified small layer k between the specified water injection well i and the specified oil well ji,j,kRepresenting the seepage resistance coefficient of the specified small layer k between the specified water injection well i and the specified oil well j;
Figure FDA0003249189690000012
to specify the number of wells in a small zone k, S represents the pair targetThe number of small layers obtained after the reservoir strata are divided.
3. The method of claim 2, wherein the seepage resistance coefficient of the designated water injection well i in the direction of the designated well j at the designated interval k is determined using the following calculation model:
Figure FDA0003249189690000021
Figure FDA0003249189690000022
wherein λ isi,j,k、Li,j,kRespectively the apparent viscosity and the effective well spacing between a designated water injection well i and a designated oil well j of a designated small layer Kro、Krw、μo、μwThe relative permeability of the oil phase, the relative permeability of the water phase, the viscosity of the crude oil and the viscosity of water are respectively.
4. The method according to any one of claims 1-3, further comprising:
determining the water absorption capacity variation coefficient of the target oil reservoir by using the following calculation model:
Figure FDA0003249189690000023
wherein V is the coefficient of variation of water absorption, WI,J,K、ZI,J,KRespectively representing a dimensionless water injection coefficient and a dimensionless effective thickness coefficient, I, K, J representing numbers corresponding to a water injection well, an oil well and a small layer, and I, K, J values are configured in a circular nesting mode; s represents the number of small layers obtained after the target oil reservoir stratum is divided, N represents the number of oil wells in the target oil reservoir, and M represents the number of water injection wells in the target oil reservoir; n, M and S represent values of I, J, K in the sequential accumulation process, M is more than or equal to 1 and less than or equal to M, N is more than or equal to 1 and less than or equal to N, and S is more than or equal to 1 and less than or equal to S;wherein, the loop nesting mode comprises the following steps: taking I as an outermost layer, J as a secondary outer layer and K as an innermost layer, wherein the cycle corresponding to each layer is that the value of the corresponding element of each layer is sequentially increased from 1 to the maximum value of the corresponding element according to the sequence of the dimensionless water injection coefficient from small to large, and the step length of each increase is 1;
and analyzing the development effect of the target oil reservoir by using the water absorption variation coefficient.
5. The method of claim 1, wherein in the case that the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient is greater than an ineffective circulation threshold, determining a seepage channel of the specified water injection well in a specified direction to a specified well in a specified stratum as an ineffective circulation channel; and determining the invalid circulation threshold according to a data set consisting of the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient before the target oil reservoir is plugged and reformed.
6. The method of claim 5, further comprising:
plugging or fracturing the target oil reservoir according to the ineffective circulation channel;
and updating the reservoir transformation coefficients of the specified small layers of the target oil reservoir between the specified water injection well and the specified oil well after plugging or fracturing treatment so as to identify the invalid circulation channel of the target oil reservoir after the current plugging by using the updated reservoir transformation coefficients.
7. A water drive invalid circulation recognition device based on a Lorentz curve is characterized by comprising the following components:
the water injection rate determining module is used for calculating single-layer unidirectional water injection rate distributed to the direction of the specified oil well by the specified water injection well in the specified small layer according to the single-well water injection rate, the plane distribution coefficient and the longitudinal distribution coefficient of the specified water injection well in the target oil reservoir; wherein the plane distribution coefficient is determined according to the seepage resistance coefficient of the specified small layer between the specified water injection well and the specified oil well and the pressure of the specified water injection well in the specified small layer; the longitudinal distribution coefficient is determined at least according to the seepage resistance coefficient, the average permeability and the reservoir transformation coefficient of the specified small layer between the specified water injection well and the specified oil well;
the water injection coefficient calculation module is used for calculating the ratio of the single-layer unidirectional water injection quantity to the total water injection quantity to be used as a dimensionless water injection coefficient of the specified water injection well in the direction of the specified oil well in the specified small layer; the total water injection amount is the accumulated value of the water injection amount of each water injection well;
the thickness coefficient calculation module is used for calculating the ratio of the average thickness of the layer position of the specified small layer between the specified water injection well and the specified oil well to the total layer position thickness as a dimensionless effective thickness coefficient of the specified small layer between the specified water injection well and the specified oil well; the total layer thickness comprises an accumulated value of the layer average thickness corresponding to each appointed small layer;
and the invalid circulation channel identification module is used for identifying whether the seepage channel of the specified water injection well in the specified small layer to the specified oil well direction is an invalid circulation channel according to the ratio of the dimensionless water injection coefficient to the dimensionless effective thickness coefficient.
8. The apparatus of claim 7, further comprising a vertical distribution coefficient calculation module configured to calculate a vertical distribution coefficient by:
Figure FDA0003249189690000031
wherein, Ci,kRepresents the longitudinal distribution coefficient, H, of a given water injection well i vertically to a given small layer ki,j,kRepresenting the average thickness of the horizon of the specified small layer K between the specified water injection well i and the specified oil well j, Ki,j,kRepresents the average permeability, M, of the given fraction k between the given water injection well i and the given oil well ji,j,kRepresenting the reservoir transformation coefficient, R, of the specified small layer k between the specified water injection well i and the specified oil well ji,j,kRepresenting the seepage resistance coefficient of the specified small layer k between the specified water injection well i and the specified oil well j;
Figure FDA0003249189690000032
to specify the number of wells in the small layer k, S represents the number of small layers obtained after the target reservoir formation is partitioned.
9. The apparatus of claim 8, further comprising a seepage resistance coefficient calculation module for calculating a seepage resistance coefficient by:
Figure FDA0003249189690000041
Figure FDA0003249189690000042
wherein λ isi,j,k、Li,j,kRespectively the apparent viscosity and the effective well spacing between a designated water injection well i and a designated oil well j of a designated small layer Kro、Krw、μo、μwThe relative permeability of the oil phase, the relative permeability of the water phase, the viscosity of the crude oil and the viscosity of water are respectively.
10. A computer-readable storage medium having computer instructions stored thereon which, when executed, implement steps comprising the method of any of claims 1-6.
CN202111040985.7A 2021-09-07 2021-09-07 Water drive invalid cycle identification method, device and medium based on Lorentz curve Active CN113738319B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111040985.7A CN113738319B (en) 2021-09-07 2021-09-07 Water drive invalid cycle identification method, device and medium based on Lorentz curve

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111040985.7A CN113738319B (en) 2021-09-07 2021-09-07 Water drive invalid cycle identification method, device and medium based on Lorentz curve

Publications (2)

Publication Number Publication Date
CN113738319A true CN113738319A (en) 2021-12-03
CN113738319B CN113738319B (en) 2022-05-03

Family

ID=78736285

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111040985.7A Active CN113738319B (en) 2021-09-07 2021-09-07 Water drive invalid cycle identification method, device and medium based on Lorentz curve

Country Status (1)

Country Link
CN (1) CN113738319B (en)

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1916358A (en) * 2006-09-06 2007-02-21 大庆油田有限责任公司 Method for coordinating relation for collecting or dissipating as well as injecting or gathering ore mass in narrow and small riverway
CN101775971A (en) * 2010-02-03 2010-07-14 中国石油天然气股份有限公司 Oil-field largest swept volume chemical-flooding oil production method
CN102041995A (en) * 2010-12-02 2011-05-04 中国海洋石油总公司 System for monitoring complicated oil deposit flooding conditions
US20120191432A1 (en) * 2011-01-26 2012-07-26 Schlumberger Technology Corporation Visualizing fluid flow in subsurface reservoirs
WO2015013697A1 (en) * 2013-07-26 2015-01-29 Schlumberger Canada Limited Well treatment
CN104879104A (en) * 2015-05-27 2015-09-02 中国石油天然气股份有限公司 Oil reservoir water injection method
US20160061020A1 (en) * 2014-08-22 2016-03-03 Chevron U.S.A. Inc. Flooding analysis tool and method thereof
CN105626010A (en) * 2016-03-16 2016-06-01 中国石油大学(华东) Method for reasonably dividing water injection layer sections in segmented water injection well
CN105696986A (en) * 2014-12-09 2016-06-22 中国海洋石油总公司 Novel combination flooding oil flooding experiment/test simulating method
US20170114619A1 (en) * 2015-10-26 2017-04-27 International Business Machines Corporation Reduced space clustering representatives and its application to long term prediction
CN110633529A (en) * 2019-09-18 2019-12-31 中国石油大学(北京) Data processing method, device and system for determining oil reservoir water drive development effect
CN110765624A (en) * 2019-10-29 2020-02-07 中国石油化工股份有限公司 Reasonable layering method for water injection oil reservoir
CN111520136A (en) * 2020-06-29 2020-08-11 东北石油大学 Method for calculating pressure behind blanking plug nozzle by considering water injection starting pressure gradient
US20210131231A1 (en) * 2019-10-31 2021-05-06 Saudi Arabian Oil Company Dynamic Calibration of Reservoir Simulation Models using Flux Conditioning
CN112901145A (en) * 2021-03-19 2021-06-04 大庆油田有限责任公司 Volume energy method for analyzing injection-production relation between oil-water wells

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1916358A (en) * 2006-09-06 2007-02-21 大庆油田有限责任公司 Method for coordinating relation for collecting or dissipating as well as injecting or gathering ore mass in narrow and small riverway
CN101775971A (en) * 2010-02-03 2010-07-14 中国石油天然气股份有限公司 Oil-field largest swept volume chemical-flooding oil production method
CN102041995A (en) * 2010-12-02 2011-05-04 中国海洋石油总公司 System for monitoring complicated oil deposit flooding conditions
US20120191432A1 (en) * 2011-01-26 2012-07-26 Schlumberger Technology Corporation Visualizing fluid flow in subsurface reservoirs
WO2015013697A1 (en) * 2013-07-26 2015-01-29 Schlumberger Canada Limited Well treatment
US20160061020A1 (en) * 2014-08-22 2016-03-03 Chevron U.S.A. Inc. Flooding analysis tool and method thereof
CN105696986A (en) * 2014-12-09 2016-06-22 中国海洋石油总公司 Novel combination flooding oil flooding experiment/test simulating method
CN104879104A (en) * 2015-05-27 2015-09-02 中国石油天然气股份有限公司 Oil reservoir water injection method
US20170114619A1 (en) * 2015-10-26 2017-04-27 International Business Machines Corporation Reduced space clustering representatives and its application to long term prediction
CN105626010A (en) * 2016-03-16 2016-06-01 中国石油大学(华东) Method for reasonably dividing water injection layer sections in segmented water injection well
CN110633529A (en) * 2019-09-18 2019-12-31 中国石油大学(北京) Data processing method, device and system for determining oil reservoir water drive development effect
CN110765624A (en) * 2019-10-29 2020-02-07 中国石油化工股份有限公司 Reasonable layering method for water injection oil reservoir
US20210131231A1 (en) * 2019-10-31 2021-05-06 Saudi Arabian Oil Company Dynamic Calibration of Reservoir Simulation Models using Flux Conditioning
CN111520136A (en) * 2020-06-29 2020-08-11 东北石油大学 Method for calculating pressure behind blanking plug nozzle by considering water injection starting pressure gradient
CN112901145A (en) * 2021-03-19 2021-06-04 大庆油田有限责任公司 Volume energy method for analyzing injection-production relation between oil-water wells

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
PENG, FJ: "Molecular dynamics simulation to estimate minimum miscibility pressure for oil with pure and impure CO2", 《JOURNAL OF PHYSICS COMMUNICATIONS》 *
ZHANG, LF等: "Research on sensitivity damage of naturally fractured carbonate reservoirs in Ordos Basin", 《ARABIAN JOURNAL OF GEOSCIENCES》 *
卢川等: "非均质底水油藏水平井合理生产压差实验研究", 《特种油气藏》 *
钟睿鸿: "注水无效循环识别及优势渗流通道定量表征方法研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅰ辑》 *

Also Published As

Publication number Publication date
CN113738319B (en) 2022-05-03

Similar Documents

Publication Publication Date Title
Friedmann et al. Assessing uncertainty in channelized reservoirs using experimental designs
CN105488583B (en) The method and device of the fine and close oil region recoverable reserves to be evaluated of prediction
Sun et al. Data-space approaches for uncertainty quantification of CO2 plume location in geological carbon storage
Tian et al. Inferring interwell connectivity using production data
CN107330578B (en) Sand body connectivity evaluation method and device
CN110593863B (en) Identification method and identification system for water consumption zone of high water-cut oil reservoir
Moharram et al. Optimal groundwater management using genetic algorithm in El-Farafra oasis, western desert, Egypt
CN106503295B (en) A kind of method and device for explaining oil field Water Flooding Layer using state-space model
CN108019204B (en) Neutron-gamma gas layer identification method based on physical property and formation water mineralization correction
CN105404726A (en) Method and device for reverting inter-well connectivity through Gaussian distribution-based capacitance model
CN112282744B (en) Unconventional oil and gas reservoir well pattern deployment optimization method and device
CN110633529B (en) Data processing method, device and system for determining oil reservoir water drive development effect
Durlofsky An approximate model for well productivity in heterogeneous porous media
CN102704924B (en) The defining method of effective dried layer and device
CN113738319B (en) Water drive invalid cycle identification method, device and medium based on Lorentz curve
Eriavbe et al. Machine learning application to permeability prediction using log & core measurements: A realistic workflow application for reservoir characterization
CN108952676A (en) A kind of shale gas reservoir heterogeneity evaluation method and its device
CN103201649A (en) Method and system for pulse neutron capture sigma inversion
Voleti et al. Normalized depths as key input and detailed QC steps for improved permeability predictions using existing machine learning techniques
CN114876454A (en) Method and system for identifying complex lithologic water flooded layer of highly deviated well
Fournier et al. Integrated rock-typing with capillary pressure curve clustering
Feizabadi et al. An integrated approach to building history-matched geomodels to understand complex long lake oil sands reservoirs, part 2: Simulation
Baig et al. Productivity Increase Estimation for Multi Stage Fracturing in Horizontal Wells for Tight Oil Reservoirs
CN113047823A (en) Method and device for determining horizontal well step strategy
Litvak et al. Uncertainty Estimation in Production Predictions Constrained by Production History and Time-Lapse Seismic in a GOM Oil Field

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