CN109630078A - A kind of layer is interior to generate CO2Transfer drive parameters design method - Google Patents

A kind of layer is interior to generate CO2Transfer drive parameters design method Download PDF

Info

Publication number
CN109630078A
CN109630078A CN201811316949.7A CN201811316949A CN109630078A CN 109630078 A CN109630078 A CN 109630078A CN 201811316949 A CN201811316949 A CN 201811316949A CN 109630078 A CN109630078 A CN 109630078A
Authority
CN
China
Prior art keywords
layer
transfer drive
release agent
inflating medium
bubble release
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
CN201811316949.7A
Other languages
Chinese (zh)
Other versions
CN109630078B (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 Oilfield Services Ltd
China National Offshore Oil Corp CNOOC
Original Assignee
China Oilfield Services Ltd
China National Offshore Oil Corp CNOOC
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 Oilfield Services Ltd, China National Offshore Oil Corp CNOOC filed Critical China Oilfield Services Ltd
Priority to CN201811316949.7A priority Critical patent/CN109630078B/en
Publication of CN109630078A publication Critical patent/CN109630078A/en
Application granted granted Critical
Publication of CN109630078B publication Critical patent/CN109630078B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16Enhanced recovery methods for obtaining hydrocarbons
    • E21B43/164Injecting CO2 or carbonated water
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/12Computing arrangements based on biological models using genetic models
    • G06N3/126Evolutionary algorithms, e.g. genetic algorithms or genetic programming

Abstract

This application discloses generate CO in a kind of layer2Transfer drive parameters design method, which comprises practical oil reservoir is reduced to the multiple zone geological model comprising transfer drive layer reconciliation blocking layer;Inflating medium and bubble release agent dosage are subjected to dynamic splitting, by the difference degree of each substratum injectivity index after analysis transfer drive, inflating medium injection rate, bubble release agent injection rate, inflating medium injection length, bubble release agent injection length and the injection slug number computation model when the injectivity index coefficient of variation minimizes are established;Establish inflating medium, bubble release agent injection pressure computation model;Model above is solved by genetic algorithm.CO is generated in the layer of the application2Transfer drive technique joins design method: profile control agent (inflating medium and bubble release agent) injection rate may be implemented and injection slug number calculates;The real-time prediction that work progress injection pressure may be implemented calculates;Measure effect characterization is intuitive.

Description

A kind of layer is interior to generate CO2Transfer drive parameters design method
Technical field
This application involves oil fields, are espespecially related to generating CO in one of reservoir waterflooding layer2Transfer drive technique ginseng Number design method.
Background technique
CO is generated in layer2Transfer drive technological parameter mainly includes inflating medium and bubble release agent dosage, slug parameter designing, construction note Enter pressure and injection discharge capacity.Scientific and reasonable transfer drive process parameters design can be polluted to avoid less permeable layer by blocking agent, most The improvement intake profile of limits expands swept volume, to improve the validity of transfer drive measure.The dosage of profile control agent is set at present Meter is mainly include the following types: first is that pressure drop indices P I (pressure index) and reservoir engineering RE (reservoir Engineering) the profile control agent dosage of decision-making technic determines method, needs the injection test at scene to test and determines use level coefficient, cost It is high;Second is that method for numerical simulation, needs the processes such as history matching, multistage plug combined prediction, takes a long time;Third is that energy is injected in front and back Power ratio method needs to be determined in advance the ratio between the water absorbing capacity of transfer drive interval before and after the processing, and flexibility is poor.And profile control agent injection pressure The design method of power mainly selective injection pressure method, discharge capacity setting method and method for numerical simulation etc..Selective implantation pressure Force method and discharge capacity setting method need ASP flooding test to determine injection pressure gradient etc., at high cost;Method for numerical simulation need Geologic modeling, History matching etc., research cycle is long.
Summary of the invention
In order to solve the above-mentioned technical problem, this application provides generate CO in a kind of layer2Transfer drive parameters design method, The method includes establishing to generate CO in layer2Transfer drive technique Chinese medicine injection rate and injection slug design a model;And it is adopted from note and is The pressure balance of system is set out, and in conjunction with calculated real time discharge volume, establishes the real-time computation model of work progress injection pressure.To The design for realizing parameters of construction technology reaches the optimal of dosing and profile modification.
In order to reach the application purpose, this application provides generate CO in a kind of layer2Transfer drive parameters design method, packet It includes:
(1) practical oil reservoir is reduced to the multiple zone geological model comprising transfer drive layer reconciliation blocking layer;
(2) inflating medium and bubble release agent dosage are subjected to dynamic splitting, pass through the difference of each substratum injectivity index after analysis transfer drive Off course degree, establish when the injectivity index coefficient of variation minimize when inflating medium injection rate, bubble release agent injection rate, inflating medium injection when Between, bubble release agent injection length and injection slug number computation model;
(3) inflating medium, bubble release agent injection pressure computation model are established;
(4) model above is solved by genetic algorithm.
In this application, reaction in stratum is injected by water injection well and generates CO2, refer in reservoir formation in the middle layer this paper. Due to the generally existing difference of reservoir injectivity index, inflating medium and bubble release agent can be selectively penetrated by stratum.Existed according to transfer drive medicament Thermally coupled flow in stratum, binder course is interior to generate CO2Transfer drive selects layer as a result, practical oil reservoir can be reduced to include " transfer drive layer " The multiple zone ideal geological model of " solution blocking layer ", is defined as transfer drive layer to the substratum for needing obviously to block in stratum, and needs Bubble release agent acidification improves formation physical property, while CO2Foam blocking acts on unconspicuous substratum and is defined as solution blocking layer.
After transfer drive measure, since profile control agent acts on the natural selection on stratum, the injectivity index of each transfer drive layer has difference The degree of irregularity of the reduction of degree, high less permeable layer relative injectivity also constantly reduces, and intake profile is effectively improved.With The continuous injection of profile control agent, the injectivity index difference of high less permeable layer can be smaller and smaller.When profile control agent total injection reaches one When determining degree, the injectivity index difference of high less permeable layer can be reduced to minimum value, continue to inject profile control agent, the difference of injectivity index It will increase instead.Therefore there are an optimal profile control agent dosage, the heterogeneity of intake profile can be improved to greatest extent.
Assuming that the well intake pressure of transfer drive well, injection displacement control are reasonable, all substratums are started, and angry The performance parameter (including concentration, viscosity, density etc.) of agent and bubble release agent is certain.Water flooding is averaged matter uniform thickness, ignores liquid and rock The elastic reaction of layer.Target first is that i represents transfer drive layer under agreement digital code, and j represents solution blocking layer.Arrange under digital code Target first is 1 to represent inflating medium, and 2 represent bubble release agent.
Single note inflating medium process
Transfer drive technique can carry out first it is single infuse inflating medium process, under a certain instantaneous time step, inflating medium is split point always Amount is equal to the inflating medium discharge capacity v at the moment1(z), the bubble release agent discharge capacity v for splitting point total amount and being equal to the moment of bubble release agent2(z)。
The seepage flow of inflating medium in the earth formation meets linear relationship for seepage.Therefore raw on unit time step on longitudinal profile The total injection rate of gas agent can use Darcy's law and be split by the injectivity index distribution before transfer drive measure and assigned to each substratum.
In formula: JiFor the injectivity index of transfer drive layer i before transfer drive measure, m3/(h·MPa);JjTo solve blocking layer j before transfer drive measure Injectivity index, m3/(h·MPa);v1It (z) is the discharge capacity of z-th of slug inflating medium, m3/h;n1For the small number of plies of transfer drive layer;n2 For the small number of plies for solving blocking layer;Q1iFor the inflating medium injection rate of transfer drive layer i, m3;kiFor the permeability of transfer drive layer i, mD;hiFor transfer drive The effective thickness of layer i, m;kjFor the permeability for solving blocking layer i, mD;hjFor the effective thickness for solving blocking layer i, m;μwFor the viscosity of water, mPa·s;rwFor water injection well wellbore radius, m;reFor drainage radius, m;A is unit conversion coefficient, a=0.0036.
Then, the inflating medium blocking radius under the time step on each substratum can be calculated according to volumetric method:
In formula, rfiFor the inflating medium blocking radius of transfer drive layer i, m;φiFor the porosity of transfer drive layer i.
For transfer drive layer, since inflating medium is the solution of certain carbonate, a large amount of aggregation can also generate weak closure Effect, therefore in transfer drive layer, caused permeability decrease be can not ignore when singly infusing inflating medium.The plugging action of inflating medium is being adjusted It drives and forms hypotonic region around layer pit shaft, so as to cause the change of stratum filtration field.According to the variation characteristic of stratum filtrational resistance, Flow region of the formation fluid in transfer drive layer can be divided into flowing within the scope of from pit shaft to blocking radius and from blocking half Diameter to supply marginal annular region in flowing.It is possible thereby to calculate the injectivity index J ' of each transfer drive layer after injection inflating mediumi:
In formula: J 'iFor the injectivity index of transfer drive layer after injection inflating medium, m3/(h·MPa);RRF1For inflating medium remnants resistance Force coefficient.
Due to ignoring the effect of solution blocking layer inflating medium injection, solution blocking layer injectivity index before and after inflating medium injection is kept It is constant:
In formula: kjFor the permeability for solving blocking layer, mD;hjFor the effective thickness for solving blocking layer, m;Jj、JjRespectively solution blocking layer is arranged Apply the injectivity index of front and back, m3/(h·MPa)。
Infuse bubble release agent process
After completing inflating medium injection process, outgassing agent solution can be injected into stratum.Under first time step, bubble release agent Split point according to each layer injectivity index after inflating medium has been infused in splitting point for each layer.
In formula: J 'iFor the injectivity index of transfer drive layer i after transfer drive measure, m3/(h·MPa);J′jFor de-plugging after transfer drive measure The injectivity index of layer j, m3/(h·MPa);v2It (z) is the discharge capacity of z-th of slug bubble release agent, m3/h;Q2iFor the outgassing of transfer drive layer i Agent injection rate, m3Q2J is the bubble release agent injection rate for solving blocking layer j, m3
Inflating medium is reacted with bubble release agent, generates CO2Foam system, and for solving blocking layer, since remaining oil saturation is higher, Meet oily defoaming effect, it is possible to calculate according to volumetric method and conciliate by the foam blocking radius under a time step on transfer drive layer The plug removal radius of blocking layer:
In formula: rfiFor the inflating medium blocking radius of transfer drive layer i, m;rfiFor the foam blocking radius of transfer drive layer i, m;rsjFor Solve the plug removal radius of blocking layer j, m.
Each layer injectivity index after calculating first time step:
In formula: θ is bubble release agent acidizing degree;RRF is foam residual resistance factor.
Then future time walks, and bubble release agent carries out dosage and splits timesharing, refers to according to each layer water suction after the variation of a upper time step Number split point.
After sometime walking, the bubble release agent dosage in transfer drive layer is greater than total dosage of inflating medium, i.e.,
The injectivity index calculation formula in transfer drive layer changes at this time, and the acid of bubble release agent also occurs in transfer drive layer at this time Change effect, transfer drive layer foam blocking radius and bubble release agent acidizing radius can be calculated according to volumetric method:
The injectivity index calculation formula of transfer drive layer reconciliation blocking layer are as follows:
After recycling by slug, the injectivity index final using each substratum is as evaluation index.Each substratum injectivity index is got over Close, intake profile is more uniform in the injecting process.Therefore in order to reflect the difference degree of the injectivity index of each permeable formation after transfer drive, Using the coefficient of variation as evaluation index.In statistics, the coefficient of variation is also known as standard rate, is defined as one group of observation data Standard deviation and mathematic expectaion ratio, that reflects observation data between discrete and difference degree.Therefore suction can be defined The aqua index coefficient of variation is as follows:
Wherein
In formula:For the mean water absorption index on stratum after transfer drive, m3/(h·MPa);VJFor the injectivity index coefficient of variation.
The injectivity index coefficient of variation is bigger after transfer drive, illustrates that the water absorbing capacity difference degree of each layer is bigger, intake profile More uneven, the measure degree that takes effect is lower;Otherwise, intake profile more tends to uniformly, and profile modification is also better.Therefore injectivity index The construction parameter when coefficient of variation reaches minimum value is that CO is generated in layer2The optimal construction parameter of transfer drive.Design process:
Above formula is substituted into respectively in the formula of the injectivity index coefficient of variation, it can be found that in formation parameter and inflating medium, outgassing In the case that agent performance parameter is certain, injectivity index coefficient of variation VJIt is inflating medium, bubble release agent discharge capacity v1(z) v2(z), angry Agent, bubble release agent inject number of days d1 d2, slug number z function, as VJ[v1(z), v2(z), d1, d2, z].Therefore by solving mesh Scalar functions VJ[v1(z), v2(z), d1, d2, z] minimum value, inflating medium, bubble release agent discharge capacity v can be obtained1 v2, inflating medium, Bubble release agent injects number of days d1 d2, slug number z best construction parameter.Since objective function is non-linear optimal problem, in order to It is enough fast and accurately to obtain globally optimal solution, it is solved herein using genetic algorithm.
Inflating medium, bubble release agent injection pressure computation model
In inflating medium and bubble release agent injection process, the pressure of injected system is balance at any time, therefore can root The pressure equation in work progress is established according to pressure balance principle and seepage theory:
Pwht=Pft+ΔP+Pr-Pht (20)
In formula: PwhtFor medicament well intake pressure, MPa;PhtFor the head of liquid of medicament in the wellbore, MPa;PftFor medicine The frictional resistance of agent in the wellbore, MPa;Δ P is displacement pressure difference, MPa;PrFor current strata pressure, MPa.
The head of liquid P of medicament in the wellborehtAre as follows:
In formula: ρtFor the density of medicament, g/cm3;H is the depth of transfer drive interval, m.
The frictional resistance P of medicament in the wellboreftAre as follows:
In formula: q is the injection rate of medicament, m3/s;ν is the kinematic viscosity of medicament, m2/s;D is tubing diameter, m;β, m are The coefficient of friction resistance, value are related with fluid flow state.
Since in transfer drive medicament injection process, Reynolds number is smaller, fluidised form is mainly rendered as laminar flow or two kinds of hydraulically smooth. When laminar flow, β=4.15, m=1;When hydraulically smooth, β=0.0246, m=0.25.
Displacement pressure differential deltap P are as follows:
In formula: v --- the discharge capacity of medicament, m3/h;
By Pht、ΔPw、ΔP、PftCalculation formula substitute into formula (20), arrange and obtain the calculation formula of well intake pressure:
By formula (24) as it can be seen that in the case where optimum capacity determines, medicament well intake pressure PwhtIt has also determined that.From And the operation pressure in optimum capacity injection can be calculated and draw operation pressure curve.
Formula (24) is solved such as Fig. 4 using genetic algorithm, basic operation step can be expressed as follows:
(1) initialization operation determines genetic manipulation methods and the crossover probabilities, mutation probability such as selection, intersection, mutation operator Etc. genetic parameters;
(2) coding strategy is selected, feasible solution set is converted chromosome knob conformational space, generates initial population;
(3) fitness function is defined, convenient for calculating the fitness value of individual;
(4) termination condition judges, if whether group's performance meets a certain index, or is completed scheduled the number of iterations, Interative computation is then terminated, the operation of (5) step is otherwise transferred to;
(5) selection duplication operation handles population using corresponding selection operator;
(6) crossing operation handles population using corresponding crossover operator;
(7) mutation operator handles population using corresponding mutation operator;
(8) parent population forms next-generation group after the operation of selection, intersection and mutation operator, then goes to step (2).
The design method of the application can also be applied to reservoir waterflooding, the life obtained using the design method of the application Gas agent dosage, bubble release agent dosage, slug number carry out generating CO in layer2After transfer drive, water injection well water injection rate, accumulation augmented injection obviously increase Add, the net oil increasing effect of well group oil well is obvious.
Compared to generating CO in conventional layer2Transfer drive technological parameter calculates, and the design method of the application has the advantages that (1) Profile control agent (inflating medium and bubble release agent) injection rate may be implemented and injection slug number calculates;(2) work progress injection pressure can be achieved The real-time prediction of power calculates;(3) comprehensively considered static and dynamic factor, calculation method more to match with process characteristic, more Science.
Other features and advantage will illustrate in the following description, also, partly become from specification It obtains it is clear that being understood and implementing the application.The purpose of the application and other advantages can be by specifications, right Specifically noted structure is achieved and obtained in claim and attached drawing.
Detailed description of the invention
Attached drawing is used to provide to further understand technical scheme, and constitutes part of specification, with this The embodiment of application is used to explain the technical solution of the application together, does not constitute the limitation to technical scheme.
Fig. 1 is " multiple zone " ideal geology model schematic.
Fig. 2 is process parameters design thinking figure.
Fig. 3 is pressure system schematic diagram.
Fig. 4 is genetic algorithm schematic diagram.
Fig. 5 is the curve that operational discharge capacity changes over time in the embodiment of the present application.
Fig. 6 is the curve that operation pressure changes over time in the embodiment of the present application.
Fig. 7 is well water flooding curve of constructing in the embodiment of the present application, and the construction time is on May 24th, 2016 to May 27.
Fig. 8 is well group oil-producing curve of constructing in the embodiment of the present application, and the construction time is on May 24th, 2016 to May 27.
Specific embodiment
For the purposes, technical schemes and advantages of the application are more clearly understood, below in conjunction with attached drawing to the application Embodiment be described in detail.It should be noted that in the absence of conflict, in the embodiment and embodiment in the application Feature can mutual any combination.
Embodiment
Below in relation to illustrating General Principle of the invention for water injection well, it should be noted that, the application be not limited to Lower water injection well.It is calculated using following water injection well data, the basic data of water injection well is as shown in Table 5-1:
The basic data of 1 water injection well of table
Substratum Thickness/m Porosity Permeability/mD Day water absorption/m/d Permeability/mD after correction Transfer drive/solution blocking layer
L30 1.3 0.247 62 0 50 Solve blocking layer
L40 14.5 0.287 1572.15 423.8 175.94 Solve blocking layer
L50 15.1 0.285 2206.60 3284.45 1309.35 Transfer drive layer
L60 13.5 0.273 954 1907.1 850.37 Transfer drive layer
L70 8.5 0.246 678.01 529.75 375.16 Solve blocking layer
L80+L90 34.6 0.256 849.14 3496.35 608.29 Transfer drive layer
L100 23.5 0.252 719.38 741.65 189.98 Transfer drive layer
L110 4.3 0.225 438 211.9 296.64 Solve blocking layer
It according to calculation formula herein, is solved using genetic algorithm, the result being finally calculated is as follows:
2 water injection well dosage of table
Each slug inflating medium of 3 water injection well of table and bubble release agent discharge capacity
Slug number Inflating medium injection rate/m3/h Bubble release agent injection rate/m3/h
1 16.14 19.49
2 16.35 21.13
3 16.37 21.25
4 16.49 18.77
5 16.66 21.56
Water injection well is according to generating CO in design result implementation level in text2Water injection rate is by 982m after transfer drive measure3/ d rises to 1732m3/ d augmented injection 107 days, accumulates augmented injection 39254m3;6 mouthfuls of producing well in well group, wherein taking effect 5 mouthfuls, accumulation has a net increase of oil 8000m3, validity period 5 months.
Although embodiment disclosed by the application is as above, the content only for ease of understanding the application and use Embodiment is not limited to the application.Technical staff in any the application fields, is taken off not departing from the application Under the premise of the spirit and scope of dew, any modification and variation, but the application can be carried out in the form and details of implementation Scope of patent protection, still should be subject to the scope of the claims as defined in the appended claims.

Claims (10)

1. generating CO in a kind of layer2Transfer drive parameters design method, which comprises
(1) practical oil reservoir is reduced to the multiple zone geological model comprising transfer drive layer reconciliation blocking layer;
(2) inflating medium and bubble release agent dosage are subjected to dynamic splitting, pass through the difference journey of each substratum injectivity index after analysis transfer drive Degree, establish when the injectivity index coefficient of variation minimize when inflating medium injection rate, bubble release agent injection rate, inflating medium injection length, Bubble release agent injection length and injection slug number computation model;
(3) inflating medium, bubble release agent injection pressure computation model are established;
(4) model above is solved by genetic algorithm.
2. according to the method described in claim 1, wherein, in step (1), being defined to the substratum for needing obviously to block in stratum For transfer drive layer, while CO2Foam blocking acts on unconspicuous substratum and is defined as solution blocking layer.
3. in step (2), carry out single note inflating medium process first according to the method described in claim 1, wherein, then into Row note bubble release agent process;Under a certain instantaneous time step, the inflating medium discharge capacity of the inflating medium split point total amount and be equal to the moment v1(z), the bubble release agent discharge capacity v for splitting point total amount and being equal to the moment of the bubble release agent2(z),
The total injection rate of the inflating medium is split by the injectivity index distribution before transfer drive measure and assigns to each substratum on unit time step:
In formula: JiFor the injectivity index of transfer drive layer i before transfer drive measure, m3/(h·MPa);JjFor the suction for solving blocking layer j before transfer drive measure Aqua index, m3/(h·MPa);v1It (z) is the discharge capacity of z-th of slug inflating medium, m3/h;n1For the small number of plies of transfer drive layer;n2For solution The small number of plies of blocking layer;QliFor the inflating medium injection rate of transfer drive layer i, m3;kiFor the permeability of transfer drive layer i, mD;hiFor transfer drive layer i's Effective thickness, m;kjFor the permeability for solving blocking layer i, mD;hjFor the effective thickness for solving blocking layer i, m;μwFor the viscosity of water, mPas; rwFor water injection well wellbore radius, m;reFor drainage radius, m;A is unit conversion coefficient a=0.0036.
4. according to the method described in claim 3, wherein, during single note inflating medium, injecting after the inflating medium respectively The injectivity index of transfer drive layer is J 'i:
Wherein the inflating medium blocking radius under the time step on each substratum is rfi:
Solution blocking layer injectivity index before and after inflating medium injection remains unchanged:
In formula, J 'iFor the injectivity index of transfer drive layer after injection inflating medium, m3/(h·MPa);JjFor the suction before and after solution blocking layer measure Aqua index, m3/(h·MPa);RRF1For inflating medium residual resistance factor;rfiFor the inflating medium blocking radius of transfer drive layer i, m;φi For the porosity of transfer drive layer i.
5. according to the method described in claim 3, wherein, after completing the inflating medium injection process, outgassing is injected into stratum Agent solution, during the note bubble release agent, under first time step, the bubble release agent is described according to having infused in splitting point for each layer Each layer injectivity index after inflating medium split point:
In formula: v2It (z) is the discharge capacity of z-th of slug bubble release agent, m3/h;Q2iFor the bubble release agent injection rate of transfer drive layer i, m3;Q2jFor solution The bubble release agent injection rate of blocking layer j, m3
Then calculate each layer injectivity index after first time step:
In formula: θ is bubble release agent acidizing degree;RRF is foam residual resistance factor,
It is wherein reacted due to the inflating medium with the bubble release agent, generates CO2Foam system, by transfer drive layer under a time step On foam blocking radius reconciliation blocking layer plug removal radius are as follows:
In formula: r 'fiFor the foam blocking radius of transfer drive layer i, m;rsjFor the plug removal radius for solving blocking layer j, m.
6. according to the method described in claim 3, in the future time step of first time step, the bubble release agent carries out dosage Timesharing is split, split point according to each layer injectivity index after the variation of a upper time step, after sometime walking, when the tune Total dosage that the bubble release agent dosage in layer is greater than inflating medium is driven, i.e.,When, the suction of the transfer drive layer reconciliation blocking layer Aqua index computation model are as follows:
The injectivity index calculation formula in the transfer drive layer changes at this time, and the acid of bubble release agent also occurs in transfer drive layer at this time Change effect, transfer drive layer foam blocking radius and bubble release agent acidizing radius are as follows:
7. according to the method described in claim 1, wherein, after slug recycles, injectivity index coefficient of variation computation model are as follows:
In formula
In formula,For the mean water absorption index on stratum after transfer drive, m3/(h·MPa);VJFor the injectivity index coefficient of variation.
8. according to the method described in claim 1, wherein, in step (3), being established according to pressure balance principle and seepage theory Calculation of pressure model in work progress:
Pwht=Pft+ΔP+Pr-Pht (20)
Wherein, PwhtFor medicament well intake pressure, MPa;PhtFor the head of liquid of medicament in the wellbore, MPa;PftExist for medicament Frictional resistance in pit shaft, MPa;Δ P is displacement pressure difference, MPa;PrFor current strata pressure, MPa.
9. the method stated according to claim 1, wherein in step (4), described solve includes chromosome coding, ideal adaptation Degree evaluation, selection duplication operation, crossing operation and mutation operator.
10. application of the method described in any one of claims 1 to 9 in reservoir waterflooding.
CN201811316949.7A 2018-11-06 2018-11-06 In-layer generation of CO2Method for designing profile control process parameters Active CN109630078B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811316949.7A CN109630078B (en) 2018-11-06 2018-11-06 In-layer generation of CO2Method for designing profile control process parameters

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811316949.7A CN109630078B (en) 2018-11-06 2018-11-06 In-layer generation of CO2Method for designing profile control process parameters

Publications (2)

Publication Number Publication Date
CN109630078A true CN109630078A (en) 2019-04-16
CN109630078B CN109630078B (en) 2021-07-09

Family

ID=66067399

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811316949.7A Active CN109630078B (en) 2018-11-06 2018-11-06 In-layer generation of CO2Method for designing profile control process parameters

Country Status (1)

Country Link
CN (1) CN109630078B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113673096A (en) * 2021-08-09 2021-11-19 中海油田服务股份有限公司 Method for calculating treatment radius of plugging removal and injection increasing agent

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2283948C2 (en) * 2004-10-25 2006-09-20 Общество с ограниченной ответственностью "Оренбурггазпром" (ООО "Оренбурггазпром") Method for gas condensate deposit development
CN103967458A (en) * 2014-02-25 2014-08-06 中国海洋石油总公司 Sand prevention section water drive method
CN105156081A (en) * 2014-06-10 2015-12-16 中国石油化工股份有限公司 Simulating and evaluating method for acidification of carbonate heavy-oil reservoir
CN106761618A (en) * 2016-12-26 2017-05-31 中国石油天然气股份有限公司 For the profile control method for determination of amount and its device of water injection well
CN107859506A (en) * 2017-11-15 2018-03-30 中国石油天然气股份有限公司 The determination method of carbon dioxide flooding layered gas-injection well gas injection parameter

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2283948C2 (en) * 2004-10-25 2006-09-20 Общество с ограниченной ответственностью "Оренбурггазпром" (ООО "Оренбурггазпром") Method for gas condensate deposit development
CN103967458A (en) * 2014-02-25 2014-08-06 中国海洋石油总公司 Sand prevention section water drive method
CN105156081A (en) * 2014-06-10 2015-12-16 中国石油化工股份有限公司 Simulating and evaluating method for acidification of carbonate heavy-oil reservoir
CN106761618A (en) * 2016-12-26 2017-05-31 中国石油天然气股份有限公司 For the profile control method for determination of amount and its device of water injection well
CN107859506A (en) * 2017-11-15 2018-03-30 中国石油天然气股份有限公司 The determination method of carbon dioxide flooding layered gas-injection well gas injection parameter

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
冯其红等: ""层间调剖注入参数优化设计"", 《油气地质与采收率》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113673096A (en) * 2021-08-09 2021-11-19 中海油田服务股份有限公司 Method for calculating treatment radius of plugging removal and injection increasing agent

Also Published As

Publication number Publication date
CN109630078B (en) 2021-07-09

Similar Documents

Publication Publication Date Title
CN101089362B (en) Improved steam oil production method
CN104847341B (en) Underground natural gas storage tank well rational productivity forecast value revision method
Forouzanfar et al. Joint optimization of number of wells, well locations and controls using a gradient-based algorithm
CN105626006A (en) CO2 drive technological limit well spacing determination method for low-permeability oil reservoir
Wang et al. Experimental investigations of hydraulic fracture growth complexity in slickwater fracturing treatments
CN100419208C (en) Multi-block equal-fluidity energy-gathering parallel synchronous oil displacement method for heterogeneous oil reservoir
Jurus et al. Modeling water flow in hydraulically-fractured shale wells
CN105937388B (en) A kind of integrated development approach of tight sandstone reservoir
CN110439515A (en) Note adopts parameter optimization method and device
CN114048695B (en) Effective shale gas seam net volume inversion method based on flowback data
GB2520057A (en) Method for Simulating and/or Controlling Fluid Injection
CN116306385A (en) Oil reservoir fracturing imbibition energy increasing numerical simulation method, system, equipment and medium
CN109630078A (en) A kind of layer is interior to generate CO2Transfer drive parameters design method
CN103104238B (en) A kind of microorganism oil displacement numerical simulation method
CN106958437A (en) A kind of wellfracturing improves recovery ratio new method
CN104675372B (en) A kind of method for injecting produced quantity distribution for polymer flooding
Rezvani et al. A novel analytical technique for determining inflow control devices flow area in CO2-EOR and CCUS projects
CN109252840B (en) Device and method for improving thin-layer heavy oil reservoir recovery ratio based on gel foam flooding
CN110188996A (en) Water-drive pool energy consumption-yield-benefit integration characterizing method
Li et al. Modeling dynamic fracture growth induced by non-Newtonian polymer injection
CN106468160A (en) A kind of determination CO2Drive method and the CO of foam stream oil ingredient2The analogy method driven
CN115935843A (en) Water injection induced dynamic fracture seepage numerical simulation method considering seepage mechanism
CN109209361A (en) A kind of Fractured extra-low permeability oil reservoirs formation parameter prediction technique
CN111027887B (en) Oil displacement mechanism determining method and system
CN109214016A (en) A kind of extra-viscous oil reservoir chemistry is cold to adopt optimization 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