CN110132971B - 一种实验室定量区分接触泥质和分散泥质的方法 - Google Patents
一种实验室定量区分接触泥质和分散泥质的方法 Download PDFInfo
- Publication number
- CN110132971B CN110132971B CN201910493722.8A CN201910493722A CN110132971B CN 110132971 B CN110132971 B CN 110132971B CN 201910493722 A CN201910493722 A CN 201910493722A CN 110132971 B CN110132971 B CN 110132971B
- Authority
- CN
- China
- Prior art keywords
- argillaceous
- contact
- dispersed
- pixel
- substances
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
- G01N2021/8405—Application to two-phase or mixed materials, e.g. gas dissolved in liquids
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (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
本发明提供了一种实验室定量区分接触泥质和分散泥质的方法,涉及石油与天然气勘探技术领域,基于薄片图像分析定性区分接触泥质和分散泥质,包括制作铸体薄片,打开铸体薄片图像,得到总像素并分别拾取颗粒部分像素和孔隙部分像素,得到总的泥质含量;对颗粒接触处的泥质进行染色;对染色部分进行识别并计数,得到接触泥质含量,进而得到分散泥质含量。该方法简单、方便,对技术人员水平、经验依赖性较小,能够准确地对接触泥质和分散泥质进行区分并对二者进行含量检测,进而减少实验误差、规避风险。
Description
技术领域
本发明涉及石油与天然气勘探技术领域,具体涉及一种实验室定量区分接触泥质和分散泥质的方法。
背景技术
在泥质砂岩的地球物理勘探中,与两个或两个以上颗粒接触的泥质为接触泥质,与一个颗粒接触的泥质为分散泥质,分散泥质和接触泥质的区分对于正确预测声速十分重要,若不对二者进行区分,使用声速模型预测声速的误差较高,会对地震和测井带来极大的风险。
目前为止,实验室定量区分接触泥质和分散泥质的方法未见文献报道,仅能够通过薄片鉴定或扫描电镜观察泥质分布和颗粒接触关系,上述方法虽简单易行、成本低廉,但对于技术人员的水平、经验依赖较大,同时,扫描电镜得到的图像虽然可以观察到颗粒表面的起伏和孔洞,但是为灰色显示,也无法定量识别泥质和颗粒。沈建国在论文用密度声波时差和自然伽玛测井曲线计算分散泥质结构泥质含量[J].测井技术,1989(2):13-17.中根据实验得到的分散泥质、结构泥质对声波波速的影响,用能反应此类影响泥质砂岩的纵波波速公式,结合密度、自然伽马测井曲线,给出了泥质砂岩中分散泥质、结构泥质的计算方法。该方法在无死孔隙时计算较好,但在有死孔隙时计算有误差,存在一定的局限性。
中国专利CN104535475B公开了一种碳酸盐岩微观结构的确定方法及装置,该方法包括:获取碳酸盐岩样本的电子计算机X射线断层扫描CT图像,并对CT图像进行二值化处理,经二值化处理的图像中碳酸盐岩样本的孔隙的像素数与碳酸盐岩样本的颗粒的像素数不同;根据孔隙的像素数计算所有孔隙的总面积,根据孔隙的像素数和颗粒的像素数计算所述碳酸盐岩样本的CT图像面积;将孔隙的总面积与碳酸盐岩样本的CT图像面积的比值,确定为孔隙度。该方案基于碳酸盐岩样本的CT图像,可以定量地计算孔隙的面积和碳酸盐岩样本的CT图像面积,进而可以准确、可靠地计算碳酸盐岩样本的孔隙度,有利于对碳酸盐岩储层的微观结构进行定量、可靠的分析。但该方法仅能区分并计算孔隙与碳酸盐岩样本的面积,而不能进一步用于接触泥质和分散泥质的定量区分,同时,CT的结果通过计算得到,数据量大,费用高昂,无法大规模应用于储层评价。
现有技术区分接触泥质和分散泥质存在对技术人员水平、经验的依赖性及误差,因此需要找寻一种能够定量区分接触泥质和分散泥质的方法,使得该方法简单、方便,能够准确地对接触泥质和分散泥质进行区分并能够对二者进行含量检测。
发明内容
本发明针对现有技术存在的问题,提供了一种实验室定量区分接触泥质和分散泥质的方法。该方法简单、方便,对技术人员水平、经验依赖性较小,能够准确地对接触泥质和分散泥质进行区分并对二者进行含量检测,进而减少实验误差、规避风险。
为实现上述目的,本发明采用的技术方案如下:
一种实验室定量区分接触泥质和分散泥质的方法,包括以下步骤:
(1)、制作铸体薄片;
(2)、打开铸体薄片图像,得到总像素并分别拾取颗粒部分像素和孔隙部分像素,得到总的泥质含量;
(3)、打开铸体薄片图像,对颗粒接触处的泥质进行染色;
(4)、对染色部分进行识别并计数,得到接触泥质像素,进而得到接触泥质含量与分散泥质含量。
进一步地,步骤(2)中所述铸体薄片图像使用图像处理软件打开,所述图像处理软件包括Image-Pro Plus、Adobe Photoshop或CoreDraw;其中,铸体薄片的制得参照SY/T6103-2004《岩石孔隙结构特征的测定、图像分析法》。
进一步地,步骤(2)中所述总像素、颗粒部分像素和孔隙部分像素使用像素计数工具得到。
进一步地,步骤(2)中所述总的泥质含量由总像素减去颗粒部分像素与孔隙部分像素得到。
进一步地,步骤(3)中所述铸体薄片图像使用画图软件打开,所述染色使用线条工具进行。
进一步地,步骤(4)中所述识别并计数使用Image-Pro-Plus软件进行;所述接触泥质含量为接触泥质像素与总像素之比;所述分散泥质含量由总的泥质含量减去接触泥质含量得到。
本发明所取得的技术效果是:
1.本发明的方法能够定量区分接触泥质和分散泥质,减少实验误差,进而规避风险。
2.本发明中的方法简单、方便,对技术人员水平、经验依赖性较小,能够大规模地应用于储层评价等。
附图说明
图1为岩心铸体薄片图片;
图2为不区分接触泥质和分散泥质的CCT预测声速与实验测量声速交会图;
图3为定量区分接触泥质和分散泥质的CCT预测声速与实验测量声速交会图;
图4为本发明与原始方法纵波速度预测误差与孔隙度交会图;
其中,1-颗粒,2-孔隙,3-接触泥质。
具体实施方式
一种实验室定量区分接触泥质和分散泥质的方法,包括以下步骤:
(1)、参照SY/T6103-2004《岩石孔隙结构特征的测定、图像分析法》制作岩心铸体薄片;
(2)、使用Image-Pro-Plus图像处理软件打开铸体薄片图像(如图1所示),使用像素计数工具得到总像素(600*800)并分别拾取颗粒1部分像素和孔隙2部分像素,由总像素减去颗粒1部分像素与孔隙2部分像素得到总的泥质含量;
(3)、使用windows自带的画图软件打开铸体薄片图像,并使用线条工具对颗粒1接触处的泥质进行染色;
(4)、使用Image-Pro-Plus图像处理软件对染色部分进行识别并计数,得到接触泥质3像素,进而得到接触泥质3含量,该值为接触泥质3像素与总像素之比,再由总的泥质含量减去接触泥质3含量得到分散泥质含量。
随机制取10组铸体薄片进行测量,将通过上述方法得到的接触泥质和分散泥质含量作为输入参数应用于CCT声速预测模型,计算出实验测量纵、横波速度(Vp:纵波速度;Vs:横波速度)如表1所示:
表1接触泥质、分散泥质含量及实验测量纵横波速度
实例 | 分散泥质含量 | 接触泥质含量 | Vp m/s | Vs m/s |
1 | 0.086 | 0.0057 | 1897 | 1219 |
2 | 0.112 | 0.0081 | 1856 | 1374 |
3 | 0.246 | 0.0039 | 2521 | 1384 |
4 | 0.245 | 0.0031 | 2505 | 1449 |
5 | 0.253 | 0.0012 | 2320 | 1525 |
6 | 0.182 | 0.0053 | 2190 | 1385 |
7 | 0.224 | 0.0043 | 2392 | 1519 |
8 | 0.229 | 0.0060 | 2436 | 1425 |
9 | 0.250 | 0.0031 | 2442 | 1491 |
10 | 0.147 | 0.0067 | 2130 | 1482 |
不对接触泥质和分散泥质进行区分(即原始方法)时,CCT预测声速与实验测量声速效果如图2所示,采用本发明的方法对接触泥质和分散泥质进行区分后得到的CCT预测声速与实验测量声速效果如图3所示;本发明与原始方法纵波速度相对误差与孔隙度交会图如图4所示。由图2、图3对比同时结合图4可知,不对接触泥质和分散泥质区分时,预测声速与实验声速误差较大,可达30%,而经过本发明的方法对接触泥质和分散泥质区分后,预测声速与实验声速误差较小,本发明中的方法相比于不区分接触泥质和分散泥质,结果预测误差下降了20%。
最后应当说明的是,以上内容仅用以说明本发明的技术方案,而非对本发明保护范围的限制,本领域的普通技术人员对本发明的技术方案进行的简单修改或者等同替换,均不脱离本发明技术方案的实质和范围。
Claims (7)
1.一种实验室定量区分接触泥质和分散泥质的方法,其特征在于:包括以下步骤:
(1)、制作铸体薄片;
(2)、打开铸体薄片图像,得到总像素并分别拾取颗粒部分像素和孔隙部分像素,得到总的泥质含量;
(3)、打开铸体薄片图像,对颗粒接触处的泥质进行染色;
(4)、对染色部分进行识别并计数,得到接触泥质像素,进而得到接触泥质含量与分散泥质含量;
步骤(2)中所述总的泥质含量由总像素减去颗粒部分像素与孔隙部分像素得到;
步骤(4)中所述接触泥质含量为接触泥质像素与总像素之比,所述分散泥质含量由总的泥质含量减去接触泥质含量得到。
2.根据权利要求1所述的方法,其特征在于:步骤(2)中所述铸体薄片图像使用图像处理软件打开。
3.根据权利要求2所述的方法,其特征在于:所述图像处理软件包括Image-Pro Plus、Adobe Photoshop或CoreDraw。
4.根据权利要求1所述的方法,其特征在于:步骤(2)中所述总像素、颗粒部分像素和孔隙部分像素使用像素计数工具得到。
5.根据权利要求1所述的方法,其特征在于:步骤(3)中所述铸体薄片图像使用画图软件打开。
6.根据权利要求1所述的方法,其特征在于:步骤(3)中所述染色使用线条工具进行。
7.根据权利要求1所述的方法,其特征在于:步骤(4)中所述识别并计数使用Image-Pro-Plus软件进行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910493722.8A CN110132971B (zh) | 2019-06-06 | 2019-06-06 | 一种实验室定量区分接触泥质和分散泥质的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910493722.8A CN110132971B (zh) | 2019-06-06 | 2019-06-06 | 一种实验室定量区分接触泥质和分散泥质的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110132971A CN110132971A (zh) | 2019-08-16 |
CN110132971B true CN110132971B (zh) | 2020-01-07 |
Family
ID=67580554
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910493722.8A Active CN110132971B (zh) | 2019-06-06 | 2019-06-06 | 一种实验室定量区分接触泥质和分散泥质的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110132971B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110132971B (zh) * | 2019-06-06 | 2020-01-07 | 克拉玛依市昂科能源科技有限公司 | 一种实验室定量区分接触泥质和分散泥质的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102175832A (zh) * | 2011-01-10 | 2011-09-07 | 中国石油天然气股份有限公司 | 一种确定典型储层最佳饱和度计算模型的方法 |
CN105181721A (zh) * | 2015-10-08 | 2015-12-23 | 中国石油大学(华东) | 基于岩心扫描的富泥披储层泥质含量计算方法 |
CN105422089A (zh) * | 2015-11-30 | 2016-03-23 | 长江大学 | 一种利用井周超声图像生成孔隙度谱的方法 |
CN110132971A (zh) * | 2019-06-06 | 2019-08-16 | 克拉玛依市昂科能源科技有限公司 | 一种实验室定量区分接触泥质和分散泥质的方法 |
-
2019
- 2019-06-06 CN CN201910493722.8A patent/CN110132971B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102175832A (zh) * | 2011-01-10 | 2011-09-07 | 中国石油天然气股份有限公司 | 一种确定典型储层最佳饱和度计算模型的方法 |
CN105181721A (zh) * | 2015-10-08 | 2015-12-23 | 中国石油大学(华东) | 基于岩心扫描的富泥披储层泥质含量计算方法 |
CN105422089A (zh) * | 2015-11-30 | 2016-03-23 | 长江大学 | 一种利用井周超声图像生成孔隙度谱的方法 |
CN110132971A (zh) * | 2019-06-06 | 2019-08-16 | 克拉玛依市昂科能源科技有限公司 | 一种实验室定量区分接触泥质和分散泥质的方法 |
Non-Patent Citations (1)
Title |
---|
铸体薄片图像的计算机自动分析;翁大丽;《河南石油》;19940228;第8卷(第1期);第46-51页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110132971A (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8583410B2 (en) | Method for obtaining consistent and integrated physical properties of porous media | |
US10422736B2 (en) | Method for determining porosity associated with organic matter in a well or formation | |
CN106198579B (zh) | 一种测量页岩中有机质含量的方法 | |
WO2017114443A1 (zh) | 一种确定储层岩溶发育程度的方法及装置 | |
US20200005013A1 (en) | Identifying geometrical properties of rock structure through digital imaging | |
CN105649615B (zh) | Ct定量、三维可视化测试储层致密油赋存状态的方法 | |
Weaver et al. | Comparison of pathologist-detected and automated computer-assisted image analysis detected sentinel lymph node micrometastases in breast cancer | |
EP3340113A1 (en) | Method and system for analyzing cuttings coming from a wellbore | |
US20150234069A1 (en) | System and Method for Quantifying Vug Porosity | |
CN110873722A (zh) | 一种岩心矿物组分鉴别方法 | |
CN108830140B (zh) | 一种基于电成像测井分形维数的火山岩岩性识别方法 | |
CN110132971B (zh) | 一种实验室定量区分接触泥质和分散泥质的方法 | |
CN110245686A (zh) | 一种计算石英相对百分含量的岩性快速识别方法 | |
CN110320569B (zh) | 一种致密砂岩储层单井裂缝发育强度定量评价方法 | |
CN106777649B (zh) | 一种裂缝型储集层孔隙结构定量评价方法 | |
CN112145165A (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
Ehrlich et al. | Image analysis of pore geometry: relationship to reservoir engineering and modeling | |
CN113027433B (zh) | 计算强非均质性砂砾岩储层渗透率的方法和装置 | |
CN115390155A (zh) | 针对非均质储层的测井解释方法、装置、电子设备及介质 | |
CN110939428A (zh) | 一种致密砂岩油气藏储层裂缝的识别方法 | |
CN116821786A (zh) | 一种低渗储层分类方法、装置及介质 | |
CN112377175B (zh) | 一种优化钻井泥浆快速识别低阻油气层的方法及系统 | |
CN108828687A (zh) | 一种基于电成像面孔率的渗透率计算方法 | |
US20220404330A1 (en) | Method for estimating hydrocarbon saturation of a rock | |
CN115700736A (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 |