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 PDFInfo
- 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
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
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 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:
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:
Δ 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.
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)
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)
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 |
-
2012
- 2012-02-27 CN CN2012100462590A patent/CN102587888A/en active Pending
Patent Citations (2)
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)
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)
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 |