CN107192725A - 一种确定泡沫油中气泡大小与数量的方法及装置 - Google Patents
一种确定泡沫油中气泡大小与数量的方法及装置 Download PDFInfo
- Publication number
- CN107192725A CN107192725A CN201710242786.1A CN201710242786A CN107192725A CN 107192725 A CN107192725 A CN 107192725A CN 201710242786 A CN201710242786 A CN 201710242786A CN 107192725 A CN107192725 A CN 107192725A
- Authority
- CN
- China
- Prior art keywords
- scan image
- pixel
- oil
- saturation
- mrow
- 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
- 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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/03—Investigating materials by wave or particle radiation by transmission
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/10—Different kinds of radiation or particles
- G01N2223/101—Different kinds of radiation or particles electromagnetic radiation
- G01N2223/1016—X-ray
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/635—Specific applications or type of materials fluids, granulates
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)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本申请实施例提供了一种确定泡沫油中气泡大小及数量的方法及装置，所述方法包括：获取洗油烘干后的第一岩心的第一CT扫描图像；获取第二岩心处于第一压力状态的第二CT扫描图像；获取所述第二岩心处于第二压力状态的第三CT扫描图像；根据所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量。本申请实施例可以更加准确的确定泡沫油中气泡的大小和数量。
Description
技术领域
本申请涉及石油开发技术领域，尤其是涉及一种确定泡沫油中气泡大小与数量的方法及装置。
背景技术
泡沫油是通过溶解气驱开采的、井口油样具有明显泡沫特征的稠油，其形成原因主要是随着油藏内压力降低，原油中释放出的气体无法很快形成连续的流动相，而是趋向于以稳定的微气泡形式在原油中以分散状态存在。这些所谓的“微气泡”，在整个体系中它们不会凝聚成一个连续的气相。产生这些微小气泡的压力即为泡点压力。而当地层压力低于热力学平衡泡点压力以下，原油中的溶解气会以扩散气的形式析出，形成一个自由移动气相，直到更低的压力，即拟泡点压力，此时原油中的微小气泡贯通，形成连续相。
泡沫油中气泡的大小和数量能直接描述原油内析出气体的聚并情况，与生产动态密切相关。目前常用的测定气泡大小和数量的方法主要是在通过玻璃刻蚀模型模拟岩心测定的，该方法无法反应油藏真实情况。
发明内容
本申请实施例的目的在于提供一种确定泡沫油气泡大小与数量的方法及装置，以更加准确的确定泡沫油中气泡的大小和数量。
为达到上述目的，本申请实施例提供了一种确定泡沫油气泡大小与数量的方法，所述方法包括：
获取洗油烘干后的第一岩心的第一CT扫描图像；所述第一CT扫描图像中的每个像素点中包含所述第一岩心对应位置处的CT值；
获取第二岩心处于第一压力状态的第二CT扫描图像；所述第二CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第一压力状态包括所述第二岩心内的压力为第一预设压力；所述第一预设压力大于等于所述泡沫油的泡点压力；所述第二岩心为所述第一岩心注满泡沫油后的岩心；
获取所述第二岩心处于第二压力状态的第三CT扫描图像；所述第三CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第二压力状态包括所述岩心内的压力为第二预设压力；所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像的成像条件相同；
根据所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量。
为达上述目的，本申请实施例还提供了一种确定泡沫油气泡大小与数量的装置，所述装置可以包括：
第一获取模块，用于获取洗油烘干后的第一岩心的第一CT扫描图像；所述第一CT扫描图像中的每个像素点中包含所述第一岩心对应位置处的CT值；
第二获取模块，用于获取第二岩心处于第一压力状态的第二CT扫描图像；所述第二CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第一压力状态包括所述第二岩心内的压力为第一预设压力；所述第一预设压力大于等于所述泡沫油的泡点压力；所述第二岩心为所述第一岩心注满泡沫油后的岩心；
第三获取模块，用于获取所述第二岩心处于第二压力状态的第三CT扫描图像；所述第三CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第二压力状态包括所述岩心内的压力为第二预设压力；所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像的成像参数相同，且成像位置相同；
确定模块，用于根据所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量。
由上述本申请实施例提供的技术方案可知，本申请实施例在测量泡沫油气泡大小及数量时，首先分别获取了烘干后岩心的第一CT扫描图像、100％包含泡沫油的岩心在压力为第一预设值时的第二CT扫描图像以及100％包含泡沫油的岩心在压力为第二预设值时的第三CT扫描图像。由于物体的成分不同，CT值就会不同，因此，可以根据这三幅CT扫描图像，确定岩心内的压力为第二预设值时，泡沫油的气泡大小和数量。本申请实施例获取泡沫油气泡大小及数量时，使用了真实岩心模拟真实油藏情况，从而提高了测量结果的准确性。
附图说明
此处所说明的附图用来提供对本申请实施例的进一步理解，构成本申请实施例的一部分，并不构成对本申请实施例的限定。在附图中：
图1为本申请实施例的一种确定泡沫油气泡大小及数量的方法示意图；
图2为本申请实施例的一种图像像素示意图；
图3为本申请实施例的一种图像含气饱和度比较结果示意图；
图4为本申请实施例的一种确定泡沫油气泡大小及数量的装置示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本申请实施例做进一步详细说明。在此，本申请实施例的示意性实施例及其说明用于解释本申请实施例，但并不作为对本申请实施例的限定。
气体从油中析出形成气泡会占据体积，对周围的油有驱动作用，大量气泡的生成产生的共同作用会对这个油藏都具有驱动作用，也就是天然的采油方式。在生成相同量的气体的情况下，气泡越分散，即越是小而多，那么对生产的正面作用越大；反之如果气体聚并严重，气泡少而大，则对生产正面作用较少并有负面作用。所以气泡的生成情况对油田的实际生产评价、预期以及针对性的调整都有重要的意义。
下面结合附图，对本申请实施例的具体实施方式作进一步的详细说明。
参考图1所示，本申请实施例所提供的一种确定泡沫油气泡大小及数量的方法可以包括以下几个步骤。该方法可以用来测量岩心内压力为第二预设压力时，泡沫油中气泡的大小和数量。所述泡沫油可以指通过溶解气驱开采的、井口油样具有明显泡沫特征的稠油。
S101，获取洗油烘干后的第一岩心的第一CT扫描图像；所述第一CT扫描图像中的每个像素点中包含所述第一岩心对应位置处的CT值。
在本实施例中，所述第一岩心可以为砂岩样品、泥岩样品等。洗油烘干的目的是为了将第一岩心中的油和水去除。所述第一CT扫描图像中各个像素点的CT值与第一岩心中对应位置的成分有关。
S102，获取第二岩心处于第一压力状态的第二CT扫描图像；所述第二CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第一压力状态包括所述第二岩心内的压力为第一预设压力；所述第一预设压力大于等于所述泡沫油的泡点压力；所述第二岩心为所述第一岩心注满泡沫油后的岩心。
在本实施例中，向所述第一岩心注入泡沫油的过程可以通过将第一岩心的一端封堵，在另一端以一定注入压力100％饱和泡沫油。需要说明的是，此处第一岩心与第二岩心之间的区别可以仅指同一岩心是否100％饱和泡沫油的区别。为了保证注入第一岩心的泡沫油中没有气泡产生，所述注入压力数值可以大于等于所注泡沫油泡点压力。向所述第一岩心样本注满泡沫油可以由实验操作人员完成。当然，也可以由其他主体完成，例如，计算机等。所述第一预设压力理论来说可以为大于等于po的任意数值。但是，具体实施过程中，为了方便操作，第一预设值只需略大于po即可。当然为了进一步简化操作，所述第一预设压力可以与所述注入压力相等。这里所述的泡沫油的泡点压力po可以是泡沫油本身的泡点压力，即在PVT筒中测量得到的泡点压力。所述第二岩心内的压力数值控制可以通过回压阀实现。
所述第二CT扫描图像可以是第二岩心内压力稳定后，得到的CT扫描图像。第二岩心内的压力是否稳定可以通过时间来判断。具体的，在调节好回压阀压力等于所述第二预设压力后，等待预设时间，然后再获取所述第二CT扫描图像。当然在其他实施例中，也可以通过其他方式来判断第二岩心内压力是否稳定。
S103，获取所述第二岩心处于第二压力状态的第三CT扫描图像；所述第三CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第二压力状态包括所述岩心内的压力为第二预设压力；所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像的成像参数条件相同。
所述第二预设值可以为待研究的压力值，即需要确定泡沫油气泡大小及数量的压力值。所述第三CT扫描图像可以是第二岩心内压力稳定后，得到的CT扫描图像。所述成像条件可以包括像素点大小、图像大小以及成像位置等。
S104，根据所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量。
在本实施例中，S104实施可以包括以下几个子步骤。
(1)根据所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像，确定所述第一CT扫描图像中各个像素点的孔隙度和所述第三CT扫描图像中各个像素点的含油饱和度。
在本申请的一个具体实施例中，所述第一CT扫描图像中各个像素点的孔隙度和所述第三CT扫描图像中各个像素点的含油饱和度可以通过以下步骤计算得到。
1)获取所述泡沫油的CT值以及空气的CT值。
由于不同油的CT值不同。因此，本实施例中，获取泡沫油的CT值可以是第二岩心中注入的那种泡沫油的CT值。不同成分的气体CT数值不同，此处的空气的CT值可以指当前实验条件下的空气的CT值。但由于空气CT的数值大致都在-1000左右，因此，可以默认空气的CT值为-1000。为了确保结果的准确性，本本申请的一个实施例中，测量空气和泡沫油的CT值的扫描条件可以与获取第一CT扫描图像、第二CT扫描图像和第三CT扫描图像相同。
2)根据所述泡沫油的CT值、所述空气的CT值、所述第一CT扫描图像、所述第二CT扫描图像以及以下公式，计算得到所述第一CT扫描图像中各个像素点的孔隙度。
式中，φi表示第一CT扫描图像中第i个像素点中的孔隙度，CToil表示泡沫油的CT值；CTair表示空气的CT值；CTdry,i表示第一CT扫描图像中第i个像素点的CT值；CTwet,i表示第二CT扫描图像中第i个像素点的CT值。
根据所述第一CT扫描图像、所述第二CT扫描图像、所述第三CT扫描图像以及以下公式，计算得到所述第三CT扫描图像中各个像素点的含油饱和度。
式中，CTi表示第三CT扫描图像中第i个像素点的CT值；So,i表示第三CT扫描图像中第i个像素点的含油饱和度；CTdry,i表示第一CT扫描图像中第i个像素点的CT值；CTwet,i表示第二CT扫描图像中第i个像素点的CT值。
需要注意的是，在本实施例中，计算所述第一CT扫描图像中各个像素点的孔隙度和所述第三CT扫描图像中各个像素点的含油饱和度之间无先后顺序，只是为了描述方便，写成如上的顺序。
(2)根据所述第一CT扫描图像中各个像素点的孔隙度和所述第三CT扫描图像中各个像素点的含油饱和度，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量。
在本申请的一个具体实施例中，可以通过以下步骤得到泡沫油的气泡大小和数量。
1)获取所述第三CT扫描图像中的各个像素点的体积。
所述述第三CT扫描图像中的各个像素点的体积可以根据成像参数确定。例如，在CT扫描中，像素体积可以通过仪器的参数获知。
2)获取所述第二岩心处于第三状态的第四CT扫描图像；所述第三状态包括所述第二岩心内的压力为泡点压力；所述第四CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第四CT扫描图像和所述第一CT扫描图像的成像条件相同。
在本实施例中，所述泡点压力可以是泡沫油在多孔介质中的泡点压力。所述第四CT扫描图像可以在第二岩心内的压力稳定后获得。
3)根据所述第一CT扫描图像、所述第二CT扫描图像和所述第四CT扫描图像，确定所述第四CT扫描图像中各个像素点的含油饱和度，并得到所述第四CT扫描图像中各像素点的含油饱和度的第一频率分布曲线。
确定所述第四CT扫描图像中的各个像素点的含油饱和度，与前面确定第三CT扫描图像中的各个像素点的含油饱和度步骤相同，具体实施可以参考前文确定第三CT扫描图像中的各个像素点的含油饱和度步骤，此处不再赘述。再得到所述第四CT扫描图像中的各个像素点的含油饱和度之后，就可以得到第四CT扫描图像中含油饱和度频率分布曲线，即第一频率分布曲线。所述第一频率分布曲线可以是一条描述不同含油饱和度数值对应的像素个数的曲线。具体的，可以将含油饱和度数值作为横坐标，分布频率作为纵坐标生成第一频率分布曲线。其中，所述分布频率可以等于含油饱和度为某一个值的像素点个数占第四CT扫描图像总像素点的百分比。
4)根据所述第一频率分布曲线，得到最小产气饱和度。
在本实施例中，当第二岩心中的泡沫油处于泡点压力时，泡沫油中气体析出并以非连续微小气泡形态存在，此时所对应的含气饱和度可以被称之为最小产气饱和度。具体的，所述含气饱和度可以等于所述第一频率分布曲线中达到峰值时，所对应的含气饱和度。
在本申请的一个具体实施例中，可以按照以下步骤最小产气饱和度。
获取所述第一频率分布曲线的峰值所对应的第一峰值含油饱和度。
将所述第一峰值含油饱和度，代入以下公式得到最小产气饱和度。
SMBGS＝1-S1 (3)
式中，SMBGS表示最小产气饱和度；S1表示第一峰值含油饱和度。
5)获取所述第二岩心处于第四状态的第五CT扫描图像；所述第四状态包括所述第二岩心内的压力为拟泡点压力；所述第五CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第五CT扫描图像和所述第一CT扫描图像的条件相同。
在本实施例中，所述拟泡点压力可以是泡沫油在多孔介质中的拟泡点压力。所述第五CT扫描图像可以在第二岩心样本内的压力稳定后获得。
6)根据所述第一CT扫描图像、所述第二CT扫描图像和所述第五CT扫描图像，确定所述第五CT扫描图像中各个像素点的含油饱和度，并得到所述第五CT扫描图像中各像素点的含油饱和度的第二频率分布曲线。
所述第二频率分布曲线与所述第一频率分布曲线生成方法类似，此处不在赘述。
7)根据所述第二频率分布曲线，确定临界含气饱和度。
当第二岩心中的泡沫油处于拟泡点压力下，泡沫油中微小气泡开始膨胀聚集贯通，并开始移动，此时的含气饱和度可以被称为临界含气饱和度。具体的，所述临界含气饱和度可以等于所述第二频率分布曲线达到峰值时，所对应的含气饱和度。
在本申请的一个具体实施例中，可以根据以下步骤得到临界含气饱和度。
所述根据所述第二频率分布曲线，确定临界含气饱和度，具体包括：
获取所述第二频率分布曲线的峰值所对应的第二峰值含油饱和度。
将所述第二峰值含油饱和度，代入以下公式得到临界含气饱和度。
SCGS＝1-S2 (4)
式中，SCGS表示临界含气饱和度；S2表示第二峰值含油饱和度。
8)根据所述第三CT扫描图像中各个像素点的含油饱和度、所述第一CT扫描图像中各个像素点的孔隙度、以及所述各个像素点的体积、最小气体饱和度以及临界含气饱和度，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量
在本实施例中，在得到泡沫油的最小气体饱和度SMBG以及临界含气饱和度SCG之后，可以按照以下步骤，得到所述第二岩心内的压力为所述第二预设压力时，所述泡沫油的气泡大小和数量。
根据所述第三CT扫描图像中各个像素点的含油饱和度，确定所述第三CT扫描图像中各个像素点的含气饱和度。
当所述第三CT扫描图像中第i个像素的含气饱和度小于所述最小产气饱和度时SMBG，则第i个像素内无气泡产生。
当所述第三CT扫描图像中第i个像素的含气饱和度大于等于所述最小产气饱和度SMBG且小于所述临界含气饱和度SCG，则第i个像素内有一个孤立气泡。
当所述第三CT扫描图像中第i个像素的含气饱和度大于等于所述临界含气饱和度SMBG，则第i个像素含有一个与周围其他像素连通的贯通气泡。
统计孤立气泡和贯通气泡的数量，得到所述泡沫油中气泡数量；
再根据以下公式计算所述泡沫油中气泡的大小；
V孤立＝ΔV像素点·φi·Sg,i (5)
V贯通＝ΔV像素点·∑φi·Sg,i (6)
式中，V孤立表示孤立气泡的大小，ΔV像素点表示第三CT扫描图像中各个像素点的体积，φi表示第一CT扫描图像中第i个像素中的孔隙度；Sg,i表示第三CT扫描图像中第i个像素的含气饱和度；V贯通表示贯通气泡的大小。
由图1所示的实施例可知，本申请实施例在测量泡沫油气泡大小及数量时，首先分别获取了烘干后岩心的第一CT扫描图像、100％包含泡沫油的岩心在压力为第一预设值时的第二CT扫描图像以及100％包含泡沫油的岩心在压力为第二预设值时的第三CT扫描图像。由于物体的成分不同，CT值就会不同，因此，可以根据这三幅CT扫描图像，确定岩心内的压力为第二预设值时，泡沫油的气泡大小和数量。本申请实施例获取泡沫油气泡大小及数量时，使用了真实岩心模拟真实油藏情况，从而提高了测量结果的准确性。
在本申请的一个实施例中，统计气泡数量时，若所述第三CT扫描图像的第i个像素中存在1个孤立气泡时，则气泡数量加1。若第三CT扫描图像的第i个像素中存在贯通气泡，则其与相邻像素的贯通气泡合起来，记为1个气泡。计算贯通气泡体积时，若第i个像素中存在贯通气泡，则将第i个像素周围所有相邻的贯通气泡按公式(6)计算得到，贯通气泡的体积。具体的，如图2所示，图中共有9个像素，分别为像素1至像素9。其中，像素1、像素4、像素7和像素9的含气饱和度是大于等于所述最小产气饱和度SMBG且小于所述临界含气饱和度SCG，则像素1、像素4、像素7和像素9中各有一个孤立气泡，气泡数记4。像素2、像素3、像素5、像素6以及像素8的含气饱和度大于临界含气饱和度SCG，且像素2、像素3、像素5、像素6以及像素8两两有公共边，所以像素2、像素3、像素5、像素6以及像素8合起来是一个贯通气泡，气泡数记1，体积根据公式(6)计算。
在本申请的另一个实施例，根据第三CT扫描图像中各个像素的含气饱和度与最小产气饱和度SMBG和临界含气饱和度SCG之间的关系，用0、1和2标记所述第三CT扫描图像中每一个像素，如图3所示。其中，0表示该像素点的含气饱和度小于SMBG；1表示该像素点的含气饱和度SMBG且小于SCG；2表示该像素点的含气饱和度大于等于SCG。统计图3中气泡个数可知，图中共有4个气泡，且这4个气泡中有三个孤立气泡，1个贯通气泡。然后根据公式(5)和公式(6)就可以得到图3中各个气泡的体积大小。
本实施例中，使用0、1和2标注了图中各个像素的含气饱和度，更加直观可以得到气泡数量和体积的统计。
本申请实施例中还提供了一种确定泡沫油气泡大小及数量的装置，如下面的实施例所述。由于该装置解决问题的原理与一种确定泡沫油气泡大小及数量的方法相似，因此该装置的实施可以参见一种确定泡沫油气泡大小及数量的方法的实施，重复之处不再赘述。
如图4所示，本申请实施例所提的一种确定泡沫油气泡大小及数量的装置可以包括以下几个模块。
第一获取模块401，用于获取洗油烘干后的第一岩心的第一CT扫描图像；所述第一CT扫描图像中的每个像素点中包含所述第一岩心对应位置处的CT值。
第二获取模块402，用于获取第二岩心处于第一压力状态的第二CT扫描图像；所述第二CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第一压力状态包括所述第二岩心内的压力为第一预设压力；所述第一预设压力大于等于所述泡沫油的泡点压力；所述第二岩心为所述第一岩心注满泡沫油后的岩心。
第三获取模块403，用于获取所述第二岩心处于第二压力状态的第三CT扫描图像；所述第三CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第二压力状态包括所述岩心内的压力为第二预设压力；所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像的成像参数相同，且成像位置相同。
确定模块404，用于根据所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量。
由上述装置的实施例可知，本申请实施例在测量泡沫油气泡大小及数量时，首先分别获取了烘干后岩心的第一CT扫描图像、100％包含泡沫油的岩心在压力为第一预设值时的第二CT扫描图像以及100％包含泡沫油的岩心在压力为第二预设值时的第三CT扫描图像。由于物体的成分不同，CT值就会不同，因此，可以根据这三幅CT扫描图像，确定岩心内的压力为第二预设值时，泡沫油的气泡大小和数量。本申请实施例获取泡沫油气泡大小及数量时，使用了真实岩心模拟真实油藏情况，从而提高了测量结果的准确性。
本申请实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中，ASIC可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用户终端中的不同的部件中。
在一个或多个示例性的设计中，本申请实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。
以上所述的具体实施例，对本申请的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本申请实施例的具体实施例而已，并不用于限定本申请的保护范围，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。
Claims (12)
1.一种确定泡沫油中气泡大小及数量的方法，其特征在于，所述方法包括以下步骤：
获取洗油烘干后的第一岩心的第一CT扫描图像；所述第一CT扫描图像中的每个像素点中包含所述第一岩心对应位置处的CT值；
获取第二岩心处于第一压力状态的第二CT扫描图像；所述第二CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第一压力状态包括所述第二岩心内的压力为第一预设压力；所述第一预设压力大于等于所述泡沫油的泡点压力；所述第二岩心为所述第一岩心注满泡沫油后的岩心；
获取所述第二岩心处于第二压力状态的第三CT扫描图像；所述第三CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第二压力状态包括所述岩心内的压力为第二预设压力；所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像的成像条件相同；
根据所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量。
2.如权利要求1所述的方法，其特征在于，所述根据所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量，具体包括：
根据所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像，确定所述第一CT扫描图像中各个像素点的孔隙度和所述第三CT扫描图像中各个像素点的含油饱和度；
根据所述第一CT扫描图像中各个像素点的孔隙度和所述第三CT扫描图像中各个像素点的含油饱和度，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量。
3.如权利要求2所述的方法，其特征在于，所述根据所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像，确定所述第一CT扫描图像中各个像素点的孔隙度和所述第三CT扫描图像中各个像素点的含油饱和度，具体包括：
获取所述泡沫油的CT值以及空气的CT值；
根据所述泡沫油的CT值、所述空气的CT值、所述第一CT扫描图像、所述第二CT扫描图像以及以下公式，计算得到所述第一CT扫描图像中各个像素点的孔隙度；
<mrow>
<msub>
<mi>&phi;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>CT</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>CT</mi>
<mrow>
<mi>d</mi>
<mi>r</mi>
<mi>y</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>CT</mi>
<mrow>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>CT</mi>
<mrow>
<mi>a</mi>
<mi>i</mi>
<mi>r</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
式中，φi表示第一CT扫描图像中第i个像素点中的孔隙度，CToil表示泡沫油的CT值；CTair表示空气的CT值；CTdry,i表示第一CT扫描图像中第i个像素点的CT值；CTwet,i表示第二CT扫描图像中第i个像素点的CT值；
根据所述第一CT扫描图像、所述第二CT扫描图像、所述第三CT扫描图像以及以下公式，计算得到所述第三CT扫描图像中各个像素点的含油饱和度；
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>o</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>CT</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>CT</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>CT</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>CT</mi>
<mrow>
<mi>d</mi>
<mi>r</mi>
<mi>y</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
式中，CTi表示第三CT扫描图像中第i个像素点的CT值；So,i表示第三CT扫描图像中第i个像素点的含油饱和度；CTdry,i表示第一CT扫描图像中第i个像素点的CT值；CTwet,i表示第二CT扫描图像中第i个像素点的CT值。
4.如权利要求3所述的方法，其特征在于，所述根据所述第一CT扫描图像中各个像素点的孔隙度和所述第三CT扫描图像中各个像素点的含油饱和度，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量，具体包括：
获取所述第三CT扫描图像中的各个像素点的体积；
获取所述第二岩心处于第三状态的第四CT扫描图像；所述第三状态包括所述第二岩心内的压力为泡点压力；所述第四CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第四CT扫描图像和所述第一CT扫描图像的成像条件相同；
根据所述第一CT扫描图像、所述第二CT扫描图像和所述第四CT扫描图像，确定所述第四CT扫描图像中各个像素点的含油饱和度，并得到所述第四CT扫描图像中各像素点的含油饱和度的第一频率分布曲线；
根据所述第一频率分布曲线，得到最小产气饱和度；
获取所述第二岩心处于第四状态的第五CT扫描图像；所述第四状态包括所述第二岩心内的压力为拟泡点压力；所述第五CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第五CT扫描图像和所述第一CT扫描图像的成像条件相同；
根据所述第一CT扫描图像、所述第二CT扫描图像和所述第五CT扫描图像，确定所述第五CT扫描图像中各个像素点的含油饱和度，并得到所述第五CT扫描图像中各像素点的含油饱和度的第二频率分布曲线；
根据所述第二频率分布曲线，确定临界含气饱和度；
根据所述第三CT扫描图像中各个像素点的含油饱和度、所述第一CT扫描图像中各个像素点的孔隙度、以及所述各个像素点的体积、最小气体饱和度以及临界含气饱和度，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量。
5.如权利要求4所述的方法，其特征在于，所述根据所述第一频率分布曲线，得到最小产气饱和度，具体包括：
获取所述第一频率分布曲线的峰值所对应的第一峰值含油饱和度；
将所述第一峰值含油饱和度，代入以下公式得到最小产气饱和度；
SMBGS＝1-S1
式中，SMBGS表示最小产气饱和度；S1表示第一峰值含油饱和度；
所述根据所述第二频率分布曲线，确定临界含气饱和度，具体包括：
获取所述第二频率分布曲线的峰值所对应的第二峰值含油饱和度；
将所述第二峰值含油饱和度，代入以下公式得到临界含气饱和度；
SCGS＝1-S2
式中，SCGS表示临界含气饱和度；S2表示第二峰值含油饱和度。
6.如权利要求5所述的方法，其特征在于，所述根据所述第三CT扫描图像中各个像素点的含油饱和度、所述第一CT扫描图像中各个像素点的孔隙度、以及所述各个像素点的体积、最小气体饱和度以及临界含气饱和度，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量，具体包括：
根据所述第三CT扫描图像中各个像素点的含油饱和度，确定所述第三CT扫描图像中各个像素点的含气饱和度；
当所述第三CT扫描图像中第i个像素的含气饱和度小于所述最小产气饱和度时，则第i个像素内无气泡产生；
当所述第三CT扫描图像中第i个像素的含气饱和度大于等于所述最小产气饱和度且小于所述临界含气饱和度，则第i个像素内有一个孤立气泡；
当所述第三CT扫描图像中第i个像素的含气饱和度大于等于所述临界含气饱和度，则第i个像素含有一个与周围其他像素连通的贯通气泡；
统计孤立气泡和贯通气泡的数量，得到所述泡沫油中气泡数量；
再根据以下公式计算所述泡沫油中气泡的大小；
V孤立＝ΔV像素点·φi·Sg,i
V贯通＝ΔV像素点·∑φi·Sg,i
式中，V孤立表示孤立气泡的大小，ΔV像素点表示第三CT扫描图像中各个像素点的体积，φi表示第一CT扫描图像中第i个像素中的孔隙度；Sg,i表示第三CT扫描图像中第i个像素的含气饱和度；V贯通表示贯通气泡的大小。
7.一种确定泡沫油中气泡大小及数量的装置，其特征在于，所述装置包括：
第一获取模块，用于获取洗油烘干后的第一岩心的第一CT扫描图像；所述第一CT扫描图像中的每个像素点中包含所述第一岩心对应位置处的CT值；
第二获取模块，用于获取第二岩心处于第一压力状态的第二CT扫描图像；所述第二CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第一压力状态包括所述第二岩心内的压力为第一预设压力；所述第一预设压力大于等于所述泡沫油的泡点压力；所述第二岩心为所述第一岩心注满泡沫油后的岩心；
第三获取模块，用于获取所述第二岩心处于第二压力状态的第三CT扫描图像；所述第三CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第二压力状态包括所述岩心内的压力为第二预设压力；所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像的成像参数相同，且成像位置相同；
确定模块，用于根据所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量。
8.如权利要求7所述的装置，其特征在于，所述确定模块具体包括：
孔隙度和含油饱和度确定子模块，用于根据所述第一CT扫描图像、所述第二CT扫描图像和所述第三CT扫描图像，确定所述第一CT扫描图像中各个像素点的孔隙度和所述第三CT扫描图像中各个像素点的含油饱和度；
气泡的大小和数量确定子模块，用于根据所述第一CT扫描图像中各个像素点的孔隙度和所述第三CT扫描图像中各个像素点的含油饱和度，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量。
9.如权利要求8所述的装置，其特征在于，所述孔隙度和含油饱和度确定子模块，具体包括：
CT值获取子单元，用于获取所述泡沫油的CT值以及空气的CT值；
孔隙度计算子单元，用于根据所述泡沫油的CT值、所述空气的CT值、所述第一CT扫描图像、所述第二CT扫描图像以及以下公式，计算得到所述第一CT扫描图像中各个像素点的孔隙度；
<mrow>
<msub>
<mi>&phi;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>CT</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>CT</mi>
<mrow>
<mi>d</mi>
<mi>r</mi>
<mi>y</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>CT</mi>
<mrow>
<mi>o</mi>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>CT</mi>
<mrow>
<mi>a</mi>
<mi>i</mi>
<mi>r</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
式中，φi表示第一CT扫描图像中第i个像素点中的孔隙度，CToil表示泡沫油的CT值；CTair表示空气的CT值；CTdry,i表示第一CT扫描图像中第i个像素点的CT值；CTwet,i表示第二CT扫描图像中第i个像素点的CT值；
含油饱和度计算子单元，用于根据所述第一CT扫描图像、所述第二CT扫描图像、所述第三CT扫描图像以及以下公式，计算得到所述第三CT扫描图像中各个像素点的含油饱和度；
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>o</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>CT</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>CT</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>CT</mi>
<mrow>
<mi>w</mi>
<mi>e</mi>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>CT</mi>
<mrow>
<mi>d</mi>
<mi>r</mi>
<mi>y</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
式中，CTi表示第三CT扫描图像中第i个像素点的CT值；So,i表示第三CT扫描图像中第i个像素点的含油饱和度；CTdry,i表示第一CT扫描图像中第i个像素点的CT值；CTwet,i表示第二CT扫描图像中第i个像素点的CT值。
10.如权利要求9所述的装置，其特征在于，所述气泡的大小和数量确定子模块，具体包括：
像素体积获取子单元，用于获取所述第三CT扫描图像中的各个像素点的体积；
第一CT图像获取子单元，用于获取所述第二岩心处于第三状态的第四CT扫描图像；所述第三状态包括所述第二岩心内的压力为泡点压力；所述第四CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第四CT扫描图像和所述第一CT扫描图像的成像条件相同；
第一频率分布曲线确定子单元，用于根据所述第一CT扫描图像、所述第二CT扫描图像和所述第四CT扫描图像，确定所述第四CT扫描图像中各个像素点的含油饱和度，并得到所述第四CT扫描图像中各像素点的含油饱和度的第一频率分布曲线；
最小产气饱和度确定子单元，用于根据所述第一频率分布曲线，得到最小产气饱和度；
第二CT图像获取子单元，用于获取所述第二岩心处于第四状态的第五CT扫描图像；所述第四状态包括所述第二岩心内的压力为拟泡点压力；所述第五CT扫描图像中的每个像素点中包含所述第二岩心对应位置处的CT值；所述第五CT扫描图像和所述第一CT扫描图像的成像条件相同；
第二频率分布曲线确定子单元，用于根据所述第一CT扫描图像、所述第二CT扫描图像和所述第五CT扫描图像，确定所述第五CT扫描图像中各个像素点的含油饱和度，并得到所述第五CT扫描图像中各像素点的含油饱和度的第二频率分布曲线；
临界含气饱和度确定子单元，用于根据所述第二频率分布曲线，确定临界含气饱和度；
气泡大小和数量确定子单元，用于根据所述第三CT扫描图像中各个像素点的含油饱和度、所述第一CT扫描图像中各个像素点的孔隙度、以及所述各个像素点的体积、最小气体饱和度以及临界含气饱和度，确定所述第二岩心内的压力为所述第二预设压力时，所述泡沫油中气泡的大小和数量。
11.如权利要求10所述的装置，其特征在于，所述最小产气饱和度确定子单元具体包括：
获取所述第一频率分布曲线的峰值所对应的第一峰值含油饱和度；
将所述第一峰值含油饱和度，代入以下公式得到最小产气饱和度；
SMBGS＝1-S1
式中，SMBGS表示最小产气饱和度；S1表示第一峰值含油饱和度；
所述临界含气饱和度确定子单元，具体包括：
获取所述第二频率分布曲线的峰值所对应的第二峰值含油饱和度；
将所述第二峰值含油饱和度，代入以下公式得到临界含气饱和度；
SCGS＝1-S2
式中，SCGS表示临界含气饱和度；S2表示第二峰值含油饱和度。
12.如权利要求11所述的装置，其特征在于，所述气泡大小和数量确定子单元具体包括：
根据所述第三CT扫描图像中各个像素点的含油饱和度，确定所述第三CT扫描图像中各个像素点的含气饱和度；
当所述第三CT扫描图像中第i个像素的含气饱和度小于所述最小产气饱和度时，则第i个像素内无气泡产生；
当所述第三CT扫描图像中第i个像素的含气饱和度大于等于所述最小产气饱和度且小于所述临界含气饱和度，则第i个像素内有一个孤立气泡；
当所述第三CT扫描图像中第i个像素的含气饱和度大于等于所述临界含气饱和度，则第i个像素含有一个与周围其他像素连通的贯通气泡；
统计孤立气泡和贯通气泡的数量，得到所述泡沫油中气泡数量；
再根据以下公式计算所述泡沫油中气泡的大小；
V孤立＝ΔV像素点·φi·Sg,i
V贯通＝ΔV像素点·∑φi·Sg,i
式中，V孤立表示孤立气泡的大小，ΔV像素点表示第三CT扫描图像中各个像素点的体积，φi表示第一CT扫描图像中第i个像素中的孔隙度；Sg,i表示第三CT扫描图像中第i个像素的含气饱和度；V贯通表示贯通气泡的大小。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710242786.1A CN107192725B (zh) | 2017-04-14 | 2017-04-14 | 一种确定泡沫油中气泡大小与数量的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710242786.1A CN107192725B (zh) | 2017-04-14 | 2017-04-14 | 一种确定泡沫油中气泡大小与数量的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107192725A true CN107192725A (zh) | 2017-09-22 |
| CN107192725B CN107192725B (zh) | 2019-10-11 |
Family
ID=59871715
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710242786.1A Active CN107192725B (zh) | 2017-04-14 | 2017-04-14 | 一种确定泡沫油中气泡大小与数量的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107192725B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108318430A (zh) * | 2017-12-18 | 2018-07-24 | 中国石油天然气股份有限公司 | 一种泡点压力测试系统及方法 |
| CN108896589A (zh) * | 2018-07-25 | 2018-11-27 | 中国科学院武汉岩土力学研究所 | 一种含气土样中气泡大小分布均匀性的评价方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101446171B1 (ko) * | 2013-08-14 | 2014-10-01 | 한국건설기술연구원 | X선 ct 영상을 이용한 시멘트 기반 재료의 기포간격계수 측정방법 |
| CN104101905A (zh) * | 2013-04-11 | 2014-10-15 | 中国石油天然气集团公司 | 一种基于岩电参数的储层分类方法 |
| CN104374682A (zh) * | 2014-11-12 | 2015-02-25 | 中国石油天然气股份有限公司 | 一种岩心ct扫描分析方法及装置 |
| CN106156439A (zh) * | 2016-07-20 | 2016-11-23 | 中国石油大学(华东) | 一种具有泡沫油现象的稠油油藏溶解气驱数值模拟方法 |
| CN106501286A (zh) * | 2016-10-19 | 2017-03-15 | 大连理工大学 | 一种应用ct测量多孔介质内气液间舍伍德数的装置及方法 |
-
2017
- 2017-04-14 CN CN201710242786.1A patent/CN107192725B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104101905A (zh) * | 2013-04-11 | 2014-10-15 | 中国石油天然气集团公司 | 一种基于岩电参数的储层分类方法 |
| KR101446171B1 (ko) * | 2013-08-14 | 2014-10-01 | 한국건설기술연구원 | X선 ct 영상을 이용한 시멘트 기반 재료의 기포간격계수 측정방법 |
| CN104374682A (zh) * | 2014-11-12 | 2015-02-25 | 中国石油天然气股份有限公司 | 一种岩心ct扫描分析方法及装置 |
| CN106156439A (zh) * | 2016-07-20 | 2016-11-23 | 中国石油大学(华东) | 一种具有泡沫油现象的稠油油藏溶解气驱数值模拟方法 |
| CN106501286A (zh) * | 2016-10-19 | 2017-03-15 | 大连理工大学 | 一种应用ct测量多孔介质内气液间舍伍德数的装置及方法 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108318430A (zh) * | 2017-12-18 | 2018-07-24 | 中国石油天然气股份有限公司 | 一种泡点压力测试系统及方法 |
| CN108896589A (zh) * | 2018-07-25 | 2018-11-27 | 中国科学院武汉岩土力学研究所 | 一种含气土样中气泡大小分布均匀性的评价方法 |
| CN108896589B (zh) * | 2018-07-25 | 2020-06-30 | 中国科学院武汉岩土力学研究所 | 一种含气土样中气泡大小分布均匀性的评价方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107192725B (zh) | 2019-10-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7693677B2 (en) | Method of characterizing the distribution of the absolute permeability of a heterogeneous sample | |
| AU2014306018B2 (en) | Image-based direct numerical simulation of petrophysical properties under simulated stress and strain conditions | |
| CA2903451C (en) | Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data | |
| CN103136102B (zh) | 一种Android平台的流畅度测试方法和装置 | |
| CN104278989B (zh) | 一种获取低孔低渗储层饱和度指数的方法 | |
| CN103674811A (zh) | 一种核磁共振孔隙度测量的校正方法、装置及系统 | |
| US11327000B2 (en) | Detecting saturation levels of a core sample using magnetic fields | |
| CN108153933A (zh) | 一种致密储层生产预测方法及装置 | |
| CN107192725A (zh) | 一种确定泡沫油中气泡大小与数量的方法及装置 | |
| CN104215652A (zh) | 确定油气饱和度的方法和装置 | |
| CN105096385A (zh) | 一种二维地震剖面三维显示方法 | |
| CN102913240B (zh) | 一种储层流体识别方法 | |
| CN107449826B (zh) | 流体分布检测方法及装置 | |
| CN106202813B (zh) | 一种计算地层原油粘度的方法及装置 | |
| CN110596772B (zh) | 一种扩散t2弛豫二维谱的校正方法和装置 | |
| Masalmeh | High oil recoveries from transition zones | |
| CN105683494A (zh) | 伪相生产模拟：在储层流动模拟中经由连续相似的阶跃函数相对渗透率受控模型来评估准多相流生产以便对多个岩石物理学实现排序的信号处理方法 | |
| CN108518220B (zh) | 基于综合流体识别指数识别储层流体的方法与装置 | |
| CN104866706A (zh) | 碳酸盐岩渗透率确定方法及装置 | |
| CN110782187A (zh) | 基于含水率的枯竭油气藏型储气库动态评价方法及系统 | |
| CN106570235B (zh) | 一种致密砂岩岩电参数获取方法及装置 | |
| CN104748906B (zh) | 碳酸盐岩毛细管压力确定方法及装置 | |
| CN109989743A (zh) | 一种确定泥浆滤液侵入深度的方法及系统 | |
| Grove et al. | Translating Perforating Laboratory Results to the Downhole Environment | |
| Shao et al. | Improved RBF-Based NMR Pore-Throat Size, Pore Typing, and Permeability Models for Middle East Carbonates |
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 |