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
visual field
moment point
displacement
imaging
displacing front
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.)
Pending
Application number
CN2012100462590A
Other languages
Chinese (zh)
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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN2012100462590A priority Critical patent/CN102587888A/en
Publication of CN102587888A publication Critical patent/CN102587888A/en
Pending legal-status Critical Current

Links

Images

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
Figure BDA0000138650620000061
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.
Figure BDA0000138650620000063
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 Pending 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 Pending 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 (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
CN111305801A (en) * 2020-03-09 2020-06-19 中国石油化工股份有限公司 Carbon dioxide flooding front edge description method

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 (5)

* 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
CN111305801A (en) * 2020-03-09 2020-06-19 中国石油化工股份有限公司 Carbon dioxide flooding front edge description method

Similar Documents

Publication Publication Date Title
Meng et al. A critical review on fundamental mechanisms of spontaneous imbibition and the impact of boundary condition, fluid viscosity and wettability
Briaud et al. SRICOS: Prediction of scour rate in cohesive soils at bridge piers
Ju et al. Improving water injectivity and enhancing oil recovery by wettability control using nanopowders
CN103257151B (en) In a kind of quantitative evaluation oil and gas secondary migration process, pore throat employs the method for rule
CN105021506A (en) Pore network model-based calculation method of three-phase relative permeability
CN104730089B (en) Flow in parallel core foam fractionation imaging detection device and method of work thereof
Lyu et al. Mechanism and influencing factors of EOR by N2 injection in fractured-vuggy carbonate reservoirs
CN109577945B (en) Experimental device and method for discriminating low-permeability-ultra-low-permeability reservoir channeling channel
CN104533370A (en) Oil deposit, crack and shaft fully-coupled simulating method of fractured horizontal well
CN105822302A (en) Oil-water distribution recognition method based on well-ground potential method
CN101942991A (en) Method for modeling predictable physical model for fractured anisotropic oil reservoir water flooding development
CN102748007B (en) Well testing analytical method and device
CN108051643A (en) Multipair dynamic monitoring displacement system of multi-functional long cores radial direction
CN104237098A (en) Method for measuring gradient of rock starting pressure
CN102587888A (en) Method for determining migration velocity of displacement front in CO2 miscible displacement physical simulation
CN103452543A (en) Fractured horizontal well pattern model and method for manufacturing same
CN107145671B (en) A kind of numerical reservoir simulation method and system
CN104179498A (en) Zone-dividing saturation simulation method
CN105137045B (en) Experimental provision and application that in a kind of synthetic determination macroscopic view flow event, foam oil phenomenon forms
CN102383783A (en) Method for analyzing flow characteristic of oil and water in seam-hole type oil reservoir
Liu et al. Experimental study of the dynamic mechanism on gas bubbles migration, fragment, coalescence and trapping in a porous media
Yu et al. Non-Darcy flow numerical simulation of XPJ low permeability reservoir
Liu et al. Displacement front behavior of near miscible CO2 flooding in decane saturated synthetic sandstone cores revealed by magnetic resonance imaging
CN106050196B (en) A kind of Offshore Heavy Oil Field oil reservoir early stage note gathers the determination method on preferable metaideophone opportunity
Ren Local capillary trapping in carbon sequestration: Parametric study and implications for leakage assessment

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