CN102587888A - Method for determining migration velocity of displacement front in CO2 miscible displacement physical simulation - Google Patents

Method for determining migration velocity of displacement front in CO2 miscible displacement physical simulation Download PDF

Info

Publication number
CN102587888A
CN102587888A CN2012100462590A CN201210046259A CN102587888A CN 102587888 A CN102587888 A CN 102587888A CN 2012100462590 A CN2012100462590 A CN 2012100462590A CN 201210046259 A CN201210046259 A CN 201210046259A CN 102587888 A CN102587888 A CN 102587888A
Authority
CN
China
Prior art keywords
displacing
displacement
visual
moment point
imaging
Prior art date
Application number
CN2012100462590A
Other languages
Chinese (zh)
Inventor
赵越超
宋永臣
刘瑜
朱宁军
王同雷
王合明
宿博
Original Assignee
大连理工大学
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 大连理工大学 filed Critical 大连理工大学
Priority to CN2012100462590A priority Critical patent/CN102587888A/en
Publication of CN102587888A publication Critical patent/CN102587888A/en

Links

Abstract

The invention discloses a method for determining the migration velocity of a displacement front in CO2 miscible displacement physical simulation, belonging to the detection field of oil pool physical simulation. The method for determining the migration velocity of a displacement front in CO2 miscible displacement physical simulation comprises the following steps of: firstly, determining the oil saturation distribution in different positions in a porous medium by use of the imaging technology; obtaining the oil saturate change curve in a sand filling model in the imaging view field; then, determining the characteristic moment point passing through the imaging view field of the displacement front; and finally, obtaining the migration velocity of the displacement front by use of an analysis method. The method can determine the migration velocity of the displacement front in the CO2 miscible displacement process and accurately obtain the flowing law of CO2 in the oil-containing porous medium so as to provide analysis basis for the nondestructive detection in the CO2 miscible displacement physical simulation process.

Description

A kind of definite CO 2The method of displacing front migration velocity in the miscible displacement of reservoir physical analogy
Technical field
The present invention relates to a kind of definite CO 2The method of displacing front migration velocity in the miscible displacement of reservoir physical analogy belongs to the reservoir physical simulation detection range.
Background technology
With CO 2Inject oil reservoir and not only can improve oil recovery factor, produce huge economic benefit, and can realize CO 2Geology is buried, and solves CO 2The environmental problem that reduces discharging, therefore, this method should be to solve human development and enhancement of environment this contradicts an effective way of problem.Pass through CO 2The technology of driving can further improve the 10%-15% of oil recovery on the basis of routine techniques, therefore further further investigation utilizes CO 2The relevant issues that improve the oil recovery factor technology are very necessary.Yet adopt the normal experiment method; Can only regard the actual formation model as one " black box "; Only through monitoring injection, the flow of extraction and the situation that pressure is analyzed the displacement of reservoir oil; This experimental technique can not accurately be grasped occurrence status and the migration process of fluid in porous media, obviously has certain limitation.
Summary of the invention
In order to overcome the problem that exists in the above-mentioned prior art, the present invention provides a kind of definite CO 2The method of displacing front migration velocity in the miscible displacement of reservoir physical analogy.This method utilizes technology such as magnetic resonance imaging (MRI) or CT imaging to CO 2Flowing law in the oil-containing porous media is held accurately, and through to saturation data analysis that imaging technique obtains, can confirm CO 2Displacing front migration velocity in the miscible displacement of reservoir process.
Technical scheme of the present invention is: a kind of definite CO 2The method of displacing front migration velocity comprises following concrete steps in the miscible displacement of reservoir physical analogy:
(1) obtains in the imaging visual field oil saturation change curve in the sandpack column
All lamella signal strength values in a certain moment in the displacement process are superimposed, and the whole three-dimensional oil-containing of filling model in the visual field that obtains forming images is distributed in the two-dimensional projection image perpendicular to the displacement direction, this projected image and CO 2Do not inject the ratio of the projected image of initial time, be the saturation distribution of projected image, it is averaged in the projection plane of the imaging visual field, oil saturation change curve in the sandpack column in the visual field that obtains to form images;
(2) confirm nine the characteristic moment point of displacing front through the imaging visual field
Through oil saturation change curve in the sandpack column in the imaging visual field is analyzed; Find with displacing front through corresponding nine the obvious characteristics moment point of imaging visual field process: characteristic moment point O moment corresponding is that the mixed phase of displacing front drives head and arrives the lower end, the visual field of forming images, and characteristic moment point A is CO 2Leading edge begins to get into imaging lower end, the visual field, and characteristic moment point B is that the displacing front afterbody gets into imaging lower end, the visual field fully, and characteristic moment point C is that the mixed phase of displacing front drives upper end, the head arrival imaging visual field, and characteristic moment point D is CO 2Leading edge begins to arrive upper end, the imaging visual field; Characteristic moment point E is that the displacing front afterbody is broken through upper end, the imaging visual field fully; Characteristic moment point F is that the displacing front afterbody is broken through whole glass sand filling model fully, and characteristic moment point G is that secondary moves back saturated end, and characteristic moment point H is that displacement test finishes;
(3) the displacing front migration velocity is obtained in analysis
Guaranteeing to inject CO 2Under the condition of constant flow, the curve linear Fitting Analysis between characteristic moment point B on the saturation curves and the characteristic moment point C is obtained the displacing front migration velocity.
The invention has the beneficial effects as follows: this definite CO 2The method of displacing front migration velocity in the miscible displacement of reservoir physical analogy; At first utilize imaging technique to confirm that the oil saturation of diverse location in the porous media distributes; Oil saturation change curve in the sandpack column in the acquisition imaging visual field; Confirm the characteristic moment point of displacing front then, utilize analytical method to obtain the displacing front migration velocity at last through the MRI imaging visual field.This method can be confirmed CO in the sandpack column 2Displacing front migration velocity in the miscible displacement of reservoir process is to CO 2Flowing law in the oil-containing porous media is held accurately, is CO 2The nondestructive testing of miscible displacement of reservoir physical analogy process provides the analysis foundation.
Description of drawings
Below in conjunction with accompanying drawing and embodiment the present invention is described further.
Fig. 1 is a kind of definite CO 2The flow chart of the method for displacing front migration velocity in the miscible displacement of reservoir physical analogy.
Fig. 2 is that the mixed phase of displacing front drives head arrival imaging lower end, visual field characteristic moment point O moment corresponding sketch map.
Fig. 3 is CO 2Leading edge begins to get into imaging lower end, visual field characteristic moment point A moment corresponding sketch map.
Fig. 4 is that the displacing front afterbody gets into imaging lower end, visual field characteristic moment point B moment corresponding sketch map fully.
Fig. 5 is that the mixed phase of displacing front drives upper end, head arrival imaging visual field characteristic moment point C moment corresponding sketch map.
Fig. 6 is CO 2Leading edge begins to arrive upper end, imaging visual field characteristic moment point D moment corresponding sketch map.
Fig. 7 is that point is broken through upper end, imaging visual field characteristic moment point E moment corresponding sketch map fully for the displacing front afterbody.
Fig. 8 is an oil saturation change curve in the sandpack column in the imaging visual field.
Among the figure: 1, back-up sand, 2, the mixed phase district, 3, CO 2Leading edge, 4, the imaging visual field; O, A, B, C, D, E, F, G, H, characteristic moment point.
The specific embodiment
Fig. 1 is a kind of definite CO 2The flow chart of the method for displacing front migration velocity in the miscible displacement of reservoir physical analogy.Introduce the specific embodiment that this method is set up below:
(1) oil saturation of confirming diverse location in the porous media distributes
With the MRI technology is example, at first obtains piece image I O0, in displacement process, obtain another width of cloth image I then Oi, oil saturation is in so any little volume elements:
S oi ( x , y , z ) = I oi ( x , y , z ) I o 0 ( x , y , z ) - - - ( 1 )
S in the formula Oi(x, y z) are the oil saturation of arbitrary volume elements in the space, I Oi(x, y z) are the signal strength signal intensity size of arbitrary volume elements in the displacement process, I O0(x, y, the signal strength signal intensity size of arbitrary corpusculum unit when being initial saturated oils z).
(2) obtain in the imaging visual field oil saturation change curve in the sandpack column
All lamella signal strength values in a certain moment in the displacement process are superimposed, and the whole three-dimensional oil-containing of filling model in the visual field that can obtain forming images is distributed in the two-dimensional projection image perpendicular to the displacement direction, this projected image and initial time (CO 2Inject) time the ratio of projected image, be the saturation distribution of projected image, it is averaged oil saturation change curve in the sandpack column in the visual field that can obtain to form images in the projection plane of the imaging visual field.
(3) confirm nine the characteristic moment point of displacing front through the imaging visual field
Analyze through oil saturation change curve in the sandpack column in the visual field that Fig. 3 is formed images; Can find to pass through the imaging visual field (shown in Fig. 2-7 with displacing front; The imaging visual field is to be positioned at certain a part of pipeline section on the sandpack column pipe) corresponding nine the obvious characteristics moment point of process: characteristic moment point O moment corresponding is that the mixed phase of displacing front drives head arrival imaging lower end, the visual field, and characteristic moment point A is CO 2Leading edge begins to get into imaging lower end, the visual field, and characteristic moment point B is that the displacing front afterbody gets into imaging lower end, the visual field fully, and characteristic moment point C is that the mixed phase of displacing front drives upper end, the head arrival imaging visual field, and characteristic moment point D is CO 2Leading edge begins to arrive upper end, the imaging visual field; Characteristic moment point E is that the displacing front afterbody is broken through upper end, the imaging visual field fully; Characteristic moment point F is that the displacing front afterbody is broken through whole glass sand filling model fully, and characteristic moment point G is that secondary moves back saturated end, and characteristic moment point H is that displacement test finishes.Thereby can displacement process be divided into eight processes and analyze:
1. OA section: along with CO 2The injection oil saturation descends, because the inhomogeneities that plug shape mixed phase drives the leading edge head, decline curve is exponential trend, and downward gradient increases gradually;
2. AB section: along with CO 2Inject oil saturation and continue to descend, same because plug shape CO 2The inhomogeneities of leading edge head, decline curve is exponential trend, and downward gradient increases gradually;
3. BC section: oil saturation continues to descend, and overall distribution of pores is more even because glass sand is filled the model porous media, CO in the displacement process 2With oily mixed phase, do not have serious fingering and channelling phenomenon to take place, thereby displacement process remain a constant speed basically, the full curve of oil-containing is exponential trend and descends;
4. CD section: break through imaging top, visual field process for mixed phase drives leading edge, because mixed phase drives the inhomogeneities of leading edge head, decline curve is exponential trend;
5. DE section: be CO 2Leading edge breaks through imaging top, visual field process, because CO 2The inhomogeneities of leading edge head, decline curve is exponential trend;
6. EF section: after oil saturation was reduced to a certain value, under constant displacement velocity condition, mixed phase drove the primary recovery rate of leading edge basically near remaining oil saturation, follow-up CO 2Little to improving the recovery ratio effect;
7. FG section: be " secondary moves back saturated " process of oil in the porous media, oil saturation further descends;
8. GH section: after final oil saturation is reduced to a certain value, follow-up CO 2Little to improving the recovery ratio effect, " secondary moves back saturation history " end be described after.
(4) the displacing front migration velocity is obtained in analysis
Analysis through to BC section on the saturation curves is visible, is guaranteeing to inject CO 2Under the condition of constant flow, though a part of CO 2Can be dissolved in the oil; And cause the displacing front distribution irregular because distribution of pores is inhomogeneous; But the displacing front propelling that in whole filling model, remains a constant speed basically in displacement process; And through to BC section linear relationship further linear fit analysis can obtain the displacing front migration velocity, detailed method is following:
Suppose in the displacement process in certain period that the oil content variable quantity is Δ V in the imaging visual field o, and available following formula is represented:
ΔV o=AφΔh (2)
Wherein A is the cross-sectional area that glass sand is filled the model porous media, cm 2φ is a porosity, mark; Δ h is the displacement of displacing front, cm.
In addition, oil saturation S oDefinition be:
S o = V o V p - - - ( 3 )
V wherein oFor glass sand in the FOV is filled the oil content in the model porous media, cm 3V pFor glass sand in the imaging visual field is filled the voids volume of model porous media, cm 3
V pCan further express with following formula:
V p=ALφ (4)
Wherein L is imaging visual field height, cm.
According to formula (2)-(4), can obtain average fltting speed expression formula of displacing front:
v ‾ = Δh Δt = ΔS o Δt L - - - ( 5 )
Δ S wherein oBe the variable quantity of oil saturation in the imaging visual field in the displacement process, mark; Δ t is the displacement process time variation amount, min.
in the formula (5) can obtain through the linear fit to BC section on the saturation curves, and then can be in the hope of the average fltting speed of displacing front.Because mixed phase drives in the process, a large amount of CO 2With oil dissolving mixed phase, this speed is much smaller than CO 2Injection rate.In addition; Can obtain displacing front and migrate to model pipe from imaging lower end, the visual field and export the used time through analyzing saturation curves BF section; According to known migration velocity, can obtain migration distance, this filling mold segment pipe range of this distance and reality is approaching basically; Thereby, further utilize the reliability of saturation curves analysis and solution displacement velocity method more than the proof.
Different characteristic when Fig. 2-7 passes through the imaging visual field for displacing front constantly.Corresponding with the obvious characteristics moment point on the oil saturation change curve in the sandpack column in the imaging visual field among Fig. 8.
Fig. 8 is oil saturation change curve in the sandpack column in the imaging visual field.Through oil saturation change curve in the sandpack column in the imaging visual field among the figure is analyzed; Can find to pass through the imaging visual field (shown in Fig. 2-7, the imaging visual field is to be positioned at certain a part of pipeline section on the sandpack column pipe) corresponding nine the obvious characteristics moment point of process with displacing front.

Claims (1)

1. definite CO 2The method of displacing front migration velocity in the miscible displacement of reservoir physical analogy is characterized in that: comprise following concrete steps:
(1) obtains in the imaging visual field oil saturation change curve in the sandpack column
All lamella signal strength values in a certain moment in the displacement process are superimposed, and the whole three-dimensional oil-containing of filling model in the visual field that obtains forming images is distributed in the two-dimensional projection image perpendicular to the displacement direction, this projected image and CO 2Do not inject the ratio of the projected image of initial time, be the saturation distribution of projected image, it is averaged in the projection plane of the imaging visual field, oil saturation change curve in the sandpack column in the visual field that obtains to form images;
(2) confirm nine the characteristic moment point of displacing front through the imaging visual field
Through oil saturation change curve in the sandpack column in the imaging visual field is analyzed; Find with displacing front through corresponding nine the obvious characteristics moment point of imaging visual field process: characteristic moment point O moment corresponding is that the mixed phase of displacing front drives head and arrives the lower end, the visual field of forming images, and characteristic moment point A is CO 2Leading edge begins to get into imaging lower end, the visual field, and characteristic moment point B is that the displacing front afterbody gets into imaging lower end, the visual field fully, and characteristic moment point C is that the mixed phase of displacing front drives upper end, the head arrival imaging visual field, and characteristic moment point D is CO 2Leading edge begins to arrive upper end, the imaging visual field; Characteristic moment point E is that the displacing front afterbody is broken through upper end, the imaging visual field fully; Characteristic moment point F is that the displacing front afterbody is broken through whole glass sand filling model fully, and characteristic moment point G is that secondary moves back saturated end, and characteristic moment point H is that displacement test finishes;
(3) the displacing front migration velocity is obtained in analysis
Guaranteeing to inject CO 2Under the condition of constant flow, the curve linear Fitting Analysis between characteristic moment point B on the saturation curves and the characteristic moment point C is obtained the displacing front migration velocity.
CN2012100462590A 2012-02-27 2012-02-27 Method for determining migration velocity of displacement front in CO2 miscible displacement physical simulation CN102587888A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012100462590A CN102587888A (en) 2012-02-27 2012-02-27 Method for determining migration velocity of displacement front in CO2 miscible displacement physical simulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012100462590A CN102587888A (en) 2012-02-27 2012-02-27 Method for determining migration velocity of displacement front in CO2 miscible displacement physical simulation

Publications (1)

Publication Number Publication Date
CN102587888A true CN102587888A (en) 2012-07-18

Family

ID=46477057

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012100462590A CN102587888A (en) 2012-02-27 2012-02-27 Method for determining migration velocity of displacement front in CO2 miscible displacement physical simulation

Country Status (1)

Country Link
CN (1) CN102587888A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104727790A (en) * 2013-12-23 2015-06-24 中国石油化工股份有限公司 Method for predicting period of re-enriching remaining oil of water flooding reservoir into new reservoir
CN105041280A (en) * 2015-06-03 2015-11-11 东北石油大学 Method and device for realizing carbon dioxide miscible-phase displacement laboratory experiment
CN106383220A (en) * 2016-10-09 2017-02-08 东北石油大学 Method and system for confirming non-Newtonian fluid injection speed in laboratory experiment

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2136868C1 (en) * 1998-04-07 1999-09-10 НГДУ "Иркеннефть" Method of developing oil deposit
CN201991509U (en) * 2011-03-16 2011-09-28 东北石油大学 Visual planar sand-filling model used for oil displacement experiments

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2136868C1 (en) * 1998-04-07 1999-09-10 НГДУ "Иркеннефть" Method of developing oil deposit
CN201991509U (en) * 2011-03-16 2011-09-28 东北石油大学 Visual planar sand-filling model used for oil displacement experiments

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YUECHAO ZHAO ET AL.: "Visualization and Measurement of CO2 Flooding in Porous Media", 《INDUSTRIAL & ENGINEERING CHEMISTRY RESEARCH》, vol. 50, no. 8, 17 March 2011 (2011-03-17) *
赵越超等: "磁共振成像在CO_2驱油实验中应用", 《大连理工大学学报》, vol. 52, no. 1, 31 January 2012 (2012-01-31) *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104727790A (en) * 2013-12-23 2015-06-24 中国石油化工股份有限公司 Method for predicting period of re-enriching remaining oil of water flooding reservoir into new reservoir
CN105041280A (en) * 2015-06-03 2015-11-11 东北石油大学 Method and device for realizing carbon dioxide miscible-phase displacement laboratory experiment
CN106383220A (en) * 2016-10-09 2017-02-08 东北石油大学 Method and system for confirming non-Newtonian fluid injection speed in laboratory experiment
CN106383220B (en) * 2016-10-09 2018-08-24 东北石油大学 The determination method and system of laboratory experiment non-newtonian fluid injection rate

Similar Documents

Publication Publication Date Title
Raeini et al. Direct simulations of two-phase flow on micro-CT images of porous media and upscaling of pore-scale forces
Krevor et al. Capillary trapping for geologic carbon dioxide storage–From pore scale physics to field scale implications
Guillen et al. Pore scale and macroscopic displacement mechanisms in emulsion flooding
Lemonnier et al. Simulation of naturally fractured reservoirs. state of the art-part 1–physical mechanisms and simulator formulation
CN103267721B (en) Method for evaluating water containing characteristic and occurrence state of compact sandstone storage layer aperture
CN104533370B (en) Pressure break horizontal well oil reservoir, crack, pit shaft coupled model method
Morel-Seytoux Two-phase flows in porous media
Rossen et al. Injection strategies to overcome gravity segregation in simultaneous gas and water injection into homogeneous reservoirs
Chierici Principles of Petroleum Reservoir Engineering: Volume 2
Elenius et al. Effects of a capillary transition zone on the stability of a diffusive boundary layer
Peters Advanced petrophysics: Dispersion, interfacial phenomena
Gong et al. Upscaling discrete fracture characterizations to dual-porosity, dual-permeability models for efficient simulation of flow with strong gravitational effects
Rossen Numerical challenges in foam simulation: a review
Maini In situ hydraulic parameters in jointed rock--Their measurement and interpretation
CN103352693B (en) A kind of output prediction system based on oily content and method thereof
Blair et al. Solution of two-phase flow problems using implicit difference equations
CN105089657A (en) Physical simulation method and physical simulation experimental device for oil and gas charging of fracture-cave carbonate reservoir
CN103556994B (en) The experiment detecting system of fractured-vuggy reservoir remaining oil distribution and detection method
Sbai et al. Numerical modeling of formation damage by two-phase particulate transport processes during CO2 injection in deep heterogeneous porous media
Song et al. Impact of permeability heterogeneity on production characteristics in water-bearing tight gas reservoirs with threshold pressure gradient
CN106223938B (en) Digitize rock core flow simulating analysis method and device
CN103926267A (en) Method for quantitatively evaluating pore throat change degree in stress sensitive process
Meng et al. A critical review on fundamental mechanisms of spontaneous imbibition and the impact of boundary condition, fluid viscosity and wettability
CN104112057B (en) A kind of large scale fractured reservoir method for numerical simulation
CN103246820A (en) Oil and gas reservoir numerical simulation calculation method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20120718