高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

致密砂岩微观非均质性定量评价方法研究

陈朝兵 付玲 陈新晶 张涛 解宇航 王泓波 朱玉双

陈朝兵, 付玲, 陈新晶, 张涛, 解宇航, 王泓波, 朱玉双. 致密砂岩微观非均质性定量评价方法研究—以鄂尔多斯盆地华庆地区延长组长6油层组为例[J]. 沉积学报, 2021, 39(5): 1086-1099. doi: 10.14027/j.issn.1000-0550.2021.037
引用本文: 陈朝兵, 付玲, 陈新晶, 张涛, 解宇航, 王泓波, 朱玉双. 致密砂岩微观非均质性定量评价方法研究—以鄂尔多斯盆地华庆地区延长组长6油层组为例[J]. 沉积学报, 2021, 39(5): 1086-1099. doi: 10.14027/j.issn.1000-0550.2021.037
CHEN ZhaoBing, FU Ling, CHEN XinJing, ZHANG Tao, XIE YuHang, WANG HongBo, ZHU YuShuang. Quantitative Evaluation Method for Micro Heterogeneity of Tight Sandstone: A case study of Chang-6 reservoir of Yanchang Formation in Huaqing area, Ordos Basin[J]. Acta Sedimentologica Sinica, 2021, 39(5): 1086-1099. doi: 10.14027/j.issn.1000-0550.2021.037
Citation: CHEN ZhaoBing, FU Ling, CHEN XinJing, ZHANG Tao, XIE YuHang, WANG HongBo, ZHU YuShuang. Quantitative Evaluation Method for Micro Heterogeneity of Tight Sandstone: A case study of Chang-6 reservoir of Yanchang Formation in Huaqing area, Ordos Basin[J]. Acta Sedimentologica Sinica, 2021, 39(5): 1086-1099. doi: 10.14027/j.issn.1000-0550.2021.037

致密砂岩微观非均质性定量评价方法研究—以鄂尔多斯盆地华庆地区延长组长6油层组为例

doi: 10.14027/j.issn.1000-0550.2021.037
基金项目: 

国家自然科学基金 41802140

国家科技重大专项 2016ZX05050006

陕西省自然科学研究基础计划项目 2019JQ-257

刘宝珺地学科学基金 DMSM2019007

详细信息
    作者简介:

    陈朝兵,男,1984年出生,讲师,石油地质与储层评价,E-mail: zbchen@xsyu.edu.cn

    通讯作者:

    朱玉双,女,教授,油气田开发,E-mail: 642617295@qq.com

  • 中图分类号: TE122.1

Quantitative Evaluation Method for Micro Heterogeneity of Tight Sandstone: A case study of Chang-6 reservoir of Yanchang Formation in Huaqing area, Ordos Basin

Funds: 

National Natural Science Foundation of China 41802140

National Science and Technology Major Project 2016ZX05050006

Natural Science Research Basic Program Project of Shaanxi Province 2019JQ-257

Liu Baojun Earth Science Foundation DMSM2019007

  • 摘要: 鄂尔多斯盆地华庆地区长6油层组为深水重力流成因的致密砂岩储层,储层整体致密、非均质性强,微观评价难度较大,注水开发效果不理想。通过粒度分析、铸体薄片、高压压汞及真实砂岩两相渗流模型等常规实验手段,从岩石骨架颗粒、成岩作用及孔喉结构三方面优选非均质性特征参数,分别提出了表征骨架颗粒非均质性的粒度标准偏差σ,表征成岩作用非均质性的充填指数C和抗溶蚀指数K以及表征孔喉结构非均质性的分形维数D。在此基础上,论证并构建了微观非均质性综合指数F,F值介于0.2~3.0,F越大,代表了储层骨架颗粒分选越差、成岩作用越强、孔喉结构越复杂,微观非均质性越强,反之则相反。通过真实砂岩模型模拟油田注水开发,结果表明,F指数能够精确表征由孔隙致密区导致的储层微观特征差异,F指数与无水期及最终期驱油效率呈较好的相关性,可用于定量评价致密砂岩油藏内部微观非均质性,预测剩余油分布规律,为深水致密砂岩油藏的勘探开发提供科学依据。
  • 图  1  鄂尔多斯盆地构造及研究区位置(据You et al.[13]修改)

    Figure  1.  Location map of the Ordos Basin structure and study area(modified from You et al.[13]

    图  2  鄂尔多斯盆地华庆地区长6储层基本特征

    (a)储层岩石类型三角图;(b)填隙物及粒间孔,HQ213井,2 056.3 m

    Figure  2.  Basic characteristics of Chang 6 reservoir in Huaqing area, Ordos Basin

    (a) triangle chart of reservoir rock types; (b) interstitial material and intergranular pore, well HQ213, 2 056.3 m

    图  3  粒度标准偏差与粒间孔发育程度、面孔率的关系图

    (a)粒度标准偏差与粒间孔关系;(b)粒度标准偏差与面孔率关系

    Figure  3.  Relationship between standard deviation of particle size, degree of intergranular pore development, and face rate

    (a) relationship between standard deviation of particle size and intergranular pore; (b) relationship between standard deviation of particle size and face rate

    图  4  充填指数C及抗溶蚀指数K示意图

    Figure  4.  Schematic diagram of filling index C and corrosion resistance index K

    图  5  充填指数C、抗溶蚀指数K与孔隙度、渗透率的关系

    Figure  5.  Relationship between filling index C, corrosion resistance index K and porosity and permeability

    图  6  华庆地区长6储层毛管压力曲线特征

    Figure  6.  Capillary pressure curve for Chang 6 reservoir in Huaqing area

    图  7  分形维数与物性关系

    (a)分形维数与孔隙度关系;(b)分形维数与渗透率关系

    Figure  7.  Relationship between fractal dimension and physical properties

    (a) relationship between fractal dimension and porosity; (b) relationship between fractal dimension and permeability

    图  8  储层致密程度变化与分形维数关系

    Figure  8.  Relationship between changes in reservoir density and fractal dimensions

    图  9  微观非均质性综合评价指数F的构建模型

    Figure  9.  Construction model of comprehensive evaluation index F of micro heterogeneity

    图  10  华庆地区长6真实砂岩模型实验典型样品渗流特征

    (a)均匀状水驱渗流特征,HQ211⁃23井,2 187.1 m;(b)网状水驱渗流特征,HQ453井,1 873.5 m;(c)指状水驱渗流特征,HQ263⁃57井,2 101.0 m

    Figure  10.  Seepage characteristics of typical samples of Chang 6 in true sandstone model experiment in Huaqing area

    (a) uniform water drive seepage characteristics, well HQ211⁃23, 2 187.1 m; (b) network water drive seepage characteristics, well HQ453, 1 873.5 m; (c) finger water drive seepage characteristics, well HQ263⁃57, 2 101.0 m

    图  11  微观非均质综合指数F与真实砂岩模型驱油效率关系

    (a)F指数与无水期驱油效率关系;(b)F指数与最终期驱油效率关系

    Figure  11.  Relationship between microscopic heterogeneous composite index F and displacement efficiency of true sandstone model

    (a) relationship between F index and oil displacement efficiency in no water production period; (b)relationship between F index and final displacement efficiency

    表  1  华庆地区长6致密砂岩粒度标准偏差与粒间孔及面孔率(据铸体薄片数据)

    Table  1.   Standard deviation of particle size and intergranular pore and face rate of Chang 6 tight sandstone in Huaqing area(according to the data of casting sheet)

    井名 深度/m 标准偏差σ 粒间孔/% 面孔率/%
    HQ224 2 254.2 0.55 1.6 2.2
    HQ224 2 259.7 0.64 2.5 3.3
    HQ239 2 094.4 0.58 3.0 4.8
    HQ239 2 096.3 1.24 0.8 1.3
    HQ239 2 099.0 1.39 0.6 0.8
    HQ239 2 103.2 1.32 0.3 1.2
    HQ239 2 112.3 0.62 0.7 1.8
    HQ240 2 164.2 0.44 2.7 3.4
    HQ240 2 166.5 1.27 0.6 1.5
    HQ264 2 105.4 0.82 1.2 1.5
    HQ266 1 965.0 0.55 2.0 2.5
    HQ266 1 975.8 0.51 1.3 1.3
    HQ464 2 084.8 0.55 3.2 4.2
    HQ472 2 070.2 0.48 1.2 1.3
    HQ483 2 017.4 1.02 0.5 1.0
    HQ488 1 872.6 0.76 2.3 3.3
    HQ54 1 994.5 1.45 0.5 1.7
    HQ54 1 999.8 1.4 1.2 1.6
    HQ70 2 059.7 1.31 0.8 1.2
    HQ70 2 064.8 1.26 1.8 2.4
    HQ70 2 072.4 0.95 1.9 2.3
    HQ70 2 082.0 1.32 1.5 1.9
    HQ70 2 088.0 1.39 0.3 0.7
    HQ70 2 097.2 1.56 0.4 0.8
    HQ72 1 844.1 0.83 0.7 2.0
    HQ156 2 071.1 0.52 1.8 2.0
    HQ156 2 074.8 0.56 0.8 1.0
    下载: 导出CSV

    表  2  华庆地区长6致密砂岩孔隙演化定量计算方法

    Table  2.   Quantitative calculation method of pore evolution of Chang 6 tight sandstone in Huaqing area

    孔隙恢复 孔隙演化定量计算公式 来源
    原始孔隙度Φ 1 Φ 1=20.91+(22.9/S0) 式中,S0:Trask分选系数,Trask=(Q1/Q31/2 Q1:粒度概率累计曲线25%处的粒径大小 Q3:粒度概率累计曲线75%处的粒径大小 Beard et al.[19]
    压实后剩余孔隙度Φ 2 Φ 2=(粒间孔面孔率+胶结物溶孔面孔率)×物性分析孔隙度/总面孔率+胶结物含量 王瑞飞等[20] 郭正权等[21] 张涵冰[22]
    胶结交代后剩余孔隙度Φ 3 Φ 3=粒间孔面孔率×物性分析孔隙度/总面孔率
    溶蚀作用增加孔隙度Φ 4 Φ 4=溶蚀孔面孔率×物性分析孔隙度/总面孔率
    误差δ δ=(计算目前孔隙度-实测物性孔隙度)/实测物性孔隙度
    下载: 导出CSV

    表  3  华庆地区长6致密砂岩储层孔隙演化

    Table  3.   Pore evolution of Chang 6 tight sandstone reservoir in Huaqing area

    井号 原始孔隙度Φ 1/% 压实后剩余 孔隙度Φ 2/% 胶结交代 剩余孔隙度Φ 3/% 溶蚀作用 增加孔隙度Φ 4/% 计算现今孔隙度/% 实测现今孔度/% 误差δ/%
    HQ120 34.2 20.3 6.3 3.7 10.0 9.7 2.9
    HQ123 33.6 19.3 7.3 4.1 11.4 10.6 7.7
    HQ60 36.5 18.6 6.3 3.9 10.2 10.3 -1.4
    HQ38 37.2 17.7 5.2 4.0 9.2 9.7 -5.6
    HQ92 34.5 18.6 5.0 4.6 9.5 9.6 -0.8
    HQ95 33.2 16.5 8.6 2.2 10.8 10.6 1.1
    HQ41 37.6 14.3 7.3 3.1 10.4 9.6 8.8
    HQ142 34.1 16.8 6.8 3.3 10.1 10.4 -2.4
    HQ143 29.6 17.2 8.1 4.6 12.7 12.4 2.4
    HQ69 38.6 16.3 5.6 3.6 9.2 8.6 7.3
    HQ149 36.9 15.2 6.7 2.5 9.2 9.4 -1.8
    HQ12 37.8 18.6 6.3 1.9 8.3 8.7 -5.1
    HQ104 32.1 19.6 8.6 2.7 11.3 10.6 6.8
    HQ140 33.2 18.6 6.2 5.7 11.9 11.7 2.2
    HQ151 37.6 17.5 5.2 4.2 9.4 8.6 9.7
    HQ180 37.3 16.4 5.6 5.0 10.6 11.2 -5.7
    HQ233 31.2 18.0 7.9 5.1 13.0 12.7 2.4
    HQ264 35.3 19.0 6.9 4.2 11.1 10.9 1.9
    平均/范围 35.0 17.7 6.7 3.8 10.5 10.2 -5.7~9.7
    下载: 导出CSV

    表  4  华庆地区长6真实砂岩模型实验数据

    Table  4.   Statistics of experimental data of Chang 6 real true sandstone model in Huaqing area

    井号 井深/m 模型渗透率/(×10-3 μm2 微观非均质性综合指数F 无水期驱油效率/% 最终期驱油效率/% 水驱油渗流特征
    HQ211-33 2 176.4 0.137 1.69 20.2 44.8 网状
    HQ211-33 2 182.2 0.203 0.87 25.8 43.5 网状
    HQ211-33 2 187.1 0.325 0.81 25.4 45.1 网状
    HQ211-33 2 187.1 0.361 0.52 34.9 47.0 均匀状
    HQ264-59 2 112.5 0.220 0.96 25.6 51.1 均匀状
    HQ411 2 109.7 0.119 1.65 22.9 44.7 网状
    HQ411 2 123.5 0.465 0.82 28.9 46.5 均匀状
    HQ413 2 180.3 0.089 1.36 19.2 49.6 网状
    HQ413 2 182.1 0.146 2.27 21.1 42.5 网状
    HQ453 1 873.5 0.095 1.33 19.7 40.7 网状
    HQ453 1 874.8 0.452 0.88 25.0 46.5 网状
    HQ156 2 059.9 0.450 1.16 21.4 42.0 均匀状
    HQ262-56 2 178.9 0.225 1.65 18.8 38.2 网状
    HQ263-57 2 094.1 0.110 2.75 17.0 34.4 指状
    HQ263-57 2 101.0 0.118 2.84 20.1 29.1 指状
    下载: 导出CSV
  • [1] 李滔,李闽,荆雪琪,等. 孔隙尺度各向异性与孔隙分布非均质性对多孔介质渗透率的影响机理[J]. 石油勘探与开发,2019,46(3):569-579.

    Li Tao, Li Min, Jing Xueqi, et al. Influence mechanism of pore-scale anisotropy and pore distribution heterogeneity on permeability of porous media[J]. Petroleum Exploration and Development, 2019, 46(3): 569-579.
    [2] 佩蒂庄 F J,波特 P E,西弗 R. 砂和砂岩[M]. 北京:科学出版社,1977.

    Pettijohn F J, Potter P E, Siever R. Sand and sandstone[M]. Beijing: Science Press, 1977.
    [3] Haldorsen H H. Reservoir characterisation procedures for numerical simulation[D]. University of Texas at Austin. PhD thesis, 1983
    [4] 裘怿楠,陈子琪. 油藏描述[M]. 北京:石油工业出版,1996:326-344.

    Qiu Yinan, Chen Ziqi. Reservoir description[M]. Beijing: Petroleum Industry Press, 1996: 326-344.
    [5] 贾承造,邹才能,李建忠,等. 中国致密油评价标准、主要类型、基本特征及资源前景[J]. 石油学报,2012,33(3):343-350.

    Jia Chengzao, Zou Caineng, Li Jianzhong, et al. Assessment criteria, main types, basic features and resource prospects of the tight oil in China[J]. Acta Petrolei Sinica, 2012, 33(3): 343-350.
    [6] 赵文光,蔡忠贤,韩中文. 应用定量方法描述储层孔隙结构的非均质性[J]. 新疆石油天然气,2006,2(3):22-25.

    Zhao Wenguang, Cai Zhongxian, Han Zhongwen. The application of a quantitative method to describe the heterogeneity of pore structures of reservoir[J]. Xinjiang Oil & Gas, 2006, 2(3): 22-25.
    [7] 白斌,朱如凯,吴松涛,等. 利用多尺度CT成像表征致密砂岩微观孔喉结构[J]. 石油勘探与开发,2013,40(3):329-333.

    Bai Bin, Zhu Rukai, Wu Songtao, et al. Multi-scale method of Nano(Micro)-CT study on microscopic pore structure of tight sandstone of Yanchang Formation, Ordos Basin[J]. Petroleum Exploration and Development, 2013, 40(3): 329-333.
    [8] 王明磊,张遂安,张福东,等. 鄂尔多斯盆地延长组长7段致密油微观赋存形式定量研究[J]. 石油勘探与开发,2015,42(6):757-762.

    Wang Minglei, Zhang Sui'an, Zhang Fudong, et al. Quantitative research on tight oil microscopic state of Chang 7 member of Triassic Yanchang Formation in Ordos Basin, NW China[J]. Petroleum Exploration and Development, 2015, 42(6): 757-762.
    [9] 孙英峰,赵毅鑫,王欣,等. 基于同步辐射装置定量表征煤孔隙结构非均质性和各向异性[J]. 石油勘探与开发,2019,46(6):1128-1137.

    Sun Yingfeng, Zhao Yixin, Wang Xin, et al. Synchrotron radiation facility-based quantitative evaluation of pore structure heterogeneity and anisotropy in coal[J]. Petroleum Exploration and Development, 2019, 46(6): 1128-1137.
    [10] 全洪慧,朱玉双,张洪军,等. 储层孔隙结构与水驱油微观渗流特征:以安塞油田王窑区长6油层组为例[J]. 石油与天然气地质,2011,32(6):952-960.

    Quan Honghui, Zhu Yushuang, Zhang Hongjun, et al. Reservoir pore structure and micro-flow characteristics of waterflooding: A case study from Chang-6 reservoir of Wangyao block in Ansai oilfield[J]. Oil & Gas Geology, 2011, 32(6): 952-960.
    [11] 杨华,梁晓伟,牛小兵,等. 陆相致密油形成地质条件及富集主控因素:以鄂尔多斯盆地三叠系延长组长7段为例[J]. 石油勘探与开发,2017,44(1):12-20.

    Yang Hua, Liang Xiaowei, Niu Xiaobing, et al. Geological conditions for continental tight oil formation and the main controlling factors for the enrichment: A case of Chang 7 member, Triassic Yanchang Formation, Ordos Basin, NW China[J]. Petroleum Exploration and Development, 2017, 44(1): 12-20.
    [12] 张国伟,孟庆任,于在平,等. 秦岭造山带的造山过程及其动力学特征[J]. 中国科学(D辑):地球科学,1996,26(3):193-200.

    Zhang Guowei, Meng Qingren, Yu Zaiping, et al. The orogenic process and dynamic characteristics of the Qinling orogenic belt[J]. Science China (Seri. D): Earth Sciences, 1996, 26(3): 193-200.
    [13] You J Y, Liu Y Q, Zhou D W,et al. Activity of hydrothermal fluid at the bottom of a lake and its influence on the development of high-quality source rocks: Triassic Yanchang Formation,southern Ordos Basin, China[J]. Australian Journal of Earth Sciences, 2020, 67(1): 115-128.
    [14] 操应长,杨田,王艳忠,等. 超临界沉积物重力流形成演化及特征[J]. 石油学报,2017,38(6):607-621.

    Cao Yingchang, Yang Tian, Wang Yanzhong, et al. Formation, evolution and sedimentary characteristics of supercritical sediment gravity-flow[J]. Acta Petrolei Sinica, 2017, 38(6): 607-621.
    [15] 曹江骏,杨友运,陈朝兵,等. 致密砂岩储层骨架砂体构型特征:以鄂尔多斯盆地合水地区延长组长6段砂体为例[J]. 沉积学报,2019,37(6):1105-1116.

    Cao Jiangjun, Yang Youyun, Chen Zhaobing, et al. Analysis of configuration characteristics for skeleton sand body with tight sandstone reservoir: A case study of Triassic Chang 6 members in Heshui area, Ordos Basin, NW China[J]. Acta Sedimentologica Sinica, 2019, 37(6): 1105-1116.
    [16] 孙军,路月,杨慧良,等. 渤海海峡BHS01孔海陆交互相沉积物粒度特征及沉积环境分析[J]. 海洋地质前沿,2019,35(8):11-19.

    Sun Jun, Lu Yue, Yang Huiliang, et al. Grain-size distribution patterns of the core sediments from BHS01 of Bohai strait and their environmental implications[J]. Marine Geology Frontiers, 2019, 35(8): 11-19.
    [17] 袁红旗,王蕾,于英华,等. 沉积学粒度分析方法综述[J]. 吉林大学学报(地球科学版),2019,49(2):380-393.

    Yuan Hongqi, Wang Lei, Yu Yinghua, et al. Review of sedimentary grain size analysis methods[J]. Journal of Jilin University (Earth Science Edition), 2019, 49(2): 380-393.
    [18] 朱筱敏. 沉积岩石学[M]. 4版. 北京:石油工业出版社,2008.

    Zhu Xiaomin. Sedimentary petrology[M]. 4th ed. Beijing: Petroleum Industry Press, 2008.
    [19] Beard D C, Weyl P K. Influence of texture on porosity and permeabil-ity of unconsolidated sand[J]. AAPG Bulletin,1973,57( 2) :349-369
    [20] 王瑞飞,陈明强. 储层沉积—成岩过程中孔隙度参数演化的定量分析:以鄂尔多斯盆地沿25区块、庄40区块为例[J]. 地质学报,2007,81(10):1432-1438.

    Wang Ruifei, Chen Mingqiang. Quantitative analysis of porosity evolution during the reservoir sedimentation-diagenesis: Taking the Yan 25 and Zhuang 40 areas in the Ordos Basín as examples[J]. Acta Geologica Sinica, 2007, 81(10): 1432-1438.
    [21] 郭正权,齐亚林,楚美娟,等. 鄂尔多斯盆地上三叠统延长组储层致密史恢复[J]. 石油实验地质,2012,34(6):594-598,603.

    Guo Zhengquan, Qi Yalin, Chu Meijuan, et al. Recovery of compact history of Yanchang reservoir in Upper Triassic, Ordos Basin[J]. Petroleum Geology & Experiment, 2012, 34(6): 594-598, 603.
    [22] 张涵冰. 深埋藏致密砂岩储层成因机理及评价:以塔中地区志留系柯坪塔格组为例[D]. 成都:西南石油大学,2017.

    Zhang Hanbing. Genetic mechanism for deep buried tight sandstone reservoir and its evalution: A case study of Silurian Kepingtage Formation in Tazhong area[D]. Chengdu: Southwest Petroleum University, 2017.
    [23] 何涛,王芳,汪伶俐. 致密砂岩储层微观孔隙结构特征:以鄂尔多斯盆地延长组长7储层为例[J]. 岩性油气藏,2013,25(4):23-26.

    He Tao, Wang Fang, Wang Lingli. Characteristics of micro-pore structure of tight sandstone reservoir: A case study from Chang 7 reservoir of Yanchang Formation in Ordos Basin[J]. Lithologic Reservoirs, 2013, 25(4): 23-26.
    [24] 赖锦,王贵文,郑懿琼,等. 低渗透碎屑岩储层孔隙结构分形维数计算方法:以川中地区须家河组储层41块岩样为例[J]. 东北石油大学学报,2013,37(1):1-7.

    Lai Jin, Wang Guiwen, Zheng Yiqiong, et al. Method for calculating the fractal dimension of the pore structure of low permeability reservoirs. A case study on the Xujiahe Formation reservoir in central Sichuan Basin[J]. Journal of Northeast Petroleum University, 2013, 37(1): 1-7.
    [25] 马新仿,张士诚,郎兆新. 孔隙结构特征参数的分形表征[J]. 油气地质与采收率,2005,12(6):34-36.

    Ma Xinfang, Zhang Shicheng, Lang Zhaoxin. Fractal characterization of characteristic parameters of pore structure[J]. Petroleum Geology and Recovery Efficiency, 2005, 12(6): 34-36.
    [26] Mandelbrot B B. The fractal geometry of nature[M]. San Franeiseo: W. H. Freeman, 1982.
    [27] Thompson A H, Katz A J, Krohn C E. The micro-geometry and transport properties of sandstones[J]. Advances in physic, 1987, 36(5), 6925-6941.
  • [1] 王乾右, 杨威, 葛云锦, 宋岩, 姜振学, 罗群, 左如斯, 李耀华, 刘聃, 张帆, 王耀华, 鲁健康.  不同沉积微相致密储层的成岩响应及其控储机理 . 沉积学报, 2021, 39(4): 841-862. doi: 10.14027/j.issn.1000-0550.2020.112
    [2] 曹江骏, 陈朝兵, 程皇辉, 朱玉杰, 罗静兰, 王茜, 马迪娜·马吾提汗null.  成岩作用对深水致密砂岩储层微观非均质性的影响 . 沉积学报, 2021, 39(4): 1031-1046. doi: 10.14027/j.issn.1000-0550.2020.040
    [3] 孙海涛, 钟大康, 王威, 王爱, 杨烁, 杜红权, 唐自成, 周志恒.  四川盆地马路背地区上三叠统须家河组致密砂岩储层成因分析 . 沉积学报, 2021, 39(5): 1057-1067. doi: 10.14027/j.issn.1000-0550.2020.121
    [4] 石太衡, 吴耿, 李荣.  硅质热泉沉积物的结构和微生物成矿成岩机制 . 沉积学报, 2020, 38(1): 113-123. doi: 10.14027/j.issn.1000-0550.2019.023
    [5] 史兵兵, 常象春, 尹伟, 毛礼鑫.  鄂尔多斯盆地镇泾地区长8致密储层成藏期临界物性厘定 . 沉积学报, 2020, 38(1): 231-243. doi: 10.14027/j.issn.1000-0550.2019.012
    [6] 陈贺贺, 朱筱敏, 陈纯芳, 尹伟, 施瑞生.  致密砂岩储层致密化与成藏史耦合关系研究——以鄂尔多斯南部镇原—泾川地区延长组长8油层组为例 . 沉积学报, 2018, 36(2): 401-414. doi: 10.14027/j.issn.1000-0550.2018.029
    [7] 曾庆鲁, 张荣虎, 王力宝, 赵继龙, 佘敏.  库车坳陷白垩系深层致密砂岩储层溶蚀作用实验模拟研究 . 沉积学报, 2018, 36(5): 946-956. doi: 10.14027/j.issn.1000-0550.2018.079
    [8] 黄少华, 秦明宽, 刘章月, 许强, 郭强.  成岩作用与油气侵位对砂岩铀成矿的约束——以准南硫磺沟地区头屯河组砂岩铀成矿为例 . 沉积学报, 2016, 34(2): 250-259. doi: 10.14027/j.cnki.cjxb.2016.02.004
    [9] 郑庆华, 柳益群.  鄂尔多斯盆地镇北地区延长组长4+5致密油层成岩作用及成岩相 . 沉积学报, 2015, 33(5): 1000-1012. doi: 10.14027/j.cnki.cjxb.2015.05.016
    [10] 杨友运, 何康宁, 任颖惠.  华庆地区长63储层内部建筑结构模式分析 . 沉积学报, 2015, 33(2): 357-363. doi: 10.14027/j.cnki.cjxb.2015.02.015
    [11]  四川盆地蜀南地区上三叠统须家河组低孔低渗储层特征及形成机理 . 沉积学报, 2013, 31(1): 167-175.
    [12] 刘占国.  四川盆地川中地区中下侏罗统砂岩储层异常致密成因机理 . 沉积学报, 2011, 29(4): 744-751.
    [13] 罗 忠.  层序界面对砂岩成岩作用及储层质量的影响——以鄂尔多斯盆地延河露头上三叠统延长组为例 . 沉积学报, 2007, 25(6): 903-914.
    [14] 漆滨汶.  山东省牛庄洼陷古近系沙河街组沙三中亚段储集层成岩作用研究 . 沉积学报, 2007, 25(1): 99-109.
    [15] 刘林玉.  鄂尔多斯盆地镇北地区长3储层微观非均质性的实验分析 . 沉积学报, 2007, 25(2): 224-229.
    [16] 孟元林.  渤海湾盆地西部凹陷南段成岩相分析与优质储层预测 . 沉积学报, 2006, 24(2): 185-192.
    [17] 刘林玉, 王震亮, 张龙.  鄂尔多斯盆地镇北地区长3砂岩的成岩作用及其对储层的影响 . 沉积学报, 2006, (5): 690-697.
    [18] 鄂尔多斯盆地镇北地区长3砂岩的成岩作用及其对储层的影响 . 沉积学报, 2006, 24(05): 690-697.
    [19] 王琪, 禚喜准, 陈国俊, 史基安, 王多云.  鄂尔多斯盆地盐池—姬源地区三叠系长4+5砂岩成岩演化特征与优质储层分布 . 沉积学报, 2005, 23(3): 397-405.
    [20] 刘宝珺, 余光明, 王成善, 张哨楠.  西藏北部西雅尔岗地区白垩—第三系红层的成岩作用及成岩环境 . 沉积学报, 1989, 7(1): 11-28.
  • 加载中
图(11) / 表 (4)
计量
  • 文章访问数:  439
  • HTML全文浏览量:  102
  • PDF下载量:  152
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-12-31
  • 修回日期:  2021-04-14
  • 刊出日期:  2021-10-10

目录

    致密砂岩微观非均质性定量评价方法研究

    doi: 10.14027/j.issn.1000-0550.2021.037
      基金项目:

      国家自然科学基金 41802140

      国家科技重大专项 2016ZX05050006

      陕西省自然科学研究基础计划项目 2019JQ-257

      刘宝珺地学科学基金 DMSM2019007

      作者简介:

      陈朝兵,男,1984年出生,讲师,石油地质与储层评价,E-mail: zbchen@xsyu.edu.cn

      通讯作者: 朱玉双,女,教授,油气田开发,E-mail: 642617295@qq.com
    • 中图分类号: TE122.1

    摘要: 鄂尔多斯盆地华庆地区长6油层组为深水重力流成因的致密砂岩储层,储层整体致密、非均质性强,微观评价难度较大,注水开发效果不理想。通过粒度分析、铸体薄片、高压压汞及真实砂岩两相渗流模型等常规实验手段,从岩石骨架颗粒、成岩作用及孔喉结构三方面优选非均质性特征参数,分别提出了表征骨架颗粒非均质性的粒度标准偏差σ,表征成岩作用非均质性的充填指数C和抗溶蚀指数K以及表征孔喉结构非均质性的分形维数D。在此基础上,论证并构建了微观非均质性综合指数F,F值介于0.2~3.0,F越大,代表了储层骨架颗粒分选越差、成岩作用越强、孔喉结构越复杂,微观非均质性越强,反之则相反。通过真实砂岩模型模拟油田注水开发,结果表明,F指数能够精确表征由孔隙致密区导致的储层微观特征差异,F指数与无水期及最终期驱油效率呈较好的相关性,可用于定量评价致密砂岩油藏内部微观非均质性,预测剩余油分布规律,为深水致密砂岩油藏的勘探开发提供科学依据。

    English Abstract

    陈朝兵, 付玲, 陈新晶, 张涛, 解宇航, 王泓波, 朱玉双. 致密砂岩微观非均质性定量评价方法研究—以鄂尔多斯盆地华庆地区延长组长6油层组为例[J]. 沉积学报, 2021, 39(5): 1086-1099. doi: 10.14027/j.issn.1000-0550.2021.037
    引用本文: 陈朝兵, 付玲, 陈新晶, 张涛, 解宇航, 王泓波, 朱玉双. 致密砂岩微观非均质性定量评价方法研究—以鄂尔多斯盆地华庆地区延长组长6油层组为例[J]. 沉积学报, 2021, 39(5): 1086-1099. doi: 10.14027/j.issn.1000-0550.2021.037
    CHEN ZhaoBing, FU Ling, CHEN XinJing, ZHANG Tao, XIE YuHang, WANG HongBo, ZHU YuShuang. Quantitative Evaluation Method for Micro Heterogeneity of Tight Sandstone: A case study of Chang-6 reservoir of Yanchang Formation in Huaqing area, Ordos Basin[J]. Acta Sedimentologica Sinica, 2021, 39(5): 1086-1099. doi: 10.14027/j.issn.1000-0550.2021.037
    Citation: CHEN ZhaoBing, FU Ling, CHEN XinJing, ZHANG Tao, XIE YuHang, WANG HongBo, ZHU YuShuang. Quantitative Evaluation Method for Micro Heterogeneity of Tight Sandstone: A case study of Chang-6 reservoir of Yanchang Formation in Huaqing area, Ordos Basin[J]. Acta Sedimentologica Sinica, 2021, 39(5): 1086-1099. doi: 10.14027/j.issn.1000-0550.2021.037
      • 储层非均质性指储层在形成过程中,经历了沉积、成岩以及后期构造作用的综合影响,使得储层空间分布及基本性质(岩性、物性、电性及含油气性)在三维空间上存在的不均匀变化[1]。学者关于储层非均质性的研究由来已久,佩蒂庄等[2]以河流沉积储层为例,考虑储层非均质性规模,将储层非均质分为5级,分别为层系规模(100 m级)、砂体规模(10 m级)、层理规模(1~10 m级)、纹层规模(10~100 mm级)、孔隙规模(10~100 μm级);Haldorsen[3]从非均质性的尺度考虑,将非均质性分为4种类型,分别为孔隙和颗粒规模的微观非均质性、岩心规模的宏观非均质性、数值模拟模型中的网块状大型非均质性和整个岩层或区域规模的巨型非均质性;裘怿楠等[4]根据中国陆相储层特征并结合油田开发生产实用性,将碎屑岩储层非均质性按照规模大小分为层内、层间、平面和微观非均质性4个级别,该分类方案也是中国各大油田技术工作者广泛采用的非均质性分类标准,其中层内、层间、平面非均质性属于宏观非均质性范畴,对油田开发过程中的注水波及系数影响较大,注水波及区域内的原油动用程度则与微观孔喉结构的非均质性密不可分,这涉及到孔隙和喉道的大小、分布特征、孔喉配位数、连通性等,即使宏观表现较为均质的储集层也不可避免地存在一定程度的微观非均质性。

        随着常规油气资源的逐渐枯竭,致密砂岩油气资源的占比越来越高。常规储层颗粒粒度粗、物性条件好、孔喉较粗、连通性好,微观非均质性较弱。与常规储层不同,致密砂岩储层物性较差,孔隙度一般小于10%,空气渗透率小于1×10-3 μm2(覆压基质渗透率小于0.1×10-3 μm2[5],颗粒粒度细、填隙物含量高、成岩强度大,微观非均质性强,其对储层内的油气聚集、渗流特征及剩余油分布的影响更大[6]。致密砂岩储层微观非均质性研究的难度要远大于常规储层,研究对象多为微米—纳米级的储集空间,研究手段涉及不同尺度的实验方法,主要包括间接测量的低温气体吸附法、核磁共振、高压压汞及恒速压汞法和直接观测的场发射扫描电镜法、CT成像扫描法等[7-10]。尽管测试方法较多,但定量评价手段较少,且各种尺度实验数据相对独立,缺乏有机结合,暂未建立起一个相对综合的微观非均质性评价指数或标准。另外,部分实验测试费用高、测试周期长,如微纳米CT、恒速压汞等,导致定量化测试手段不能快速、广泛的运用于油田科研生产,因此亟需一种资料获取容易、快捷简便的微观非均质性定量评价方法。

        基于此,本文通过粒度分析、铸体薄片、高压压汞及真实砂岩两相渗流模型等实验手段,对鄂尔多斯盆地华庆地区长6致密砂岩储层微观非均质性展开研究,分别从岩石骨架颗粒非均质性、成岩作用非均质性及孔喉结构非均质性三方面展开研究,优选影响微观非均质性的核心参数,构建微观非均质性综合指数F,并探讨微观非均质性与渗流特征之间的关系,为评价致密砂岩微观特征及致密油开发效果提供依据。

      • 鄂尔多斯盆地地处中国中部,为稳定的大型多旋回克拉通盆地,现今构造形态为一平缓的西倾单斜,坡度通常小于1°。盆地内可划分为六个二级构造单元,分别为北部的伊蒙隆起、西部的天环凹陷和西缘逆冲推覆带、东部的晋西绕折带、南部的渭北隆起带及中部的伊陕斜坡[11]。晚三叠世延长期,鄂尔多斯盆地发育大型内陆河流—三角洲—湖泊沉积体系,主要物源供给方向为盆地北部阴山古陆和西南部秦祁造山带[12],延长组地层总厚度约1 300 m,自下向上划分为长10—长1十个油层组[11]。长6期,研究区处于盆地西南部及东北部物源的交汇处[13]图1),受盆地周边地震、火山及湖底暗流等触发因素的影响,研究区发育重力流沉积,西南部坡折带较陡,主要发育浊流及少量砂质碎屑流砂体[14];东北部坡折带较缓,发育砂质碎屑流及少量浊流、泥质碎屑流沉积体[15]。长6油层组厚度约80~110 m,自下向上可划分为长63、长62、长61三个小层,主力含油层为长63小层,不同重力流砂体相互叠置,储层整体非均质性强,开发效果不理想。

        图  1  鄂尔多斯盆地构造及研究区位置(据You et al.[13]修改)

        Figure 1.  Location map of the Ordos Basin structure and study area(modified from You et al.[13]

        华庆地区长6储层砂岩以灰色细粒岩屑长石砂岩和长石岩屑砂岩为主,其次为少量长石砂岩(图2a),普遍含泥质杂基,碎屑具有“低石英,高长石”的特点,颗粒分选中等,磨圆度中等—好,呈次棱角状。砂岩整体成分成熟度较低,结构成熟度中等。胶结物以伊利石、绿泥石、碳酸盐为主,含少量硅质、长石质。孔隙类型以粒间孔(1.59%)为主(图2b),其次为长石溶孔(0.49%)和岩屑溶孔(0.15%),含少量晶间孔及微裂隙(0.17%),平均面孔率2.4%。孔喉类型以小孔—微喉型为主,并有少量的小孔—微细喉型,平均孔隙度10.2%,平均渗透率0.35×10-3 μm2,储层整体致密,成岩作用较强,为典型的致密砂岩储层。

        图  2  鄂尔多斯盆地华庆地区长6储层基本特征

        Figure 2.  Basic characteristics of Chang 6 reservoir in Huaqing area, Ordos Basin

        与常规牵引流致密砂岩不同,华庆地区长6深水重力流致密砂岩的碎屑颗粒粒度分选及成岩作用差异更为突出。由于深水环境水动力十分复杂,重力流砂体的规模、分布受到周围火山、地震、波浪等外部触发因素的影响[13],因此不同期次沉积的重力流砂体内的原始机械杂基及碎屑云母含量差别较大,其中长6砂岩内的杂基含量介于0.8%~18.6%,平均含量为5.9%。碎屑云母含量介于0.1%~21.3%,平均含量为7.4%,导致深水致密砂岩内的泥质、云母等塑性组分明显多于常规牵引流致密砂岩,这势必影响到沉积物的原始粒度分选及成岩压实、胶结等作用的强弱。因此,选取华庆地区长6深水重力流致密砂岩作为微观非均质性评价的研究对象更加具有代表性。

      • 粒度是沉积物颗粒大小的表征,其组成与分布受控于物源供给、搬运条件、沉积动力及地形地貌等因素,包含了丰富的物源、水动力条件、沉积环境等环境变化信息[16]。因此,通过对沉积物粒度的分析,可以有效的追溯由物源、水动力条件及沉积环境差异造成的岩石骨架颗粒的非均质性。

        沉积学中引入粒度分析的方法已超过一个世纪。将粒度分析仪所测粒度数据成图,常见的有直方图、频率曲线图、累积曲线图,用于辅助分析沉积环境并获得粒度参数,其中累积曲线是将不同粒级颗粒质量累加的百分比图,可直接从曲线上获得某些累积百分比处的颗粒直径,进而计算粒度参数,是目前运用最广的粒度辅助图形[17-18]。常见的粒度参数有平均粒径(Mz)、偏度(Sk)和峰度(Kg)、标准偏差σ(或分选系数),各项参数均有各自不同的意义,其中平均粒径(Mz)反映粒度分布的集中趋势,表示沉积物粒度的粗细情况;偏度(Sk)用来衡量频率曲线的不对称程度,即非正态特征情况;峰度(Kg)是指分布曲线与正态分布曲线相比时,其峰的宽窄尖锐程度,是衡量频率曲线尖峰凸起程度的参数;标准偏差σ是表示沉积物分选程度的参数,包括颗粒大小的均匀程度或沉积物围绕集中趋势的离差,较好的分选性对应较小的分选系数,当分选系数为0时,沉积物颗粒是绝对均匀的[17]。因此,标准偏差σ是表征岩石颗粒非均质性的重要参数,计算分选系数的公式不止一种,本次研究选用福克和沃德提出的标准偏差公式[18]

        σ = Φ 84 - Φ 16 4 + Φ 95 - Φ 5 6.6 (1)

        式中: Φ 16~ Φ 84代表粒级分布的中央部分, Φ 95 Φ 5代表对水动力条件反映最灵敏的粗、细尾部的分选情况。

        通过对华庆地区13口井/27块长6砂岩铸体薄片数据分析统计(表1),研究区致密砂岩粒度标准偏差σ介于0.44~1.56,平均0.94;粒度标准偏差σ与粒间孔呈现较强相关性(相关系数R为0.652 6),与面孔率也呈现一定相关性(相关系数R为0.551 9),但相关性弱于前者(图3),表明原始碎屑颗粒的分选程度对现今致密砂岩残余粒间孔产生一定影响,尽管后期压实、胶结等成岩作用破坏了大部分原始粒间孔隙,但分选性好的碎屑颗粒所搭建的岩石骨架更易保存粒间孔隙;研究区长6储层普遍发育长石溶孔,导致薄片下的整体面孔率与粒度标准偏差σ的相关性弱于粒间孔。因此,粒度标准偏差σ可用于表征致密砂岩原始碎屑颗粒的非均质性,它代表了原始沉积环境对致密砂岩微观特征的影响强弱。

        表 1  华庆地区长6致密砂岩粒度标准偏差与粒间孔及面孔率(据铸体薄片数据)

        Table 1.  Standard deviation of particle size and intergranular pore and face rate of Chang 6 tight sandstone in Huaqing area(according to the data of casting sheet)

        井名 深度/m 标准偏差σ 粒间孔/% 面孔率/%
        HQ224 2 254.2 0.55 1.6 2.2
        HQ224 2 259.7 0.64 2.5 3.3
        HQ239 2 094.4 0.58 3.0 4.8
        HQ239 2 096.3 1.24 0.8 1.3
        HQ239 2 099.0 1.39 0.6 0.8
        HQ239 2 103.2 1.32 0.3 1.2
        HQ239 2 112.3 0.62 0.7 1.8
        HQ240 2 164.2 0.44 2.7 3.4
        HQ240 2 166.5 1.27 0.6 1.5
        HQ264 2 105.4 0.82 1.2 1.5
        HQ266 1 965.0 0.55 2.0 2.5
        HQ266 1 975.8 0.51 1.3 1.3
        HQ464 2 084.8 0.55 3.2 4.2
        HQ472 2 070.2 0.48 1.2 1.3
        HQ483 2 017.4 1.02 0.5 1.0
        HQ488 1 872.6 0.76 2.3 3.3
        HQ54 1 994.5 1.45 0.5 1.7
        HQ54 1 999.8 1.4 1.2 1.6
        HQ70 2 059.7 1.31 0.8 1.2
        HQ70 2 064.8 1.26 1.8 2.4
        HQ70 2 072.4 0.95 1.9 2.3
        HQ70 2 082.0 1.32 1.5 1.9
        HQ70 2 088.0 1.39 0.3 0.7
        HQ70 2 097.2 1.56 0.4 0.8
        HQ72 1 844.1 0.83 0.7 2.0
        HQ156 2 071.1 0.52 1.8 2.0
        HQ156 2 074.8 0.56 0.8 1.0

        图  3  粒度标准偏差与粒间孔发育程度、面孔率的关系图

        Figure 3.  Relationship between standard deviation of particle size, degree of intergranular pore development, and face rate

      • 成岩作用是砂岩储层致密化的重要因素[19-20],沉积物沉积后经历了漫长的深埋成岩过程,包括上覆地层的重力压实作用、各类自生矿物的充填胶结作用、不稳定矿物的蚀变溶解作用以及不同矿物之间的交代作用,是影响微观非均质性的主要因素之一。

      • 成岩作用是造成储层孔隙丧失与再生的重要因素,目前关于碎屑岩储层孔隙演化的定量研究已相对成熟,根据前人对不同成岩作用残余孔隙恢复的理论计算公式[20-22]表2),对研究区长6致密砂岩储层的孔隙演化进行了定量计算(表3),结果表明,华庆地区长6致密砂岩储层原始孔隙度(Φ 1)为35.0%,压实作用后剩余孔隙度(Φ 2)为17.7%,压实损失率为49.5%;经历胶结、交代作用后剩余孔隙度(Φ 3)仅为6.7%,胶结交代损失率为31.5%;溶蚀作用增加的次生孔隙度(Φ 4)为3.8%,溶蚀增孔率为10.8%;计算现今孔隙度为10.5%,实测现今孔隙度为10.2%,误差范围介于-5.7%~9.7%。因此,压实作用和胶结交代作用是导致孔隙度丧失的主要原因,溶蚀作用在一定程度上改善了储层的储集性能。

        表 2  华庆地区长6致密砂岩孔隙演化定量计算方法

        Table 2.  Quantitative calculation method of pore evolution of Chang 6 tight sandstone in Huaqing area

        孔隙恢复 孔隙演化定量计算公式 来源
        原始孔隙度Φ 1 Φ 1=20.91+(22.9/S0) 式中,S0:Trask分选系数,Trask=(Q1/Q31/2 Q1:粒度概率累计曲线25%处的粒径大小 Q3:粒度概率累计曲线75%处的粒径大小 Beard et al.[19]
        压实后剩余孔隙度Φ 2 Φ 2=(粒间孔面孔率+胶结物溶孔面孔率)×物性分析孔隙度/总面孔率+胶结物含量 王瑞飞等[20] 郭正权等[21] 张涵冰[22]
        胶结交代后剩余孔隙度Φ 3 Φ 3=粒间孔面孔率×物性分析孔隙度/总面孔率
        溶蚀作用增加孔隙度Φ 4 Φ 4=溶蚀孔面孔率×物性分析孔隙度/总面孔率
        误差δ δ=(计算目前孔隙度-实测物性孔隙度)/实测物性孔隙度

        表 3  华庆地区长6致密砂岩储层孔隙演化

        Table 3.  Pore evolution of Chang 6 tight sandstone reservoir in Huaqing area

        井号 原始孔隙度Φ 1/% 压实后剩余 孔隙度Φ 2/% 胶结交代 剩余孔隙度Φ 3/% 溶蚀作用 增加孔隙度Φ 4/% 计算现今孔隙度/% 实测现今孔度/% 误差δ/%
        HQ120 34.2 20.3 6.3 3.7 10.0 9.7 2.9
        HQ123 33.6 19.3 7.3 4.1 11.4 10.6 7.7
        HQ60 36.5 18.6 6.3 3.9 10.2 10.3 -1.4
        HQ38 37.2 17.7 5.2 4.0 9.2 9.7 -5.6
        HQ92 34.5 18.6 5.0 4.6 9.5 9.6 -0.8
        HQ95 33.2 16.5 8.6 2.2 10.8 10.6 1.1
        HQ41 37.6 14.3 7.3 3.1 10.4 9.6 8.8
        HQ142 34.1 16.8 6.8 3.3 10.1 10.4 -2.4
        HQ143 29.6 17.2 8.1 4.6 12.7 12.4 2.4
        HQ69 38.6 16.3 5.6 3.6 9.2 8.6 7.3
        HQ149 36.9 15.2 6.7 2.5 9.2 9.4 -1.8
        HQ12 37.8 18.6 6.3 1.9 8.3 8.7 -5.1
        HQ104 32.1 19.6 8.6 2.7 11.3 10.6 6.8
        HQ140 33.2 18.6 6.2 5.7 11.9 11.7 2.2
        HQ151 37.6 17.5 5.2 4.2 9.4 8.6 9.7
        HQ180 37.3 16.4 5.6 5.0 10.6 11.2 -5.7
        HQ233 31.2 18.0 7.9 5.1 13.0 12.7 2.4
        HQ264 35.3 19.0 6.9 4.2 11.1 10.9 1.9
        平均/范围 35.0 17.7 6.7 3.8 10.5 10.2 -5.7~9.7
      • 由于压实作用对于致密砂岩储层的影响相对复杂,不仅与上覆地层重力压实相关,还受到同期孔隙内流体异常压力或所处区域构造应力的影响[20-21]。同时考虑到华庆地区长6致密砂岩储层埋深稳定,无大型断裂构造控制,镜下观察到的岩石骨架颗粒均以点—线接触为主,因此本次研究将压实作用对不同样品孔隙度的破坏程度视为同等。基于此,本次研究重点考虑填隙物充填胶结作用及溶蚀作用对于致密砂岩微观非均质性的影响(交代作用为矿物成分的置换,对孔隙度影响不大,不作考虑)。

        以油田常见的的铸体薄片资料入手,提出以填隙物充填指数C(简称充填指数C)和易溶颗粒抗溶蚀指数K(简称抗溶蚀指数K)两个参数来评价成岩作用的非均质性(图4),其中填隙物充填指数C=填隙物含量/(填隙物含量+剩余粒间孔含量),填隙物包括杂基和胶结物,因此填隙物充填指数C代表了原始沉积机械杂基及成岩胶结物对孔隙的充填程度,数值介于0~1,数值越大,充填越致密,剩余粒间孔越不发育;由于溶蚀指数(溶蚀指数=易溶颗粒溶孔含量/(易溶颗粒含量+易溶颗粒溶孔含量))数值通常很小,不方便使用,因此采用“抗溶蚀指数”的概念,抗溶蚀指数K=易溶颗粒含量/(易溶颗粒含量+易溶颗粒溶孔含量),研究区易溶颗粒主要为长石和岩屑,因此抗溶蚀指数代表的是长石和岩屑的溶蚀程度,该数值介于0~1,数值越大,长石和岩屑的溶蚀程度越低,溶孔越不发育。

        图  4  充填指数C及抗溶蚀指数K示意图

        Figure 4.  Schematic diagram of filling index C and corrosion resistance index K

        根据铸体薄片及实测物性数据统计结果(图5),充填指数C和抗溶蚀指数K与孔隙度具有较好的相关性,相关系数分别为0.674 7和0.597 3,与渗透率的相关性较差,相关系数分别为0.353 5和0.280 9,表明充填指数C和抗溶蚀指数K能够较准确的表征现今致密砂岩的储集空间大小,其中充填指数C对孔隙度的控制作用(相关性)强于抗溶蚀指数K,这与华庆地区长6储集空间以粒间孔为主,溶孔为辅的特征相匹配。而充填指数C和抗溶蚀指数K与渗透率的低相关性反映了致密砂岩的渗透率受多种因素控制、单一参数无法准确描述渗透率变化规律的事实[20]。因此,充填指数C和抗溶蚀指数K可用于定量评价成岩作用导致的孔隙度变化的非均质性。

        图  5  充填指数C、抗溶蚀指数K与孔隙度、渗透率的关系

        Figure 5.  Relationship between filling index C, corrosion resistance index K and porosity and permeability

      • 上述充填指数C和抗溶蚀指数K能够定量表征致密砂岩所经历的充填胶结及溶蚀强度,但无法反映成岩作用导致的孔隙空间的复杂程度,即不能表征孔喉结构的非均质性。对于致密砂岩而言,孔喉结构非均质性是影响微观渗流特征及剩余油的重要因素,是微观非均质性评价的核心[1,5]

        高压压汞实验是研究致密砂岩孔喉结构的常用方法,其实验压力最高可达220 MPa,能够反映3.7 nm以上孔径分布、结构等信息。根据高压压汞实验结果,华庆地区长6砂岩的排驱压力介于0.06~25.00 MPa,平均3.90 MPa;中值压力介于0.90~36.88 MPa,平均10.14 MPa;中值半径介于0.01~0.81 mm,平均0.14 mm;分选系数介于0.04~3.87,平均1.56;最大汞饱和度介于28.60%~84.30%,平均58.60%;退汞效率介于15.30%~64.50%,平均为26.80%。毛管压力曲线形态差异较大(图6),孔喉结构非均质性整体较强。

        图  6  华庆地区长6储层毛管压力曲线特征

        Figure 6.  Capillary pressure curve for Chang 6 reservoir in Huaqing area

        尽管高压压汞能够获取表征孔喉结构特征的多个实验参数,但各参数表征的意义不尽相同,缺少一个相对聚焦、能够评价孔喉大小、分选及复杂程度的综合指数[23-24]。因此,本次研究通过几何学中的分形维数来研究致密砂岩孔喉结构的非均质性。分形维数是目前定量表征孔喉结构非均质性及复杂程度简单而有效的途径[25],最早由美籍法国数学家Mandelbrot[26]提出,用于表征复杂形体的不规则性,数值介于2~3;之后,Thomvson et al.[27]认为沉积岩也具有分形特征,分形维数越接近3,表明孔喉分选越差、孔喉结构越复杂,微观非均质性越强;反之,则相反。根据分形几何理论,储层毛细管压力遵从拉普拉斯(Laplace)方程(即P c=2σcosθ/r,P c为孔径相应的毛细管压力;σ为液体的表面张力;θ为润湿接触角)的条件下,可得:

        l g S w = D - 3 l g P c + 3 - D l g P m i n (2)

        式中:D为孔隙结构分形维数;P c为任意孔隙半径所对应的毛细管压力;P min为入口毛细管压力,即最大孔径对应的毛细管压力;S w为压力为P c时的润湿相的饱和度。若样品孔隙结构具有分形特征,则lgS w和lgP c在双对数坐标系下为显著相关的直线,通过直线的斜率及截距就可以求出分形维数D

        结果表明,分形维数与压汞实测物性呈现出一定相关性(图7),其中渗透率与分形维数相关性较明显,相关系数R为0.630 8,孔隙度与分形维数相关性较差,相关系数R仅为0.472 2。如图8所示,由样品1—样品4,随着致密砂岩孔隙度和渗透率的降低,填隙物充填程度增加,大孔喉逐渐被填隙物分隔、瓦解成为数量众多的次级小孔隙,孔径分布也由微米级向亚微米—纳米级转变,分形维数趋近于3,孔隙结构也趋于复杂。因此,分形维数较好地描述、表征了致密砂岩微观孔喉结构非均质性的强弱。

        图  7  分形维数与物性关系

        Figure 7.  Relationship between fractal dimension and physical properties

        图  8  储层致密程度变化与分形维数关系

        Figure 8.  Relationship between changes in reservoir density and fractal dimensions

      • 对标准偏差σ、充填指数C、抗溶蚀指数K及分形维数D等4个关键参数进行反复论证,构建了微观非均质性综合指数F,用于评价致密砂岩储层的微观特征差异,具体公式如图9所示。微观非均质性综合指数F是一个相对全面的多因素复合参数,其包含了骨架颗粒粒度分选、成岩作用及孔隙结构复杂程度等相关信息。从上述单因素与面孔率及孔渗相关性结果来看,粒度标准偏差σ与粒间孔、面孔率的相关系数分别为0.652 6和0.551 9,充填指数C和抗溶蚀指数K与孔隙度的相关系数分别为0.674 7和0.597 3,分形维数D与孔隙度、渗透率的相关系数分别为0.472 2和0.630 8,从单因素所占权重来看,4个特征参数基本相当,并未表现出明显的权重优势。同时,考虑到不同实验原理及实验条件的差异,对特征参数所占权重定量化存在一定难度,因此本研究认为粒度标准偏差σ、充填指数C、抗溶蚀指数K和分形维数D对微观非均质性综合指数F的影响权重相当。通过计算,华庆地区长6致密砂岩的F值介于0.2~3.0,F值越大,代表储层骨架颗粒分选越差、成岩作用越强、孔喉结构越复杂,微观非均质性越强,反之则相反。

        图  9  微观非均质性综合评价指数F的构建模型

        Figure 9.  Construction model of comprehensive evaluation index F of micro heterogeneity

        油藏的许多宏观生产规律及储集层的渗流特征,均是在储集层微观结构及孔隙尺度下的各相流体运移的综合反映,即岩石的微观非均质性是根本,宏观特征是表象[5]。微观非均质性与油田生产的结合是目前微观非均质评价的难点,油水生产动态往往受生产制度、注水量、宏观非均质性(油层厚度、隔夹层发育)等多重因素控制,因此微观非均质性评价方法和标准的生产验证往往无法操作,但微观非均质性对肉眼及显微尺度下的储层油水运动及分布的影响可在实验室条件下得到很好验证。本次研究通过真实砂岩模型油水两相渗流实验,观察分析微观非均质对渗流特征及采收率的影响,论证微观非均质性综合评价指数F与生产特征之间的关系。

        真实砂岩模型(1997,西北大学地质系获国家专利)是由实际岩心经抽提、烘干、切片、磨平等工序之后,粘贴在两片玻璃之间制作而成的。模型长宽约为2.5 cm×2.5 cm,承压上限为0.2 MPa,加压耐温上限为100 ℃,实验模拟油、水的性质与实际油藏接近,红色为油,蓝色为水,模型入口及出口处均连接有玻璃量管,可精确记录油水体积变化。真实砂岩模型的优点是能够直接观察流体在实际油层岩石孔隙空间中的渗流特征。实验过程中,对模型依次进行饱和水、饱和油(油驱水)、水驱油过程,其中水驱油过程分为无水期驱油(模型出口未见水之前)和最终驱油(模型出口不出油时)过程,通过模型入口及出口玻璃量管记录的油水体积变化,计算无水期及最终期驱油效率。本次研究共制作真实砂岩微观模型15块,实验开始前对模型进行模拟水渗透率测定。

        通过油水渗流实验的实时观察,各模型无论是饱和油还是水驱油过程,均存在团块状的局部孔隙致密区,油水均无法进入,各模型孔隙致密区的发育面积存在差异,体现了不同模型之间的微观非均质性差异,而这种差异在油气储层中是广泛存在的。孔隙致密区对于水驱油过程的影响程度要强于对饱和油过程的影响程度,体现在水驱油过程中,注入水的渗流路径并未按照之前油进入模型的通道推进,而是呈现出3种不同的渗流路径,大致可分为均匀型、网状型及指状型(图10),其中均匀型渗流对应的模型孔隙致密区面积小且分布较分散,注入水均匀推进,驱油效率相对较高(图10a);网状型渗流对应的模型孔隙致密区呈较大团块状,注入水绕过孔隙致密区,沿着高渗通道呈网状驱油,形成团块状残余油,驱油效率相对较低(图10b);指状型渗流通常既受到孔隙致密区的影响,又可能受到诸如微层理及微裂缝的影响,微观非均质性最强,形成大面积的残余油,驱油效率最低(图10c)。因此,微观非均质性控制了流体的渗流特征,进而影响到剩余油的分布和油田注水开发效果。

        图  10  华庆地区长6真实砂岩模型实验典型样品渗流特征

        Figure 10.  Seepage characteristics of typical samples of Chang 6 in true sandstone model experiment in Huaqing area

        通过建立无水期及最终期驱油效率与微观非均质综合指数F的关系可以看出(表4图11),F指数与无水期及最终期驱油效率呈较好的相关性,相关系数R分别为0.736和0.784,这是因为F指数大小是以孔隙致密区发育规模的形式体现,而孔隙致密区往往是沉积水动力波动导致的细粒物质发育区,或是孔隙充填程度高且溶蚀作用欠发育区,因此F指数越大,孔隙致密区发育面积越大,储层微观非均质性越强,导致注入水无法均匀推进,形成网状型或指状型的渗流通道,注入水快速突进到油井导致水淹现象,形成绕流残余油,驱油效率较低。反之,F指数越小,储层微观非均质性越弱,孔隙致密区不发育,注入水以均匀推进为主,驱油效率较高。综上,F指数能够对模型中由孔隙致密区导致的储层微观特征差异定量化,可用于评价致密砂岩油藏内部微观非均质性,从而预测油田注入水的地下渗流特征及剩余油分布规律。

        表 4  华庆地区长6真实砂岩模型实验数据

        Table 4.  Statistics of experimental data of Chang 6 real true sandstone model in Huaqing area

        井号 井深/m 模型渗透率/(×10-3 μm2 微观非均质性综合指数F 无水期驱油效率/% 最终期驱油效率/% 水驱油渗流特征
        HQ211-33 2 176.4 0.137 1.69 20.2 44.8 网状
        HQ211-33 2 182.2 0.203 0.87 25.8 43.5 网状
        HQ211-33 2 187.1 0.325 0.81 25.4 45.1 网状
        HQ211-33 2 187.1 0.361 0.52 34.9 47.0 均匀状
        HQ264-59 2 112.5 0.220 0.96 25.6 51.1 均匀状
        HQ411 2 109.7 0.119 1.65 22.9 44.7 网状
        HQ411 2 123.5 0.465 0.82 28.9 46.5 均匀状
        HQ413 2 180.3 0.089 1.36 19.2 49.6 网状
        HQ413 2 182.1 0.146 2.27 21.1 42.5 网状
        HQ453 1 873.5 0.095 1.33 19.7 40.7 网状
        HQ453 1 874.8 0.452 0.88 25.0 46.5 网状
        HQ156 2 059.9 0.450 1.16 21.4 42.0 均匀状
        HQ262-56 2 178.9 0.225 1.65 18.8 38.2 网状
        HQ263-57 2 094.1 0.110 2.75 17.0 34.4 指状
        HQ263-57 2 101.0 0.118 2.84 20.1 29.1 指状

        图  11  微观非均质综合指数F与真实砂岩模型驱油效率关系

        Figure 11.  Relationship between microscopic heterogeneous composite index F and displacement efficiency of true sandstone model

        基于上述评价方法及计算结果,提出了华庆地区长6深水致密砂岩储层微观非均质性评价标准,F指数<1.0,储层微观特征表现为Ⅰ类弱非均质性,1.0<F指数<2.0,储层微观特征表现为Ⅱ类中等非均质性,F指数>2.0,储层微观特征表现为Ⅲ类强非均质性。

      • (1) 粒度标准偏差σ较好地表征了岩石骨架颗粒的非均质性,反映原始水动力变化导致的颗粒分选差异。粒度标准偏差σ与粒间孔呈现较强相关性,表明原始碎屑颗粒的分选程度对现今致密砂岩残余粒间孔产生一定影响,尽管后期压实、胶结等成岩作用破坏了大部分原始粒间孔隙,但分选性好的碎屑颗粒所搭建的岩石骨架更易保存粒间孔隙。

        (2) 充填指数C和抗溶蚀指数K从填隙物充填及颗粒溶蚀程度上定量表征了成岩作用的非均质性,二者与实测孔隙度具有较好的相关性,表明充填指数C和抗溶蚀指数K能够较准确的表征现今致密砂岩的残余储集空间大小,但对渗透率的表征程度较低,反映致密砂岩渗透率受多重因素控制,单一参数难以准确描述渗透率的变化规律。

        (3) 分形维数D是描述孔喉大小、孔径分布、孔喉配置关系的综合结构指数,D指数与渗透率的相关性更为明显,表明渗透率与微观孔隙结构的复杂程度关系密切,当储层致密程度增加时,孔径分布由微米级向亚微米—纳米级转变,分形维数趋于3,孔隙结构也趋于复杂。

        (4) 孔隙致密区广泛存在于致密砂岩储层中,采用微观非均质性综合指数F可以精确表征致密砂岩微观非均质性,F指数<1.0,储层为弱非均质性,1.0<F指数<2.0,储层为中等非均质性,F指数>2.0,储层为强非均质性。

    参考文献 (27)

    目录

      /

      返回文章
      返回