CN109490963B - 裂缝储层岩石物理建模方法及系统 - Google Patents
裂缝储层岩石物理建模方法及系统 Download PDFInfo
- Publication number
- CN109490963B CN109490963B CN201710824095.2A CN201710824095A CN109490963B CN 109490963 B CN109490963 B CN 109490963B CN 201710824095 A CN201710824095 A CN 201710824095A CN 109490963 B CN109490963 B CN 109490963B
- Authority
- CN
- China
- Prior art keywords
- fracture
- seismic data
- parameters
- obtaining
- threshold value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000011435 rock Substances 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 238000003384 imaging method Methods 0.000 claims abstract description 13
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 8
- 238000009499 grossing Methods 0.000 claims abstract description 8
- 238000012937 correction Methods 0.000 claims description 18
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 230000010287 polarization Effects 0.000 claims description 5
- 230000004044 response Effects 0.000 abstract description 13
- 238000010586 diagram Methods 0.000 description 26
- 239000012530 fluid Substances 0.000 description 15
- 238000004458 analytical method Methods 0.000 description 9
- 238000002310 reflectometry Methods 0.000 description 8
- 229920006395 saturated elastomer Polymers 0.000 description 6
- 239000013598 vector Substances 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000012512 characterization method Methods 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
公开了一种裂缝储层岩石物理建模方法及系统。该方法可以包括:根据测井曲线,对速度曲线与密度曲线进行平滑处理,获取骨架参数C0;根据地层成像测井资料,获得目的层段的裂缝参数,进而建立裂缝参数初始模型;基于骨架参数与裂缝参数初始模型,获得裂缝介质的有效弹性常数矩阵Ceff,进而获得正演的地震数据;设定门槛值,如果正演的地震数据与实际的地震数据误差大于门槛值,根据误差修正模型参数,直到正演的地震数据与实际的地震数据误差小于或等于门槛值,获得最终的裂缝储层岩石物理模型。本发明采用地震属性约束建立的裂缝储层岩石物理模型,能得到与地震尺度一致裂缝参数,能更好的模拟裂缝储层参数变化的地震响应变化,指导地震属性裂缝预测。
Description
技术领域
本发明涉及地球物理技术领域,更具体地,涉及一种裂缝储层岩石物理建模方法及系统。
背景技术
裂缝储层进行建模研究主要包括裂缝性储层表征及建模。裂缝表征常用的方法包括岩心分析法、野外露头分析法、测井分析法、地震预测分析法等。此发明主要提出的裂缝储层建模方法是基于测井分析法得到裂缝性表征并通过等效模型进行建模,反演得到的裂缝基本参数是否能最大程度表征裂缝成为该技术的关键。
在常规测井曲线方法中,不同测井资料由于各自数据采集原理的差异,反应的裂缝属性不同。邓攀等(2002)结合三孔隙度曲线特征对火山岩储层的构造裂缝进行了识别并计算了裂缝孔隙度大小,在此基础上进行建模,然而该方法的适用性比较差。李琼等(2006)通过神经网络进行裂缝识别;唐洪等(2012)利用常规测井资料,采用曲线元的方法计算存在概率,表征出老区块裂缝的发育特征。这些常规测井资料虽然能够预测并表征裂缝的各项参数,但干扰信息较多同时受限于纵向分辨率,预测和表征的精度相对于成像测井资料较差。成像测井资料一方面能够获得高精度、高分辨率的井壁扫描图像,另一方面在垂向上具有连续性和直观性的优势,成为目前识别裂缝分布和产状最为有效的测井手段。
在对裂缝储层岩石物理建模方面,印兴耀(2014)在岩石基质中利用等效介质理论添加含流体裂隙,然后结合Hudson理论和各向异性Gassmann理论构建岩石物理模型。选择理论模型骨架参数模拟分析不同裂缝参数介质地震波场响应。曾勇坚(2014)对基于致密砂岩裂缝储层的岩石物理建模作了相关分析(曾勇坚,印兴耀,曹丹平.基于致密砂岩裂缝储层的岩石物理建模分析.2014年中国地球科学联合学术年会—专题20:岩石物理与非常规油气勘探开发论文集),首先是通过利用Voigt-Reuss-Hill平均计算岩石基质的体积模量和剪切模量,采用DEM理论计算岩石骨架的相关参数,然后利用Thomsen裂缝理论和Brown-Korringa理论建模。上述裂缝储层岩石物理建模方法都是应用测井解释裂缝参数正演模拟不同裂缝介质的地震响应,但是测井解释的裂缝参数与地震频带反映的裂缝参数是两个不同尺度裂缝参数,存在不相容性,测井尺度解释结果建立的储层裂缝模型地震响应往往与实际地震记录存在残差,裂缝参数和含流体性改变后的地震响应也不能有效指导地震振幅变化预测裂缝参数和含流体性。因此,有必要开发一种裂缝储层岩石物理建模方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种裂缝储层岩石物理建模方法及系统,采用地震属性约束建立的裂缝储层岩石物理模型,能得到与地震尺度一致裂缝参数,能更好的模拟裂缝储层参数变化的地震响应变化,指导地震属性裂缝预测。
根据本发明的一方面,提出了一种裂缝储层岩石物理建模方法。所述方法可以包括:根据测井曲线,对速度曲线与密度曲线进行平滑处理,获取骨架参数C0;根据地层成像测井资料,获得目的层段的裂缝参数,进而建立裂缝参数初始模型;基于所述骨架参数与所述裂缝参数初始模型,获得裂缝介质的有效弹性常数矩阵Ceff,进而获得正演的地震数据;设定门槛值,如果所述正演的地震数据与实际的地震数据误差小于或等于所述门槛值,获得最终的裂缝储层岩石物理模型。
优选地,还包括:如果所述正演的地震数据与实际的地震数据误差大于所述门槛值,修正裂缝参数初始模型,直至所述正演的地震数据与实际的地震数据误差小于或等于所述门槛值,获得最终的裂缝储层岩石物理模型。
优选地,根据平滑处理后的速度曲线与密度曲线,获取所述骨架参数C0。
优选地,所述裂缝参数包括:目的层段裂缝纵横比与裂缝密度。
优选地,还包括:基于所述最终的裂缝储层岩石物理模型,修改裂缝参数,获得该裂缝参数对应的裂缝储层岩石物理模型;正演分析不同的裂缝储层岩石物理模型的振幅、旅行时与偏振随偏移距变化属性,指导地震属性裂缝预测。
优选地,所述有效弹性常数矩阵为:
ceff=c0+c1+c2 (4)
其中,c1表示骨架参数的一阶更正,c2表示骨架参数的二阶更正。
优选地,所述骨架参数的一阶更正为:
优选地,所述骨架参数的二阶更正为:
根据本发明的另一方面,提出了一种裂缝储层岩石物理建模系统,可以包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:根据测井曲线,对速度曲线与密度曲线进行平滑处理,获取骨架参数C0;根据地层成像测井资料,获得目的层段的裂缝参数,进而建立裂缝参数初始模型;基于所述骨架参数与所述裂缝参数初始模型,获得裂缝介质的有效弹性常数矩阵Ceff,进而获得正演的地震数据;设定门槛值,如果所述正演的地震数据与实际的地震数据误差小于或等于所述门槛值,获得最终的裂缝储层岩石物理模型。
优选地,还包括:如果所述正演的地震数据与实际的地震数据误差大于所述门槛值,修正裂缝参数初始模型,直至所述正演的地震数据与实际的地震数据误差小于或等于所述门槛值,获得最终的裂缝储层岩石物理模型。
本发明的有益效果:通过迭代修正裂缝参数初始模型,最终得到与井旁地震道误差最小对应的储层裂缝参数,建立地震频带尺度的裂缝储层岩石物理模型;再修改裂缝参数和裂缝中含流体性,正演模拟不同裂缝参数时,地震振幅属性的变化特征,指导地震资料裂缝解释。
本发明的方法和装置具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的裂缝储层岩石物理建模方法的步骤的流程图。
图2a、图2b和图2c分别示出了根据本发明的一个实施例的测井资料包括纵波速度、横波速度、密度的示意图。
图3a、图3b和图3c分别示出了根据本发明的一个实施例的裂缝的倾向、走向、倾角统计的示意图。
图4a、图4b和图4c分别示出了根据本发明的一个实施例的成像测井得到的裂缝半径的示意图。
图5示出了根据本发明的一个实施例的裂缝储层岩石物理模型的示意图。
图6a、图6b、图6c分别示出了根据本发明的一个实施例的反射率法正演得到的方位角为0°的垂直分量、径向分量与切向分量的示意图;图6d、图6e、图6f分别示出了根据本发明的一个实施例的反射率法正演得到的方位角为30°的垂直分量、径向分量与切向分量的示意图;图6g、图6h、图6i分别示出了根据本发明的一个实施例的反射率法正演得到的方位角为60°的垂直分量、径向分量与切向分量的示意图;图6j、图6k、图6m分别示出了根据本发明的一个实施例的反射率法正演得到的方位角为90°的垂直分量、径向分量与切向分量的示意图。
图7a和图7b分别示出了根据本发明的一个实施例的归一化的正演数据和实际数据的示意图。
图8a、图8b、图8c、图8d、图8e和图8f示出了根据本发明的一个实施例的第一次迭代至第六次迭代的归一化正演数据和实际数据的残差的示意图。
图9a和图9b分别示出了根据本发明的一个实施例的裂缝层饱含水时的顶层与底层的PP波椭圆拟合的示意图。
图10a和图10b分别示出了根据本发明的一个实施例的裂缝层饱含水时的顶层与底层的PS波椭圆拟合的示意图。
图11a和图11b分别示出了根据本发明的一个实施例的裂缝层不含流体时的顶层与底层的PP波椭圆拟合的示意图。
图12a和图12b分别示出了根据本发明的一个实施例的裂缝层不含流体时的顶层与底层的PS波椭圆拟合的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
针对测井和地震裂缝尺度不相容问题,本发明提出一种裂缝储层岩石物理建模方法。本发明应用测井数据建立裂缝参数初始模型,并通过正演地震波场与实际地震记录误差优化裂缝参数,最终建立地震频带尺度的裂缝储层岩石物理模型,并进一步正演分析储层裂缝参数和含流体性变化的地震波场响应变化,指导地震资料预测裂缝。
图1示出了根据本发明的裂缝储层岩石物理建模方法的步骤的流程图。
在该实施例中,根据本发明的裂缝储层岩石物理建模方法可以包括:
步骤101,根据测井曲线,对速度曲线与密度曲线进行平滑处理,获取骨架参数C0;在一个示例中,根据平滑处理后的速度曲线与密度曲线,获取骨架参数C0。
具体地,测井资料与相关地质信息可以反映地下复杂岩性的地质参数,设有一组经过环境影响校正的实际测井值向量a为公式(1):
a=(a1,a2,a3,…,am)T (1)
其中,a1,a2,a3,…,am分别表示不同深度的测井值,m表示不同深度的点的数量,且有一组待求的未知储层参数及矿物相对体积向量x为公式(2)
x=(x1,x2,x3,…,xm)T (2)
其中,x1,x2,x3,…,xm分别表示不同深度的点的矿物的体积含量。
则根据所采用的测井解释模型,各测井响应方程均可表示为公式(3):
a1=fi(x,z),(i=1,2,…,m) (3)
其中,z为一组区域性解释参数向量,不同的测井响应方程有不同的解释参数。这样,根据测井资料可以得到地层厚度及背景介质的弹性性质建立的裂缝模型,进而对速度曲线与密度曲线进行平滑处理,获取骨架参数C0。
步骤102,根据地层成像测井资料,获得目的层段的裂缝参数,进而建立裂缝参数初始模型;在一个示例中,裂缝参数包括:目的层段裂缝纵横比与裂缝密度。
步骤103,基于骨架参数与裂缝参数初始模型,获得裂缝介质的有效弹性常数矩阵Ceff,进而获得正演的地震数据;
在一个示例中,有效弹性常数矩阵为:
ceff=c0+c1+c2 (4)
其中,c1表示骨架参数的一阶更正,c2表示骨架参数的二阶更正。
在一个示例中,骨架参数的一阶更正为:
在一个示例中,骨架参数的二阶更正为:
步骤104,设定门槛值,如果正演的地震数据与实际的地震数据误差小于或等于门槛值,获得最终的裂缝储层岩石物理模型。
在一个示例中,还包括:如果正演的地震数据与实际的地震数据误差大于门槛值,修正裂缝参数初始模型,直至正演的地震数据与实际的地震数据误差小于门槛值,获得最终的裂缝储层岩石物理模型。
在一个示例中,还包括:基于最终的裂缝储层岩石物理模型,修改裂缝参数和裂缝中含流体性,获得该裂缝参数对应的裂缝储层岩石物理模型;正演分析不同的裂缝储层岩石物理模型的振幅、旅行时与偏振随偏移距变化属性,指导地震资料裂缝解释。
具体地,Hudson模型是基于对含有薄硬币形状的椭球裂缝或包含物的弹性固体中的平均波长的散射理论分析得到的,在已知岩石的成分和孔隙空间,Hudson模型可用来估算其等效弹性模量和衰减。Hudson模型具有如下假设:连续性假设;长波近似,即ka<1(k为波数,a为裂隙的特征尺度);在地震波长范围内裂纹位置分布均匀;裂纹稀疏且不连通,裂纹总体积占岩石总体积的比例很少;裂纹是具有小的纵横比的扁球体,即近似为薄币形。
裂隙半径和裂隙间距必须远小于波长。由于裂隙彼此之间是隔离的,流体不能相互流动,因此Hudson模型模拟的是非常高频情况下的饱和岩石属性,即获取有效弹性常数矩阵为公式(4)。
对于一个缝隙垂直方向沿3轴方向的单一的裂缝组,裂缝介质显示横向各向同性,且骨架参数的一阶更正为公式(5),骨架参数的二阶更正为公式(6)。
设定门槛值,如果正演的地震数据与实际的地震数据误差小于或等于门槛值,获得最终的裂缝储层岩石物理模型;如果正演的地震数据与实际的地震数据误差大于门槛值,修正裂缝参数初始模型,进而重复上述步骤,直至正演的地震数据与实际的地震数据误差小于或等于门槛值。基于最终的裂缝储层岩石物理模型,修改裂缝参数和裂缝中含流体性,获得该裂缝参数对应的裂缝储层岩石物理模型;正演分析不同的裂缝储层岩石物理模型的振幅、旅行时与偏振随偏移距变化属性,指导地震资料裂缝解释。
本方法采用地震属性约束建立的裂缝储层岩石物理模型,能得到与地震尺度一致裂缝参数,能更好的模拟裂缝储层参数变化的地震响应变化,指导地震属性裂缝预测。
应用示例
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
图2a、图2b和图2c分别示出了根据本发明的一个实施例的测井资料包括纵波速度、横波速度、密度的示意图。
测井资料与相关地质信息可以反映地下复杂岩性的地质参数,设有一组经过环境影响校正的实际测井值向量a为公式(1),且有一组待求的未知储层参数及矿物相对体积向量x为公式(2)。则根据所采用的测井解释模型,各测井响应方程均可表示为公式(3)。这样,根据测井资料可以得到地层厚度及背景介质的弹性性质建立的裂缝模型,如图2a-图2c所示,由图可知,岩性以岩屑石英砂岩、粉砂岩、泥质粉砂岩、页岩为主,裂缝孔隙度相对较高,裂缝有效性较好。进而对速度曲线与密度曲线进行平滑处理,获取骨架参数C0,其中,模型的参数如表1所示。
表1
层位 | 密度(g/cm<sup>3</sup>) | 纵波速度(km/s) | 横波速度(km/s) | 厚度(m) |
盖层 | 2.6 | 4.7 | 2.5 | 1000 |
储层 | 2.7 | 5.52 | 2.9 | 300 |
半空间 | 2.8 | 6 | 3.15 | ∞ |
图3a、图3b和图3c分别示出了根据本发明的一个实施例的裂缝的倾向、走向、倾角统计的示意图。
图4a、图4b和图4c分别示出了根据本发明的一个实施例的成像测井得到的裂缝半径的示意图,其中,图4a的裂缝半径为1.5m,图4b和图4c的裂缝半径均为2m。
图5示出了根据本发明的一个实施例的裂缝储层岩石物理模型的示意图。
图6a、图6b、图6c分别示出了根据本发明的一个实施例的反射率法正演得到的方位角为0°的垂直分量、径向分量与切向分量的示意图;图6d、图6e、图6f分别示出了根据本发明的一个实施例的反射率法正演得到的方位角为30°的垂直分量、径向分量与切向分量的示意图;图6g、图6h、图6i分别示出了根据本发明的一个实施例的反射率法正演得到的方位角为60°的垂直分量、径向分量与切向分量的示意图;图6j、图6k、图6m分别示出了根据本发明的一个实施例的反射率法正演得到的方位角为90°的垂直分量、径向分量与切向分量的示意图。
根据地层成像测井资料,获得目的层段的裂缝参数,进而建立裂缝参数初始模型,其中,裂缝参数包括:目的层段裂缝的条数、纵横比与裂缝密度,如图图3a-图3c、图4a-图4c所示,分析裂缝的倾角发现,高导缝在70-90度的占到了接近50%,单90度的就达到了20%,所以可以设定裂缝为垂直裂缝,从方位角和倾向也能看出,裂缝主要以垂直裂缝为主。初始的裂缝参数以及得到的弹性矩阵如表2所示。
表2
Hudson模型是基于对含有薄硬币形状的椭球裂缝或包含物的弹性固体中的平均波长的散射理论分析得到的,在已知岩石的成分和孔隙空间,Hudson模型可用来估算其等效弹性模量和衰减。Hudson模型具有如下假设:连续性假设;长波近似,即ka<1(k为波数,a为裂隙的特征尺度);在地震波长范围内裂纹位置分布均匀;裂纹稀疏且不连通,裂纹总体积占岩石总体积的比例很少;裂纹是具有小的纵横比的扁球体,即近似为薄币形。
裂隙半径和裂隙间距必须远小于波长。由于裂隙彼此之间是隔离的,流体不能相互流动,因此Hudson模型模拟的是非常高频情况下的饱和岩石属性,即获取有效弹性常数矩阵为:
对于一个缝隙垂直方向沿3轴方向的单一的裂缝组,裂缝介质显示横向各向同性,且骨架参数的一阶更正为公式(5),骨架参数的二阶更正为公式(6)。
图7a和图7b分别示出了根据本发明的一个实施例的归一化的正演数据和实际数据的示意图。
设定门槛值,如果正演的地震数据与实际的地震数据误差小于或等于门槛值,获得最终的裂缝储层岩石物理模型,如图7a所示;如果正演的地震数据与实际的地震数据误差大于门槛值,修正裂缝参数初始模型,进而重复上述步骤,直至正演的地震数据与实际的地震数据误差小于或等于门槛值。基于最终的裂缝储层岩石物理模型,修改裂缝参数和裂缝中含流体性,获得该裂缝参数对应的裂缝储层岩石物理模型;正演分析不同的裂缝储层岩石物理模型的振幅、旅行时与偏振随偏移距变化属性,指导地震资料裂缝解释。
图8a、图8b、图8c、图8d、图8e和图8f分别示出了根据本发明的一个实施例的第一次迭代至第六次迭代的归一化正演数据和实际数据的残差的示意图,图8a-图8f的残差逐渐减小,说明迭代修正的裂缝模型正演结果与实际地震记录越吻合,即裂缝模型越接近真实裂缝。
图9a和图9b分别示出了根据本发明的一个实施例的裂缝层饱含水时的顶层与底层的PP波椭圆拟合的示意图。
图10a和图10b分别示出了根据本发明的一个实施例的裂缝层饱含水时的顶层与底层的PS波椭圆拟合的示意图。
图11a和图11b分别示出了根据本发明的一个实施例的裂缝层不含流体时的顶层与底层的PP波椭圆拟合的示意图。
图12a和图12b分别示出了根据本发明的一个实施例的裂缝层不含流体时的顶层与底层的PS波椭圆拟合的示意图。
图9a-图12b分别是修改裂缝饱和流体特征后,正演的含裂缝层顶底的PP和PS波响应的椭圆特征。椭圆特征的变化结果对利用实际地震资料的振幅属性预测裂缝参数就有指导意义。
综上所述,本发明采用地震属性约束建立的裂缝储层岩石物理模型,能得到与地震尺度一致裂缝参数,能更好的模拟裂缝储层参数变化的地震响应变化,指导地震属性裂缝预测。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
根据本发明的实施例,提供了一种裂缝储层岩石物理建模系统,可以包括:存储器,存储有计算机可执行指令;处理器,处理器运行存储器中的计算机可执行指令,执行以下步骤:根据测井曲线,对速度曲线与密度曲线进行平滑处理,获取骨架参数C0;根据地层成像测井资料,获得目的层段的裂缝参数,进而建立裂缝参数初始模型;基于骨架参数与裂缝参数初始模型,获得裂缝介质的有效弹性常数矩阵Ceff,进而获得正演的地震数据;设定门槛值,如果正演的地震数据与实际的地震数据误差小于或等于门槛值,获得最终的裂缝储层岩石物理模型。
在一个示例中,还包括:如果正演的地震数据与实际的地震数据误差大于门槛值,修正裂缝参数初始模型,直至正演的地震数据与实际的地震数据误差小于或等于门槛值,获得最终的裂缝储层岩石物理模型。
本发明采用地震属性约束建立的裂缝储层岩石物理模型,能得到与地震尺度一致裂缝参数,能更好的模拟裂缝储层参数变化的地震响应变化,指导地震属性裂缝预测。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (6)
1.一种裂缝储层岩石物理建模方法,包括:
根据测井曲线,对速度曲线与密度曲线进行平滑处理,获取骨架参数c0;
根据地层成像测井资料,获得目的层段的裂缝参数,进而建立裂缝参数初始模型;
基于所述骨架参数与所述裂缝参数初始模型,获得裂缝介质的有效弹性常数矩阵ceff,进而获得正演的地震数据;
设定门槛值,如果所述正演的地震数据与实际的地震数据误差小于或等于所述门槛值,获得最终的裂缝储层岩石物理模型;
其中,还包括:
如果所述正演的地震数据与实际的地震数据误差大于所述门槛值,修正裂缝参数初始模型,直至所述正演的地震数据与实际的地震数据误差小于或等于所述门槛值,获得最终的裂缝储层岩石物理模型;
其中,根据平滑处理后的速度曲线与密度曲线,获取所述骨架参数c0;
其中,所述有效弹性常数矩阵为:
ceff=c0+c1+c2 (4)
其中,c1表示骨架参数的一阶更正,c2表示骨架参数的二阶更正。
2.根据权利要求1所述的裂缝储层岩石物理建模方法,其中,所述裂缝参数包括:目的层段裂缝纵横比与裂缝密度。
3.根据权利要求1所述的裂缝储层岩石物理建模方法,其中,还包括:
基于所述最终的裂缝储层岩石物理模型,修改裂缝参数,获得该裂缝参数对应的裂缝储层岩石物理模型;
正演分析不同的裂缝储层岩石物理模型的振幅、旅行时与偏振随偏移距变化属性,预测地震资料属性裂缝。
6.一种裂缝储层岩石物理建模系统,其特征在于,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
根据测井曲线,对速度曲线与密度曲线进行平滑处理,获取骨架参数c0;
根据地层成像测井资料,获得目的层段的裂缝参数,进而建立裂缝参数初始模型;
基于所述骨架参数与所述裂缝参数初始模型,获得裂缝介质的有效弹性常数矩阵ceff,进而获得正演的地震数据;
设定门槛值,如果所述正演的地震数据与实际的地震数据误差小于或等于所述门槛值,获得最终的裂缝储层岩石物理模型;
其中,还包括:
如果所述正演的地震数据与实际的地震数据误差大于所述门槛值,修正裂缝参数初始模型,直至所述正演的地震数据与实际的地震数据误差小于或等于所述门槛值,获得最终的裂缝储层岩石物理模型;
其中,根据平滑处理后的速度曲线与密度曲线,获取所述骨架参数c0;
其中,所述有效弹性常数矩阵为:
ceff=c0+c1+c2 (4)
其中,c1表示骨架参数的一阶更正,c2表示骨架参数的二阶更正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710824095.2A CN109490963B (zh) | 2017-09-13 | 2017-09-13 | 裂缝储层岩石物理建模方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710824095.2A CN109490963B (zh) | 2017-09-13 | 2017-09-13 | 裂缝储层岩石物理建模方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109490963A CN109490963A (zh) | 2019-03-19 |
CN109490963B true CN109490963B (zh) | 2021-04-23 |
Family
ID=65689093
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710824095.2A Active CN109490963B (zh) | 2017-09-13 | 2017-09-13 | 裂缝储层岩石物理建模方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109490963B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111781637B (zh) * | 2019-04-04 | 2023-06-30 | 中国石油天然气集团有限公司 | 近海底水合物储层建模方法及装置 |
CN110988988B (zh) * | 2019-11-25 | 2020-10-23 | 中国矿业大学(北京) | 基于垂直裂缝介质的地震波波场模拟方法及装置 |
CN111767647B (zh) * | 2020-06-19 | 2021-09-14 | 中国石油大学(华东) | 一种定量求取页岩中构造裂缝规模的方法 |
CN111751883B (zh) * | 2020-07-15 | 2021-07-30 | 中南大学 | 一种基于各向异性流体因子的致密裂缝储层地震预测方法 |
CN112198551A (zh) * | 2020-09-12 | 2021-01-08 | 北京恺标技术发展有限公司 | 一种通用的定量化评价储层裂缝强度的新方法 |
CN113156499B (zh) * | 2021-02-26 | 2023-05-26 | 中国石油天然气股份有限公司 | 台盆区针对裂缝型储层的地震资料叠后定量预测方法 |
CN115113266B (zh) * | 2021-03-19 | 2024-08-27 | 中国石油化工股份有限公司 | 火山岩深度域相控速度建模方法、存储介质和计算机设备 |
CN113189677B (zh) * | 2021-04-22 | 2022-05-31 | 西南石油大学 | 一种三维油藏物性参数模型自动更新方法 |
CN115576007B (zh) * | 2022-11-22 | 2023-03-14 | 西南石油大学 | 一种基于乱序矩阵的半确定性裂缝建模方法及系统 |
CN117166996B (zh) * | 2023-07-27 | 2024-03-22 | 中国地质大学(北京) | 地质参数门槛值的确定方法、装置、设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102445709A (zh) * | 2010-10-14 | 2012-05-09 | 中国石油大学(北京) | 与孔隙结构相关的全频带速度预测模型 |
CN104252007A (zh) * | 2013-06-26 | 2014-12-31 | 中国石油化工股份有限公司 | 一种相容性岩石物理建模方法 |
CN104360389A (zh) * | 2014-11-12 | 2015-02-18 | 中国石油大学(华东) | 一种致密砂岩储层岩石弹性模量计算方法 |
CN106597545A (zh) * | 2016-12-07 | 2017-04-26 | 中国矿业大学(北京) | 一种水平裂缝地震叠前反演方法和装置 |
-
2017
- 2017-09-13 CN CN201710824095.2A patent/CN109490963B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102445709A (zh) * | 2010-10-14 | 2012-05-09 | 中国石油大学(北京) | 与孔隙结构相关的全频带速度预测模型 |
CN104252007A (zh) * | 2013-06-26 | 2014-12-31 | 中国石油化工股份有限公司 | 一种相容性岩石物理建模方法 |
CN104360389A (zh) * | 2014-11-12 | 2015-02-18 | 中国石油大学(华东) | 一种致密砂岩储层岩石弹性模量计算方法 |
CN106597545A (zh) * | 2016-12-07 | 2017-04-26 | 中国矿业大学(北京) | 一种水平裂缝地震叠前反演方法和装置 |
Non-Patent Citations (3)
Title |
---|
储层地震岩石物理建模研究现状与进展;印兴耀 等;《石油物探》;20160531;第55卷(第3期);第309-325页 * |
地震约束油藏地质模型的更新;李德鹏;《中国优秀硕士学位论文全文数据库•基础科学辑》;20150615(第6期);A011-24 * |
基于岩石物理模型的页岩孔隙结构反演及横波速度预测;逄硕 等;《吉林大学学报(地球科学版)》;20170331;第47卷(第2期);第606-615页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109490963A (zh) | 2019-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109490963B (zh) | 裂缝储层岩石物理建模方法及系统 | |
CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
US7062072B2 (en) | Methods of producing images of underground formations surrounding a borehole | |
Hu et al. | Simultaneous multifrequency inversion of full-waveform seismic data | |
CN107203005B (zh) | 一种定量化计算裂缝描述参数的方法 | |
US9852373B2 (en) | Properties link for simultaneous joint inversion | |
US10422922B2 (en) | Method for predicting rock strength by inverting petrophysical properties | |
CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
US20060285437A1 (en) | Radial profiling of formation mobility using horizontal and vertical shear slowness profiles | |
WO2008070596A1 (en) | Identification of fracture clusters in rock formations | |
WO2016008105A1 (zh) | 一种基于柯西分布的叠后波阻抗反演方法 | |
US10705241B2 (en) | Determining sea water resistivity | |
CN103869362A (zh) | 体曲率获取方法和设备 | |
CN110095811B (zh) | 膏岩层速度模型构建处理方法及装置 | |
US11215035B2 (en) | Method to predict reservoir formation permeability using combined acoustic and multi-frequency dielectric measurements | |
Luo et al. | Probabilistic seismic petrophysical inversion with statistical double-porosity Biot-Rayleigh model | |
Shen et al. | Data-driven interpretation of ultradeep azimuthal propagation resistivity measurements: Transdimensional stochastic inversion and uncertainty quantification | |
Eikrem et al. | Bayesian estimation of reservoir properties—effects of uncertainty quantification of 4D seismic data | |
CN109521470B (zh) | 分析地质构造对地震反演裂缝密度影响的方法 | |
CN114325832A (zh) | 一种裂缝参数和弹性参数同步反演方法及系统 | |
US11867856B2 (en) | Method and system for reflection-based travel time inversion using segment dynamic image warping | |
Zhang | Ensemble methods of data assimilation in porous media flow for non-Gaussian prior probability density | |
CN109541692B (zh) | 一种基于地震资料共振成像中的时频分析方法 | |
CN114859407A (zh) | 火山岩储层声学特征参数的确定方法及装置 | |
Shen et al. | Statistical Inversion of Electromagnetic Logging Data |
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 |