简世凯^1,2,,
符力耘^1,2,,,
王志伟³,
韩同城^1,2,
刘建林¹
1. 中国石油大学(华东)深层油气重点实验室, 青岛 266580
2. 中国石油大学(华东)地球科学与技术学院, 青岛 266580
3. 中国科学院油气资源研究重点实验室, 中国科学院地质与地球物理研究所, 北京 100029

基金项目: 国家重大专项课题"页岩气勘探地球物理技术研究"（2017ZX05036-005）和中国石油大学（华东）研究生创新工程项目（YCX2019002）联合资助


详细信息

作者简介: 简世凯, 男, 1991年生, 博士研究生, 主要研究数字岩心与裂缝储层反演.E-mail:jianshiakai@163.com

通讯作者: 符力耘, 男, 1964年生, 教授, 博士生导师, 主要从事地震学、岩石物理与勘探地震学相关研究.E-mail:lfu@upc.edu.cn

中图分类号: P631

收稿日期:2020-02-14
修回日期:2020-05-23
上线日期:2020-07-25

Elastic equivalent numerical modeling based on the dynamic method of Longmaxi Formation shale digital core

JIAN ShiKai^1,2,,
FU LiYun^1,2,,,
WANG ZhiWei³,
HAN TongCheng^1,2,
LIU JianLin¹
1. Key Laboratory of Deep Oil and Gas, China University of Petroleum(East China), Qingdao 266580, China
2. School of Geosciences, China University of Petroleum(East China), Qingdao 266580, China
3. Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China


More Information

Corresponding author: FU LiYun,E-mail:lfu@upc.edu.cn

MSC: P631

--> Received Date: 14 February 2020
Revised Date: 23 May 2020
Available Online: 25 July 2020


摘要
摘要:数字岩心是计算岩石弹性性质的一类常用方法.龙马溪组页岩具有多矿物构成、复杂微结构和强非均质等特征，常规岩石物理的弹性等效解析建模局限性较大，目前流行的静态数值等效建模方法的精度有限.本文基于高分辨率的页岩数字岩心数据，采用多阈值分割方法将数字岩心分解为黏土、石英、孔隙、TOC、长石类、黄铁矿类等六种矿物类型；利用矿物组分等效模量法计算各类矿物的弹性模量；采用二元函数分水岭方法表征不同压力下的岩石孔隙变形和颗粒接触关系变化；通过取向分布函数（ODF）定量分析矿物颗粒展布造成的各向异性特征.最后基于Biot孔弹方程，采用不分裂卷积完全匹配层（CPML）旋转交错网格有限差分法模拟不同压力下弹性波在数字岩心中的传播.以未加压的数字岩心为参考模型，计算不同压力下弹性波走时的平均时间差，进而估算各压力点的数字岩心等效速度.与该岩心样品的超声实验测量速度比较，动态法数值计算结果略偏高，据此校正数值计算过程中表征岩石微结构及颗粒接触关系随压力变化的二元函数，有效改善动态法弹性等效数值建模精度.
关键词: 龙马溪组页岩/
数字岩心/
动态法弹性等效数值建模/
Biot孔弹方程/
有限差分模拟

Abstract:It is a widespread method to analyze rock elastic properties based on digital cores. As the result of many minerals, complex microstructures and strong heterogeneity in Longmaxi Formation shale, the accuracy of elastic equivalent analysis modeling by conventional rock physics method is limited, and the current popular static numerical equivalent modeling method also has strong limitations. We apply multi-threshold segmentation method to divide the high-resolution digital core images of shale into six types including pore, TOC, clay, quartz, feldspar and pyrite, whose elastic modulus were calculated by the mineral components equivalent modulus method. The changes of pore deformation and grain-to-grain contact in rock mass under different pressures were characterized by binary function watershed method. And we use orientation distribution function (ODF) to quantitatively characterize the anisotropy of mineral grain distribution. Finally, Biot poroelastic equation, applying the rotated staggered grid finite-difference technique with unsplit convolution perfectly matched layer (CPML), is used to simulate the propagation of elastic waves in digital core models under different pressures. Taking the uncompressed digital core model as the criterion, the average time delay of elastic wave propagation under different pressures is calculated, which contributes to estimate the equivalent velocity of the digital core models in each pressure. The numerical equivalent velocity is slightly higher than the velocity measured by ultrasonic experiment on the same core, so the binary function characterizing change of the rock mass microstructure and grain-to-grain contact dependent-pressure can be corrected to effectively improve the accuracy of dynamic elastic equivalent numerical modeling.
Key words:Longmaxi Formation shale/
Digital core/
Elastic equivalent numerical modeling based on the dynamic method/
Biot poroelastic equation/
Finite-difference simulation



PDF全文下载地址:

http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_15525