CN110988004B - 一种粒间压溶成因石英胶结物含量的评价方法 - Google Patents

一种粒间压溶成因石英胶结物含量的评价方法 Download PDF

Info

Publication number
CN110988004B
CN110988004B CN201911288542.2A CN201911288542A CN110988004B CN 110988004 B CN110988004 B CN 110988004B CN 201911288542 A CN201911288542 A CN 201911288542A CN 110988004 B CN110988004 B CN 110988004B
Authority
CN
China
Prior art keywords
sandstone
content
quartz
sheet
quartz cement
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
Application number
CN201911288542.2A
Other languages
English (en)
Other versions
CN110988004A (zh
Inventor
王文广
林承焰
林建力
张宪国
董春梅
任丽华
黄导武
段冬平
黄鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201911288542.2A priority Critical patent/CN110988004B/zh
Publication of CN110988004A publication Critical patent/CN110988004A/zh
Application granted granted Critical
Publication of CN110988004B publication Critical patent/CN110988004B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating 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/22Investigating 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 measuring secondary emission from the material
    • G01N23/225Investigating 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 measuring secondary emission from the material using electron or ion
    • G01N23/2251Investigating 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 measuring secondary emission from the material using electron or ion using incident electron beams, e.g. scanning electron microscopy [SEM]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications

Abstract

一种粒间压溶成因石英胶结物含量的评价方法,属于砂岩成岩数值模拟方法研究技术领域。收集整理研究目标层位砂岩的铸体薄片、阴极发光薄片和扫描电镜图像资料;砂岩样品分类为部分A和部分B,部分A砂岩样品包括铸体薄片、扫描电镜和阴极发光薄片图像资料,部分B砂岩样品包含铸体薄片和扫描电镜图像资料;在铸体薄片和阴极发光薄片图像上定量统计出总体石英胶结物含量Tqc;在阴极发光薄片图像上定量统计出压溶成因石英胶结物含量PSqc;基于部分A砂岩样品的总体石英胶结物含量和压溶成因石英胶结物含量,评价出压溶成因石英胶结物的比例系数;结合部分B砂岩样品的铸体薄片上的总体石英胶结物含量,评价出所有砂岩样品的压溶成因石英胶结物含量EPSqc。

Description

一种粒间压溶成因石英胶结物含量的评价方法
技术领域
本发明提出一种粒间压溶成因石英胶结物含量的评价方法,属性砂岩成岩数值模拟方法研究技术领域。
背景技术
储层质量研究起始于19世纪80年代,主要依靠经验关系或一代地球化学模型;在这之后,成岩作用动力学模型逐渐被开发出来,如压实作用、石英胶结作用、蒙脱石伊利石化、斜长石钠长石化、纤维状伊利石形成等模型。目前,斜长石钠长石化和纤维状伊利石形成的动力学模型并没有被整合进入储层质量预测模型,当前储层质量预测模型主要包括压实作用、石英胶结作用和纤维状伊利石形成动力学模型。
致密砂岩储层是近年油气勘探的选区目标。成岩数值模拟方法是评价出有利目标区域的重要方法手段,那么如何获得准确的成岩参数是致密砂岩成岩数值模拟研究的关键。在砂岩成岩作用研究中,时常会遇到取芯层段上分析测试资料不匹配问题;这限制了分析测试资料的应用效率,进而影响着成岩参数的评价。
在砂岩成岩作用过程中,压实作用和胶结作用是造成孔隙度减少的2个主要因素。相对于胶结作用堵塞粒间体积,压实作用是直接减少粒间体积。压实作用包括机械压实作用和化学压实作用。在压实作用条件下砂岩孔隙度降低是由机械压实作用和化学压实作用共同控制的。机械压实作用是在上覆水体或沉积层的重荷作用下,沉积物中发生骨架颗粒重新排列、塑性颗粒变形、孔隙水排出和粒间体积减少的过程;胶结作用是指矿物溶解物质在孔隙空间中发生的沉淀作用。压溶作用在石英矿物含量较高的石英砂岩、长石砂岩和岩屑砂岩中较为普遍。尤其是当砂岩中石英矿物颗粒多被伊利石包裹或具有较少颗粒间支撑碳酸盐胶结物时,化学压实作用较为常见。
化学压实作用是三个物理过程的连续,包括颗粒接触处溶蚀作用、颗粒边界上扩散作用和沉淀作用(或胶结作用)。释放的物质(如,硅质)然后再次局部地沉淀在自由碎屑颗粒表面。由Sibley and Blatt(1976)和Houseknecht(1991)使用的化学压实作用评价方法被用来定量分析砂岩样品中粒间压溶作用的效果。这个定量方法是在阴极发光薄片图像上通过石英颗粒-石英颗粒重叠区域来定量地评价石英胶结物含量。沿着颗粒接触位置投影出颗粒边界,并且呈现出压溶作用发生的位置,这些部位被定量统计为压溶石英胶结物含量。尽管这个方法具有主观性,但它仍然是砂岩中可定义粒间压溶作用的一种常用方法。
关于砂岩机械压实作用模型的研究,粒间体积-深度曲线是一种砂岩机械压实作用参数;而关于砂岩化学压实作用模型的研究,一些模型考虑沉淀作用是速率限制过程,但是由于不考虑源项(溶解作用),因此在埋藏深度低于1500m或温度低于70℃时模型必须被抑制;否则模型将过早开始化学压实作用;一些模型认为颗粒边界上扩散作用是速率限制过程,但在埋藏深度大于4500m时必须加以抑制;否则,该模型将预测负孔隙度。当前化学压实作用模型的关键参数评价方案仍然有待深入研究。基于压实作用数值模拟研究需求及存在问题,现有的技术不能满足当前低渗透-致密砂岩储层的有利区目标评价的需求。
发明内容
为了克服现有技术的不足,本发明提供一种粒间压溶成因石英胶结物含量的评价方法。
一种粒间压溶成因石英胶结物含量的评价方法,含有以下步骤;收集整理研究目标层位砂岩的铸体薄片、阴极发光薄片和扫描电镜图像资料;砂岩样品分类为部分A和部分B,
部分A砂岩样品包括铸体薄片、扫描电镜和阴极发光薄片图像资料,
部分B砂岩样品包含铸体薄片和扫描电镜图像资料;
在铸体薄片和阴极发光薄片图像上定量统计出总体石英胶结物含量Tqc;在阴极发光薄片图像上定量统计出压溶成因石英胶结物含量PSqc;基于部分A砂岩样品的总体石英胶结物含量和压溶成因石英胶结物含量,评价出压溶成因石英胶结物的比例系数;结合部分B砂岩样品的铸体薄片上的总体石英胶结物含量,评价出所有砂岩样品的压溶成因石英胶结物含量EPSqc。
还包括以下步骤:
步骤1)、获取研究区图像资料和取芯样品,收集整理研究区内目标层位砂岩的铸体薄片、阴极发光薄片和扫描电镜等图像分析测试资料,并且收集钻井取芯井段上的取芯砂岩样品;
步骤2)、分类统计砂岩样品,这些砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩;将研究区内收集到的砂岩样品分类为部分A和部分B;
部分A砂岩样品包括制作铸体薄片、阴极发光薄片和扫描电镜样品块;
部分B砂岩样品包含制作铸体薄片和扫描电镜样品块;
步骤3)、定量地统计铸体薄片上总体石英胶结物含量,在Axio Vision SE 64图像分析系统下采用图像分析软件拍摄铸体薄片的图像资料,统计出部分A砂岩样品和部分B砂岩样品的总体石英胶结物含量;
步骤4)、定量地统计阴极发光薄片上总体石英胶结物和压溶成因石英胶结物含量;在CL8200 Mk5-2阴极发光系统下采用图像分析软件拍摄部分A砂岩样品的阴极发光薄片图像资料,恢复出石英颗粒初始形态,识别出石英颗粒-石英颗粒的重叠部分,作为石英压溶作用部分;从而确定出压溶成因石英胶结物含量;并且采用图像识别方法在阴极发光薄片图像资料上获取总体石英胶结物含量;
步骤5)、确定铸体薄片上总体石英胶结物含量的校正系数;在阴极发光薄片图像资料上,采用图像识别方法识别出来的石英胶结物含量,具有更高的可靠性;当对比分析铸体薄片上总体石英胶结物含量和阴极发光薄片上总体石英胶结物含量时,采用阴极发光薄片上统计的石英胶结物含量来校正铸体薄片上统计的总体石英胶结物含量;
步骤6)、计算出压溶成因石英胶结物的比例系数,部分A砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩,统计出部分A砂岩样品的总体石英胶结物含量Tqc和压溶成因石英胶结物含量PSqc;评价出部分A砂岩样品的压溶成因石英胶结物的比例系数;
步骤7)、获得所有砂岩样品的压溶成因石英胶结物含量;根据部分A砂岩样品的压溶成因石英胶结物的比例系数,结合部分B砂岩样品铸体薄片上的总体石英胶结物含量Tqc,评价出部分B砂岩样品的压溶成因石英胶结物含量;联合部分A砂岩样品和部分B砂岩样品,获得所有砂岩样品的压溶成因石英胶结物含量EPSqc。
步骤3)中总体石英胶结物含量的定量统计,步骤如下;
总体砂岩样品包括部分A砂岩样品和部分B砂岩样品;统计所有砂岩样品铸体薄片上的总体石英胶结物含量;
对于粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩不同粒径岩性的砂岩样品,首先使用Axio Vision SE 64图像分析系统下图像分析软件对每个铸体薄片拍摄20张显微照片,拍摄照片的目镜是10×倍;然后,使用CorelDraw软件在电脑屏幕上画出和识别出每张铸体薄片图像上的石英胶结物;并且通过Image-Pro Plus软件获取每张铸体薄片图像上石英胶结物的总面积;求取铸体薄片的所有显微照片中石英胶结物含量的平均值,得到这个铸体薄片的总体石英胶结物含量,见公式(1);
Tqc=average(figure(Tqci))i=1,2,3…20…公式(1)
式中:Tqc是在铸体薄片上统计的总体石英胶结物含量;figure(Tqci)是铸体薄片的第i张显微照片上统计的石英胶结物含量。
步骤4)中总体石英胶结物和压溶成因石英胶结物含量的定量统计,步骤如下;
在CL8200 Mk5-2阴极发光系统下采用图像分析软件拍摄部分A砂岩样品的阴极发光薄片图像资料,恢复出石英颗粒初始形态,识别出石英颗粒-石英颗粒的重叠部分,作为石英压溶作用部分;从而定量统计出压溶成因石英胶结物含量;评价步骤如下;
在阴极发光薄片图像资料上识别出石英颗粒-石英颗粒的重叠部分,从而确定出压溶成因石英胶结物含量;
使用阴极发光仪、研究级正立多功能偏光显微镜、摄像头和计算机在加速电压10kV和束电流250uA条件下拍摄阴极发光薄片图像,曝光持续时间10s(放大倍数50×)、30s(放大倍数100×)和60s(放大倍数200×);在阴极发光系统下采用图像分析软件拍摄30张阴极发光薄片图像照片;在阴极发光薄片图像资料上使用一个100(10×10)正方形网格,手动恢复出石英颗粒初始形态,识别出石英颗粒-石英颗粒重叠区域,从而确定出压溶成因石英胶结物含量;然后,对阴极发光薄片的所有显微照片统计出来的压溶成因石英胶结物含量取平均值,这样评价出来的压溶成因石英胶结物含量更具有统计规律,见公式(2);
PSqc=average(figure(PSqci))i=1,2,3…20…公式(2)
式中:PSqc是阴极发光薄片上统计的压溶成因石英胶结物含量;figure(PSqci)是阴极发光薄片上的第i张显微照片上统计的压溶成因石英胶结物含量;
划分出一个100(10×10)正方形网格,在网格空间内恢复出石英颗粒的初始形态,然后识别出石英颗粒-石英颗粒的重叠部分,进而统计出压溶成因石英胶结物含量;
在阴极发光薄片图像资料上,基于石英颗粒和石英胶结物具有不同发光的特征,采用图像分割技术、可视化技术和图像提取方法来刻画出石英胶结物和计算出石英胶结物含量;对多张阴极发光薄片图像上统计的石英胶结物含量取平均值,见公式(3);
CLqc=average(figure(CLqci))i=1,2,3…20…公式(3)
式中:CLqc是阴极发光薄片上统计的总体石英胶结物含量;figure(CLqc)是阴极发光薄片上的第i张显微照片上图像识别方法统计的总体石英胶结物含量。
步骤5)确定铸体薄片上总体石英胶结物含量的校正系数,步骤如下:
在阴极发光薄片图像资料上,根据石英颗粒和石英胶结物具有不同的发光特征,采用图像识别方法识别出来的石英胶结物含量;当对比分析铸体薄片上总体石英胶结物含量和阴极发光薄片上总体石英胶结物含量时,采用阴极发光薄片上统计出来的总体石英胶结物含量CLqc来校正铸体薄片上统计出来的总体石英胶结物含量Tqc;这里给出一个校正系数TCqc,见公式(4);
Figure BDA0002315097910000061
式中:Tqc是在铸体薄片上统计的总体石英胶结物含量;CLqc是在阴极发光薄片上统计的总体石英胶结物含量;TCqc是铸体薄片上统计的总体石英胶结物含量的校正系数。
步骤6)中压溶成因石英胶结物的比例系数,步骤如下;
基于包含粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩的部分A砂岩样品的铸体薄片、扫描电镜和阴极发光薄片图像资料;在铸体薄片上定量地统计出总体石英胶结物含量Tqc;在阴极发光薄片图像上识别出石英颗粒-石英颗粒的重叠区域,从而获得压溶成因石英胶结物含量PSqc;综合考虑铸体薄片图像上统计的总体石英胶结物含量、铸体薄片上总体石英胶结物含量的校正系数和阴极发光薄片图像上统计的压溶成因石英胶结物含量,建立压溶成因石英胶结物的比例系数评价公式,见公式5;
Figure BDA0002315097910000062
式中:Tqc是铸体薄片上统计的总体石英胶结物含量;PSqc是阴极发光薄片上统计的压溶成因石英胶结物含量;TCqc是铸体薄片上统计的总体石英胶结物含量的校正系数;Rpsqc是指压溶成因石英胶结物的比例系数。
步骤7)中获得所有砂岩样品的压溶成因石英胶结物含量,采用压溶成因石英胶结物的比例系数来评价出更多的压溶成因石英胶结物含量;步骤如下;
根据部分A砂岩样品的铸体薄片图像资料的所有石英胶结物含量和阴极发光薄片图像资料的压溶成因石英胶结物含量而计算的压溶成因石英胶结物的比例系数;结合部分B砂岩样品的铸体薄片上定量统计的总体石英胶结物含量Tqc,评价出部分B砂岩样品的压溶成因石英胶结物含量,见公式(6);
EPSqc=Rpsqc×Tqc×TCqc 公式(6)
式中:Tqc是铸体薄片上定量统计的总体石英胶结物含量;Rpsqc是指压溶成因石英胶结物含量的比例系数;TCqc是铸体薄片上统计的总体石英胶结物含量的校正系数;EPSqc是部分B砂岩样品的评价的压溶成因石英胶结物含量;联合部分A砂岩样品和部分B砂岩样品,就获得了所有砂岩样品的压溶成因石英胶结物含量EPSqc。
结合步骤1)、步骤2)、步骤3)、步骤4)、步骤5)、步骤6)和步骤7)来评价出粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩的压溶成因石英胶结物含量;步骤如下;
步骤1)、获取研究区图像资料和取芯样品,收集整理研究区内目标层位砂岩的铸体薄片、阴极发光薄片和扫描电镜等图像分析测试资料,并且收集钻井取芯井段上的取芯砂岩样品;
步骤2)、分类统计砂岩样品,这些砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩;将研究区内收集到的砂岩样品分类为部分A和部分B;部分A砂岩样品包括铸体薄片、阴极发光薄片和扫描电镜图像资料,或者制作铸体薄片、阴极发光薄片和扫描电镜样品块;部分B砂岩样品包含铸体薄片和阴极发光薄片图像资料,或者制作铸体薄片和扫描电镜样品块;
步骤3)、定量地统计铸体薄片上总体石英胶结物含量,在Axio VisionSE64图像分析系统下采用图像分析软件拍摄铸体薄片的图像资料,统计出部分A砂岩样品和部分B砂岩样品的总体石英胶结物含量;
步骤4)、定量地统计阴极发光薄片上总体石英胶结物和压溶成因石英胶结物含量;在CL8200 Mk5-2阴极发光系统下采用图像分析软件拍摄部分A砂岩样品的阴极发光薄片图像资料,恢复出石英颗粒初始形态,识别出石英颗粒-石英颗粒的重叠部分,作为石英压溶作用部分;从而确定出压溶成因石英胶结物含量;并且采用图像识别方法在阴极发光薄片图像资料上获取总体石英胶结物含量;
步骤5)、确定铸体薄片上总体石英胶结物含量的校正系数;在阴极发光薄片图像资料上,采用图像识别方法识别出来的石英胶结物含量,具有更高的可靠性;当对比分析铸体薄片上总体石英胶结物含量和阴极发光薄片上总体石英胶结物含量时,采用阴极发光薄片上统计的石英胶结物含量来校正铸体薄片上统计的总体石英胶结物含量;
步骤6)、计算出压溶成因石英胶结物的比例系数,部分A砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩,统计出部分A砂岩样品的总体石英胶结物含量Tqc和压溶成因石英胶结物含量PSqc;在基础上,评价出部分A砂岩样品的压溶成因石英胶结物的比例系数;
步骤7)、获得所有砂岩样品的压溶成因石英胶结物含量;根据部分A砂岩样品的压溶成因石英胶结物含量的比例系数,结合部分B砂岩样品铸体薄片上的总体石英胶结物含量Tqc,评价出部分B砂岩样品的压溶成因石英胶结物含量;联合部分A砂岩样品和部分B砂岩样品,获得了所有砂岩样品的压溶成因石英胶结物含量EPSqc。
本发明具有操作可行、计算简单方便的特征,具备以下几方面创新性:
1)在铸体薄片图像资料上定量统计总体石英胶结物含量;在Axio Vision SE 64图像分析系统下采用图像分析软件拍摄铸体薄片图像资料,拍摄20张铸体薄片图像照片,使用CorelDRAW软件在电脑上画出和识别出每张铸体薄片图像上总体石英胶结物含量;最后,对这个铸体薄片的所有图像资料定量统计的石英胶结物含量取平均值,从而确定这个铸体薄片的总体石英胶结物含量。
2)在阴极发光薄片图像资料上定量统计压溶成因石英胶结物含量;在阴极发光薄片图像资料上使用一个100(10×10)正方形网格,手动恢复出石英颗粒初始形态,确定石英颗粒-石英颗粒重叠区域,从而确定出压溶成因石英胶结物含量。这可以在阴极发光薄片图像上定量统计的压溶成因石英胶结物含量。
3)根据粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩的总体石英胶结物和粒间压溶成因石英胶结物含量,评价出压溶成因石英胶结物含量的比例系数;结合部分B砂岩样品上所有石英胶结物含量,评价出部分B砂岩样品上压溶成因石英胶结物含量,进而获得所有砂岩样品的压溶成因石英胶结物含量。这样,通过对比压溶成因石英胶结物含量与石英胶结作用模型模拟结果确保化学压实作用的准确性。本发明能很好地满足致密砂岩压实作用数值模拟研究的需求,对成岩学研究具有一定指导意义。
本发明是在致密砂岩的压实作用数值模拟、有利区目标评价的背景下,确定化学压实作用的关键参数,即压溶成因石英胶结物含量。
本发明提出了一种粒间压溶成因石英胶结物含量的评价方法,这主要包括所有砂岩样品的铸体薄片的总体石英胶结物含量和部分砂岩样品的阴极发光薄片的压溶成因石英胶结物含量的定量统计;利用部分A砂岩样品的总体石英胶结物含量和压溶成因石英胶结物含量来评价压溶成因石英胶结物含量的比例系数;然后,利用该比例系数,根据部分B砂岩样品的总体石英胶结物含量反过来评价压溶成因石英胶结物含量。本发明满足了致密砂岩的压实作用数值模拟研究的需求,补充了这方面的空白;可以为重现沉积时期和压实作用条件下砂岩孔隙度演化研究提供关键参数。
附图说明
当结合图考虑时,通过参照下面的详细描述,能够更完整更好地理解本发明以及容易得知其中许多伴随的优点,但此处所说明的图用来提供对本发明的进一步理解,构成本发明的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定,如图其中:
图1研究区目标层位砂岩的粒间压溶成因石英胶结物含量评价的技术路线图;
图2铸体薄片图像资料上总体石英胶结物含量的评价图;
图3阴极发光薄片图像资料上压溶成因石英胶结物含量的评价图;
图4粗砂岩的压溶成因石英胶结物含量的剖面图;
图5中砂岩的压溶成因石英胶结物含量的剖面图;
图6砂质砾岩的压溶成因石英胶结物含量的剖面图;
图7细砂岩的压溶成因石英胶结物含量的剖面图;
图8粉砂岩的压溶成因石英胶结物含量的剖面图;
图9不同粒径岩性砂岩的实测压溶成因石英胶结物含量与石英胶结作用模拟结果对比分析图;
图10研究区目标层位砂岩的温度演化史;
图11机械压实作用下目标层位砂岩孔隙度演化史;
图12化学压实作用下目标层位砂岩的孔隙度响应特征;
图13XX盆地XX凹陷目标层位砂岩快速沉降阶段末期石英胶结物平面分布图;
图14XX盆地XX凹陷目标层位砂岩抬升剥蚀阶段末期石英胶结物平面分布图;
图15XX盆地XX凹陷目标层位砂岩区域稳定沉降阶段末期石英胶结物平面分布图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
下面结合图及具体实施例对本发明的应用原理做进一步描述。
实施例1:如图1、图2、图3、图4、图5、图6、图7、图8、图9、图10、图11、图12、图13、图14及图15所示,一种粒间压溶成因石英胶结物含量的评价方法,收集整体研究区内铸体薄片、阴极发光薄片和扫描电镜等图像分析测试资料;将砂岩样品分类为包括铸体薄片、扫描电镜和阴极发光薄片的部分A砂岩样品和包括铸体薄片和扫描电镜的部分B砂岩样品;在铸体薄片图像资料上定量统计出总体石英胶结物含量,在阴极发光薄片图像资料上定量统计出压溶成因石英胶结物含量;基于部分A砂岩样品的总体石英胶结物和压溶成因石英胶结物含量评价出压溶成因石英胶结物的比例系数;结合铸体薄片上统计获得的总体石英胶结物含量,评价出所有砂岩样品上压溶成因石英胶结物含量。本实施方案提供了一套详细的压溶成因石英胶结物的定量评价方法(图1),包括以下步骤:
(1)获取研究区图像资料和取芯样品:以XX盆地XX凹陷目标层位4口探井为基础,收集整理研究区内目标层位砂岩的铸体薄片、阴极发光薄片和扫描电镜等图像分析测试资料,并且收集钻井取芯井段上的岩芯砂岩样品。
在研究区范围内,统计钻井取芯段数量和钻井取芯段长度;并且收集已经取芯段砂岩样品的铸体薄片、扫描电镜和阴极发光薄片图像资料,见表1:铸体薄片、阴极发光薄片和扫描电镜图片数量的统计表;同时也收集取芯段砂岩样品40块次。
表1
Figure BDA0002315097910000111
(2)分类统计砂岩样品:这些砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩;将研究区内已经收集到的砂岩样品分类为部分A和部分B;部分A砂岩样品包括铸体薄片、扫描电镜和阴极发光薄片图像资料,或者制作铸体薄片、扫描电镜和阴极发光薄片;部分B砂岩样品包括铸体薄片和扫描电镜图像资料,或者制作了铸体薄片和扫描电镜样品块。
(3)定量地统计出总体石英胶结物含量:在Axio Vision SE 64图像分析系统下采用图像分析软件拍摄铸体薄片的图像资料,统计出部分A砂岩样品和部分B砂岩样品的总体石英胶结物含量;
在铸体薄片上统计分析部分A砂岩样品和部分B砂岩样品的总体石英胶结物含量。其具体评价流程为,见图2:
在铸体薄片图像上定量统计出粗砂岩、中砂岩、细砂岩、砂质砾岩和粉砂岩的总体石英胶结物含量。为了满足一定的地质统计学规律,在铸体薄片上进行至少300个点数来检查收集的岩石成分数据的准确性,其标准偏差为6%或更低。
对于不同粒径岩性的砂岩样品,首先使用Axio Vision SE 64图像分析系统下图像分析软件对每个铸体薄片拍摄20张显微照片,拍摄照片时目镜是10×倍。然后,使用CorelDRAW软件在电脑屏幕上画出和识别出每张显微照片上的石英胶结物;并且通过Image-Pro Plus软件来获取每张显微照片上石英胶结物的总面积。最后,求取这个铸体薄片的所有显微照片的石英胶结物含量的平均值,从而得到这个铸体薄片的总体石英胶结物含量。
(4)定量统计压溶成因石英胶结物含量:在CL8200 Mk5-2阴极发光系统下采用图像分析软件拍摄部分A砂岩样品的阴极发光薄片图像资料,恢复出石英颗粒初始形态,定量统计出石英颗粒-石英颗粒的重叠部分,作为石英压溶作用部分;从而定量统计出压溶成因石英胶结物含量;
使用阴极发光仪、研究级正立多功能偏光显微镜、摄像头和计算机在加速电压10kV和束电流250uA条件下拍摄阴极发光薄片图像资料,曝光持续时间10s(放大倍数50×)、30s(放大倍数100×)和60s(放大倍数200×)。在阴极发光系统下采用图像分析软件拍摄30张阴极发光薄片图像照片;在阴极发光薄片图像资料上使用一个100(10×10)正方形网格,手动恢复出石英颗粒初始形态,确定石英颗粒-石英颗粒重叠区域,从而确定出压溶成因石英胶结物含量。
在阴极发光薄片图像资料上,基于石英颗粒和石英胶结物具有不同发光的特征,采用图像分割技术、可视化技术和图像提取方法来刻画出石英胶结物和计算出石英胶结物含量;为了满足统计学规律,对多张阴极发光薄片图像上统计的石英胶结物含量取平均值。
(5)、确定铸体薄片上总体石英胶结物含量的校正系数;采用图像识别方法在阴极发光薄片图像资料上定量统计出来的石英胶结物含量,具有更高的可靠性;因而,当对比分析铸体薄片上总体石英胶结物含量和阴极发光薄片上总体石英胶结物含量时,采用阴极发光薄片上统计出来的石英胶结物含量来校正铸体薄片上统计出来的总体石英胶结物含量;
(6)计算出压溶成因石英胶结物的比例系数,部分A砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩,统计出包含粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩的部分A砂岩样品的总体石英胶结物含量Tqc和压溶成因石英胶结物含量PSqc(见表2);在此基础上,评价出部分A砂岩样品的压溶成因石英胶结物含量的比例系数见表2:部分A砂岩样品的总体石英胶结物含量和压溶成因石英胶结物含量的统计表;
表2
Figure BDA0002315097910000121
Figure BDA0002315097910000131
(7)获取所有砂岩样品的压溶成因石英胶结物含量:根据部分A砂岩样品的压溶成因石英胶结物含量的比例系数,结合铸体薄片上定量统计的总体石英胶结物含量Tqc,评价出部分B砂岩样品的压溶成因石英胶结物含量;联合部分A砂岩样品和部分B砂岩样品,从而获得了所有砂岩样品的压溶成因石英胶结物含量EPSqc。
本发明提出的一种粒间压溶成因石英胶结物含量的评价方法,联合XX盆地XX凹陷目标层段砂岩的步骤(1)、步骤(2)、步骤(3)、步骤(4)、步骤(5)、步骤(6)和步骤(7)评价过程,评价出粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩的压溶成因石英胶结物含量;图4是粗砂岩的压溶成因石英胶结物含量的剖面图;图5是中砂岩的压溶成因石英胶结物含量的剖面图;图6是砂质砾岩的压溶成因石英胶结物含量的剖面图;图7是细砂岩的压溶成因石英胶结物含量的剖面图;图8是粉砂岩的压溶成因石英胶结物含量的剖面图;对比分析实测压溶成因石英胶结物含量与石英胶结作用模型模拟结果确保了化学压实作用的准确性(图9)。
在此基础上,重现XX盆地XX凹陷目标层段砂岩的温度演化史、机械压实作用下砂岩孔隙度演化史和化学压实作用下砂岩孔隙度减少量;图10是研究区目标层位砂岩的温度演化史;图11是机械压实作用下目标层位砂岩孔隙度演化史;图12是化学压实作用下目标层位砂岩的孔隙度响应特征。在此基础上,重现了目标层位砂岩的化学压实作用过程。
本发明创新地提出了一种粒间压溶成因石英胶结物含量的评价方法,为压实作用数值模拟研究提供关键参数;这有利于重现致密砂岩沉积时期和压实作用条件下孔隙度演化过程。这一定意义上提供了低渗透-致密砂岩储层精细评价的需求,为后续低渗透-致密砂岩储层的有利目标优选提供技术支撑。
实施例2:如图1、图2、图3、图4、图5、图6、图7、图8、图9、图10、图11、图12、图13、图14及图15所示,一种粒间压溶成因石英胶结物含量的评价方法,含有以下步骤;收集整理研究目标层位砂岩的铸体薄片、阴极发光薄片和扫描电镜图像资料;砂岩样品分类为部分A和部分B,部分A砂岩样品包括铸体薄片、扫描电镜和阴极发光薄片图像资料,部分B砂岩样品包含铸体薄片和扫描电镜图像资料;在铸体薄片和阴极发光薄片图像上定量统计出总体石英胶结物含量Tqc;在阴极发光薄片图像上定量统计出压溶成因石英胶结物含量PSqc;基于部分A砂岩样品的总体石英胶结物含量和压溶成因石英胶结物含量,评价出压溶成因石英胶结物的比例系数;结合部分B砂岩样品的铸体薄片上的总体石英胶结物含量,评价出所有砂岩样品的压溶成因石英胶结物含量EPSqc。
还包括以下步骤:
步骤1)、获取研究区图像资料和取芯样品,收集整理研究区内目标层位砂岩的铸体薄片、阴极发光薄片和扫描电镜等图像分析测试资料,并且收集钻井取芯井段上的取芯砂岩样品;
步骤2)、分类统计砂岩样品,这些砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩;将研究区内收集到的砂岩样品分类为部分A和部分B;部分A砂岩样品包括铸体薄片、阴极发光薄片和扫描电镜图像资料,或者制作铸体薄片、阴极发光薄片和扫描电镜样品块;部分B砂岩样品包含铸体薄片和阴极发光薄片图像资料,或者制作铸体薄片和扫描电镜样品块;
步骤3)、定量地统计铸体薄片上总体石英胶结物含量,在Axio Vision SE 64图像分析系统下采用图像分析软件拍摄铸体薄片的图像资料,统计出部分A砂岩样品和部分B砂岩样品的总体石英胶结物含量;
步骤4)、定量地统计阴极发光薄片上总体石英胶结物和压溶成因石英胶结物含量;在CL8200 Mk5-2阴极发光系统下采用图像分析软件拍摄部分A砂岩样品的阴极发光薄片图像资料,恢复出石英颗粒初始形态,识别出石英颗粒-石英颗粒的重叠部分,作为石英压溶作用部分;从而确定出压溶成因石英胶结物含量;并且采用图像识别方法在阴极发光薄片图像资料上获取总体石英胶结物含量;
步骤5)、确定铸体薄片上总体石英胶结物含量的校正系数;在阴极发光薄片图像资料上,采用图像识别方法识别出来的石英胶结物含量,具有更高的可靠性;当对比分析铸体薄片上总体石英胶结物含量和阴极发光薄片上总体石英胶结物含量时,采用阴极发光薄片上统计的石英胶结物含量来校正铸体薄片上统计的总体石英胶结物含量;
步骤6)、计算出压溶成因石英胶结物的比例系数,部分A砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩,统计出部分A砂岩样品的总体石英胶结物含量Tqc和压溶成因石英胶结物含量PSqc;在基础上,评价出部分A砂岩样品的压溶成因石英胶结物的比例系数;
步骤7)、获得所有砂岩样品的压溶成因石英胶结物含量;根据部分A砂岩样品的压溶成因石英胶结物的比例系数,结合部分B砂岩样品铸体薄片上的总体石英胶结物含量Tqc,评价出部分B砂岩样品的压溶成因石英胶结物含量;联合部分A砂岩样品和部分B砂岩样品,获得所有砂岩样品的压溶成因石英胶结物含量EPSqc。
步骤3)中总体石英胶结物含量的定量统计,主要的步骤如下;
总体砂岩样品包括部分A砂岩样品和部分B砂岩样品;统计所有砂岩样品铸体薄片上的总体石英胶结物含量;
对于粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩不同粒径岩性的砂岩样品,首先使用Axio Vision SE 64图像分析系统下图像分析软件对每个铸体薄片拍摄20张显微照片,拍摄照片的目镜是10×倍;然后,使用CorelDraw软件在电脑屏幕上画出和识别出每张铸体薄片图像上的石英胶结物;并且通过Image-Pro Plus软件获取每张铸体薄片图像上石英胶结物的总面积;为了满足统计学规律,求取这个铸体薄片的所有显微照片中石英胶结物含量的平均值,从而得到这个铸体薄片的总体石英胶结物含量,见公式(1);这样,对铸体薄片的所有显微照片统计出来的总体石英胶结物含量更具有统计规律;
Tqc=average(figure(Tqci))i=1,2,3…20…公式(1)
式中:Tqc是在铸体薄片上统计的总体石英胶结物含量;figure(Tqci)是铸体薄片的第i张显微照片上统计的石英胶结物含量。
步骤4)中总体石英胶结物和压溶成因石英胶结物含量的定量统计,其主要的技术特征如下;
在CL8200 Mk5-2阴极发光系统下采用图像分析软件拍摄部分A砂岩样品的阴极发光薄片图像资料,恢复出石英颗粒初始形态,识别出石英颗粒-石英颗粒的重叠部分,作为石英压溶作用部分;从而定量统计出压溶成因石英胶结物含量;其步骤如下;
在阴极发光薄片图像资料上识别出石英颗粒-石英颗粒的重叠部分,从而确定出压溶成因石英胶结物含量;
使用阴极发光仪、研究级正立多功能偏光显微镜、摄像头和计算机在加速电压10kV和束电流250uA条件下拍摄阴极发光薄片图像,曝光持续时间10s(放大倍数50×)、30s(放大倍数100×)和60s(放大倍数200×);在阴极发光系统下采用图像分析软件拍摄30张阴极发光薄片图像照片;在阴极发光薄片图像资料上使用一个100(10×10)正方形网格,手动恢复出石英颗粒初始形态,识别出石英颗粒-石英颗粒重叠区域,从而确定出压溶成因石英胶结物含量;然后,对阴极发光薄片的所有显微照片统计出来的压溶成因石英胶结物含量取平均值,这样评价出来的压溶成因石英胶结物含量更具有统计规律,见公式(2);
PSqc=average(figure(PSqci))i=1,2,3…20…公式(2)
式中:PSqc是阴极发光薄片上统计的压溶成因石英胶结物含量;figure(PSqci)是阴极发光薄片上的第i张显微照片上统计的压溶成因石英胶结物含量;
划分出一个100(10×10)正方形网格,在网格空间内恢复出石英颗粒的初始形态,然后识别出石英颗粒-石英颗粒的重叠部分,进而统计出压溶成因石英胶结物含量;
在阴极发光薄片图像资料上,基于石英颗粒和石英胶结物具有不同发光的特征,采用图像分割技术、可视化技术和图像提取方法来刻画出石英胶结物和计算出石英胶结物含量;为了满足统计学规律,对多张阴极发光薄片图像上统计的石英胶结物含量取平均值,见公式(3);
CLqc=average(figure(CLqci))i=1,2,3…20…公式(3)
式中:CLqc是阴极发光薄片上统计的总体石英胶结物含量;figure(CLqc)是阴极发光薄片上的第i张显微照片上图像识别方法统计的总体石英胶结物含量。
步骤5)确定铸体薄片上总体石英胶结物含量的校正系数,主要的步骤如下:
在阴极发光薄片图像资料上,由于石英颗粒和石英胶结物具有不同的发光特征;采用图像识别方法识别出来的石英胶结物含量,具有更高的可靠性;当对比分析铸体薄片上总体石英胶结物含量和阴极发光薄片上总体石英胶结物含量时,采用阴极发光薄片上统计出来的总体石英胶结物含量CLqc来校正铸体薄片上统计出来的总体石英胶结物含量Tqc;这里给出一个校正系数TCqc,见公式(4);
Figure BDA0002315097910000181
式中:Tqc是在铸体薄片上统计的总体石英胶结物含量;CLqc是在阴极发光薄片上统计的总体石英胶结物含量;TCqc是铸体薄片上统计的总体石英胶结物含量的校正系数。
步骤6)中压溶成因石英胶结物的比例系数,其主要的技术特征如下;
基于包含粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩的部分A砂岩样品的铸体薄片、扫描电镜和阴极发光薄片图像资料;在铸体薄片上定量地统计出总体石英胶结物含量Tqc;在阴极发光薄片图像上识别出石英颗粒-石英颗粒的重叠区域,从而获得压溶成因石英胶结物含量PSqc;综合考虑铸体薄片图像上统计的总体石英胶结物含量、铸体薄片上总体石英胶结物含量的校正系数和阴极发光薄片图像上统计的压溶成因石英胶结物含量,建立压溶成因石英胶结物的比例系数评价公式,见公式5;
Figure BDA0002315097910000182
式中:Tqc是铸体薄片上统计的总体石英胶结物含量;PSqc是阴极发光薄片上统计的压溶成因石英胶结物含量;TCqc是铸体薄片上统计的总体石英胶结物含量的校正系数;Rpsqc是指压溶成因石英胶结物含量的比例系数。
步骤7)中获得所有砂岩样品的压溶成因石英胶结物含量,主要步骤如下;
由于铸体薄片制作便宜、方便且技术较为成熟,因而铸体薄片图像资料比阴极发光薄片图像资料更丰富;并且铸体薄片图像资料的定量统计资料比阴极发光薄片图像资料的定量统计资料更为丰富;因而,采用压溶成因石英胶结物的比例系数来评价出更多的压溶成因石英胶结物含量;
在收集到的资料中,具有铸体薄片图像资料的粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩的样品数据更多;而具有阴极发光薄片图像资料的粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩的样品数据有限;也就是说,可以获得总体石英胶结物含量的样品数据点数远大于可以获得压溶成因石英胶结物含量的样品数据点数;
根据部分A砂岩样品的铸体薄片图像资料的所有石英胶结物含量和阴极发光薄片图像资料的压溶成因石英胶结物含量而计算的压溶成因石英胶结物含量的比例系数;结合部分B砂岩样品的铸体薄片上定量统计的总体石英胶结物含量Tqc,评价出部分B砂岩样品的压溶成因石英胶结物含量,见公式(6);
EPSqc=Rpsqc×Tqc×TCqc 公式(6)
式中:Tqc是铸体薄片上定量统计的总体石英胶结物含量;Rpsqc是指压溶成因石英胶结物含量的比例系数;TCqc是铸体薄片上统计的总体石英胶结物含量的校正系数;EPSqc是部分B砂岩样品的评价的压溶成因石英胶结物含量;联合部分A砂岩样品和部分B砂岩样品,就获得了所有砂岩样品的压溶成因石英胶结物含量EPSqc。
结合步骤1)、步骤2)、步骤3)、步骤4)、步骤5)、步骤6)和步骤7)来评价出粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩的压溶成因石英胶结物含量;其主要的技术特征如下;
步骤1)、获取研究区图像资料和取芯样品,收集整理研究区内目标层位砂岩的铸体薄片、阴极发光薄片和扫描电镜等图像分析测试资料,并且收集钻井取芯井段上的取芯砂岩样品;
步骤2)、分类统计砂岩样品,这些砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩;将研究区内收集到的砂岩样品分类为部分A和部分B;部分A砂岩样品包括铸体薄片、阴极发光薄片和扫描电镜图像资料,或者制作铸体薄片、阴极发光薄片和扫描电镜样品块;部分B砂岩样品包含铸体薄片和阴极发光薄片图像资料,或者制作铸体薄片和扫描电镜样品块;
步骤3)、定量地统计铸体薄片上总体石英胶结物含量,在Axio Vision SE 64图像分析系统下采用图像分析软件拍摄铸体薄片的图像资料,统计出部分A砂岩样品和部分B砂岩样品的总体石英胶结物含量;
步骤4)、定量地统计阴极发光薄片上总体石英胶结物和压溶成因石英胶结物含量;在CL8200 Mk5-2阴极发光系统下采用图像分析软件拍摄部分A砂岩样品的阴极发光薄片图像资料,恢复出石英颗粒初始形态,识别出石英颗粒-石英颗粒的重叠部分,作为石英压溶作用部分;从而确定出压溶成因石英胶结物含量;并且采用图像识别方法在阴极发光薄片图像资料上获取总体石英胶结物含量;
步骤5)、确定铸体薄片上总体石英胶结物含量的校正系数;在阴极发光薄片图像资料上,采用图像识别方法识别出来的石英胶结物含量,具有更高的可靠性;当对比分析铸体薄片上总体石英胶结物含量和阴极发光薄片上总体石英胶结物含量时,采用阴极发光薄片上统计的石英胶结物含量来校正铸体薄片上统计的总体石英胶结物含量;
步骤6)、计算出压溶成因石英胶结物的比例系数,部分A砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩,统计出部分A砂岩样品的总体石英胶结物含量Tqc和压溶成因石英胶结物含量PSqc;在基础上,评价出部分A砂岩样品的压溶成因石英胶结物的比例系数;
步骤7)、获得所有砂岩样品的压溶成因石英胶结物含量;根据部分A砂岩样品的压溶成因石英胶结物含量的比例系数,结合部分B砂岩样品铸体薄片上的总体石英胶结物含量Tqc,评价出部分B砂岩样品的压溶成因石英胶结物含量;联合部分A砂岩样品和部分B砂岩样品,就获得了所有砂岩样品的压溶成因石英胶结物含量EPSqc。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种粒间压溶成因石英胶结物含量的评价方法,其特征在于含有以下步骤;收集整理研究目标层位砂岩的铸体薄片、阴极发光薄片和扫描电镜图像资料;砂岩样品分类为部分A和部分B,部分A砂岩样品包括铸体薄片、扫描电镜和阴极发光薄片图像资料,部分B砂岩样品包含铸体薄片和扫描电镜图像资料;在铸体薄片和阴极发光薄片图像上定量统计出总体石英胶结物含量Tqc;在阴极发光薄片图像上定量统计出压溶成因石英胶结物含量PSqc;基于部分A砂岩样品的总体石英胶结物含量和压溶成因石英胶结物含量,根据部分A砂岩样品的压溶成因石英胶结物的比例系数,结合部分B砂岩样品铸体薄片上的总体石英胶结物含量Tqc,评价出部分B砂岩样品的压溶成因石英胶结物含量;联合部分A砂岩样品和部分B砂岩样品,获得所有砂岩样品的压溶成因石英胶结物含量EPSqc。
2.如权利要求1所述的一种粒间压溶成因石英胶结物含量的评价方法,其特征在于还包括以下步骤:
步骤1)、获取研究区图像资料和取芯样品,收集整理研究区内目标层位砂岩的铸体薄片、阴极发光薄片和扫描电镜图像分析测试资料,并且收集钻井取芯井段上的取芯砂岩样品;
步骤2)、分类统计砂岩样品,这些砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩;将研究区内收集到的砂岩样品分类为部分A和部分B;
部分A砂岩样品包括制作铸体薄片、阴极发光薄片和扫描电镜样品块;
部分B砂岩样品包含制作铸体薄片和扫描电镜样品块;
步骤3)、定量地统计铸体薄片上总体石英胶结物含量,在Axio Vision SE 64图像分析系统下采用图像分析软件拍摄铸体薄片的图像资料,统计出部分A砂岩样品和部分B砂岩样品的总体石英胶结物含量;
步骤4)、定量地统计阴极发光薄片上总体石英胶结物和压溶成因石英胶结物含量;在CL8200 Mk5-2阴极发光系统下采用图像分析软件拍摄部分A砂岩样品的阴极发光薄片图像资料,恢复出石英颗粒初始形态,识别出石英颗粒-石英颗粒的重叠部分,作为石英压溶作用部分;从而确定出压溶成因石英胶结物含量;并且采用图像识别方法在阴极发光薄片图像资料上获取总体石英胶结物含量;
步骤5)、确定铸体薄片上总体石英胶结物含量的校正系数;在阴极发光薄片图像资料上,采用图像识别方法识别出来的石英胶结物含量,具有更高的可靠性;当对比分析铸体薄片上总体石英胶结物含量和阴极发光薄片上总体石英胶结物含量时,采用阴极发光薄片上统计的石英胶结物含量来校正铸体薄片上统计的总体石英胶结物含量;
步骤6)、计算出压溶成因石英胶结物的比例系数,部分A砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩,统计出部分A砂岩样品的总体石英胶结物含量Tqc和压溶成因石英胶结物含量PSqc;评价出部分A砂岩样品的压溶成因石英胶结物的比例系数;
步骤7)、获得所有砂岩样品的压溶成因石英胶结物含量;根据部分A砂岩样品的压溶成因石英胶结物的比例系数,结合部分B砂岩样品铸体薄片上的总体石英胶结物含量Tqc,评价出部分B砂岩样品的压溶成因石英胶结物含量;联合部分A砂岩样品和部分B砂岩样品,获得所有砂岩样品的压溶成因石英胶结物含量EPSqc。
3.如权利要求2所述的一种粒间压溶成因石英胶结物含量的评价方法,其特征在于步骤3)中总体石英胶结物含量的定量统计,步骤如下;
总体砂岩样品包括部分A砂岩样品和部分B砂岩样品;统计所有砂岩样品铸体薄片上的总体石英胶结物含量;
对于粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩不同粒径岩性的砂岩样品,首先使用Axio Vision SE 64图像分析系统下图像分析软件对每个铸体薄片拍摄20张显微照片,拍摄照片的目镜是10×倍;然后,使用CorelDraw软件在电脑屏幕上画出和识别出每张铸体薄片图像上的石英胶结物;并且通过Image-Pro Plus软件获取每张铸体薄片图像上石英胶结物的总面积;求取铸体薄片的所有显微照片中石英胶结物含量的平均值,得到这个铸体薄片的总体石英胶结物含量,见公式(1);
Tqc=average(figure(Tqci)) i=1,2,3…20…公式(1)
式中:Tqc是在铸体薄片上统计的总体石英胶结物含量;figure(Tqci)是铸体薄片的第i张显微照片上统计的石英胶结物含量。
4.如权利要求2所述的一种粒间压溶成因石英胶结物含量的评价方法,其特征在于步骤4)中总体石英胶结物和压溶成因石英胶结物含量的定量统计,步骤如下;
在CL8200 Mk5-2阴极发光系统下采用图像分析软件拍摄部分A砂岩样品的阴极发光薄片图像资料,恢复出石英颗粒初始形态,识别出石英颗粒-石英颗粒的重叠部分,作为石英压溶作用部分;从而定量统计出压溶成因石英胶结物含量;评价步骤如下;
在阴极发光薄片图像资料上识别出石英颗粒-石英颗粒的重叠部分,从而确定出压溶成因石英胶结物含量;
使用阴极发光仪、研究级正立多功能偏光显微镜、摄像头和计算机在加速电压10kV和束电流250uA条件下拍摄阴极发光薄片图像,曝光持续时间10s(放大倍数50×)、30s(放大倍数100×)和60s(放大倍数200×);在阴极发光系统下采用图像分析软件拍摄30张阴极发光薄片图像照片;在阴极发光薄片图像资料上使用一个100(10×10)正方形网格,手动恢复出石英颗粒初始形态,识别出石英颗粒-石英颗粒重叠区域,从而确定出压溶成因石英胶结物含量;然后,对阴极发光薄片的所有显微照片统计出来的压溶成因石英胶结物含量取平均值,这样评价出来的压溶成因石英胶结物含量更具有统计规律,见公式(2);
PSqc=average(figure(PSqci)) i=1,2,3…20…公式(2)
式中:PSqc是阴极发光薄片上统计的压溶成因石英胶结物含量;figure(PSqci)是阴极发光薄片上的第i张显微照片上统计的压溶成因石英胶结物含量;
划分出一个100(10×10)正方形网格,在网格空间内恢复出石英颗粒的初始形态,然后识别出石英颗粒-石英颗粒的重叠部分,进而统计出压溶成因石英胶结物含量;
在阴极发光薄片图像资料上,基于石英颗粒和石英胶结物具有不同发光的特征,采用图像分割技术、可视化技术和图像提取方法来刻画出石英胶结物和计算出石英胶结物含量;对多张阴极发光薄片图像上统计的石英胶结物含量取平均值,见公式(3);
CLqc=average(figure(CLqci)) i=1,2,3…20…公式(3)
式中:CLqc是阴极发光薄片上统计的总体石英胶结物含量;figure(CLqc)是阴极发光薄片上的第i张显微照片上图像识别方法统计的总体石英胶结物含量。
5.如权利要求2所述的一种粒间压溶成因石英胶结物含量的评价方法,其特征在于步骤5)确定铸体薄片上总体石英胶结物含量的校正系数,步骤如下:
在阴极发光薄片图像资料上,根据石英颗粒和石英胶结物具有不同的发光特征,采用图像识别方法识别出来的石英胶结物含量;当对比分析铸体薄片上总体石英胶结物含量和阴极发光薄片上总体石英胶结物含量时,采用阴极发光薄片上统计出来的总体石英胶结物含量CLqc来校正铸体薄片上统计出来的总体石英胶结物含量Tqc;这里给出一个校正系数TCqc,见公式(4);
Figure FDA0003432345270000041
式中:Tqc是在铸体薄片上统计的总体石英胶结物含量;CLqc是在阴极发光薄片上统计的总体石英胶结物含量;TCqc是铸体薄片上统计的总体石英胶结物含量的校正系数。
6.如权利要求2所述的一种粒间压溶成因石英胶结物含量的评价方法,其特征在于步骤6)中压溶成因石英胶结物的比例系数,步骤如下;
基于包含粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩的部分A砂岩样品的铸体薄片、扫描电镜和阴极发光薄片图像资料;在铸体薄片上定量地统计出总体石英胶结物含量Tqc;在阴极发光薄片图像上识别出石英颗粒-石英颗粒的重叠区域,从而获得压溶成因石英胶结物含量PSqc;综合考虑铸体薄片图像上统计的总体石英胶结物含量、铸体薄片上总体石英胶结物含量的校正系数和阴极发光薄片图像上统计的压溶成因石英胶结物含量,建立压溶成因石英胶结物的比例系数评价公式,见公式5;
Figure FDA0003432345270000051
式中:Tqc是铸体薄片上统计的总体石英胶结物含量;PSqc是阴极发光薄片上统计的压溶成因石英胶结物含量;TCqc是铸体薄片上统计的总体石英胶结物含量的校正系数;Rpsqc是指压溶成因石英胶结物的比例系数。
7.如权利要求2所述的一种粒间压溶成因石英胶结物含量的评价方法,其特征在于步骤7)中获得所有砂岩样品的压溶成因石英胶结物含量,采用压溶成因石英胶结物的比例系数来评价出更多的压溶成因石英胶结物含量;步骤如下;
根据部分A砂岩样品的铸体薄片图像资料的所有石英胶结物含量和阴极发光薄片图像资料的压溶成因石英胶结物含量而计算的压溶成因石英胶结物的比例系数;结合部分B砂岩样品的铸体薄片上定量统计的总体石英胶结物含量Tqc,评价出部分B砂岩样品的压溶成因石英胶结物含量,见公式(6);
EPSqc=Rpsqc×Tqc×TCqc 公式(6)
式中:Tqc是铸体薄片上定量统计的总体石英胶结物含量;Rpsqc是指压溶成因石英胶结物含量的比例系数;TCqc是铸体薄片上统计的总体石英胶结物含量的校正系数;EPSqc是部分B砂岩样品的评价的压溶成因石英胶结物含量;联合部分A砂岩样品和部分B砂岩样品,就获得了所有砂岩样品的压溶成因石英胶结物含量EPSqc。
8.如权利要求3所述的一种粒间压溶成因石英胶结物含量的评价方法,其特征在于结合步骤1)、步骤2)、步骤3)、步骤4)、步骤5)、步骤6)和步骤7)来评价出粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩的压溶成因石英胶结物含量;步骤如下;
步骤1)、获取研究区图像资料和取芯样品,收集整理研究区内目标层位砂岩的铸体薄片、阴极发光薄片和扫描电镜图像分析测试资料,并且收集钻井取芯井段上的取芯砂岩样品;
步骤2)、分类统计砂岩样品,这些砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩;将研究区内收集到的砂岩样品分类为部分A和部分B;部分A砂岩样品包括铸体薄片、阴极发光薄片和扫描电镜图像资料,或者制作铸体薄片、阴极发光薄片和扫描电镜样品块;部分B砂岩样品包含铸体薄片和阴极发光薄片图像资料,或者制作铸体薄片和扫描电镜样品块;
步骤3)、定量地统计铸体薄片上总体石英胶结物含量,在Axio Vision SE 64图像分析系统下采用图像分析软件拍摄铸体薄片的图像资料,统计出部分A砂岩样品和部分B砂岩样品的总体石英胶结物含量;
步骤4)、定量地统计阴极发光薄片上总体石英胶结物和压溶成因石英胶结物含量;在CL8200 Mk5-2阴极发光系统下采用图像分析软件拍摄部分A砂岩样品的阴极发光薄片图像资料,恢复出石英颗粒初始形态,识别出石英颗粒-石英颗粒的重叠部分,作为石英压溶作用部分;从而确定出压溶成因石英胶结物含量;并且采用图像识别方法在阴极发光薄片图像资料上获取总体石英胶结物含量;
步骤5)、确定铸体薄片上总体石英胶结物含量的校正系数;在阴极发光薄片图像资料上,采用图像识别方法识别出来的石英胶结物含量,具有更高的可靠性;当对比分析铸体薄片上总体石英胶结物含量和阴极发光薄片上总体石英胶结物含量时,采用阴极发光薄片上统计的石英胶结物含量来校正铸体薄片上统计的总体石英胶结物含量;
步骤6)、计算出压溶成因石英胶结物的比例系数,部分A砂岩样品包括粗砂岩、中砂岩、砂质砾岩、细砂岩和粉砂岩,统计出部分A砂岩样品的总体石英胶结物含量Tqc和压溶成因石英胶结物含量PSqc;在基础上,评价出部分A砂岩样品的压溶成因石英胶结物的比例系数;
步骤7)、获得所有砂岩样品的压溶成因石英胶结物含量;根据部分A砂岩样品的压溶成因石英胶结物含量的比例系数,结合部分B砂岩样品铸体薄片上的总体石英胶结物含量Tqc,评价出部分B砂岩样品的压溶成因石英胶结物含量;联合部分A砂岩样品和部分B砂岩样品,获得了所有砂岩样品的压溶成因石英胶结物含量EPSqc。
CN201911288542.2A 2019-12-12 2019-12-12 一种粒间压溶成因石英胶结物含量的评价方法 Active CN110988004B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911288542.2A CN110988004B (zh) 2019-12-12 2019-12-12 一种粒间压溶成因石英胶结物含量的评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911288542.2A CN110988004B (zh) 2019-12-12 2019-12-12 一种粒间压溶成因石英胶结物含量的评价方法

Publications (2)

Publication Number Publication Date
CN110988004A CN110988004A (zh) 2020-04-10
CN110988004B true CN110988004B (zh) 2022-04-01

Family

ID=70093610

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911288542.2A Active CN110988004B (zh) 2019-12-12 2019-12-12 一种粒间压溶成因石英胶结物含量的评价方法

Country Status (1)

Country Link
CN (1) CN110988004B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105334149A (zh) * 2015-11-24 2016-02-17 中国石油大学(北京) 一种致密储层微观孔隙结构评价与储层分类方法
CN105334150A (zh) * 2015-11-24 2016-02-17 中国石油大学(北京) 前陆盆地砂岩储层成岩过程与孔隙演化过程的确定方法
WO2016161914A1 (zh) * 2015-04-07 2016-10-13 四川行之智汇知识产权运营有限公司 一种利用地质和测井信息预测储层成岩相的方法
CN106597548A (zh) * 2016-12-02 2017-04-26 中国石油大学(华东) 一种多因素定量评价地质时期3d孔隙度的方法
CN109298464A (zh) * 2018-10-17 2019-02-01 中国石油天然气股份有限公司 致密砂岩储层成岩相测井识别方法及装置
CN109709301A (zh) * 2018-11-30 2019-05-03 中国石油天然气股份有限公司 一种裂缝孔隙型致密砂岩储层分布确定方法、装置及系统
CN109856029A (zh) * 2019-02-01 2019-06-07 中海石油(中国)有限公司上海分公司 一种基于图像分析的孔隙度评价方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107228816B (zh) * 2017-06-01 2018-11-13 中国石油大学(华东) 一种泥页岩中不同类型孔隙演化评价方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016161914A1 (zh) * 2015-04-07 2016-10-13 四川行之智汇知识产权运营有限公司 一种利用地质和测井信息预测储层成岩相的方法
CN105334149A (zh) * 2015-11-24 2016-02-17 中国石油大学(北京) 一种致密储层微观孔隙结构评价与储层分类方法
CN105334150A (zh) * 2015-11-24 2016-02-17 中国石油大学(北京) 前陆盆地砂岩储层成岩过程与孔隙演化过程的确定方法
CN106597548A (zh) * 2016-12-02 2017-04-26 中国石油大学(华东) 一种多因素定量评价地质时期3d孔隙度的方法
CN109298464A (zh) * 2018-10-17 2019-02-01 中国石油天然气股份有限公司 致密砂岩储层成岩相测井识别方法及装置
CN109709301A (zh) * 2018-11-30 2019-05-03 中国石油天然气股份有限公司 一种裂缝孔隙型致密砂岩储层分布确定方法、装置及系统
CN109856029A (zh) * 2019-02-01 2019-06-07 中海石油(中国)有限公司上海分公司 一种基于图像分析的孔隙度评价方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Reservoir quality, natural fractures, and gas productivity of upper Triassic Xujiahe tight gas sandstones in western Sichuan Basin, China;Dali Yue etc.;《Marine and Petroleum Geology》;20171212;第89卷;第370-386页 *
应用阴极发光技术探讨含矿岩系成岩作用及其与成矿关系--以云南金顶铅锌矿床为例;宋志敏等;《现代地质》;19920630;第6卷(第2期);第139-149页 *
砂岩储层中石英胶结作用的形成机理及评价;陈振林等;《2011 International Conference on Information, Services and Management Engineering》;20111231;第1451-1457页 *
砂岩压实作用研究现状及进展;林承焰等;《沉积学报》;20191204;第1-15页 *
胶莱盆地莱阳凹陷莱阳组储层成岩作用及孔隙演化;周勇等;《吉林大学学报》;20130331;第43卷(第2期);第340-350页 *

Also Published As

Publication number Publication date
CN110988004A (zh) 2020-04-10

Similar Documents

Publication Publication Date Title
Liu et al. Multi-scale fractal analysis of pores in shale rocks
CN107228816B (zh) 一种泥页岩中不同类型孔隙演化评价方法
WO2017114443A1 (zh) 一种确定储层岩溶发育程度的方法及装置
CN110021030B (zh) 一种岩土体材料数字图像的分割阈值确定方法
Åkesson et al. The influence of foliation on the fragility of granitic rocks, image analysis and quantitative microscopy
CN104374682B (zh) 一种岩心ct扫描分析方法及装置
CN109856029B (zh) 一种基于图像分析的孔隙度评价方法
CN105649615A (zh) Ct定量、三维可视化测试储层致密油赋存状态的方法
CN103278436A (zh) 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法
Liu et al. A new method for threshold determination of gray image
CN109540929B (zh) 一种测定盆地砂岩型铀矿成矿年龄的方法
CN103063558A (zh) 一种机制砂颗粒形状评价方法
CN111521539B (zh) 一种致密砂砾岩储层中不同类型孔隙结构定量表征方法
CN103308433A (zh) 一种基于孔隙演化分析评价致密砂岩储层成岩相的方法
GB2524810A (en) Method of analysing a drill core sample
Li et al. Morphology-based indices and recommended sampling sizes for using image-based methods to quantify degradations of compacted aggregate materials
CN109444190B (zh) 利用数字岩心分析提高复杂油气层含水饱和度评价精度的方法
Barbera et al. A multi-technique approach for the determination of the porous structure of building stone
CN103149130B (zh) 砾岩岩心颗粒结构中粒度的分析方法
CN110988004B (zh) 一种粒间压溶成因石英胶结物含量的评价方法
Lewis et al. An automated system for the statistical analysis of sediment texture and structure at the micro scale
Cole et al. Faunal and sedimentological analysis of a latest Silurian stromatoporoid biostrome from the central Appalachian Basin
CN112348330A (zh) 确定蚀变火山碎屑砂砾岩储层定量分类的方法
Dąbski Application of the Handysurf E-35B electronic profilometer for the study of weathering micro-relief in glacier forelands in SE Iceland
CN115931948B (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