CN117892654A - 基于ct扫描和多相物理场耦合的砂岩数字岩心多孔介质渗流特征的模拟方法 - Google Patents
基于ct扫描和多相物理场耦合的砂岩数字岩心多孔介质渗流特征的模拟方法 Download PDFInfo
- Publication number
- CN117892654A CN117892654A CN202410036071.0A CN202410036071A CN117892654A CN 117892654 A CN117892654 A CN 117892654A CN 202410036071 A CN202410036071 A CN 202410036071A CN 117892654 A CN117892654 A CN 117892654A
- Authority
- CN
- China
- Prior art keywords
- simulation method
- seepage
- sandstone
- image
- simulation
- 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
- 238000004088 simulation Methods 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000002591 computed tomography Methods 0.000 title claims abstract description 13
- 230000008878 coupling Effects 0.000 title claims abstract description 9
- 238000010168 coupling process Methods 0.000 title claims abstract description 8
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 8
- 239000011148 porous material Substances 0.000 claims abstract description 59
- 239000012530 fluid Substances 0.000 claims abstract description 36
- 239000011435 rock Substances 0.000 claims abstract description 31
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 24
- 230000004044 response Effects 0.000 claims abstract description 13
- 238000011160 research Methods 0.000 claims abstract description 9
- 238000012805 post-processing Methods 0.000 claims abstract description 5
- 230000011218 segmentation Effects 0.000 claims description 19
- 238000012545 processing Methods 0.000 claims description 12
- 230000035699 permeability Effects 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 5
- 230000007246 mechanism Effects 0.000 abstract description 5
- 238000002474 experimental method Methods 0.000 abstract description 4
- 230000007547 defect Effects 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000009736 wetting Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
-
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
- G06T17/205—Re-meshing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Pathology (AREA)
- Fluid Mechanics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Immunology (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Computer Graphics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Algebra (AREA)
- Dispersion Chemistry (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于CT扫描和多相物理场耦合的砂岩数字岩心多孔介质渗流特征的模拟方法。所述模拟方法包括如下步骤：S1、对岩心进行CT扫描得到CT图像；S2、采用Avizo软件对CT图像进行后处理，构建三维数字化岩心模型；S3、提取渗流孔隙结构侯进行网格精细划分，得到重构后的三维数字化岩心模型；S4、将重构后的三维数字化岩心模型进行渗流特征模拟研究。本发明模拟方法得到的模拟结果接近流体在砂岩地层中渗流情况，本发明建立了流体在多孔介质中其粘度与含水饱和度响应方程，弥补了在物理实验中探索砂岩多尺度复杂性质的不足，为流体在砂岩孔隙中渗流机理提供了参考，在砂岩油(气)藏流体力学技术中具有很好的应用价值。
Description
技术领域
本发明涉及一种基于CT扫描和多相物理场耦合的砂岩数字岩心多孔介质渗流特征的模拟方法，属于砂岩数字岩心渗流模拟技术领域。
背景技术
数字岩心技术是以真实岩心的CT扫描图像为基础，并基于逆向重构软件建立出岩心物理结构的三维模型，常用的CT扫描技术可达微纳米级别，现已广泛应用于地质、化学、医学等领域。砂岩此类多孔介质的内部孔隙结构由细小孔隙相和岩相组成，利用阈值分割的方法，连通性分析提取其中的孔隙，以便于后续对砂岩多尺度特性进行研究，有利于岩石物理化学渗流理论研究开发探索。
研究多孔介质中孔隙结构的形态、连通性、渗透率、孔隙度等宏观性质是探究多孔岩心的渗流能力以预测储层流体。孔隙是指岩石内部未被固体物质占据的空间结构，而喉道就是链接相邻孔隙的部位。孔隙反映了岩石的储集能力，喉道反映了岩石的渗流特性。低渗砂岩是一种具有复杂的微观孔隙结构的低孔、低渗的岩石，且它的孔喉细小、连通性差、几何形状多变。
由于砂岩多尺度特征、复杂的微观孔隙结构使得常规室内实验无法对其性质变化以及渗流表现出的非线性等现象进行准确的描述。但探究岩石微观参数、渗流机理特征是必不可缺的，鉴于此，当前先进的微观成像技术和计算机科学的发展，结合三维可视化数字软件，反映岩石实际的孔隙空间结构分布，建立三维数字化岩心模型并与多物理场耦合模拟进行微观渗流分析，势在必行。
发明内容
本发明的目的是提供一种基于CT扫描和多相物理场耦合的砂岩数字岩心多孔介质渗流特征的模拟方法，可探索驱替过程中流体进入砂岩后的渗流机理和自身粘度与含水饱和度之间的关系，以适应砂岩自身复杂的结构与性质，解决现有的常规室内试验存在着成本高、时间空间占用大、不可控因素较多，并且难以得到岩心微观层面的现象、机理以及无法解释相关的渗流特性等问题。
本发明提供的基于CT扫描和多相物理场耦合的砂岩数字岩心多孔介质渗流特征的模拟方法，包括如下步骤：
S1、对岩心进行CT扫描得到CT图像；
S2、采用Avizo软件对所述CT图像进行后处理，构建三维数字化岩心模型；
S3、提取渗流孔隙结构侯进行网格精细划分，得到重构后的所述三维数字化岩心模型；
S4、将重构后的所述三维数字化岩心模型进行渗流特征模拟研究。
本发明模拟方法中，步骤S2中，所述后处理包括依次进行的图像优化、调整对比度、阈值设置、二值化分割和分割操作。
本发明模拟方法中，所述图像优化的方式为灰度处理，目的是使原始图像更加清晰、明亮，孔隙和骨架细节部分更加突出。整体效果更好，以便于后续进行滤波、阈值分割处理等，使数字岩心更加接近真实模型。
本发明模拟方法中，所述二值化分割的方式为阈值分割处理，处理后的图像以蓝色显示灰度值为1的孔隙，灰色显示灰度值为0的骨架。
本发明模拟方法中，步骤S3之前，所述方法还包括对所述三维数字化岩心模型进行孔隙分析的步骤：删除非连通的孔隙，选择标记孤立的孔隙与连通的孔隙，并将模型进行球棍化处理，以更直观地观察微观孔隙与喉道，建立基于速度场压力场流线模型。
本发明模拟方法中，步骤S3中，为了保证数字岩心能够准确地表征岩样特征，考虑到计算机储存和运算能力，提取部分体积，选择合理尺寸的网格精度，既能有效准确地表征数字岩心特征，又能适应现有的计算机硬件水平。
步骤S4中，在Comsol软件中进行所述渗流特征模拟研究，基于所述三维数字化岩心模型，在Comsol软件生成一个没有内部几何误差的自由四面体3D网格模型。确定边界条件、初始条件等参数其中包括初始状态下两相流体的位置、壁面处的润湿性设置、入口端和出口端的位置。基于达西方程与N-S方程相关理论，对流体在多孔介质中微观流动特征进行求解；再根据流体性质拟合建立粘度响应含水率方程并绘制曲线图，得到渗流过程中粘度随含水饱和度变化情况。
其中，Navier-Stokes方程为：
连续模型基于连续性假设，流体遵循质量守恒方程、动量守恒方程，多孔介质中的纳米流体流动仍然满足N-S方程，其描述方程为Navier-Stokes方程：
式中，ρ为流体宏观密度；U为流体速度；μ为流体动力粘度；p是压力
其中，粘度响应含水率方程为；
μoeff＝aSw 4+bSw 3+cSw 2+dSw+μo
式中，Sw为含水率；a,b,c,d为流体系数；μ0为流体初始粘度；μoeff为流体响应粘度
其中，达西公式为；
其中，v为表观流速；k为渗透率；μ为动力粘度；为压力梯度
本明模拟方法中，步骤S4中，基于所述渗流特征模拟研究，得到岩心的渗透率和孔隙度。
本发明模拟方法得到的模拟结果接近流体在砂岩地层中渗流情况，本发明建立了流体在多孔介质中其粘度与含水饱和度响应方程，弥补了在物理实验中探索砂岩多尺度复杂性质的不足，为流体在砂岩孔隙中渗流机理提供了参考，在砂岩油(气)藏流体力学技术中具有很好的应用价值。
附图说明
图1所示为三维数字岩心灰度均衡化图。
图2所示为三维数字岩心二值化图。
图3所示为三维孔隙结构空间图。
图4所示为部分四面体3D网格图。
图5所示为速度、压力场流线模型图。
图6所示为粘度与含水率响应曲线图。
图7所示为本发明模拟方法的流程图。
具体实施方式
下述实施例中所使用的实验方法如无特殊说明，均为常规方法。
下述实施例中所用的材料、试剂等，如无特殊说明，均可从商业途径得到。
下述实施例所用低渗砂岩相关参数如下：长度2.92cm，直径2.5cm，干重31.8048g，湿重33.8079g，孔隙体积2.0031mL，孔隙度13.9％，气测渗透率9mD，水测渗透率3.2mD。
如图7所示，为本发明模拟方法的流程图：
通过Avizo软件对得到CT图像进行后处理，包括灰度处理、阈值分割等。考虑到CT灰度均衡化处理目的是使原始图像更加清晰、明亮，孔隙和骨架细节部分更加突出。整体效果更好，以便于后续进行滤波、阈值分割处理等，使数字岩心更加接近真实模型。
因其在CT扫描中放射源的固定方向，导致岩样在多方向上其边缘区域没有信息，对灰度均衡化处理的图像进行图像类型转换采集并将其图像外部信息修正为0，采集目的是使岩样图像更为平滑精确。采用阈值分割法进行二值化分割分离岩样中灰度图像的孔隙空间与岩石骨架，阈值分割后的二值图像中可简化为0为孔隙，1为骨架。将阈值分割后的二维均衡化灰度图像累积叠加，从而重构出岩石内部真实的孔隙结构模型，即数字岩心。我们对模型进行了孔隙分析，删除了非连通的孔隙，选择标记孤立的孔隙与连通的孔隙，并将模型进行球棍化处理，更直观地观察微观孔隙与喉道，建立基于速度场压力场流线模型。
在Avizo软件中进行网格划分，保证数字岩心能够准确地表征岩样特征，考虑到计算机储存和运算能力，提取部分体积，选择合理尺寸的网格精度，既能有效准确地表征数字岩心特征，又能适应现有的计算机硬件水平。
将Avizo软件所重构的三维数字化岩心导入至Comsol多物理场耦合软件中。通过Comsol软件研究数字岩心的渗流特征，对模型结果进行验证求解，构建流体在多孔介质中自组装粘度变化方程。
基于三维数字化岩心模型，在Comsol软件生成一个没有内部几何误差的自由四面体3D网格模型。确定边界条件、初始条件等参数其中包括初始状态下两相流体的位置、壁面处的润湿性设置、入口端和出口端的位置。基于达西方程与N-S方程相关理论，对流体在多孔介质中微观流动特征进行求解。再根据流体性质拟合建立粘度响应含水率方程并绘制曲线图，得到渗流过程中粘度随含水饱和度变化情况。
N-S方程为：
连续模型基于连续性假设，流体遵循质量守恒方程、动量守恒方程，多孔介质中的纳米流体流动仍然满足N-S方程，其描述方程为Navier-Stokes方程：
式中，ρ为流体宏观密度；U为流体速度；μ为流体动力粘度；p是压力其中，粘度响应含水率方程为；
μoeff＝aSw 4+bSw 3+cSw 2+dSw+μo
式中，Sw为含水率；a,b,c,d为流体系数；μ0为流体初始粘度；μoeff为流体响应粘度。
其中，达西公式为；
其中，v为表观流速；k为渗透率；μ为动力粘度；为压力梯度
图1给出的是本发明CT初步扫描后的灰度均衡化展示图。如图所示，通过灰度均衡化进行图像优化、细节突出，随后利用Convert image type、Corp editor以及Volume edit功能对灰度图像进行剪裁提取有效信息灰度图像。鉴于Volume edit功能，需要额外输出.mask文件以便后续孔隙分析使用。
图2给出的是本发明利用阈值分割功能处理后展示图。如图所示，通过阈值分割进行二值化处理，调节阈值强度范围，使灰度图像中孔隙结构与岩石骨架区分。
在Avizo软件中通过Volume fraction功能选择.mask文件作为计算分母，二值化处理后图像数据为分子进行孔隙度计算。
图3给出的是通过孔隙结构分析的结果展示图。如图所示，通过Axisconnectivity功能分析阈值处理后的图像，选择适中的相邻关系强度范围进行Z轴上的孔隙连通性分析。随后利用And not image功能排除连通孔隙进行计算得到孤立孔隙并标记。对连通孔隙模型采用Extract subvolume功能提取部分体积进行孔隙分布区分以及孔喉球棍模型建立：通过Separate objects功能，选择较小的参数因子2、输出类型为Connectedobject。随后利用Generate pore network model生成孔隙网格模型，调整孔隙和喉道参数。通过二值图像构建三维数字岩心，经过计算，分割后的岩样孔隙度为14.9％，真实岩心孔隙度为13.9％，误差在1.0％，计算孔隙度略大于真实孔隙度，原因可能是在阈值分割图像处理过程中，主观地选择阈值强度范围，导致小部分骨架被选择视为孔隙。
图4和图5给出的是通过网格划分以及速度压力场流线模型的展示图。利用绝对渗透率测试功能建立速度场、压力场流线模型。在Avizo中利用Meshing模块对连通的孔隙模型进行网格划分，构建部分自由四面体3D孔隙网格模型作为导入Comsol软件的数据。
利用Comsol软件，设置全局参数：导入材料并设置两相流体密度、粘度。
在孔隙网格模型中设置入口和出口条件、边界、润湿壁条件：选择网格模型四周设为边界、两侧分别设为入口和出口、润湿壁选择所有的边界、润湿接触角导入参数(以Pi/3为基准)设置为水湿。
设置瞬态求解参数，进行渗流模拟，分析流体在多孔介质空间中的流动，定义变量与积分，孔隙度por：intop1(1)/V_tot、渗透率perm：u_out*mu_f/dPdL其中u_out为表观出口速度，计算孔隙网格的孔隙度与渗透率。
图6给出的是本发明通过渗流模拟分析的粘度变化结果展示图。如图所示，根据模拟结果，设置流体粘度的函数并拟合构建粘度-含水饱和度响应方程，求解并绘制响应关系曲线图。可以看出，基于粘度-含水饱和度响应方程，在含水率低于70％时，流体尺寸增大，粘度升高，粘度与含水率正相关。在含水率高于70％时，自由水相增多，油水界面膜强度减小，粘度与含水率负相关。
Claims (8)
1.一种基于CT扫描和多相物理场耦合的砂岩数字岩心多孔介质渗流特征的模拟方法，包括如下步骤：
S1、对岩心进行CT扫描得到CT图像；
S2、采用Avizo软件对所述CT图像进行后处理，构建三维数字化岩心模型；
S3、提取渗流孔隙结构侯进行网格精细划分，得到重构后的所述三维数字化岩心模型；
S4、将重构后的所述三维数字化岩心模型进行渗流特征模拟研究。
2.根据权利要求1所述的模拟方法，其特征在于：步骤S2中，所述后处理包括依次进行的图像优化、调整对比度、阈值设置、二值化分割和分割操作。
3.根据权利要求2所述的模拟方法，其特征在于：所述图像优化的方式为灰度处理，以使图像效果增强，加强孔隙和骨架的图像对比程度。
4.根据权利要求2所述的模拟方法，其特征在于：所述二值化分割的方式为阈值分割处理，处理后的图像以蓝色显示灰度值为1的孔隙，灰色显示灰度值为0的骨架。
5.根据权利要求1-4中任一项所述的模拟方法，其特征在于：步骤S3之前，所述方法还包括对所述三维数字化岩心模型进行孔隙分析的步骤：删除非连通的孔隙，选择标记孤立的孔隙与连通的孔隙，并将模型进行球棍化处理，以更直观地观察微观孔隙与喉道，建立基于速度场压力场流线模型。
6.根据权利要求1-5中任一项所述的模拟方法，其特征在于：步骤S4中，在Comsol软件中进行所述渗流特征模拟研究。
7.根据权利要求1-6中任一项所述的模拟方法，其特征在于：步骤S4中，基于所述渗流特征模拟研究，得到岩心的渗透率和孔隙度。
8.根据权利要求1-6中任一项所述的模拟方法，其特征在于：步骤S4中，基于所述渗流特征模拟研究，建立流体在多孔介质中其粘度与含水饱和度的响应方程。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410036071.0A CN117892654A (zh) | 2024-01-10 | 2024-01-10 | 基于ct扫描和多相物理场耦合的砂岩数字岩心多孔介质渗流特征的模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410036071.0A CN117892654A (zh) | 2024-01-10 | 2024-01-10 | 基于ct扫描和多相物理场耦合的砂岩数字岩心多孔介质渗流特征的模拟方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117892654A true CN117892654A (zh) | 2024-04-16 |
Family
ID=90651740
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410036071.0A Pending CN117892654A (zh) | 2024-01-10 | 2024-01-10 | 基于ct扫描和多相物理场耦合的砂岩数字岩心多孔介质渗流特征的模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117892654A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118153404A (zh) * | 2024-05-11 | 2024-06-07 | 中国地质大学（北京） | 气体互驱模型的构建方法和装置、计算机设备和存储介质 |
| CN118150441A (zh) * | 2024-05-11 | 2024-06-07 | 中国地质大学（北京） | 一种基于流体流动特征评价岩石微观孔隙结构的方法 |
| CN118410684A (zh) * | 2024-06-28 | 2024-07-30 | 长安大学 | 用于岩土类材料渗透性评估模型的构建方法及应用 |
-
2024
- 2024-01-10 CN CN202410036071.0A patent/CN117892654A/zh active Pending
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118153404A (zh) * | 2024-05-11 | 2024-06-07 | 中国地质大学（北京） | 气体互驱模型的构建方法和装置、计算机设备和存储介质 |
| CN118150441A (zh) * | 2024-05-11 | 2024-06-07 | 中国地质大学（北京） | 一种基于流体流动特征评价岩石微观孔隙结构的方法 |
| CN118150441B (zh) * | 2024-05-11 | 2024-07-30 | 中国地质大学（北京） | 一种基于流体流动特征评价岩石微观孔隙结构的方法 |
| CN118410684A (zh) * | 2024-06-28 | 2024-07-30 | 长安大学 | 用于岩土类材料渗透性评估模型的构建方法及应用 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN117892654A (zh) | 基于ct扫描和多相物理场耦合的砂岩数字岩心多孔介质渗流特征的模拟方法 | |
| Gerke et al. | Improving watershed-based pore-network extraction method using maximum inscribed ball pore-body positioning | |
| Baychev et al. | Reliability of algorithms interpreting topological and geometric properties of porous media for pore network modelling | |
| Garing et al. | Pore-scale capillary pressure analysis using multi-scale X-ray micromotography | |
| Jing et al. | Coal cleat reconstruction using micro-computed tomography imaging | |
| Tomutsa et al. | Analysis of chalk petrophysical properties by means of submicron-scale pore imaging and modeling | |
| Vogel | Morphological determination of pore connectivity as a function of pore size using serial sections | |
| CN108802073B (zh) | 一种基于数字岩心的岩石电性参数获取方法及装置 | |
| CN103698803B (zh) | 一种岩石孔隙结构表征方法及装置 | |
| Rabbani et al. | Evolution of pore-scale morphology of oil shale during pyrolysis: a quantitative analysis | |
| Thompson et al. | Application of a new grain-based reconstruction algorithm to microtomography images for quantitative characterization and flow modeling | |
| CN110146525A (zh) | 一种基于分形理论及ct扫描的煤体孔渗参数预测方法 | |
| Dalton et al. | Methods to measure contact angles in scCO2-brine-sandstone systems | |
| Jiang et al. | Permeability estimation of porous media by using an improved capillary bundle model based on micro-CT derived pore geometries | |
| CN104331579A (zh) | 一种低渗透储层原油边界层的模拟方法 | |
| CN109142875B (zh) | 一种利用数字岩心获取致密砂岩岩石电学特性的方法 | |
| Van Marcke et al. | An improved pore network model for the computation of the saturated permeability of porous rock | |
| Qin et al. | Pore-scale modeling of pore structure properties and wettability effect on permeability of low-rank coal | |
| Wang et al. | Coupled pressure-driven flow and spontaneous imbibition in shale oil reservoirs | |
| Wang et al. | Microscale numerical simulation of non-Darcy flow of coalbed methane | |
| Jackson et al. | Deep learning of multi-resolution X-Ray micro-CT images for multi-scale modelling | |
| CN110609046B (zh) | 一种基于微ct水驱油图像的毛管力预测方法 | |
| Jing et al. | Coupling of pipe network modelling and domain decomposition for flow in mineralised coal cores | |
| Korost et al. | Computation of reservoir properties based on 3D-structure of porous media | |
| CN114705606B (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 |