CN110609046A - 一种基于微ct水驱油图像的毛管力预测方法 - Google Patents
一种基于微ct水驱油图像的毛管力预测方法 Download PDFInfo
- Publication number
- CN110609046A CN110609046A CN201911006948.7A CN201911006948A CN110609046A CN 110609046 A CN110609046 A CN 110609046A CN 201911006948 A CN201911006948 A CN 201911006948A CN 110609046 A CN110609046 A CN 110609046A
- Authority
- CN
- China
- Prior art keywords
- oil
- water
- image
- capillary force
- micro
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N13/02—Investigating surface tension of liquids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Theoretical Computer Science (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于微CT水驱油图像的毛管力预测方法。该方法首先通过微CT成像技术,获取岩心水驱油过程中的微观油水分布图像;通过滤波、分割等图像处理技术,提取孔隙空间中的油水分布;利用Generalized Marching Cubes算法在油水相界面生成光滑的面网格模型,通过二次型曲面逼近面网格,从而求解油水界面的平均曲率;通过Young‑Laplace等式,求得微观毛管力的分布。该方法与传统岩心室内试验相比,对岩样破坏性小,实现了孔隙尺度毛管力的可视化表征。
Description
技术领域
本发明涉及油气田开发、地下污染物运移、二氧化碳地质封存等工程领域,尤其涉及一种基于微CT水驱油图像的毛管力预测方法。
背景技术
由于岩石孔隙尺寸极小,地下多相流体的输运过程往往受到毛细管力的控制。因此,准确预测毛管力的大小及分布,对理解孔隙尺度流动机理,合理制定油气田开发、地下污染物治理、二氧化碳地质封存方案等至关重要。
本发明基于岩石水驱油过程中获取的油水分布微CT图像,结合ImageJ软件,通过图像处理及分析技术,提出了一种求解微观水驱油过程毛细管力分布的方法,为孔隙尺度多相渗流过程毛细管力的求解及可视化提供了一种新方法。
发明内容
本发明基于岩石水驱油过程中获取的油水分布微CT图像,结合ImageJ开源软件,通过图像处理及分析技术,获得孔隙尺度岩石水驱油过程中毛管力的可视化分布。为了实现以上目的,本发明将采取以下实施方法:
(一)开展微观水驱油实验,通过微CT扫描技术,获取岩心水驱油过程中的高分辨率油水分布微CT图像。
(二)基于ImageJ软件,利用中值滤波器对原始图像进行滤波处理,以提高图像质量;利用分水岭分割算法,将获取的微CT图像分割为油相、水相和岩石颗粒,进而提取孔隙空间中的油水分布。
(三)通过Generalized Marching Cubes算法,在油水相界面处生成光滑的面网格,通过二次型曲面近似替代面网格,求解油水相界面的平均曲率。
(四)将求解得到的平均曲率带入Young-Laplace等式,即可求得毛细管力参数。
传统上,毛管力的预测方法主要包括:多孔板法、压汞法和离心分离法。与传统方法相比,本发明的益处在于:
与多孔板法相比,本发明的试验周期短;与压汞法相比,本发明避免了对岩样的破坏,可有效地预测渗吸及排水过程中的毛管力;与离心分离法相比,该方法避免了样品边界的影响,且可模拟储层原始温压条件。此外,本发明实现了孔隙尺度毛细管力分布的可视化表征。
附图说明
为了更清晰地说明本发明的实施过程,下面将结合附图和具体实施步骤对本发明作进一步说明。
图1为本发明实施的工作流程图。
图2为岩石微观水驱油实验的装置图。
图3为滤波前后残余油状态下油水分布微CT图像。(a)为原始图像;(b)中值滤波后的图像。其中,灰色代表骨架,黑色代表水相,中间灰度值代表油相。
图4为分水岭算法分割前后的油水分布图像。(a)分割前油水分布的原始图像;(b)分割后的油水分布二维图像。其中黑色代表油相,灰白色代表水相,灰色代表岩石骨架。
图5为提取的油水相界面面网格。(a)油水分布的二维图像;(b)提取的单个油簇;(c)包含油水界面的多区域网格;(d)提取的油水相界面面网格。
图6为油簇毛管力分布云图。
图7为毛管力分布直方图。
具体实施方式
为了详细说明本发明实现的方法、达到的目的和具有的功效,下面将结合附图对本发明作进一步说明。需要指出,下文描述的实施例仅是本发明的部分实施例,并不是所有实施例。相关领域的其他研究人员、技术人员在本发明所提方法的基础上,没有其他创新性劳动的前提下获得的所有其他实施例,都应在本申请的保护范围之内。
图1是本申请一种基于微CT技术的毛管力预测方法实施例的流程图,主要包括以下步骤。
S1:开展岩心水驱油实验,通过微CT扫描技术获取各驱替阶段的油水分布图像。具体过程为:首先将岩心烘干,减少岩心中原存流体对实验结果的影响;将岩心饱和地层水,通过微CT扫描获取流体的初始分布图像;开展油驱水实验建立束缚水饱和度,通过微CT扫描技术获取束缚水饱和度状态下的油水分布图像;开展水驱油实验模拟油藏注水开发过程,待驱替结束,通过微CT扫描技术获取残余油状态下的油水分布图像。本发明实施例采用的岩心为砂岩,采用的水相为地层水,油相为混有9%体积碘代正丁烷的煤油。岩石及流体物性参数如表1所示。实验装置如图2所示。
表1岩石及流体物性参数
S2:利用中值滤波对原始图像进行滤波处理,以提高图像质量。通过对比原始图像和滤波后的图像可知,中值滤波器在滤除原始图像中的噪音、孤立小岩石颗粒及小液滴的同时,可以有效保护相边界,如图3所示。
S3:基于滤波后的图像,利用分水岭分割算法,将图像分为油相、水相及岩石颗粒。通过体渲染,可以对油水在孔隙中的分布进行三维可视化,如图4所示。
S4:基于三维可视化的油水重建图像,选取典型残余油簇,通过GeneralizedMarching Cubes算法生成包含油、水的多区域网格模型,通过布尔运算提取只包含油水相界面的面网格模型,如图5所示。
S5:基于提取的油水相界面网格模型,通过二次曲面拟合,完成相界面曲率的求解:
ax2+by2+cz2+2exy+2fyz+2gzx+2lx+2my+2nz+d=0 (1)
其中,a,b,c,d,e,f,g,l,m,n为拟合常系数。当二次型曲面无限逼近油水相界面面网格时,各参数取得最优解。此时,二次型曲面的特征值及特征向量分别为油水相界面的主曲率和主曲率方向。因此,油水相界面的平均曲率为两主曲率的平均值。
k=k1+k2 (2)
其中,k1和k2分别为油水相界面的两主曲率,k为油水相界面的平均曲率。
S6:将求得的平均曲率代入Young-Laplace公式,完成毛管力的求取:
Pc=2σk (3)
其中,Pc为毛管力,σ为油水表面张力系数。图6所示为所选残余油簇的毛管力分布云图。
通过选定其他目标区域的油簇,重复S4-S6步骤,可完成其毛管力的求取与可视化。经过统计分析,可求得其毛管力分布直方图。图7所示为选取的5个残余油簇的毛管力分布直方图,从直方图可知毛细管力主要分布在(-10kPa,10kPa)之间。
以上用于对本发明的基本方法、达到的目的和功效进行说明的实施例仅是本发明的部分实施例,并不能限制本发明的使用范围。凡是根据本发明提出的方法、原理,对以上实施例所做的任何修改、修饰,均在本发明技术方案之内。
Claims (5)
1.一种基于微CT水驱油图像的毛管力预测方法,其特征在于,该方法包括:
步骤S1,开展岩心水驱油实验,通过微CT扫描技术获取各驱替阶段的油水分布图像;
步骤S2,通过滤波、分割等图像处理技术,获取岩心孔隙中的油水分布;
步骤S3,利用Generalized Marching Cubes算法,在油水相界面处生成光滑的面网格,通过二次型曲面近似替代油、水相界面的面网格,求解油水相界面的平均曲率;
步骤S4,通过Young-Laplace等式,求得毛管力的微观分布。
2.按照权利要求1所述的一种基于微CT水驱油图像的毛管力预测方法,其特征在于:所述步骤(1)中,通过微CT扫描技术获取岩心水驱油过程中的油水分布图像,类似具有清晰界面流动的多相渗流图像也可作为毛管力计算的基础数据。
3.按照权利要求1所述的一种基于微CT水驱油图像的毛管力预测方法,其特征在于:所述步骤(2)中,通过滤波提高图像质量,通过图像分割技术提取孔隙空间中的油、水分布。
4.按照权利要求1所述的一种基于微CT水驱油图像的毛管力预测方法,其特征在于:所述步骤(3)中,通过Generalized Marching Cubes算法,生成包含油、水的多区域面网格模型,通过布尔运算提取只包含油、水相界面的面网格,通过二次型曲面近似替代面网格,求得油水相界面的平均曲率。
5.按照权利要求1所述的一种基于微CT水驱油图像的毛管力预测方法,其特征在于:所述步骤(4)中,在获得油水相界面的平均曲率后,通过Young-Laplace等式,求得毛管力的微观分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911006948.7A CN110609046B (zh) | 2019-10-22 | 2019-10-22 | 一种基于微ct水驱油图像的毛管力预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911006948.7A CN110609046B (zh) | 2019-10-22 | 2019-10-22 | 一种基于微ct水驱油图像的毛管力预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110609046A true CN110609046A (zh) | 2019-12-24 |
CN110609046B CN110609046B (zh) | 2022-12-13 |
Family
ID=68893348
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911006948.7A Active CN110609046B (zh) | 2019-10-22 | 2019-10-22 | 一种基于微ct水驱油图像的毛管力预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110609046B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111739149A (zh) * | 2020-06-15 | 2020-10-02 | 中国石油大学(华东) | 一种岩石ct扫描图像的油水分布连续性修复方法 |
CN117849079A (zh) * | 2024-03-04 | 2024-04-09 | 中国电建集团华东勘测设计研究院有限公司 | 一种基于ct成像的土体蒸发微观表征方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102628354A (zh) * | 2012-04-12 | 2012-08-08 | 东北石油大学 | 孔隙微米级油水分布识别量化方法 |
CN109540765A (zh) * | 2019-01-08 | 2019-03-29 | 西南石油大学 | 一种基于岩心微ct图像的孔渗参数预测方法 |
-
2019
- 2019-10-22 CN CN201911006948.7A patent/CN110609046B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102628354A (zh) * | 2012-04-12 | 2012-08-08 | 东北石油大学 | 孔隙微米级油水分布识别量化方法 |
CN109540765A (zh) * | 2019-01-08 | 2019-03-29 | 西南石油大学 | 一种基于岩心微ct图像的孔渗参数预测方法 |
Non-Patent Citations (2)
Title |
---|
吕祥锋等: "煤层气-水两相流渗透率测定实验研究", 《水资源与水工程学报》 * |
宋睿: "基于微尺度重建模型的岩石热—流—固耦合细观机理研究", 《中国博士学位论文全文数据库 工程科技I辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111739149A (zh) * | 2020-06-15 | 2020-10-02 | 中国石油大学(华东) | 一种岩石ct扫描图像的油水分布连续性修复方法 |
CN111739149B (zh) * | 2020-06-15 | 2023-09-01 | 中国石油大学(华东) | 一种岩石ct扫描图像的油水分布连续性修复方法 |
CN117849079A (zh) * | 2024-03-04 | 2024-04-09 | 中国电建集团华东勘测设计研究院有限公司 | 一种基于ct成像的土体蒸发微观表征方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110609046B (zh) | 2022-12-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Berg et al. | Industrial applications of digital rock technology | |
Tahmasebi et al. | Rapid multiscale modeling of flow in porous media | |
US11590708B2 (en) | Three-dimensional fluid micromodels | |
Herring et al. | Topological persistence for relating microstructure and capillary fluid trapping in sandstones | |
Bultreys et al. | Multi-scale, micro-computed tomography-based pore network models to simulate drainage in heterogeneous rocks | |
Jouini et al. | Numerical estimation of carbonate rock properties using multiscale images | |
Zhang et al. | Characterization of coal porosity and permeability evolution by demineralisation using image processing techniques: A micro-computed tomography study | |
Alyafei et al. | A sensitivity study of the effect of image resolution on predicted petrophysical properties | |
Berg et al. | Onset of oil mobilization and nonwetting-phase cluster-size distribution | |
CN103206208A (zh) | 一种微观不同赋存状态剩余油宏观量化的方法 | |
AU2013274604A1 (en) | Digital rock analysis with multiphase flow REV determination | |
CN110609046B (zh) | 一种基于微ct水驱油图像的毛管力预测方法 | |
Jiang et al. | Permeability estimation of porous media by using an improved capillary bundle model based on micro-CT derived pore geometries | |
Liu et al. | Pore‐scale visualization and quantitative analysis of the spontaneous imbibition based on experiments and micro‐CT technology in low‐permeability mixed‐wettability rock | |
CN104737037A (zh) | 使用多点模拟的储层建模 | |
Cantisano et al. | Relative permeability in a shale formation in Colombia using digital rock physics | |
Hajizadeh et al. | An algorithm for 3D pore space reconstruction from a 2D image using sequential simulation and gradual deformation with the probability perturbation sampler | |
CN113916916A (zh) | 一种用于页岩数字岩心三维渗流-颗粒流动耦合的模拟方法 | |
Esser et al. | Network model of porous media–Review of old ideas with new methods | |
CN117892654A (zh) | 基于ct扫描和多相物理场耦合的砂岩数字岩心多孔介质渗流特征的模拟方法 | |
Lei | Reconstruction and Analysis of Tight Sandstone Digital Rock Combined with X‐Ray CT Scanning and Multiple‐Point Geostatistics Algorithm | |
Youssef et al. | High resolution μ-CT combined to numerical models to assess electrical properties of bimodal carbonates | |
CN112816388A (zh) | 基于ct和数字岩心三维重构的油砂渗流性能测试方法 | |
Hu et al. | Pore-network model for geo-materials | |
CN115147539A (zh) | 一种基于指示克里金的孔隙结构剩余油表征方法 |
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 |