CN108169095A - 中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法 - Google Patents
中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法 Download PDFInfo
- Publication number
- CN108169095A CN108169095A CN201711329867.1A CN201711329867A CN108169095A CN 108169095 A CN108169095 A CN 108169095A CN 201711329867 A CN201711329867 A CN 201711329867A CN 108169095 A CN108169095 A CN 108169095A
- Authority
- CN
- China
- Prior art keywords
- sandstone
- particle
- mid
- feldspar
- range
- 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
- 238000009933 burial Methods 0.000 title claims abstract description 75
- 239000011148 porous material Substances 0.000 title claims abstract description 61
- XYQHCDPZBXIAGW-UHFFFAOYSA-N Andesine Natural products COC(=O)C1=Cc2ccc3c(CCN(C)C)cc(OC)c(O)c3c2C(=O)O1 XYQHCDPZBXIAGW-UHFFFAOYSA-N 0.000 title claims abstract description 58
- 229910052658 andesine Inorganic materials 0.000 title claims abstract description 58
- 238000003556 assay Methods 0.000 title claims abstract description 23
- 239000002245 particle Substances 0.000 claims abstract description 331
- 239000010433 feldspar Substances 0.000 claims abstract description 307
- 239000000523 sample Substances 0.000 claims abstract description 167
- 239000013049 sediment Substances 0.000 claims abstract description 146
- 239000013074 reference sample Substances 0.000 claims abstract description 115
- 239000000203 mixture Substances 0.000 claims abstract description 70
- 238000005266 casting Methods 0.000 claims abstract description 44
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000011084 recovery Methods 0.000 claims abstract description 7
- 239000004575 stone Substances 0.000 claims description 22
- 239000003973 paint Substances 0.000 claims description 18
- 239000011435 rock Substances 0.000 claims description 17
- 238000003776 cleavage reaction Methods 0.000 claims description 16
- 230000007017 scission Effects 0.000 claims description 16
- 239000013078 crystal Substances 0.000 claims description 15
- 239000004576 sand Substances 0.000 claims description 11
- 239000004568 cement Substances 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 abstract description 27
- 238000005056 compaction Methods 0.000 abstract description 7
- 238000005259 measurement Methods 0.000 abstract description 7
- 230000006378 damage Effects 0.000 abstract description 5
- 239000009671 shengli Substances 0.000 description 21
- 239000010438 granite Substances 0.000 description 6
- 239000010453 quartz Substances 0.000 description 5
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 5
- 238000004090 dissolution Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000004512 die casting Methods 0.000 description 3
- 239000003822 epoxy resin Substances 0.000 description 3
- 229920000647 polyepoxide Polymers 0.000 description 3
- 238000003703 image analysis method Methods 0.000 description 2
- 239000004615 ingredient Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 102100030762 Apolipoprotein L1 Human genes 0.000 description 1
- 101100237637 Bos taurus APOL gene Proteins 0.000 description 1
- 101100323521 Homo sapiens APOL1 gene Proteins 0.000 description 1
- 230000002378 acidificating effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000004106 carminic acid Substances 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000000879 optical micrograph Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000004062 sedimentation 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
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Soil Conditioners And Soil-Stabilizing Materials (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明提出一种中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，包括以下步骤：确定与中深层砂岩样品相似的地表现代沉积物样品作为地表现代沉积物参考样品，分别磨制中深层砂岩样品和地表现代沉积物参考样品的铸体薄片，并对其中的每个长石颗粒进行形态恢复和溶解孔隙圈绘，求取中深层砂岩样品和地表现代沉积物参考样品中长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量，并将二者作差，进而求得中深层砂岩埋藏成岩中长石颗粒的溶解孔隙面孔率。该方法消除了地表长石溶解作用形成的次生孔隙对测定的影响，并恢复了压实作用和胶结作用破坏的长石溶解形成的次生孔隙，其测定结果更准确。
Description
技术领域
本发明属于石油天然气勘探与开发领域，尤其涉及一种中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法。
背景技术
随着油气需求的不断提高，以及中浅层油气勘探程度和难度的不断增加，油气勘探的目标正逐步转向埋深大于2500m的中深层。由于埋藏成岩过程中压实作用和胶结作用会破坏孔隙，溶解作用会形成次生孔隙，因而，随着埋藏深度增加，砂岩中原生孔隙含量逐渐减小、次生孔隙含量逐渐增加，具体表现为：在埋深2500-3500m的中层砂岩中，储集空间以次生孔隙为主；在埋深大于3500m的深层砂岩中，原生孔隙可能消失殆尽。
长石是我国陆相砂岩的重要组类成矿物，埋藏成岩过程中长石溶解形成的次生孔隙被认为是中深层砂岩油气储层重要的次生孔隙型之一，对中深层砂岩埋藏成岩中长石溶解形成的次生孔隙面孔率进行定量测定，对中深层油气勘探至关重要。长石溶解次生孔隙面孔率指长石溶解形成的次生孔隙面积占岩石铸体薄片总面积的百分含量，是影响中深层砂岩优质油气储层评价和预测的关键因素之一，决定了中深层砂岩油气勘探是以寻找埋藏成岩过程中长石溶解形成大量次生孔隙发育区为目标，还是以寻找原生孔隙发育的有利沉积砂体为目标。
目前，国内外学者主要依据Schmidt(1979)建立的次生孔隙识别标志，利用肉眼估计法、点计数法、图像分析法等方法，定量确定中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率。
然而，这些方法均存在不足：(1)地表现代碎屑沉积物中发育大量长石溶解形成的次生孔隙，中深层砂岩中发育的全部形态的长石溶解次生孔隙在地表现代碎屑沉积物中均有发育，然而，在现有的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法中，并未考虑现代地表碎屑沉积物中长石溶解形成的次生孔隙含量，这导致其测定结果出现严重错误；(2)中深层砂岩经历了较强的压实作用和胶结作用，长石溶解形成的部分次生孔隙在埋藏成岩过程中会因压实作用和胶结作用而被破坏，现有的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法不能确定埋藏成岩过程中因压实作用和胶结作用破坏的长石溶解次生孔隙含量，导致其测定结果出现误差。
因此，急需发展一种更准确的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法。
发明内容
本发明针对上述的现有的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法中存在的不足，提出一种中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，其消除了地表长石溶解作用形成的次生孔隙对中深层砂岩埋藏成岩过程中长石溶解作用形成的次生孔隙面孔率测定的影响，并恢复了压实作用和胶结作用破坏的长石溶解形成的次生孔隙，其测定结果更准确，更符合中深层真实的地质情况。
为了达到上述目的，本发明采用的技术方案为：
一种中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，包括以下步骤：
(1)通过野外踏勘，确定与中深层砂岩样品在岩石组构特征和沉积条件上相似的地表现代沉积物样品作为地表现代沉积物参考样品；
(2)分别磨制中深层砂岩样品和地表现代沉积物参考样品的铸体薄片，其中，磨制地表现代沉积物参考样品的铸体薄片时，需针对地表现代沉积物参考样品中的不同粒径范围的样品分别磨制对应的铸体薄片；
(3)分别获取每个铸体薄片在透射光显微镜下对应的图像，并对图像中的每个长石颗粒的形态进行恢复，得到其溶解前形态，测量获得每个长石颗粒对应的溶解前面积和粒径；针对每个长石颗粒，圈绘长石颗粒边缘与其溶解前形态的轮廓之间的区域作为粒缘溶解孔隙，圈绘长石颗粒内部溶解区域作为粒内溶解孔隙，测量获得每个长石颗粒对应的粒缘溶解孔隙面积和粒内溶解孔隙面积；
(4)针对地表现代沉积物参考样品，求取地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；针对中深层砂岩样品，依据与步骤(2)中地表现代沉积物参考样品相同的粒径范围划分标准将长石颗粒划分为不同粒径范围，求取中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；
(5)将求得的中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量与地表现代沉积物参考样品中对应粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量作差，获得中深层砂岩中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量，进而求得中深层砂岩中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积，以及中深层砂岩埋藏成岩中长石颗粒的溶解孔隙面孔率。
作为优选，步骤(1)中判定地表现代沉积物样品与中深层砂岩样品在岩石组构特征和沉积条件上相似的标准为：当地表现代沉积物样品与所述中深层砂岩样品在岩石成分、长石颗粒的粒径范围、长石颗粒的磨圆特征和物源区母岩类型上均相同，且地表现代沉积物样品的沉积物搬运距离小于或等于所述中深层砂岩样品的最大延伸距离时，判定为相似。
作为优选，步骤(2)中，磨制地表现代沉积物参考样品的铸体薄片时，将地表现代沉积物参考样品划分为0.1-0.25mm、0.25-0.5mm、0.5-1.0mm和1.0-2.0mm四个粒径范围。
作为优选，步骤(3)中，恢复长石颗粒形态的具体步骤为：针对双晶或解理发育的长石颗粒，沿其双晶或解理的延伸方向将长石颗粒恢复成长条形，作为其溶解前形态；针对无双晶、无解理且具有未溶解的平直边缘的长石颗粒，沿其平直边缘的延伸方向将长石颗粒恢复成长条形，作为其溶解前形态；针对无双晶、无解理且无平直边缘的长石颗粒，沿其长轴方向将长石颗粒恢复成长条形，作为其溶解前形态。
作为优选，步骤(3)中，将长石颗粒恢复成长条形后，当图像中长石颗粒的长条形轮廓间呈线接触或缝合接触状态时，需将其恢复成点接触状态后再进行圈绘；当图像中长石颗粒的长条形轮廓内填充有胶结物时，需将胶结物剔除后再进行圈绘。
作为优选，步骤(4)中，求取中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量的具体步骤为：
利用式(1)-式(4)分别计算中深层砂岩样品中不同粒径范围长石颗粒的溶解前总面积，式(1)-式(4)的表达式分别为：
式(1)-式(4)中，Sa为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的长石颗粒的溶解前总面积，μm2；Sb为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的长石颗粒的溶解前总面积，μm2；Sc为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的长石颗粒的溶解前总面积，μm2；Sd为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的长石颗粒的溶解前总面积，μm2；Si为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的第i个长石颗粒的溶解前面积，μm2；Sj为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的第j个长石颗粒的溶解前面积，μm2；Sp为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的第p个长石颗粒的溶解前面积，μm2；Sq为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的第q个长石颗粒的溶解前面积，μm2；m1为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的长石颗粒的个数；m2为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的长石颗粒的个数；m3为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的长石颗粒的个数；m4为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的长石颗粒的个数；
利用式(5)-式(8)分别计算中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积，式(5)-式(8)的表达式分别为：
式(5)-式(8)中，A1为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的长石颗粒的溶解孔隙总面积，μm2；A2为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的长石颗粒的溶解孔隙总面积，μm2；A3为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的长石颗粒的溶解孔隙总面积，μm2；A4为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的长石颗粒的溶解孔隙总面积，μm2；S1yi为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的第i个长石颗粒的粒缘溶解孔隙面积，μm2；S2yj为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的第j个长石颗粒的粒缘溶解孔隙面积，μm2；S3yp为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的第p个长石颗粒的粒缘溶解孔隙面积，μm2；S4yq为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的第q个长石颗粒的粒缘溶解孔隙面积，μm2；S1ni为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的第i个长石颗粒的粒内溶解孔隙面积，μm2；S2nj为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的第j个长石颗粒的粒内溶解孔隙面积，μm2；S3np为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的第p个长石颗粒的粒内溶解孔隙面积，μm2；S4nq为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的第q个长石颗粒的粒内溶解孔隙面积，μm2；
利用式(9)-式(12)分别计算中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量，式(9)-式(12)的表达式分别为：
P1＝A1/Sa×100％ (9)
P2＝A2/Sb×100％ (10)
P3＝A3/Sc×100％ (11)
P4＝A4/Sd×100％ (12)
式(9)-式(12)中，P1为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；P2为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；P3为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；P4为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量。
作为优选，步骤(4)中，求取地表现代沉积物参考样品中不同粒径范围的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量的具体步骤为：
利用式(13)-式(16)分别计算地表现代沉积物参考样品中不同粒径范围的长石颗粒的溶解前总面积，式(13)-式(16)的表达式分别为：
式(13)-式(16)中，S'a为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的长石颗粒的溶解前总面积，μm2；S'b为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的长石颗粒的溶解前总面积，μm2；S'c为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的长石颗粒的溶解前总面积，μm2；S'd为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的长石颗粒的溶解前总面积，μm2；S'i为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的第i个长石颗粒的溶解前面积，μm2；S'j为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的第j个长石颗粒的溶解前面积，μm2；S'p为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的第p个长石颗粒的溶解前面积，μm2；S'q为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的第q个长石颗粒的溶解前面积，μm2；n1为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的长石颗粒的个数；n2为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的长石颗粒的个数；n3为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的长石颗粒的个数；n4为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的长石颗粒的个数；
利用式(17)-式(20)分别计算地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解孔隙总面积，式(17)-式(20)的表达式分别为：
式(5)-式(8)中，A'1为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的长石颗粒的溶解孔隙总面积，μm2；A'2为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的长石颗粒的溶解孔隙总面积，μm2；A'3为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的长石颗粒的溶解孔隙总面积，μm2；A'4为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的长石颗粒的溶解孔隙总面积，μm2；S'1yi为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的第i个长石颗粒的粒缘溶解孔隙面积，μm2；S'2yj为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的第j个长石颗粒的粒缘溶解孔隙面积，μm2；S'3yp为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的第p个长石颗粒的粒缘溶解孔隙面积，μm2；S'4yq为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的第q个长石颗粒的粒缘溶解孔隙面积，μm2；S'1ni为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的第i个长石颗粒的粒内溶解孔隙面积，μm2；S'2nj为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的第j个长石颗粒的粒内溶解孔隙面积，μm2；S'3np为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的第p个长石颗粒的粒内溶解孔隙面积，μm2；S'4nq为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的第q个长石颗粒的粒内溶解孔隙面积，μm2；
利用式(21)-式(24)分别计算地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量，式(21)-式(24)的表达式分别为：
P'1＝A'1/S'a×100％ (21)
P'2＝A'2/S'b×100％ (22)
P'3＝A'3/S'c×100％ (23)
P'4＝A'4/S'd×100％ (24)
式(21)-式(24)中，P'1为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；P'2为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；P'3为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；P'4为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量。
作为优选，步骤(5)中，计算中深层砂岩中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量的公式为：
ΔP1＝P1-P'1 (25)
ΔP2＝P2-P'2 (26)
ΔP3＝P3-P'3 (27)
ΔP4＝P4-P'4 (28)
式(25)-式(28)中，ΔP1为中深层砂岩中恢复后粒径为0.1-0.25mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量；ΔP2为中深层砂岩中恢复后粒径为0.25-0.5mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量；ΔP3为中深层砂岩中恢复后粒径为0.5-1.0mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量；ΔP4为中深层砂岩中恢复后粒径为1.0-2.0mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量。
作为优选，步骤(5)中，求取中深层砂岩中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积的公式为：
St1＝ΔP1×Sa (29)
St2＝ΔP2×Sb (30)
St3＝ΔP3×Sc (31)
St4＝ΔP4×Sd (32)
式(29)-式(32)中，St1为中深层砂岩中恢复后粒径为0.1-0.25mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积，μm2；St2为中深层砂岩中恢复后粒径为0.25-0.5mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积，μm2；St3为中深层砂岩中恢复后粒径为0.5-1.0mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积，μm2；St4为中深层砂岩中恢复后粒径为1.0-2.0mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积，μm2。
作为优选，步骤(5)中，求取中深层砂岩埋藏成岩中长石颗粒的溶解孔隙面孔率的公式为：
Pt＝(St1+St2+St3+St4)/(Sa+Sb+Sc+Sd) (33)
P＝Pt*Vf (34)
式(33)-式(34)中，Pt为中深层砂岩埋藏成岩中长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；Vf为中深层砂岩样品铸体薄片中长石颗粒面积占铸体薄片面积的百分含量，通过点计数法求得；P为中深层砂岩埋藏成岩中长石颗粒的溶解孔隙面孔率。
与现有技术相比，本发明的优点和积极效果在于：
1、本发明提供的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，首次通过将中深层砂岩长石溶解形成的次生孔隙量与地表现代沉积物中长石溶解孔隙量作差的方法，消除了地表长石溶解作用形成的次生孔隙对中深层砂岩埋藏成岩过程中长石溶解作用形成的次生孔隙面孔率测定的影响，修正了现有方法的理论错误；
2、本发明提供的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，首次恢复了压实作用和胶结作用破坏的长石溶解形成的次生孔隙，其测定结果更准确，更符合中深层真实的地质情况；
3、本发明提供的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，首次考虑了长石颗粒粒径对长石溶解形成的次生孔隙的影响，大大提高了中深层砂岩埋藏成岩过程中长石溶解形成的次生孔隙面孔率的测定精度；
4、本发明提供的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法测定的面孔率与原有的图像分析方法获得的面孔率，二者对勘探方向的指导意义完全不同，原有的图像分析方法获得的面孔率指导勘探以寻找酸性流体优势运移方向为主，采用本发明提供的方法测定的面孔率指导勘探以寻找有利的沉积砂体为主。
附图说明
图1为本发明实施例所提供的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法的流程图；
图2为本发明实施例所提供的在透射光显微镜下铸体薄片中长石颗粒图像的示意图；其中，a为正交光下具有卡斯巴双晶的长石颗粒，b为正交光下具有聚片双晶的长石颗粒，c为完全解理的长石颗粒，d为正交光下一级灰白干涉色的长石颗粒；
图3为本发明实施例所提供的长石颗粒形态恢复和溶解孔隙圈绘的示意图；
图4为本发明实施例所提供的胜利油田某地区沙四段3343.9m砂岩的岩石学特征(单偏光)；其中，a为透射光显微镜下砂岩铸体薄片中花岗岩岩屑的特征，b为透射光显微镜下砂岩铸体薄片的特征(图中Q指石英，F指长石，R指岩屑)；
图5为本发明实施例所提供的青岛某山区河道沙的岩石学特征；其中，a为地表现代沉积物的野外取样特征，b为透射光显微镜下地表现代沉积物铸体薄片的特征(单偏光)(图中Q指石英，F指长石)；
图6为本发明实施例所提供的青岛某山区地表现代沉积物取样的物源区花岗岩的岩石学特征；其中，a为花岗岩野外露头特征，b为透射光显微镜下花岗岩铸体薄片的特征(正交光)(图中Q指石英，F指长石)；
图7为本发明实施例所提供的胜利油田某地区沙四段3343.9m砂岩样品的长石颗粒的形态恢复和溶解孔隙圈绘示意图；其中，a为样品铸体薄片在透射光显微镜下的原貌特征，b为长石颗粒的形态恢复结果，c为长石颗粒点接触状态的恢复结果，d为长石颗粒的溶解孔隙圈绘结果；
图8为本发明实施例所提供的青岛某山区地表现代沉积物参考样品的长石颗粒的形态恢复和溶解孔隙圈绘示意图；其中，a为样品铸体薄片在透射光显微镜下的原貌特征，b为长石颗粒的形态恢复和溶解孔隙圈绘结果。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明实施例提供了一种中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，其流程图如图1所示，包括以下步骤：
(1)通过野外踏勘，确定与中深层砂岩样品在岩石组构特征和沉积条件上相似的地表现代沉积物样品作为地表现代沉积物参考样品。
在本步骤中，需要说明的是，岩石组构特征和沉积条件可依据样品的成分、结构、构造、颜色等特征确定。为了便于确定地表现代沉积物参考样品，步骤(1)中判定地表现代沉积物样品与中深层砂岩样品在岩石组构特征和沉积条件上相似的一种优选标准为：当地表现代沉积物样品与所述中深层砂岩样品在岩石成分、长石颗粒的粒径范围、长石颗粒的磨圆特征和物源区母岩类型上均相同，且地表现代沉积物样品的沉积物搬运距离小于或等于所述中深层砂岩样品的最大延伸距离时，判定为相似。
(2)分别磨制中深层砂岩样品和地表现代沉积物参考样品的铸体薄片，其中，磨制地表现代沉积物参考样品的铸体薄片时，需针对地表现代沉积物参考样品中的不同粒径范围的样品分别磨制对应的铸体薄片。
在本步骤中，需要说明的是，制备铸体薄片的具体步骤为本领域技术人员所熟知，例如可采用ZT-2型高压铸体仪充注蓝色环氧树脂，制作铸体砂岩样品，采用TX-PG-250调速磨片机对铸体砂岩样品进行磨制，获得铸体薄片。
此外，还需要说明的是，磨制地表现代沉积物参考样品的铸体薄片时，对不同粒径范围的样品分别磨制对应的铸体薄片，是为了保证所选择的地表现代沉积物参考样品与中深层砂岩样品在粒级分类方面具有可比性。由于在本领域常规沉积学参考书中定义了砂岩的粒级分类为：细砂为0.1-0.25mm，中砂为0.25-0.5mm，粗砂为0.5-1mm，巨砂为1.0-2.0mm。因而，作为一种优选，磨制地表现代沉积物参考样品的铸体薄片时，将地表现代沉积物参考样品划分为0.1-0.25mm、0.25-0.5mm、0.5-1.0mm和1.0-2.0mm四个粒径范围。
(3)分别获取每个铸体薄片在透射光显微镜下对应的图像，并对图像中的每个长石颗粒的形态进行恢复，得到其溶解前形态，测量获得每个长石颗粒对应的溶解前面积和粒径；针对每个长石颗粒，圈绘长石颗粒边缘与其溶解前形态的轮廓之间的区域作为粒缘溶解孔隙，圈绘长石颗粒内部溶解区域作为粒内溶解孔隙，测量获得每个长石颗粒对应的粒缘溶解孔隙面积和粒内溶解孔隙面积。
在本步骤中，需要说明的是，如图2所示，在透射光显微镜下铸体薄片中长石颗粒的图像有四种。如图3所示，作为一种优选，恢复长石颗粒形态的具体步骤为：针对双晶或解理发育的长石颗粒，沿其双晶或解理的延伸方向将长石颗粒恢复成长条形，作为其溶解前形态；针对无双晶、无解理且具有未溶解的平直边缘的长石颗粒，沿其平直边缘的延伸方向将长石颗粒恢复成长条形，作为其溶解前形态；针对无双晶、无解理且无平直边缘的长石颗粒，沿其长轴方向将长石颗粒恢复成长条形，作为其溶解前形态。恢复后，以长条形轮廓的面积为长石颗粒的溶解前面积，以恢复后长石颗粒溶解前形态的长轴作为粒径。
此外，还需要说明的是，为了消除压实作用对溶解孔隙的破坏，将长石颗粒恢复成长条形后，当图像中长石颗粒的长条形轮廓间呈线接触或缝合接触状态时，需将其恢复成点接触状态后再进行圈绘。为了消除胶结作用对溶解孔隙的破坏，当图像中长石颗粒的长条形轮廓内填充有胶结物时，需将胶结物剔除后再进行圈绘。
(4)针对地表现代沉积物参考样品，求取地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；针对中深层砂岩样品，依据与步骤(2)中地表现代沉积物参考样品相同的粒径范围划分标准将长石颗粒划分为不同粒径范围，求取中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量。
在本步骤中，分别求取了地表现代沉积物参考样品和中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量，以便于后续作差。
其中，求取中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量的具体步骤为：
利用式(1)-式(4)分别计算中深层砂岩样品中不同粒径范围长石颗粒的溶解前总面积，式(1)-式(4)的表达式分别为：
利用式(5)-式(8)分别计算中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积，式(5)-式(8)的表达式分别为：
利用式(9)-式(12)分别计算中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量，式(9)-式(12)的表达式分别为：
P1＝A1/Sa×100％ (9)
P2＝A2/Sb×100％ (10)
P3＝A3/Sc×100％ (11)
P4＝A4/Sd×100％ (12)
求取地表现代沉积物参考样品中不同粒径范围的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量的具体步骤为：
利用式(13)-式(16)分别计算地表现代沉积物参考样品中不同粒径范围的长石颗粒的溶解前总面积，式(13)-式(16)的表达式分别为：
利用式(17)-式(20)分别计算地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解孔隙总面积，式(17)-式(20)的表达式分别为：
利用式(21)-式(24)分别计算地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量，式(21)-式(24)的表达式分别为：
P'1＝A'1/S'a×100％ (21)
P'2＝A'2/S'b×100％ (22)
P'3＝A'3/S'c×100％ (23)
P'4＝A'4/S'd×100％ (24)。
(5)将求得的中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量与地表现代沉积物参考样品中对应粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量作差，获得中深层砂岩中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量，进而求得中深层砂岩中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积，以及中深层砂岩埋藏成岩中长石颗粒的溶解孔隙面孔率。
在本步骤中，通过将求得的中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量与地表现代沉积物参考样品中对应粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量作差，能够消除地表长石溶解作用形成的次生孔隙对中深层砂岩埋藏成岩过程中长石溶解作用形成的次生孔隙面孔率测定的影响，使测定结果更准确。
其中，计算中深层砂岩中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量的公式为：
ΔP1＝P1-P'1 (25)
ΔP2＝P2-P'2 (26)
ΔP3＝P3-P'3 (27)
ΔP4＝P4-P'4 (28)
求取中深层砂岩中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积的公式为：
St1＝ΔP1×Sa (29)
St2＝ΔP2×Sb (30)
St3＝ΔP3×Sc (31)
St4＝ΔP4×Sd (32)
求取中深层砂岩埋藏成岩中长石颗粒的溶解孔隙面孔率的公式为：
Pt＝(St1+St2+St3+St4)/(Sa+Sb+Sc+Sd) (33)
P＝Pt*Vf (34)
为了更清楚详细地介绍本发明实施例所提供的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，下面将结合具体实施例进行描述。
选取胜利油田某地区沙四段3343.9m砂岩作为中深层砂岩样品，测定其埋藏成岩中长石溶解次生孔隙面孔率，包括以下步骤：
第一步：确定地表现代沉积物参考样品
(1)胜利油田某地区沙四段3343.9m砂岩的岩石学特征如图4所示，由图4a可见，沙四段3343.9m砂岩的物源区母岩类型主要为花岗岩，且其沉积单元的最大延伸距离为3.2km；由图4b可见，沙四段3343.9m砂岩为由长石和石英颗粒组成的长石砂岩，长石颗粒粒径为0.1.0-2.0mm，磨圆为次棱角状。
(2)通过对青岛某山区(N37°，E120°)的野外踏勘，获得一地表现代沉积物样品，其岩石学特征如图5和图6所示。由图5可见，该地表现代沉积物为由长石、石英组成的松散沙，长石颗粒粒径为0.1.0-2.0mm，磨圆为次棱角状。由图6可见，该地表现代沉积物样品的物源区母岩类型为花岗岩，且从物源到地表现代沉积物的沉积物搬运距离不超过2.5km。
(3)由于该地表现代沉积物样品与胜利油田某地区沙四段3343.9m砂岩样品在岩石成分、长石颗粒的粒径范围、长石颗粒的磨圆特征和物源区母岩类型上均相同，且地表现代沉积物样品的沉积物搬运距离小于胜利油田某地区沙四段3343.9m砂岩样品的最大延伸距离，因而，该地表现代沉积物样品与胜利油田某地区沙四段3343.9m砂岩样品在岩石组构特征和沉积条件上相似，作为地表现代沉积物参考样品。
第二步：分别磨制中深层砂岩样品和地表现代沉积物参考样品的铸体薄片
(1)采用ZT-2型高压铸体仪对胜利油田某地区沙四段3343.9m砂岩样品充注蓝色环氧树脂，制作铸体砂岩样品，采用TX-PG-250调速磨片机对铸体砂岩样品进行磨制，获得铸体薄片。
(2)分别用150目、65目、35目、20目、10目的标准样品筛对从青岛某山区获得的地表现代沉积物参考样品进行筛选，筛选出0.1-0.25mm、0.25-0.5mm、0.5-1.0mm和1.0-2.0mm四个粒径范围的样品，采用ZT-2型高压铸体仪分别对筛选出的不同粒径范围的样品充注蓝色环氧树脂，制作铸体砂岩样品，采用TX-PG-250调速磨片机对铸体砂岩样品进行磨制，获得铸体薄片。
第三步：分别对中深层砂岩样品和地表现代沉积物参考样品中的长石颗粒进行形态恢复和溶解孔隙圈绘
(1)利用蔡司Axioscope A1 APOL.数字透反偏光显微镜及摄像系统分别获取每个铸体薄片在透射光显微镜下对应的图像，并对图像中的每个长石颗粒的形态进行恢复，得到其溶解前形态。恢复长石颗粒形态的具体步骤为：针对双晶或解理发育的长石颗粒，沿其双晶或解理的延伸方向将长石颗粒恢复成长条形，作为其溶解前形态；针对无双晶、无解理且具有未溶解的平直边缘的长石颗粒，沿其平直边缘的延伸方向将长石颗粒恢复成长条形，作为其溶解前形态；针对无双晶、无解理且无平直边缘的长石颗粒，沿其长轴方向将长石颗粒恢复成长条形，作为其溶解前形态。上述长石颗粒的形态恢复过程在Axio VisionSoftware Rel.软件中完成，恢复后，以长条形轮廓的面积为长石颗粒的溶解前面积，以恢复后长石颗粒溶解前形态的长轴作为粒径。
(2)针对每个长石颗粒，圈绘长石颗粒边缘与其溶解前形态的轮廓之间的区域作为粒缘溶解孔隙，圈绘长石颗粒内部溶解区域作为粒内溶解孔隙，测量获得每个长石颗粒对应的粒缘溶解孔隙面积和粒内溶解孔隙面积。其中，将长石颗粒恢复成长条形后，当图像中长石颗粒的长条形轮廓间呈线接触或缝合接触状态时，需将其恢复成点接触状态后再进行圈绘；当图像中长石颗粒的长条形轮廓内填充有胶结物时，需将胶结物剔除后再进行圈绘。上述长石颗粒的溶解孔隙圈绘过程，在Axio Vision Software Rel.软件中完成。
具体的，针对胜利油田某地区沙四段3343.9m砂岩样品，其长石颗粒的形态恢复和溶解孔隙圈绘的部分结果如图7所示；针对从青岛某山区获得的地表现代沉积物参考样品，其长石颗粒的形态恢复和溶解孔隙圈绘的部分结果如图8所示。
第四步：分别求取中深层砂岩样品和地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量
(1)针对胜利油田某地区沙四段3343.9m砂岩样品，依据形态恢复后的粒径将长石颗粒划分为0.1-0.25mm、0.25-0.5mm、0.5-1.0mm和1.0-2.0mm四个粒径范围，其中，0.1-0.25mm粒径范围内共计33个长石颗粒，0.25-0.5mm粒径范围内共计55个长石颗粒，0.5-1.0mm粒径范围内共计17个长石颗粒，1.0-2.0mm粒径范围内共计9个长石颗粒。统计每个长石颗粒的溶解前面积、粒缘溶解孔隙面积和粒内溶解孔隙面积，以粒径为0.1-0.25mm范围内的长石颗粒为例，表1给出了统计数据。
表1胜利油田某地区沙四段3343.9m砂岩样品中恢复后粒径为0.1-0.25mm范围内的长石颗粒的溶解前面积、粒缘溶解孔隙面积和粒内溶解孔隙面积数据
| 长石颗粒编号 | Si/μm2 | S1yi/μm2 | S1ni/μm2 |
| 1 | 21323.1 | 5463.33 | 653.14 |
| 2 | 12313.61 | 3457.8 | - |
| 3 | 20923.53 | 4391.41 | - |
| 4 | 21300.05 | 4372.2 | - |
| 5 | 25476.3 | 8959.55 | - |
| 6 | 27915.97 | 7430.42 | - |
| 7 | 31400.66 | 7161.49 | - |
| 8 | 63585.09 | 19421.3 | - |
| 9 | 28953.31 | 4648.82 | 4710.29 |
| 10 | 27804.55 | 6327.77 | - |
| 11 | 10008.41 | 3642.22 | - |
| 12 | 44129.21 | 5086.81 | - |
| 13 | 27278.2 | 10361.87 | - |
| 14 | 30701.42 | 11679.68 | - |
| 15 | 32948.99 | 12959.06 | - |
| 16 | 21907.08 | 11380 | 3584.59 |
| 17 | 14895.43 | 5624.69 | - |
| 18 | 28807.31 | 9789.41 | - |
| 19 | 42442.57 | 15191.27 | 8294.88 |
| 20 | 9101.7 | 5336.54 | - |
| 21 | 10388.77 | 3876.58 | - |
| 22 | 44782.35 | 15590.83 | - |
| 23 | 37859.06 | 12851.49 | - |
| 24 | 50068.94 | 20039.87 | - |
| 25 | 19006.37 | 6731.18 | - |
| 26 | 22617.85 | 6965.55 | - |
| 27 | 28115.75 | 4022.57 | - |
| 28 | 24550.38 | 6708.13 | - |
| 29 | 33525.29 | 6954.02 | - |
| 30 | 49404.27 | 14330.66 | - |
| 31 | 7111.54 | 3830.47 | - |
| 32 | 35273.4 | 8018.25 | - |
| 33 | 23766.61 | 9505.11 | - |
对其他粒径范围内的长石颗粒做同样的统计，利用上述式(1)-式(4)分别计算胜利油田某地区沙四段3343.9m砂岩样品中不同粒径范围长石颗粒的溶解前总面积，利用上述式(5)-式(8)分别计算胜利油田某地区沙四段3343.9m砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积，利用上述式(9)-式(12)分别计算胜利油田某地区沙四段3343.9m砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量，计算结果如表2所示。
表2胜利油田某地区沙四段3343.9m砂岩样品中不同粒径范围长石颗粒的溶解前总面积、溶解孔隙总面积以及溶解孔隙总面积占其溶解前总面积的百分含量
(2)针对青岛某山区地表现代沉积物参考样品，其0.1-0.25mm粒径范围内的样品中共计64个长石颗粒，0.25-0.5mm粒径范围内的样品中共计45个长石颗粒，0.5-1.0mm粒径范围内的样品中共计53个长石颗粒，1.0-2.0mm粒径范围内的样品中共计11个长石颗粒。统计每个长石颗粒的溶解前面积、粒缘溶解孔隙面积和粒内溶解孔隙面积，以0.1-0.25mm粒径范围内的样品为例，表3给出了统计数据。
表3青岛某山区地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的长石颗粒的溶解前面积、粒缘溶解孔隙面积和粒内溶解孔隙面积数据
对其他粒径范围内的长石颗粒做同样的统计，利用上述式(13)-式(16)分别计算青岛某山区地表现代沉积物参考样品中不同粒径范围的长石颗粒的溶解前总面积，利用上述式(17)-式(20)分别计算青岛某山区地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解孔隙总面积，利用上述式(21)-式(24)分别计算青岛某山区地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量，计算结果如表4所示。
表4青岛某山区地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解前总面积、溶解孔隙总面积以及溶解孔隙总面积占其溶解前总面积的百分含量
第五步：求取中深层砂岩埋藏成岩中长石颗粒的溶解孔隙面孔率
(1)利用上述式(25)-式(28)，计算胜利油田某地区沙四段3343.9m砂岩样品中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量，其结果如表5所示。
表5胜利油田某地区沙四段3343.9m砂岩样品中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量
(2)利用式(29)-式(32)，求取胜利油田某地区沙四段3343.9m砂岩样品中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积，其结果如表6所示。
表6胜利油田某地区沙四段3343.9m砂岩样品中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积
(3)利用上述式(33)求得胜利油田某地区沙四段3343.9m砂岩样品埋藏成岩中长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量Pt为10.60％；通过点计数法求得胜利油田某地区沙四段3343.9m砂岩样品中长石颗粒溶解前总面积占铸体薄片面积的百分含量Vf为20％，进而利用上述式(34)求得胜利油田某地区沙四段3343.9m砂岩样品埋藏成岩中长石颗粒的溶解孔隙面孔率为2.12％。
Claims (10)
1.一种中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，其特征在于，包括以下步骤：
(1)通过野外踏勘，确定与中深层砂岩样品在岩石组构特征和沉积条件上相似的地表现代沉积物样品作为地表现代沉积物参考样品；
(2)分别磨制中深层砂岩样品和地表现代沉积物参考样品的铸体薄片，其中，磨制地表现代沉积物参考样品的铸体薄片时，需针对地表现代沉积物参考样品中的不同粒径范围的样品分别磨制对应的铸体薄片；
(3)分别获取每个铸体薄片在透射光显微镜下对应的图像，并对图像中的每个长石颗粒的形态进行恢复，得到其溶解前形态，测量获得每个长石颗粒对应的溶解前面积和粒径；针对每个长石颗粒，圈绘长石颗粒边缘与其溶解前形态的轮廓之间的区域作为粒缘溶解孔隙，圈绘长石颗粒内部溶解区域作为粒内溶解孔隙，测量获得每个长石颗粒对应的粒缘溶解孔隙面积和粒内溶解孔隙面积；
(4)针对地表现代沉积物参考样品，求取地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；针对中深层砂岩样品，依据与步骤(2)中地表现代沉积物参考样品相同的粒径范围划分标准将长石颗粒划分为不同粒径范围，求取中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；
(5)将求得的中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量与地表现代沉积物参考样品中对应粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量作差，获得中深层砂岩中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量，进而求得中深层砂岩中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积，以及中深层砂岩埋藏成岩中长石颗粒的溶解孔隙面孔率。
2.根据权利要求1所述的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，其特征在于，步骤(1)中判定地表现代沉积物样品与中深层砂岩样品在岩石组构特征和沉积条件上相似的标准为：当地表现代沉积物样品与所述中深层砂岩样品在岩石成分、长石颗粒的粒径范围、长石颗粒的磨圆特征和物源区母岩类型上均相同，且地表现代沉积物样品的沉积物搬运距离小于或等于所述中深层砂岩样品的最大延伸距离时，判定为相似。
3.根据权利要求1所述的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，其特征在于：步骤(2)中，磨制地表现代沉积物参考样品的铸体薄片时，将地表现代沉积物参考样品划分为0.1-0.25mm、0.25-0.5mm、0.5-1.0mm和1.0-2.0mm四个粒径范围。
4.根据权利要求1所述的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，其特征在于，步骤(3)中，恢复长石颗粒形态的具体步骤为：针对双晶或解理发育的长石颗粒，沿其双晶或解理的延伸方向将长石颗粒恢复成长条形，作为其溶解前形态；针对无双晶、无解理且具有未溶解的平直边缘的长石颗粒，沿其平直边缘的延伸方向将长石颗粒恢复成长条形，作为其溶解前形态；针对无双晶、无解理且无平直边缘的长石颗粒，沿其长轴方向将长石颗粒恢复成长条形，作为其溶解前形态。
5.根据权利要求4所述的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，其特征在于：步骤(3)中，将长石颗粒恢复成长条形后，当图像中长石颗粒的长条形轮廓间呈线接触或缝合接触状态时，需将其恢复成点接触状态后再进行圈绘；当图像中长石颗粒的长条形轮廓内填充有胶结物时，需将胶结物剔除后再进行圈绘。
6.根据权利要求3所述的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，其特征在于，步骤(4)中，求取中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量的具体步骤为：
利用式(1)-式(4)分别计算中深层砂岩样品中不同粒径范围长石颗粒的溶解前总面积，式(1)-式(4)的表达式分别为：
式(1)-式(4)中，Sa为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的长石颗粒的溶解前总面积，μm2；Sb为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的长石颗粒的溶解前总面积，μm2；Sc为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的长石颗粒的溶解前总面积，μm2；Sd为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的长石颗粒的溶解前总面积，μm2；Si为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的第i个长石颗粒的溶解前面积，μm2；Sj为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的第j个长石颗粒的溶解前面积，μm2；Sp为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的第p个长石颗粒的溶解前面积，μm2；Sq为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的第q个长石颗粒的溶解前面积，μm2；m1为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的长石颗粒的个数；m2为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的长石颗粒的个数；m3为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的长石颗粒的个数；m4为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的长石颗粒的个数；
利用式(5)-式(8)分别计算中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积，式(5)-式(8)的表达式分别为：
式(5)-式(8)中，A1为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的长石颗粒的溶解孔隙总面积，μm2；A2为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的长石颗粒的溶解孔隙总面积，μm2；A3为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的长石颗粒的溶解孔隙总面积，μm2；A4为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的长石颗粒的溶解孔隙总面积，μm2；S1yi为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的第i个长石颗粒的粒缘溶解孔隙面积，μm2；S2yj为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的第j个长石颗粒的粒缘溶解孔隙面积，μm2；S3yp为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的第p个长石颗粒的粒缘溶解孔隙面积，μm2；S4yq为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的第q个长石颗粒的粒缘溶解孔隙面积，μm2；S1ni为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的第i个长石颗粒的粒内溶解孔隙面积，μm2；S2nj为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的第j个长石颗粒的粒内溶解孔隙面积，μm2；S3np为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的第p个长石颗粒的粒内溶解孔隙面积，μm2；S4nq为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的第q个长石颗粒的粒内溶解孔隙面积，μm2；
利用式(9)-式(12)分别计算中深层砂岩样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量，式(9)-式(12)的表达式分别为：
P1＝A1/Sa×100％ (9)
P2＝A2/Sb×100％ (10)
P3＝A3/Sc×100％ (11)
P4＝A4/Sd×100％ (12)
式(9)-式(12)中，P1为中深层砂岩样品中恢复后粒径为0.1-0.25mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；P2为中深层砂岩样品中恢复后粒径为0.25-0.5mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；P3为中深层砂岩样品中恢复后粒径为0.5-1.0mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；P4为中深层砂岩样品中恢复后粒径为1.0-2.0mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量。
7.根据权利要求6所述的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，其特征在于，步骤(4)中，求取地表现代沉积物参考样品中不同粒径范围的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量的具体步骤为：
利用式(13)-式(16)分别计算地表现代沉积物参考样品中不同粒径范围的长石颗粒的溶解前总面积，式(13)-式(16)的表达式分别为：
式(13)-式(16)中，S'a为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的长石颗粒的溶解前总面积，μm2；S'b为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的长石颗粒的溶解前总面积，μm2；S'c为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的长石颗粒的溶解前总面积，μm2；S'd为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的长石颗粒的溶解前总面积，μm2；S'i为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的第i个长石颗粒的溶解前面积，μm2；S'j为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的第j个长石颗粒的溶解前面积，μm2；S'p为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的第p个长石颗粒的溶解前面积，μm2；S'q为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的第q个长石颗粒的溶解前面积，μm2；n1为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的长石颗粒的个数；n2为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的长石颗粒的个数；n3为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的长石颗粒的个数；n4为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的长石颗粒的个数；
利用式(17)-式(20)分别计算地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解孔隙总面积，式(17)-式(20)的表达式分别为：
式(5)-式(8)中，A'1为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的长石颗粒的溶解孔隙总面积，μm2；A'2为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的长石颗粒的溶解孔隙总面积，μm2；A'3为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的长石颗粒的溶解孔隙总面积，μm2；A'4为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的长石颗粒的溶解孔隙总面积，μm2；S'1yi为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的第i个长石颗粒的粒缘溶解孔隙面积，μm2；S'2yj为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的第j个长石颗粒的粒缘溶解孔隙面积，μm2；S'3yp为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的第p个长石颗粒的粒缘溶解孔隙面积，μm2；S'4yq为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的第q个长石颗粒的粒缘溶解孔隙面积，μm2；S'1ni为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的第i个长石颗粒的粒内溶解孔隙面积，μm2；S'2nj为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的第j个长石颗粒的粒内溶解孔隙面积，μm2；S'3np为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的第p个长石颗粒的粒内溶解孔隙面积，μm2；S'4nq为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的第q个长石颗粒的粒内溶解孔隙面积，μm2；
利用式(21)-式(24)分别计算地表现代沉积物参考样品中不同粒径范围长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量，式(21)-式(24)的表达式分别为：
P'1＝A'1/S'a×100％ (21)
P'2＝A'2/S'b×100％ (22)
P'3＝A'3/S'c×100％ (23)
P'4＝A'4/S'd×100％ (24)
式(21)-式(24)中，P'1为地表现代沉积物参考样品中粒径为0.1-0.25mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；P'2为地表现代沉积物参考样品中粒径为0.25-0.5mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；P'3为地表现代沉积物参考样品中粒径为0.5-1.0mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；P'4为地表现代沉积物参考样品中粒径为1.0-2.0mm范围内的长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量。
8.根据权利要求7所述的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，其特征在于，步骤(5)中，计算中深层砂岩中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量的公式为：
ΔP1＝P1-P'1 (25)
ΔP2＝P2-P'2 (26)
ΔP3＝P3-P'3 (27)
ΔP4＝P4-P'4 (28)
式(25)-式(28)中，ΔP1为中深层砂岩中恢复后粒径为0.1-0.25mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量；ΔP2为中深层砂岩中恢复后粒径为0.25-0.5mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量；ΔP3为中深层砂岩中恢复后粒径为0.5-1.0mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量；ΔP4为中深层砂岩中恢复后粒径为1.0-2.0mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积占其溶解前总面积的百分含量。
9.根据权利要求8所述的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，其特征在于，步骤(5)中，求取中深层砂岩中不同粒径范围长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积的公式为：
St1＝ΔP1×Sa (29)
St2＝ΔP2×Sb (30)
St3＝ΔP3×Sc (31)
St4＝ΔP4×Sd (32)
式(29)-式(32)中，St1为中深层砂岩中恢复后粒径为0.1-0.25mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积，μm2；St2为中深层砂岩中恢复后粒径为0.25-0.5mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积，μm2；St3为中深层砂岩中恢复后粒径为0.5-1.0mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积，μm2；St4为中深层砂岩中恢复后粒径为1.0-2.0mm范围内的长石颗粒在埋藏成岩过程中产生的溶解孔隙总面积，μm2。
10.根据权利要求9所述的中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法，其特征在于，步骤(5)中，求取中深层砂岩埋藏成岩中长石颗粒的溶解孔隙面孔率的公式为：
Pt＝(St1+St2+St3+St4)/(Sa+Sb+Sc+Sd) (33)
P＝Pt*Vf (34)
式(33)-式(34)中，Pt为中深层砂岩埋藏成岩中长石颗粒的溶解孔隙总面积占其溶解前总面积的百分含量；Vf为中深层砂岩样品铸体薄片中长石颗粒面积占铸体薄片面积的百分含量，通过点计数法求得；P为中深层砂岩埋藏成岩中长石颗粒的溶解孔隙面孔率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711329867.1A CN108169095B (zh) | 2017-12-13 | 2017-12-13 | 中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711329867.1A CN108169095B (zh) | 2017-12-13 | 2017-12-13 | 中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108169095A true CN108169095A (zh) | 2018-06-15 |
| CN108169095B CN108169095B (zh) | 2020-01-10 |
Family
ID=62525847
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711329867.1A Active CN108169095B (zh) | 2017-12-13 | 2017-12-13 | 中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108169095B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109187161A (zh) * | 2018-08-30 | 2019-01-11 | 中国石油大学(北京) | 一种碎屑岩中长石溶蚀程度的定量评价方法 |
| CN114441405A (zh) * | 2021-12-22 | 2022-05-06 | 中国地质大学（北京） | 基于压实和胶结减孔趋势的次生增孔幅度定量评价方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102562048A (zh) * | 2010-12-30 | 2012-07-11 | 长江大学 | 一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法 |
| CN102778421A (zh) * | 2012-07-10 | 2012-11-14 | 中国石油大学(华东) | 地质历史时期砂岩储层渗透率演化恢复方法 |
| CN103344541A (zh) * | 2013-07-08 | 2013-10-09 | 中国石油大学(华东) | 一种泥页岩总孔隙度测定方法 |
| KR20170006663A (ko) * | 2015-07-09 | 2017-01-18 | 주식회사 엘지화학 | 이차 전지 양극 내부 기공 분포 분석 방법 및 이를 위한 고분자 |
| CN106596375A (zh) * | 2016-12-19 | 2017-04-26 | 中国石油大学(华东) | 一种地质历史时期储层孔隙度恢复的方法 |
-
2017
- 2017-12-13 CN CN201711329867.1A patent/CN108169095B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102562048A (zh) * | 2010-12-30 | 2012-07-11 | 长江大学 | 一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法 |
| CN102778421A (zh) * | 2012-07-10 | 2012-11-14 | 中国石油大学(华东) | 地质历史时期砂岩储层渗透率演化恢复方法 |
| CN103344541A (zh) * | 2013-07-08 | 2013-10-09 | 中国石油大学(华东) | 一种泥页岩总孔隙度测定方法 |
| KR20170006663A (ko) * | 2015-07-09 | 2017-01-18 | 주식회사 엘지화학 | 이차 전지 양극 내부 기공 분포 분석 방법 및 이를 위한 고분자 |
| CN106596375A (zh) * | 2016-12-19 | 2017-04-26 | 中国石油大学(华东) | 一种地质历史时期储层孔隙度恢复的方法 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109187161A (zh) * | 2018-08-30 | 2019-01-11 | 中国石油大学(北京) | 一种碎屑岩中长石溶蚀程度的定量评价方法 |
| CN114441405A (zh) * | 2021-12-22 | 2022-05-06 | 中国地质大学（北京） | 基于压实和胶结减孔趋势的次生增孔幅度定量评价方法 |
| CN114441405B (zh) * | 2021-12-22 | 2023-06-06 | 中国地质大学（北京） | 基于压实和胶结减孔趋势的次生增孔幅度定量评价方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108169095B (zh) | 2020-01-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104747183B (zh) | 一种碳酸盐岩储层综合分类方法 | |
| CN109856029B (zh) | 一种基于图像分析的孔隙度评价方法 | |
| CN104564041B (zh) | 基于开发渗透率下限的低渗透碎屑岩储层有效性评价方法 | |
| CN107589470B (zh) | 小区域范围内岩溶古地貌恢复方法 | |
| CN102748016B (zh) | 地质历史时期砂岩储层孔隙度演化恢复方法 | |
| CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
| CN107861171A (zh) | 一种陆相湖盆沉积微相平面分布图编制方法 | |
| CN109324345A (zh) | 叠合盆地油气成藏期岩石孔隙度恢复方法 | |
| CN104089964B (zh) | 基于测井米氏旋回分析方法的测年方法 | |
| US11441420B2 (en) | System and method for using geological analysis for the designing of stimulation operations | |
| CN108072748A (zh) | 一种基于优势岩性的碳酸盐岩岩相表征的方法 | |
| Chandra et al. | Impact of depositional and diagenetic features on petrophysical and rock mechanical properties in Arab-D reservoir equivalent upper Jubaila Formation, Saudi Arabia | |
| BR112015021666B1 (pt) | método implantado por computador para atualizar uma estimativa para uma propriedade de material de um volume e sistema de computador para atualizar uma estimativa para uma propriedade de material de um volume | |
| CN105869060A (zh) | 一种细粒岩微细纹层分类方法 | |
| CN108169095A (zh) | 中深层砂岩埋藏成岩中长石溶解次生孔隙面孔率测定方法 | |
| CN105158802A (zh) | 坳陷湖盆重力流沉积物测井定量识别方法 | |
| CN104991286A (zh) | 一种基于沉积模式的沉积相表征方法 | |
| CN108629459A (zh) | 储层含烃孔隙的检测方法及装置 | |
| CN114897767A (zh) | 一种致密混积岩储层储集空间多尺度表征与储层分类方法 | |
| CN107479100A (zh) | 一种陡倾斜深部隐伏铁矿找矿方法 | |
| CN110095488A (zh) | 一种示踪油气活动的方法 | |
| CN107944671A (zh) | 一种确定地层束缚水饱和度的方法 | |
| CN110244357A (zh) | 一种构造破坏型古油藏的判识方法 | |
| CN108979631A (zh) | 一种白云岩潜山油藏低渗层精细表征方法和装置 | |
| Zhang et al. | Physical modeling of complex lithofacies volcanic edifices and seismic characteristics |
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 |