CN109444015A - 多视域、多尺度下的页岩储层微孔系统鉴定方法 - Google Patents
多视域、多尺度下的页岩储层微孔系统鉴定方法 Download PDFInfo
- Publication number
- CN109444015A CN109444015A CN201811283717.6A CN201811283717A CN109444015A CN 109444015 A CN109444015 A CN 109444015A CN 201811283717 A CN201811283717 A CN 201811283717A CN 109444015 A CN109444015 A CN 109444015A
- Authority
- CN
- China
- Prior art keywords
- shale
- shale reservoir
- hole
- kens
- ken
- 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
- 238000000034 method Methods 0.000 title claims abstract description 28
- 239000011148 porous material Substances 0.000 claims abstract description 52
- 238000011160 research Methods 0.000 claims description 27
- 238000012545 processing Methods 0.000 claims description 6
- 230000003321 amplification Effects 0.000 claims description 5
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 230000000877 morphologic effect Effects 0.000 claims description 3
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 12
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 6
- 239000003345 natural gas Substances 0.000 abstract description 3
- 239000007789 gas Substances 0.000 description 15
- 238000004458 analytical method Methods 0.000 description 10
- XKRFYHLGVUSROY-UHFFFAOYSA-N argon Substances [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 description 7
- 229910052786 argon Inorganic materials 0.000 description 7
- -1 argon ion Chemical class 0.000 description 7
- 239000011435 rock Substances 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 5
- 238000005498 polishing Methods 0.000 description 5
- 239000003079 shale oil Substances 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 239000005416 organic matter Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 238000001000 micrograph Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- BVPWJMCABCPUQY-UHFFFAOYSA-N 4-amino-5-chloro-2-methoxy-N-[1-(phenylmethyl)-4-piperidinyl]benzamide Chemical group COC1=CC(N)=C(Cl)C=C1C(=O)NC1CCN(CC=2C=CC=CC=2)CC1 BVPWJMCABCPUQY-UHFFFAOYSA-N 0.000 description 1
- 244000137852 Petrea volubilis Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 201000009310 astigmatism Diseases 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000005422 blasting Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003703 image analysis method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000001755 magnetron sputter deposition Methods 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 239000003973 paint Substances 0.000 description 1
- 238000010422 painting Methods 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 230000003685 thermal hair damage Effects 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
Classifications
-
- 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/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- 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/22—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 measuring secondary emission from the material
- G01N23/225—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 measuring secondary emission from the material using electron or ion
- G01N23/2251—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 measuring secondary emission from the material using electron or ion using incident electron beams, e.g. scanning electron microscopy [SEM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- 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/10056—Microscopic image
- G06T2207/10061—Microscopic image from scanning electron microscope
-
- 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/30108—Industrial image inspection
- G06T2207/30132—Masonry; Concrete
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Quality & Reliability (AREA)
- Dispersion Chemistry (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及一种天然气勘探技术,特别涉及一种多视域、多尺度下的页岩储层微孔系统鉴定方法。该方法通过对页岩样品的表面实施观察,并采用发射环境扫描电镜采用低倍镜扫描其表面,选定多种孔隙结构的观察区域,对该观察区域采用相对于与低倍镜的高倍镜实施扫描,将孔隙结构实施连续拍摄,制作成多视域页岩储层孔隙结构图像,针对多视域页岩储层孔隙结构图像计算盒维数D(r),从而得到页岩储层孔隙的分形维数。
Description
技术领域
本发明涉及天然气勘探技术,多视域下的页岩储层微孔系统多尺度鉴定及统计技术,主要针对页岩气储层微孔结构进行多尺度观察、分类并与分形维数相结合进行统计计算,是页岩油气储层的孔隙结构鉴定和地质参数的确定方法。
尤其涉及对页岩微孔隙结构的研究,把页岩的微观孔喉和不同类型孔隙的结构特征进行定性和定量的描述、评价。
背景技术
页岩气是指主体位于暗色泥页岩或高碳泥页岩中,主要以吸附状态或游离状态为存在方式的天然气聚集(张金川,2004)。泥页岩中的成分和结构因素,尤其是在细粒泥页岩中,其显著的特点是孔喉结构较为复杂且细小,以微米级孔隙和中孔隙为主,这同时也决定了页岩气赋存状态以吸附作用为主(于炳松,2012)。正是因为非常规页岩储层孔隙结构的特殊性质和天然气赋存状态、方式的特殊性,带来了一系列常规油气储层的评价方法和体系难以适用于非常规页岩油气储层的难题。
针对页岩气储层微孔尺度分类及统计,目前国内还未见有较好的方法和技术能够对微孔结构特征进行多视域、多尺度鉴定、统计和分类,而页岩储层微孔结构体系描述对于页岩储层研究则是一项关键技术。
目前对于页岩扫描电镜实验分析方法主要有以下两类:
1)新鲜断面分析方法
由于页岩的沉积环境、埋藏深度等条件的不同,页岩油气储层的物性发育特征有着很大的区别,其孔隙尺度从微米级到纳米级广泛分布。对于微米级孔隙较为发育的页岩通过采用常规的制样方法能够在扫描电镜下直接观察。常规扫面电镜岩样的制备方法步骤如下:①样品的预处理,通常情况下样品的预处理随分析目的而异,一般可按洗油、磨制、酸化、净化和干燥五个程序进行处理,并选取表面平整、大小合适的自然新鲜断面;②将已进行预处理过的样品粘接到样品台上,因地质分析上观察的样品绝大多数是不导电的绝缘体,通常使用导电胶进行固定;③地质岩石样品大多数是高电阻物质,在电子束的轰击下将迅速产生电荷累积形成荷电效应,及扫描电镜观察图像上出现不规则亮区、云雾状象散、图像歪斜漂移等现象,将严重影响图像质量,甚至无法进行分析。因此通常使用镀金仪对样品镀膜后观察,这样能够防止荷电,减少电子束轰击引起的热损伤,增强样品表面的机械稳定性(图1)。
2)氩离子抛光—SEM分析方法
泥页岩结构致密,孔隙结构微小且复杂,自然断面表面粗糙,由于部分样品疏松,还易出现碎屑脱落覆盖孔隙,难以观察纳米级孔隙的大小、形态和结构特征等。而通过利用氩离子抛光技术对样品表面进行处理后,能够获得一个平滑的观察面,且能出去样品表面凹凸不平部分,不会对样品造成机械损害。实验步骤为:①首先将页岩样品用切片机切割成大小合适的小块(约10mm×10mm×3mm),在选定好需要抛光的截面后,用砂纸对样品进行初步打磨,然后固定在抛光仪上利用氩离子进行抛光处理;②将经过氩离子抛光处理后的页岩样品利用导电胶固定在样品台上;③对固定好的样品进行镀金膜处理。
因氩离子抛光—SEM分析方法的样品观察表面平整,无散落,尽可能的减少了样品处理过程中人为因素的影响,能够对页岩孔隙结构大小、形状和结构特点等信息进行直观的观察(图2)。
目前国内外对于研究页岩储层微孔的扫描电镜传统成像分析技术虽被广泛运用,但技术方法并不够成熟,也有其局限性和视域选取随机性,且在观察方式上不够便利。
首先,基于扫描电镜(SEM)的观测结果图像具有视域观察上的局限性。当在低倍扫描电镜下观察时,大孔发育特征显示较为清晰,而微孔及纳米级孔隙发育形态特征无法清楚、准确的识别,对于孔隙结构的鉴定描述并不能反映出整个页岩微孔体系情况;当在高倍扫描电镜下观察时,微孔发育特征能够清晰展示,但无法对较大孔径孔喉结构进行全面、完整的观察(图3)。
其次,对观测图像位置选取具有随机性,样品表面的选取,制作样品的人工损坏以及电镜下观察的可能无法代表非均质性较强的整个页岩样品,且受到矿物颗粒大小影响,在固定的观察视域内,若颗粒越大,所能观察到的颗粒数量越少,相对应所能观察到的孔隙也就越少,这样便难以反映出页岩微孔结构体系的特征,且对于孔隙结构的辨识和统计时研究者根据自己的关注点和研究重点而忽略对其它细节的把握,易对真实结果产生误差。
分形理论为描述页岩微孔结构特征的不规则性研究提供了科学的方法。分形理论是由Mandelbrot在1975年提出并建立,是一种较为复杂的起源于数学几何学的方法和理论。分形是一种具有自相似性和标度不变形的数学集合,即是在某形态上某些组成部位与整体的相似性,可被认为是一种工具,用于描述自然界中不规则物体或形体,可以模拟出各种自然现象。
传统的实验方法和几何理论共同结合研究也页岩孔隙参数,但并不适合用于研究较为复杂和不规则的微孔隙结构,并且需要其它更加精细的非物理实验方法并与各类高精密仪器相结合来寻找孔隙特征参数,这样得出的特征参数才可以准确的用来描述页岩微孔隙结构。页岩孔隙结构的研究也可以通过建立模型的方法来研究,通过以真实孔隙结构空间为研究基础,并借助于较为成熟的数字图像处理方法,经过计算处理后得到页岩微孔结构的特征及分布,对页岩孔隙结构的研究开创了新的思路。
目前,国内外学者也有采用分形理论来研究多孔介质孔隙的特征:谢晓永(2006)等人采用吸附法研究了岩石孔隙的分形特征,杨峰(2013)等人通过高压压汞法和氮气吸附法分析了页岩孔隙结构的分形特征,安士凯、杨宇(2011)等人研究煤层孔隙的分形特征等。
综合来看,虽然已有学者通过运用分形理论来研究岩石孔隙的分形特征,也在研究中取得了较好的效果,但目前尚未有把页岩孔隙分形特征与扫描电镜下具体孔隙结构形态相结合统计研究,并采用多尺度、多视域新型图像分析方法的先例。
发明内容
页岩油气储层作为多孔介质,孔隙结构具有极强的复杂性和不规则性,正是因为其强烈的复杂性和不均质性,对页岩的整体结构和气体运移赋存的稳定性起着至关重要的作用,且页岩储层发育着丰富的纳米级孔隙,微孔隙的表面形态和结构特征对页岩气的吸附和运移有着极为重要的影响,而其中的纳米级微孔结构难以准确测量。研究多孔介质中多尺度微孔隙结构参数的分布规律和特征,是进一步理解页岩气在地层中的赋存状态、解吸扩散与渗流机理的重要基础。
多视域、多尺度观察页岩微孔隙特征技术对页岩微孔结构体系的研究提供了一种更为精细、准确的研究方法。在此图像分析技术的基础上,通过运用分形理论,将页岩孔隙结构的研究与数字图像处理方法相结合,经过计算处理后得到页岩储层微孔结构及分布特征模型,实现对孔隙的定性和定量分析表征,将为页岩储层微孔结构类型、统计分析、油气渗流能力以及进一步的油气资源潜力评价等研究项目提供新的研究思路和方法。
为实现多视域下的页岩储层微孔系统多尺度鉴定及统计,本发明通过使用场发射环境扫描电镜(FESEM)提供了如下技术方案,步骤如下:
a、对页岩样品进行制备,根据页岩样品观察面特征,在扫描电镜低倍数(5x—50000x)下确定一个研究所需要的观察区域。
b、根据对页岩孔隙结构观察的需要,在步骤a中确定观察区域后,选择合适扫描电镜图像放大倍数(如10000x—300000x),其放大倍数要求确保能够观察到所需最小孔隙结构。根据放大倍数的具体数值,将该矩形观察视域进行网格化处理,分成n×m个网格统计视域,图像倍数越大,n、m值越大,结果越精细(图4)。
c、通过步骤b的网格化处理,在扫描电镜下设定相应的倍数分别对每个网格区域进行观察并连续拍摄孔隙结构图像。
d、依次将步骤c所拍摄网格照片无缝拼接成一幅多视域页岩孔隙形态结构照片,其图像视域与步骤a中所需研究矩形区域相同,即呈现多视域下的页岩储层微孔结构形态。
e、根据步骤d所得到的多视域页岩储层孔隙结构图像,按照研究需要,对不同类型孔隙结构按照形态特征进行鉴定划分、分类,并针对不同类型孔隙完成轮廓形态勾画、上色区分。
f、计算盒维数D(r)。根据e步骤在图像中识别出目标孔隙后,将目标孔隙提取出来,根据盒维数数学定义,利用可全部包含最小孔隙的最小网格,对该图像进行网格化处理,统计包含目标孔隙的网格数N(r)。依次改变网格边长r,并统计包含有目标孔隙的网格数,获得多个边长对应的网格数N(r)。将多组数据代入盒维数公式:
计算得出多个盒维数D(r)。
g、通过步骤f所获得的数值,将多个边长r对应的网格数N(r)代入双对数坐标系中,横坐标为log,纵坐标为log N(r)。可获得多个点(log,log N(r)),并根据线性拟合,作图得到的斜率值为盒维数D(r)的数学期望值,即为页岩储层孔隙的分形维数。
步骤b和步骤c通过使用新型图像分析技术,通过在低倍数下网格化估算矩形观察视域,在高倍数下拍摄孔隙结构照片,实现对页岩储层微孔系统的多视域观察。
步骤d将多视域连续拍摄图片无缝拼接成一幅多视域图片,将研究区域内所需微孔结构系统的呈现出,同时实现了对页岩储层微孔系统的多尺度鉴定。
步骤d通过对多视域、多尺度页岩储层微孔系统照片进行孔隙形态结构识别、鉴定和划分,对各微孔类型能够进行精确、直观观察,全面了解微孔结构的分布情况。
将目标孔隙结构与分形维数相结合,根据盒维数、网格化处理和曲线拟合计算出孔隙分形维数值,该方法可通过循环反复的计算研究逐渐将数据理论研究认识向具体的实际地质影响因素进行逼近,真实反映出地质研究规律。
其中,步骤c中对区域进行观察后呈S型连续拍摄孔隙结构图像。
选择S型连续拍摄可使每副需要后续拼接的图片实施重叠,以此便于后续的拼接,缩短了鉴定时间,提高了鉴定效率。
步骤a中选定的观察区域为矩形观察区域是由于拍摄的照片基本上为矩形,当然,也可借助特定工具将其选定为圆形或三角型区域,可根据需要进行设定。
本发明与现有技术相比较,其发明效果如下:
通过运用新型图像分析技术,即基于多视域下的页岩储层微孔系统多尺度鉴定技术:①尽可能的解决了传统方法在扫描电镜下对岩石样品观察位置选取上的随机性;②通过多视域图像无缝链接,能为研究者提供更为直接、全面的观察方式,减小了在视域观察上的局限性;③提供了一种多尺度孔隙结构观察方式,能够在一幅图像上同时呈现多尺度、多量级的孔隙结构形态。多视域、多尺度观察页岩微孔隙特征技术更为精确的对页岩储层孔隙结构体系做出了统计评价,为页岩油气开发研究提供了基础。
进一步来说,通过把新型图像分析技术与分形维数相结合研究,借助于数字图像处理方法,建立了一种基于“形态-数据-统计”从定性到定量的综合分析,能够更精确的定量描述出孔隙结构的复杂程度和非均质性。
下面结合附图对本发明作进一步描述。
附图说明
图1页岩储层新鲜断面扫描电镜图;
图2氩离子抛光—SEM分析方法页岩表面扫描电镜图;
图3高倍数下扫描电镜页岩微孔结构图;
图4低倍数下扫描电镜视域网格化(n=8,m=8)图;
图5为本发明实施例中N150样品多视域下的页岩储层微孔结构形态12x12拼接图;
图6为本发明实施例中N150样品孔隙划分刻画图;
图7为本发明实施例中粒间(内)孔分形维数双对数坐标线性拟合图;
图8为本发明实施例中微裂缝分形维数双对数坐标线性拟合图;
图9为本发明实施例中有机质孔分形维数双对数坐标线性拟合图。
具体实施方式
案例所选样品为川东南志留系龙马溪组页岩,具备良好的区域地质基础资料,基本地球化学测试资料等,其页岩TOC含量较高,孔隙类型多样复杂,包括粒间(内)孔、溶蚀孔、晶间孔、微裂缝以及有机质孔等多种,采用本发明对其页岩储层微孔系统进行鉴定和统计。
“油气藏地质及开发工程”国家重点实验室提供了场发射环境扫面电子显微镜(Quanta250 FEG),全自动磁控离子溅射镀膜仪(Emitech K55X)等配套测试系统为本方法在实施过程中提供各类实验支持。
针对川东南志留系龙马溪组页岩N150号岩样微孔结构分析为例,采用本发明技术方案的具体实施步骤为:
对N150页岩氩离子抛光样品进行前期处理制备后,在扫描电镜仪器下观察样品抛光面特征,并在扫描电镜低倍数下确定一个研究所需要的观察区域。
根据对页岩孔隙结构观察的需要,确定观察区域后,选择合适扫描电镜图像放大倍数,其放大倍数要求确保能够观察到所需最小孔隙结构,为能够研究此样品中有机质纳米级孔隙,选择放大倍数为30000倍。根据所需要研究的观察区域,将该区域确定为12x12个网格统计区域。
通过网格化处理,在扫描电镜下设定30000倍分别对每个网格区域进行观察并呈S型连续拍摄孔隙结构图像。
将所拍摄N150网格照片无缝拼接成一幅多视域页岩孔隙形态结构照片,即呈现多视域下的页岩储层微孔结构形态(图5)。
拼接获得的多视域页岩储层孔隙结构图像,按照研究需要,对不同类型孔隙结构按照形态特征进行鉴定划分、分类,并针对不同类型孔隙完成轮廓形态勾画、上色区分(图6),对此样品选择研究区域后根据所需研究孔隙类型进行刻画,主要为粒间(内)孔(红色)、微裂缝(蓝色)和有机质孔(黄色)。
计算盒维数D(r)。根据经过对不同类型孔隙进行刻画的图像,分别将三种类型孔隙提取出来,根据盒维数数学定义,利用可全部包含最小孔隙的最小网格,对该图像进行网格化处理,统计包含目标孔隙的网格数N(r)。依次改变网格边长r,并统计包含有目标孔隙的网格数,获得多个边长对应的网格数N(r)。将多组数据代入盒维数公式:计算得出多个盒维数D(r)。
通过计算所获得的不同类型孔隙的多个盒维数值,将多个边长r对应的网格数N(r)代入双对数坐标系中,横坐标为log,纵坐标为log N(r)。可获得多个点(log,log N(r)),并根据线性拟合,作图得到的斜率值为盒维数D(r)的数学期望值(图7、图8、图9),即为页岩储层孔隙的分形维数,从而得出不同类型孔隙分形维数值以及拟合曲线。
本实施例中,通过对高分辨率页岩储层扫描电镜图像进行分析,以一个新的方法和思路确定页岩储层地质参数的方法,提高了对页岩储层地质参数的计算精度,解决了现有技术中对页岩储层地质参数计算不准确的技术问题,以及对页岩储层评价提供了新的研究思路。
Claims (3)
1.一种多视域、多尺度下的页岩储层微孔系统鉴定方法,其特征在于,包括以下步骤:
a、制备页岩样品,根据页岩样品观察面特征,采用发射环境扫描电镜,使用5x—50000x的倍数对页岩样品实施扫描,确定该页岩样品上具有多种孔隙结构的观察区域;
b、根据对页岩孔隙结构观察的需要,在步骤a中确定观察区域后,选择10000x—300000x的倍数进行拍摄,其放大倍数要求确保能够观察到研究所需最小孔隙结构;根据高倍数的具体数值,将该矩形观察视域进行网格化估算,分成n×m个网格统计视域;
c、通过步骤b的网格化估算,在扫描电镜下设定相应的高倍数对该观察区域进行观察并连续拍摄孔隙结构图像;
d、依次将步骤c所拍摄网格照片无缝拼接成一幅多视域、多尺度页岩孔隙形态结构照片,其图像视域与步骤a中所需研究矩形区域相同,即呈现多视域下的页岩储层微孔结构形态;
e、根据步骤d所得到的多视域页岩储层孔隙结构图像,按照研究需要,对不同类型孔隙结构按照形态特征进行鉴定划分、分类,并针对不同类型孔隙完成轮廓形态勾画、充填颜色区分;
f、计算盒维数D(r),根据e步骤在图像中识别出目标孔隙后,将目标孔隙提取出来,根据盒维数数学定义,利用可全部包含最小孔隙的最小网格,对该图像进行网格化处理,统计包含目标孔隙的网格数N(r)。依次改变网格边长r,并统计包含有目标孔隙的网格数,获得多个边长对应的网格数N(r)。将多组数据代入盒维数公式:计算得出多个盒维数D(r);
g、通过步骤f所获得的数值,将多个边长r对应的网格数N(r)代入双对数坐标系中,横坐标为log,纵坐标为log N(r),可获得多个点(log,logN(r)),并根据线性拟合,作图得到的斜率值为盒维数D(r)的数学期望值,即为页岩储层孔隙的分形维数。
2.根据权利要求1所述的多视域、多尺度下的页岩储层微孔系统鉴定方法,其特征在于:所述的步骤c中对区域进行观察后呈S型连续拍摄孔隙结构图像。
3.根据权利要求1或2所述的多视域、多尺度下的页岩储层微孔系统鉴定方法,其特征在于:所述的步骤a中选定的观察区域为矩形观察区域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811283717.6A CN109444015A (zh) | 2018-10-31 | 2018-10-31 | 多视域、多尺度下的页岩储层微孔系统鉴定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811283717.6A CN109444015A (zh) | 2018-10-31 | 2018-10-31 | 多视域、多尺度下的页岩储层微孔系统鉴定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109444015A true CN109444015A (zh) | 2019-03-08 |
Family
ID=65549211
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811283717.6A Pending CN109444015A (zh) | 2018-10-31 | 2018-10-31 | 多视域、多尺度下的页岩储层微孔系统鉴定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109444015A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110411911A (zh) * | 2019-08-28 | 2019-11-05 | 四川大学 | 一种基于分形维数的散粒体几何组成特征评价方法 |
CN111007233A (zh) * | 2019-12-25 | 2020-04-14 | 西南石油大学 | 一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法 |
CN112147169A (zh) * | 2019-06-27 | 2020-12-29 | 中国石油化工股份有限公司 | 储层显微组分的定量方法 |
CN112305000A (zh) * | 2020-10-22 | 2021-02-02 | 西南石油大学 | 一种基于页岩有机质孔形态表征页岩储集能力的方法 |
CN113192119A (zh) * | 2021-05-27 | 2021-07-30 | 宜宾学院 | 一种多尺度孔隙面孔率的定量统计方法 |
CN113466102A (zh) * | 2021-06-25 | 2021-10-01 | 清华大学 | 一种2.n维油藏芯片及其制作方法 |
CN114049252A (zh) * | 2021-09-27 | 2022-02-15 | 中国科学院自动化研究所 | 一种面向序列切片的扫描电镜三维影像获取系统及方法 |
CN114371187A (zh) * | 2021-12-17 | 2022-04-19 | 成都理工大学 | 页岩孔隙定量结构复杂性测定方法、系统、介质、设备 |
CN114965538A (zh) * | 2022-05-25 | 2022-08-30 | 大庆油田有限责任公司 | 陆相页岩有机纳米孔隙识别方法 |
CN116448648A (zh) * | 2023-06-20 | 2023-07-18 | 浙江大学海南研究院 | 一种土体微观孔隙结构的显微观测区间推荐方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101413871A (zh) * | 2008-12-05 | 2009-04-22 | 西北有色金属研究院 | 一种多孔材料孔结构分形维数的计算方法 |
CN101655913A (zh) * | 2009-09-17 | 2010-02-24 | 上海交通大学 | 分形维数的计算机生成图像被动检测方法 |
CN106053316A (zh) * | 2016-06-08 | 2016-10-26 | 燕山大学 | 土体微观结构孔隙分布及分形程序与计算方法 |
CN107067379A (zh) * | 2017-03-16 | 2017-08-18 | 中国科学院地质与地球物理研究所 | 基于三维fib‑sem图像的页岩孔隙定量表征方法 |
US20180120215A1 (en) * | 2016-11-01 | 2018-05-03 | Korea Institute Of Geoscience And Mineral Resources | Method for porosity measurement using sem images of rock samples reacted with a gadolinium compound |
-
2018
- 2018-10-31 CN CN201811283717.6A patent/CN109444015A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101413871A (zh) * | 2008-12-05 | 2009-04-22 | 西北有色金属研究院 | 一种多孔材料孔结构分形维数的计算方法 |
CN101655913A (zh) * | 2009-09-17 | 2010-02-24 | 上海交通大学 | 分形维数的计算机生成图像被动检测方法 |
CN106053316A (zh) * | 2016-06-08 | 2016-10-26 | 燕山大学 | 土体微观结构孔隙分布及分形程序与计算方法 |
US20180120215A1 (en) * | 2016-11-01 | 2018-05-03 | Korea Institute Of Geoscience And Mineral Resources | Method for porosity measurement using sem images of rock samples reacted with a gadolinium compound |
CN107067379A (zh) * | 2017-03-16 | 2017-08-18 | 中国科学院地质与地球物理研究所 | 基于三维fib‑sem图像的页岩孔隙定量表征方法 |
Non-Patent Citations (2)
Title |
---|
HUCHENG DENG,XIAOFEI HU,HUAZHOU ANDY LI,BIN LUO,WEI WANG: "Improved pore-structure characterization in shale formations with FESEM technique", 《JOURNAL OF NATURAL GAS SCIENCE AND ENGINEERING》 * |
张鹏飞,卢双舫,李俊乾等: "基于扫描电镜的页岩微观孔隙结构定量表征", 《中国石油大学学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112147169A (zh) * | 2019-06-27 | 2020-12-29 | 中国石油化工股份有限公司 | 储层显微组分的定量方法 |
CN112147169B (zh) * | 2019-06-27 | 2024-05-14 | 中国石油化工股份有限公司 | 储层显微组分的定量方法 |
CN110411911A (zh) * | 2019-08-28 | 2019-11-05 | 四川大学 | 一种基于分形维数的散粒体几何组成特征评价方法 |
CN111007233A (zh) * | 2019-12-25 | 2020-04-14 | 西南石油大学 | 一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法 |
CN111007233B (zh) * | 2019-12-25 | 2022-03-11 | 西南石油大学 | 一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法 |
CN112305000A (zh) * | 2020-10-22 | 2021-02-02 | 西南石油大学 | 一种基于页岩有机质孔形态表征页岩储集能力的方法 |
CN113192119B (zh) * | 2021-05-27 | 2023-01-06 | 宜宾学院 | 一种多尺度孔隙面孔率的定量统计方法 |
CN113192119A (zh) * | 2021-05-27 | 2021-07-30 | 宜宾学院 | 一种多尺度孔隙面孔率的定量统计方法 |
CN113466102A (zh) * | 2021-06-25 | 2021-10-01 | 清华大学 | 一种2.n维油藏芯片及其制作方法 |
CN114049252A (zh) * | 2021-09-27 | 2022-02-15 | 中国科学院自动化研究所 | 一种面向序列切片的扫描电镜三维影像获取系统及方法 |
CN114371187B (zh) * | 2021-12-17 | 2023-11-03 | 成都理工大学 | 页岩孔隙定量结构复杂性测定方法、系统、介质、设备 |
CN114371187A (zh) * | 2021-12-17 | 2022-04-19 | 成都理工大学 | 页岩孔隙定量结构复杂性测定方法、系统、介质、设备 |
CN114965538A (zh) * | 2022-05-25 | 2022-08-30 | 大庆油田有限责任公司 | 陆相页岩有机纳米孔隙识别方法 |
CN116448648A (zh) * | 2023-06-20 | 2023-07-18 | 浙江大学海南研究院 | 一种土体微观孔隙结构的显微观测区间推荐方法及系统 |
CN116448648B (zh) * | 2023-06-20 | 2023-08-11 | 浙江大学海南研究院 | 一种土体微观孔隙结构的显微观测区间推荐方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109444015A (zh) | 多视域、多尺度下的页岩储层微孔系统鉴定方法 | |
CN105115874B (zh) | 多源信息融合的多组分三维数字岩心构建方法 | |
CN107067379B (zh) | 基于三维fib-sem图像的页岩孔隙定量表征方法 | |
Ruzyla | Characterization of pore space by quantitative image analysis | |
Sutton et al. | Virtual paleontology—an overview | |
CN104237103A (zh) | 一种孔隙连通性定量表征方法及装置 | |
CN106198579A (zh) | 一种测量页岩中有机质含量的方法 | |
CN113609696A (zh) | 基于图像融合的多尺度多组分数字岩心构建方法及系统 | |
Arasan et al. | Effect of particle size and shape on the grain-size distribution using Image analysis | |
CN109115819B (zh) | 一种基于fib-him对页岩有机微孔和介孔进行表征的方法及系统 | |
CN102062725B (zh) | 一种馆藏壁画表面颜料特征的无损分析方法 | |
CN111624146A (zh) | 一种定量表征致密储层全尺度孔径分布特征的方法 | |
Biswal et al. | Modeling of multiscale porous media | |
CN112164103B (zh) | 多数据交叉的野外露头信息采集及建模方法 | |
CN108956392A (zh) | 一种潮滩沉积物类型的无人机识别方法 | |
Toth et al. | Automated analysis of heterogeneous carbon nanostructures by high-resolution electron microscopy and on-line image processing | |
CN107038700A (zh) | 基于gssim的土壤湿度空间分布特征定量分析方法 | |
Babadagli | Scanline method to determine the fractal nature of 2-D fracture networks | |
Jia et al. | Effects of point density on DEM accuracy of airborne LiDAR | |
CN114813728A (zh) | 一种基于测量Dpar值的物源分析方法 | |
Randle | Recent developments in electron backscatter diffraction | |
CN114460053A (zh) | 基于量子点荧光效应的光学元件亚表面缺陷三维重构方法 | |
Cooper et al. | Quantification of detrital, authigenic and porosity components of the Fontainebleau Sandstone: A comparison of conventional optical and combined scanning electron microscope‐based methods of modal analyses | |
CN108931497B (zh) | 激光共聚焦三维重构对烃源岩进行评价的方法和装置 | |
Saghi et al. | Three-dimensional metrology and fractal analysis of dendritic nanostructures |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190308 |
|
RJ01 | Rejection of invention patent application after publication |