CN105866144B - 一种利用ct扫描分辨水合物与水两相的方法 - Google Patents
一种利用ct扫描分辨水合物与水两相的方法 Download PDFInfo
- Publication number
- CN105866144B CN105866144B CN201610362372.8A CN201610362372A CN105866144B CN 105866144 B CN105866144 B CN 105866144B CN 201610362372 A CN201610362372 A CN 201610362372A CN 105866144 B CN105866144 B CN 105866144B
- Authority
- CN
- China
- Prior art keywords
- hydrate
- phase
- intensity profile
- water
- scan
- 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.)
- Active
Links
Classifications
-
- 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
- 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]
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Radiology & Medical Imaging (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
一种利用CT技术实现水合物与水两相分辨的方法,属于化学工程和石油工程技术领域。这种CT技术实现水合物与水两相分辨的方法,首先利用成像技术获得容器中水合物生成前的灰度分布,然后获得容器中水合物生成后的灰度分布,再次获得水合物生成后所在区域从而实现水合物与水的两相分辨,最后利用公式定量获得水合物的分解过程。该方法可以较准确的确定水合物生成后所在区域,从而直观判断出水合物相与水相各自位置,因此分辨效果非常好,且计算出水合物分解过程中的饱和度变化,从而实现定量分析水合物的分解过程。
Description
技术领域
本发明属于化学工程和石油工程技术领域,涉及一种利用CT技术实现水合物与水两相分辨的方法。
背景技术
水合物作为一种具有商业开发前景的新型能源,实现水合物的有效开采是一个重要的研究课题。然而由于水合物与水的密度非常接近,在利用CT扫描的可视化技术获取水合物的生成过程时,水合物与水两相边界模糊,在灰度图像中无法将两项准确区分开。造影剂的出现可以改变各相的成像效果,但是不同种类造影剂的最佳分辨浓度有很大争议。因此利用可视化技术CT扫描水合物的过程中,仍旧无法将水合物与水两相分辨开。
发明内容
为了克服上述现有技术所存在的问题,本发明提供了一种利用CT技术实现水合物与水两相分辨的方法。该方法利用CT技术对水合物相与水相在加入碘化钾的容器中的存在位置进行准确的把握,通过CT扫描成像获取水合物生成前后各相物质的灰度分布情况判断出水合物所在区域,从而将水合物与水两相分辨开。通过成像技术对水合物相的饱和度进行计算,可以定量分析得出水合物的分解过程。
本发明提供了一种利用CT技术实现水合物与水两相分辨的方法,该方法可以较准确的确定水合物生成后所在区域,从而直观判断出水合物相与水相各自位置,分辨效果好,而且计算出水合物分解过程中的饱和度变化,从而实现定量分析水合物的分解过程。
本发明的一种利用CT扫描分辨水合物与水两相的方法,采用CT技术实现水合物与水两相分辨,方法首先利用成像技术获得容器中水合物生成前的灰度分布,然后获得容器中水合物生成后的灰度分布,再次获得水合物生成后所在区域从而实现水合物与水的两相分辨,最后利用公式定量获得水合物的分解过程。
具体步骤如下:
步骤一:获得生成前状态的水合物灰度分布
将玻璃砂作为多孔介质填入容器中,配制碘化钾的质量分数为5.5-6.5%的水合物溶液,将该水溶液置于容器中得到水合物生成前的状态;通过CT扫描获取水合物生成前的灰度分布。
步骤二:获得生成后状态的水合物灰度分布
将有水溶液的容器置于合适的工况下,使得温度、压力满足水合物的生成条件,进行水合物的生成过程。当水合物生成结束时,通过CT扫描获取水合物生成后的灰度分布。
步骤三:获得水合物生成后所在区域
对比步骤一和步骤二的水合物生成前后的灰度分布情况,根据CT扫描成像结果,判断水合物生成后所在区域;水合物相密度越小灰度值越小,从而将容器中水合物和水两相分辨开。
将容器置于合适工况下进行水合物的分解过程,通过CT扫描获取不同分解时刻下容器的灰度分布,用公式(1)至(3)计算相应时刻下水合物的饱和度:
其中,Vh为灰度分布中水合物相的像素个数,Vslice为所有相物质的像素个数,Vpore为孔隙部分的像素个数。通过计算各个时刻水合物相的饱和度,可以定量获得水合物的分解过程。
本发明的利用CT扫描分辨水合物与水两相的方法,首先利用成像技术获得容器中水合物生成前的灰度分布,然后获得容器中水合物生成后的灰度分布,再次获得水合物生成后所在区域从而实现水合物与水的两相分辨,最后利用公式定量获得水合物的分解过程。该方法可以较准确的确定水合物生成后所在区域,从而直观判断出水合物相与水相各自位置,因此分辨效果非常好,且计算出水合物分解过程中的饱和度变化,从而实现定量分析水合物的分解过程。
附图说明
图1(a)是容器中水合物生成前的灰度分布图。
图1(b)是容器中水合物生成后的灰度分布图。
图2(a)是分解时刻为0分钟时水合物的灰度分布图。
图2(b)是分解时刻为22分钟时水合物的灰度分布图。
图2(c)是分解时刻为57分钟时水合物的灰度分布图。
图2(d)是分解时刻为100分钟时水合物的灰度分布图。
图2(e)是分解时刻为139分钟时水合物的灰度分布图。
图2(f)是分解时刻为199分钟时水合物的灰度分布图。
图3是水合物饱和度随分解时间的变化图。
具体实施方式
下面结合技术方案和附图详细叙述本发明的实施例。
实施例1
(1)获得容器中水合物生成前的灰度分布
以四氢呋喃水合物与水溶液的两相分辨为例,首先将玻璃砂作为多孔介质填入试管中,然后按照四氢呋喃和去离子水的质量比为3:2配制四氢呋喃水溶液,加入碘化钾的质量分数为6%,再次将四氢呋喃水溶液倒入试管中,得到四氢呋喃水合物生成前的状态。通过CT扫描获取四氢呋喃水合物生成前的灰度分布,如图1(a)所示;
(2)获得容器中水合物生成后的灰度分布
将盛装有四氢呋喃水溶液的试管置于水浴中,设置温度为-6℃,压力为常压,进行四氢呋喃水合物的生成过程。静置24小时,保证四氢呋喃水合物生成后,通过CT扫描技术获取水合物生成后的灰度分布,如图1(b)所示;
(3)获得水合物生成后所在区域
对比图1(a)和(b)所示的四氢呋喃水合物生成前后的灰度分布情况,根据CT技术成像原理,相密度越小灰度值越小,判断出灰度值最小的黑色区域即为四氢呋喃水合物生成后所在区域,而灰度值居中的灰色区域即为四氢呋喃水溶液所在区域,灰度值最大的白色区域即为玻璃砂,从而将试管中四氢呋喃水合物和水溶液两相分辨开;
(4)获得水合物的分解过程
将试管置于10℃、常压下进行四氢呋喃水合物的分解过程,通过CT扫描获取不同分解时刻下试管内的灰度分布,图2(a)为分解0分钟,图2(b)为分解22分钟,图2(c)为分解57分钟,图2(d)为分解100分钟,图2(e)为分解139分钟,图2(f)为分解199分钟,使用公式(1)至(3)计算相应时刻下四氢呋喃水合物的饱和度:
其中,Vh为灰度分布中水合物相的像素个数,Vslice为所有相物质的像素个数,Vpore为孔隙部分的像素个数。通过计算各个分解时刻四氢呋喃水合物的饱和度,如图3所示,可以定量获得水合物的分解过程。
Claims (3)
1.一种利用CT扫描分辨水合物与水两相的方法,其特征在于以下步骤:
步骤一:获得生成前状态的水合物灰度分布
将玻璃砂作为多孔介质填入容器中,配制碘化钾的质量分数为5.5-6.5%的水合物溶液,将该水合物溶液置于容器中得到水合物生成前的状态,通过CT扫描获取水合物生成前的灰度分布;
步骤二:获得生成后状态的水合物灰度分布
将有水合物溶液的容器置于合适的工况下,使得温度、压力满足水合物的生成条件,进行水合物的生成过程;当水合物生成结束时,通过CT扫描获取水合物生成后的灰度分布;
步骤三:获得水合物生成后所在区域
对比步骤一和步骤二的水合物生成前后的灰度分布情况,根据CT扫描成像结果,判断水合物生成后所在区域;水合物相密度越小灰度值越小,从而将容器中水合物和水两相分辨开;
将容器置于合适工况下进行水合物的分解过程,通过CT扫描获取不同分解时刻下容器的灰度分布,用公式(1)至(3)计算相应时刻下水合物的饱和度:
其中,Vh为灰度分布中水合物相的像素个数,Vslice为所有相物质的像素个数,Vpore为孔隙部分的像素个数,通过计算各个时刻水合物相的饱和度,可以定量获得水合物的分解过程。
2.根据权利要求1所述的一种利用CT扫描分辨水合物与水两相的方法,其特征还在于:配制碘化钾的质量分数为6%的水合物溶液。
3.根据权利要求1或2所述的一种利用CT扫描分辨水合物与水两相的方法,其特征还在于:所述的水合物为四氢呋喃水合物。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610362372.8A CN105866144B (zh) | 2016-05-26 | 2016-05-26 | 一种利用ct扫描分辨水合物与水两相的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610362372.8A CN105866144B (zh) | 2016-05-26 | 2016-05-26 | 一种利用ct扫描分辨水合物与水两相的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105866144A CN105866144A (zh) | 2016-08-17 |
CN105866144B true CN105866144B (zh) | 2019-01-18 |
Family
ID=56642027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610362372.8A Active CN105866144B (zh) | 2016-05-26 | 2016-05-26 | 一种利用ct扫描分辨水合物与水两相的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105866144B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107290358B (zh) * | 2017-06-20 | 2019-11-08 | 大连理工大学 | 一种应用ct测量多孔介质内co2-盐水界面面积变化的方法 |
CN109374436B (zh) * | 2018-09-19 | 2021-05-11 | 大连理工大学 | 一种基于ct图像的水合物沉积物剪切带识别方法 |
CN114034725A (zh) * | 2021-12-13 | 2022-02-11 | 中国石油大学(华东) | 一种四氢呋喃水合物生成过程观测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4916170A (en) * | 1985-12-27 | 1990-04-10 | Nippon Oil Co., Ltd. | Process for making skin marker |
US5920649A (en) * | 1996-07-30 | 1999-07-06 | Sanyo Electric Co., Ltd. | Braille recognition system with ternary-coding |
CN201043955Y (zh) * | 2006-12-28 | 2008-04-02 | 同方威视技术股份有限公司 | 双视角扫描辐射成像装置 |
CN103063687A (zh) * | 2013-01-06 | 2013-04-24 | 中国石油大学(华东) | 多孔介质中剩余油微观分布图像采集实验装置 |
CN104483334A (zh) * | 2014-12-22 | 2015-04-01 | 中国矿业大学 | 一种可视化表征水泥基材料水分空间分布和含量的方法 |
CN105021505A (zh) * | 2015-07-08 | 2015-11-04 | 上海大学 | 一种凝胶在岩石微通道中的滞留性的可视化评价方法 |
-
2016
- 2016-05-26 CN CN201610362372.8A patent/CN105866144B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4916170A (en) * | 1985-12-27 | 1990-04-10 | Nippon Oil Co., Ltd. | Process for making skin marker |
US5920649A (en) * | 1996-07-30 | 1999-07-06 | Sanyo Electric Co., Ltd. | Braille recognition system with ternary-coding |
CN201043955Y (zh) * | 2006-12-28 | 2008-04-02 | 同方威视技术股份有限公司 | 双视角扫描辐射成像装置 |
CN103063687A (zh) * | 2013-01-06 | 2013-04-24 | 中国石油大学(华东) | 多孔介质中剩余油微观分布图像采集实验装置 |
CN104483334A (zh) * | 2014-12-22 | 2015-04-01 | 中国矿业大学 | 一种可视化表征水泥基材料水分空间分布和含量的方法 |
CN105021505A (zh) * | 2015-07-08 | 2015-11-04 | 上海大学 | 一种凝胶在岩石微通道中的滞留性的可视化评价方法 |
Non-Patent Citations (1)
Title |
---|
Direct observations of three dimensional growth of hydrates hosted in porous media;Prasad Kerkar等;《APPLIED PHYSICS LETTERS》;20091231(第95期);第024102-1至024102-3页 |
Also Published As
Publication number | Publication date |
---|---|
CN105866144A (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ramandi et al. | Micro-CT image calibration to improve fracture aperture measurement | |
CN105866144B (zh) | 一种利用ct扫描分辨水合物与水两相的方法 | |
Bultreys et al. | Imaging and image-based fluid transport modeling at the pore scale in geological materials: A practical introduction to the current state-of-the-art | |
CN105372168B (zh) | 一种混凝土孔隙三维结构的构建方法 | |
Pini et al. | Moving across scales: a quantitative assessment of X-ray CT to measure the porosity of rocks | |
Boone et al. | 3D mapping of water in oolithic limestone at atmospheric and vacuum saturation using X-ray micro-CT differential imaging | |
KR101707440B1 (ko) | 가도리니움 화합물을 반응시킨 암석 시료의 전자현미경 이미지를 활용한 공극률 측정방법 | |
CN104751473B (zh) | 碳酸岩多尺度孔洞特征的确定方法及装置 | |
CN106198579A (zh) | 一种测量页岩中有机质含量的方法 | |
Ketcham et al. | Accurate measurement of small features in X‐ray CT data volumes, demonstrated using gold grains | |
CA2896465A1 (en) | Method for producing a three-dimensional characteristic model of a porous material sample for analysis of permeability characteristics | |
Cheng et al. | Visualization study on fluid distribution and end effects in core flow experiments with low-field mri method | |
RU2467315C1 (ru) | Способ определения пространственного распределения и концентрации глины в образце керна | |
AU2016201034B2 (en) | Method of quantitative calculation and assessment of bituman content in ancient carbonate reservoirs | |
Jaques et al. | Review of porosity uncertainty estimation methods in computed tomography dataset | |
Payton et al. | The influence of grain shape and size on the relationship between porosity and permeability in sandstone: A digital approach | |
Jackson et al. | Deep learning of multiresolution X-ray micro-computed-tomography images for multiscale modeling | |
KR20110054092A (ko) | 엑스레이 씨티촬영을 통한 미세 토사의 간극비 측정방법 | |
CN113167713B (zh) | 数字表征岩石渗透率的方法 | |
Armstrong et al. | Estimation of curvature from micro-CT liquid-liquid displacement studies with pore scale resolution | |
CN116993736A (zh) | 一种骨质疏松区域智能识别系统 | |
Zhu et al. | Multi-scale characterization of organic matter pore space in deep marine shale combined with mathematical morphology | |
CN104568921B (zh) | 钢中单一稀土夹杂硫化物含量的测定方法 | |
Tripathi et al. | Exporting petrophysical properties of sandstones from thin section image analysis | |
US20210063298A1 (en) | Method and device for characterizing microbial carbonate pores, and server |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |