CN105319603B - 致密砂岩储层复杂网状裂缝的预测方法 - Google Patents
致密砂岩储层复杂网状裂缝的预测方法 Download PDFInfo
- Publication number
- CN105319603B CN105319603B CN201510752304.8A CN201510752304A CN105319603B CN 105319603 B CN105319603 B CN 105319603B CN 201510752304 A CN201510752304 A CN 201510752304A CN 105319603 B CN105319603 B CN 105319603B
- Authority
- CN
- China
- Prior art keywords
- crack
- rock
- stress
- fracture
- density
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明属于石油勘探领域,具体地,涉及一种致密砂岩储层复杂网状裂缝的预测方法。预测方法步骤如下:建立构造地质模型及裂缝发育模型,测试古、今地应力大小及方向,完整岩石力学参数实验,测试岩石力学参数及裂缝应力敏感性,建立裂缝性岩石多级复合破裂准则,进行岩石变形物理测试获取峰值强度,建立单轴状态应力‑应变与裂缝体密度关系模型,建立三轴状态应力‑应变与裂缝体密度、产状关系模型,建立单轴状态应力‑应变与裂缝体密度关系模型,计算模拟现今条件下裂缝参数,验证裂缝定量预测结果的可靠性。本发明能准确地获取致密砂岩储层复杂裂缝参数并进行定量表征,适合于任何脆性储层为主的裂缝定量预测工作,减少了勘探开发的风险和成本。
Description
技术领域
本发明属于石油勘探领域,具体地,涉及一种致密砂岩储层复杂网状裂缝的预测方法。
背景技术
随着油气资源勘探开发逐渐由东部向西部、由常规储层向非常规储层的转变,寻找裂缝性油气藏已成为热点,如何预测裂缝的空间分布,并定量表征裂缝参数是石油地质研究的前沿问题。与其它储层相比,致密砂岩储层埋藏深度大、地层压力系数高、成岩作用强烈、构造运动期次多、网状裂缝非常发育,构造裂缝作为重要的储集空间和渗流通道,其形成和发育直接受控于构造应力场和储层岩石物理性质。因此,国内外学者多采用构造应力场法来预测裂缝分布,另外还有已知井点约束法、分形几何法、曲率分析法、测井分析法、地震分析法、应变量分析法、离散裂缝网络建模技术(DFN)等。实践证明,基于地质力学理论进行裂缝空间预测是必然的趋势,但仍存在缺陷和不足,储层地质结构非均质性强、应力条件复杂多变、构造活动期次多都影响着裂缝的发育和分布,同时局部构造和沉积特征也会改变应局部应力的再分配,使得不同组系裂缝之间相互叠加和改造,最终形成了复杂的网状系统;另外基于均质地质力学模型的裂缝成因机制分析,只考虑了各种弹性能量的相互转换,忽略了裂缝面摩擦耗能和塑形变形能,从而无法建立更合适的岩石复合破裂准则和更准确的裂缝参数力学模型,浪费了大量的计算资源,使得裂缝建模和预测遇到了瓶颈效应。
发明内容
为克服现有技术存在的缺陷,本发明提供一种致密砂岩储层复杂网状裂缝的预测方法,综合考虑储层力学非均质性和多期叠加机制,基于能量守恒定律,考虑弹性应变能、塑性耗散能、裂缝表面能和摩擦耗能,在复合破裂准则优选的前提下,建立有效的岩石破裂力学模型,在古应力场数值模拟的基础上,计算裂缝密度参数,并给出空间展布特征。
为实现上述目的,本发明采用下述方案:
致密砂岩储层多期裂缝的预测方法,步骤如下:
步骤1:建立构造地质模型及裂缝发育模型
步骤2:测试古、今地应力大小及方向
步骤3:完整岩石力学参数实验
步骤4:测试岩石力学参数及裂缝应力敏感性
步骤5:建立裂缝性岩石多级复合破裂准则
步骤6:进行岩石变形物理测试获取峰值强度
步骤7:建立单轴状态应力-应变与裂缝体密度关系模型
步骤8:建立三轴状态应力-应变与裂缝体密度、产状关系模型
步骤7:建立单轴状态应力-应变与裂缝体密度关系模型
步骤10:计算模拟现今条件下裂缝参数
步骤11:验证裂缝定量预测结果的可靠性
相对于现有技术,本发明具有如下有益效果:解决了准确获取致密砂岩储层复杂裂缝参数并进行定量表征的问题,适合于任何脆性储层为主的裂缝定量预测工作;有效预测多期应力场作用下致密砂岩储层裂缝的产状、密度,为研究裂缝性储层的有利发育区提供了可靠依据,为裂缝性储层的压裂改造设计和开发方案优化提供了保障,减少了勘探开发的风险和成本。
附图说明
图1是致密砂岩储层多期裂缝的预测方法的流程示意图。
具体实施方式
致密砂岩储层多期裂缝的预测方法,步骤如下:
步骤1:建立构造地质模型及裂缝发育模型,具体方法如下:
(1)、收集三维地震解释深度域的断层、地层数据,按砂层组或砂层级别建立现今三维构造地质模型;
(2)、分析构造演化史,确定关键造缝期,基于现今三维构造地质模型和断层相关褶皱理论和方法,建立研究区深度-声波时差关系图,获得压实规律,通过弯曲去褶皱、恢复到基准面、拼板恢复法恢复造缝时期的古构造形态,建立关键期古构造地质模型;
(3)、依据野外露头观察、岩心观察、成像测井解释和岩心CT扫描,统计裂缝产状、密度、开度、充填度、孔隙度、渗透率、力学性质参数,分析裂缝优势充填方向和共轭裂缝代表的最大主应力方向,结合构造演化史、裂缝充填物包裹体测试划分裂缝发育期次,确定造缝关键期,建立不同产状裂缝的充填程度模版,划分裂缝充填分布区域,综合建立裂缝地质发育模式;充填等级分为未充填、半充填和全充填。
步骤2:测试古、今地应力大小及方向,具体方法如下:
(1)、对钻井岩心按照不同层位进行采样,每个层位对应一组岩石样品,每组包含四块,岩石样品平行于岩心长轴方向,在实验室内加工成标准圆柱柱塞样,其长度和直径分别50mm和25mm;
(2)、选取标准柱塞样,通过岩石破裂声发射法测试凯赛尔点,获得不同期次构造运动的古最大主应力值及范围;
(3)、选取标准柱塞样,测试岩心声波速度的各向异性,确定现今最大应力、最小主应力方向,并测试不同方向的应变;
(4)、收集油田现场井壁崩落资料、储层压裂改造资料,分析获得现今地应力大小及方位;
(5)、选取每口井的标准柱塞样,在实验室中进行粘滞剩余磁性测量,提取地理北方向的信息,用于岩心裂缝定向和岩心声速定向,并校正以上测试分析获得的地应力方向;
步骤3:完整岩石力学参数实验,具体方法如下:
(1)、对钻井岩心按照不同岩性进行水平间隔采样,间隔为45°,每个岩性对应一组岩石样品,每组包含六块,岩样垂直于岩心长轴方向,在实验室内加工成标准圆柱柱塞样,其长度和直径分别50mm和25mm;
(2)、每组选取一块岩心在MTS真三轴岩石力学仪上进行单轴加载实验,测试获得岩石的单轴抗压强度、泊松比、弹性模量;
(3)、每组选取一块岩心在直剪岩石力学仪上进行直接剪切实验,测试获得岩石的单轴抗剪强度;
(4)、每组选取一块岩心在进行巴西劈裂实验,测试获得岩石的单轴抗张强度;
(5)、首先根据地层经历的构造演化史,确定岩样的历史古埋深,计算上覆地层围压范围,每隔5-10Mpa设置围压等级,每组选取三块岩心在MTS真三轴岩石力学仪上进行真三轴加载实验,测试获得岩石的三轴抗压强度、泊松比、弹性模量、内摩擦角、内聚力;
(6)、选择测井解释模型,对杨氏模量、泊松比、密度动态岩石力学参数进行测井解释,并对岩石力学实验结果进行动静态校正,以获得更真实的地层力学强度参数;
步骤4:测试岩石力学参数及裂缝应力敏感性,具体方法如下:
(1)、按照步骤3的标准,钻取岩心中含有不同产状、不同充填程度、不同数量裂缝的岩样,加工成标准柱塞样,其长度和直径分别50mm和25mm;;
(2)、按照步骤3进行岩石力学实验,获得单轴、直剪、三轴围压条件下的裂缝性岩石力学强度参数以及裂缝面摩擦系数、抗剪强度;
步骤5:建立裂缝性岩石多级复合破裂准则,具体方法如下:
(1)、根据完整岩石单轴、直剪、三轴岩石力学实验结果及应力-应变曲线,首先分析单轴破裂形态是否为剪切、拉张或张剪性质,量取剪裂角,再分析围压状态下岩石破裂的力学性质,绘出莫尔应力圆,求取包络线及斜率,斜率即内摩擦系数;
(2)、在三向压缩状态下,脆性致密砂岩破裂判据适用两段式莫尔-库仑准则:
τ=C+fσ或
式中:σ1为最大主应力,单位为MPa;σ3为最小主应力,单位为MPa;τ为剪应力,单位为MPa;α为岩石剪切破裂角,单位为弧度;为岩石内摩擦角,单位为弧度;C为岩石内聚力,单位为MPa;f为岩石内摩擦系数;莫尔-库仑准则合理给出了岩石发生剪切破裂所需的应力和剪切破裂方向,可以用于判断裂缝是否产生及裂缝的产状;
(3)、拉张应力状态下采用改进的格里菲斯准则,当(σ1+3σ3)>0时,破裂准则为:
当(σ1+3σ3)≤0时,破裂准则为:
θ=0,σ3=-σT
式中:σT为岩石的抗拉强度,单位为MPa;θ为张性破裂角,单位为度;
(4)、对于裂缝性岩体,采用各向异性拉-剪复合破裂准则:
或τ=C+f·σ
C=K·Cj+(1-K)·Cr
f=K·fj+(1-K)·fr
式中:Ft为岩石发生拉张破裂应力条件,单位为MPa;Fs为岩石无法发生拉张破裂后发生剪切破裂的应力条件,单位为Pa;Cj、fj为裂缝面内聚力和内摩擦系数;Cr、fr为基岩内聚力和内摩擦系数;K为裂缝充填率,其值位于0-1之间,以上参数通过野外和单井裂缝统计结果、岩石力学实验获得;
步骤6:进行岩石变形物理测试获取峰值强度,具体方法如下:
通过岩石物理测试定量获得研究区岩石的变形阶段,基于砂岩单轴压缩全应力应变曲线,设置好加载等级,加载等级为单轴抗压强度的10%、20%、30%、40%、50%、60%、70%、80%、90%、95%,同时采用计算机CT扫描仪器对不同受力阶段的裂缝扩展结果进行扫描成像,统计单元体内裂缝的产状、开度、体密度、长度及应变能密度,同时划分砂岩变形-破裂阶段:岩石压密阶段-岩石体积扩容阶段-新增微裂缝阶段-岩石破坏阶段-微裂缝集结产生宏观缝阶段,确定微观裂缝大量集结连通时的峰值强度σp;
步骤7:建立单轴状态应力-应变与裂缝体密度关系模型,具体方法如下:
(1)、基于岩石物理测试结果及应力-应变曲线计算总输入能量,以弹性胡克定律和能量守恒定律为前提,总输入能量=弹性应变能+裂缝摩擦耗能+残余能,其中裂缝面摩擦能占总能量的1/4-1/5,得到用主应力表示的弹性应变能Ue或应变能密度
(2)、选取平行六面表征单元体计算裂缝摩擦耗能,其边长等分别为L1、L2、L3,单位为m;沿σ1方向单元体边长为L1,沿σ2方向单元体边长为L2,沿σ3方向单元体边长为L3,单轴压缩情况下,裂纹尖端表面的剪切滑动位移为:
L=(1-μ)/GKΠ(8r/π)1/2
式中:L为剪切滑移位移,单位为m;G为剪切模量,单位为Gpa;μ为岩石泊松比;KΠ为二类强度因子,因此,岩石从破坏开始产生裂缝并受到摩擦阻力停止前的释放能主要包括裂缝表面能和摩擦耗能:
Ua=σ2acos2θtanθ/(sin2θE)+(1-u2)KΠπhσ2acos2θsin2θ/(4E)
式中:θ为岩石剪裂角,单位为弧度;σ为裂缝面正应力,单位为Mpa;a为裂缝贯穿后的长度,单位为m;
(3)、在各种能量分别表征的基础上,计算单轴压缩状态下裂缝体密度和应变能密度之间的关系:
式中:为新增裂缝表面积的摩擦能密度,单位J/m3;为裂缝产生后剩余的应变能密度,单位J/m3;J为产生单位面积裂缝所需要的总能量,单位为J/m2;为要产生裂缝必须克服的弹性应变能密度,单位J/m3;Dvf为裂缝体密度,单位为m2/m3;为系数,直接通过实验数据拟合回归法求得。
步骤8:建立三轴状态应力-应变与裂缝体密度、产状关系模型,具体方法如下:
(1)、基于单轴压缩状态下应力-应变与裂缝密度参数的定量关系表达式,通过区分不同应力状态,计算应变能密度和裂缝体密度之间的定量化关系,关键是确定公式系数,三轴挤压应力状态下,裂缝体密度计算公式:
Dvf=[σ1 2+σ2 2+σ3 2-2μ(σ1+σ2+σ3)-σp 2+2μ(σ2+σ3)σp]/(2EJ)
裂缝线密度计算公式:
(2)、有张应力存在时,裂缝体密度计算公式:
Dvf=[(σ1ε1+σ2ε2+σ3ε3)E-σt 2]/(2EJ)
裂缝线密度计算公式:
当(σ1+3σ3)>0时,裂缝线密度:
当(σ1+3σ3)≤0时,θ=0,裂缝线密度:
Dlf=Dvf
(3)、基于多级复合破裂准则,获得裂缝的走向和倾角产状数据,在大地三维坐标系中,裂缝倾角是裂缝表面与XY平面的夹角,也就是平面lx+my+nz=0与平面y=0之间的夹角αdip(0°≤αdip≤90°),其计算式为:
步骤9:求取裂缝叠加因子并定量计算模拟多期古裂缝参数,具体方法如下:
(1)、将以上得到的裂缝古密度计算公式,基于应力场模拟平台,在步骤1-4获得的古构造地质模型、岩石力学参数、应力状态的基础上,确定模型力学边界,进行力的加载,模拟古构造应力场三维分布,由程序自动判别破裂状态,给出应变能密度,并计算出第一期裂缝古密度、倾角、走向参数;
(2)、以第一期裂缝参数计算结果为基础,继续加载晚期构造作用力,产生新的构造应力场和相应裂缝参数,这时采用裂缝干扰因子,包括裂缝走向夹角、充填程度、应力状态,分两种情况进行叠加计算;
(3)、拉张应力条件下,如果新产生裂缝与原始裂缝夹角超过30°则会发生切割或终止现象,如果夹角小于30°,且裂缝半充填或未充填时,则两者叠加继续延伸,裂缝规模变大,如果夹角小于30°,且裂缝全部充填,则低角度斜交或平行原始裂缝;
(4)、压应力为主的条件下,如果新产生裂缝与原始裂缝夹角超过45°则会发生切割或终止现象,如果夹角小于45°、大于15°,且裂缝半充填或未充填时,则两者叠加继续延伸,裂缝规模变大,如果夹角小于45°,且裂缝全部充填,则低角度斜交或平行原始裂缝;
(5)、晚期裂缝与早期裂缝夹角β≤30°或15°≤β≤45°时,两期叠加后裂缝体密度计算公式为:
晚期裂缝走向与早期裂缝走向夹角β>30°或β>30°时,两期叠加后裂缝体密度计算公式为:
式中:σy1、σy2、σy3分别为早期最大、中间、最小主应力,单位为MPa;σx1、σx2、σx3分别为晚期最大、中间、最小主应力,单位为MPa;εx1、εx2、εx3分别为晚期主应变;
(6)、根据步骤1统计出的裂缝充填规律和平面分布结果,进行多期次裂缝参数叠加,通过前两期构造应力场的模拟计算,运行裂缝计算程序,产生两套单期裂缝参数,接着把后一期的古应力场及裂缝参数数据导出,叠加到前一期应力场和裂缝参数结果中,判断每一节点的充填程度、应力状态及后期最大主应力与早期裂缝走向夹角大小,调用裂缝叠加干扰因子进行判别,最终得到每一节点的叠加裂缝参数值;
步骤10:计算模拟现今条件下裂缝参数,具体方法如下:
经过多期裂缝参数叠加之后,再模拟计算现今应力场分布,通过井点应力值的校正,导出每一节点的应力值,叠加到古裂缝参数中,基于“现今应力场不能产生裂缝,仅对裂缝有效性进行改造”的思路,最终获得现今裂缝的密度参数,采用软件三维显示功能,实现裂缝倾角、走向、密度参数的三维空间显示。
步骤11:验证裂缝定量预测结果的可靠性,具体方法如下:
针对以上得到的裂缝三维分布预测结果,通过单井岩心统计、成像测井解释和CT扫描裂缝结果进行精度验证,如模拟结果与实际数据吻合度吻合度大于90%,则认为计算模拟结果可靠,否则回到步骤7,重新分析岩石力学及应力测试结果,并重新修正应力、应变与裂缝参数关系模型。
Claims (1)
1.一种致密砂岩储层多期裂缝的预测方法,其特征在于,步骤如下:
步骤1:建立构造地质模型及裂缝发育模型,具体方法如下:
(1)、收集三维地震解释深度域的断层、地层数据,按砂层组或砂层级别建立现今三维构造地质模型;
(2)、分析构造演化史,确定关键造缝期,基于现今三维构造地质模型和断层相关褶皱理论和方法,建立研究区深度-声波时差关系图,获得压实规律,通过弯曲去褶皱、恢复到基准面、拼板恢复法恢复造缝时期的古构造形态,建立关键期古构造地质模型;
(3)、依据野外露头观察、岩心观察、成像测井解释和岩心CT扫描,统计裂缝产状、密度、开度、充填度、孔隙度、渗透率、力学性质参数,分析裂缝优势充填方向和共轭裂缝代表的最大主应力方向,结合构造演化史、裂缝充填物包裹体测试划分裂缝发育期次,确定造缝关键期,建立不同产状裂缝的充填程度模版,划分裂缝充填分布区域,综合建立裂缝地质发育模式;充填等级分为未充填、半充填和全充填;
步骤2:测试古、今地应力大小及方向,具体方法如下:
(1)、对钻井岩心按照不同层位进行采样,每个层位对应一组岩石样品,每组包含四块,岩石样品平行于岩心长轴方向,在实验室内加工成标准圆柱柱塞样,其长度和直径分别50mm和25mm;
(2)、选取标准柱塞样,通过岩石破裂声发射法测试凯赛尔点,获得不同期次构造运动的古最大主应力值及范围;
(3)、选取标准柱塞样,测试岩心声波速度的各向异性,确定现今最大应力、最小主应力方向,并测试不同方向的应变;
(4)、收集油田现场井壁崩落资料、储层压裂改造资料,分析获得现今地应力大小及方位;
(5)、选取每口井的标准柱塞样,在实验室中进行粘滞剩余磁性测量,提取地理北方向的信息,用于岩心裂缝定向和岩心声速定向,并校正以上测试分析获得的地应力方向;
步骤3:完整岩石力学参数实验,具体方法如下:
(1)、对钻井岩心按照不同岩性进行水平间隔采样,间隔为45°,每个岩性对应一组岩石样品,每组包含六块,岩样垂直于岩心长轴方向,在实验室内加工成标准圆柱柱塞样,其长度和直径分别50mm和25mm;
(2)、每组选取一块岩心在MTS真三轴岩石力学仪上进行单轴加载实验,测试获得岩石的单轴抗压强度、泊松比、弹性模量;
(3)、每组选取一块岩心在直剪岩石力学仪上进行直接剪切实验,测试获得岩石的单轴抗剪强度;
(4)、每组选取一块岩心进行巴西劈裂实验,测试获得岩石的单轴抗张强度;
(5)、首先根据地层经历的构造演化史,确定岩样的历史古埋深,计算上覆地层围压范围,每隔5-10Mpa设置围压等级,每组选取三块岩心在MTS真三轴岩石力学仪上进行真三轴加载实验,测试获得岩石的三轴抗压强度、泊松比、弹性模量、内摩擦角、内聚力;
(6)、选择测井解释模型,对杨氏模量、泊松比、密度动态岩石力学参数进行测井解释,并对岩石力学实验结果进行动静态校正,以获得更真实的地层力学强度参数;
步骤4:测试岩石力学参数及裂缝应力敏感性,步骤4具体方法如下:
(1)、按照步骤3的标准,钻取岩心中含有不同产状、不同充填程度、不同数量裂缝的岩样,加工成标准柱塞样,其长度和直径分别50mm和25mm;
(2)、按照步骤3进行岩石力学实验,获得单轴、直剪、三轴围压条件下的裂缝性岩石力学强度参数以及裂缝面摩擦系数、抗剪强度;
步骤5:建立裂缝性岩石多级复合破裂准则,具体方法如下:
(1)、根据完整岩石单轴、直剪、三轴岩石力学实验结果及应力-应变曲线,首先分析单轴破裂形态是否为剪切、拉张或张剪性质,量取剪裂角,再分析围压状态下岩石破裂的力学性质,绘出莫尔应力圆,求取包络线及斜率,斜率即内摩擦系数;
(2)、在三向压缩状态下,脆性致密砂岩破裂判据适用两段式莫尔-库仑准则:
τ=C+fσ或
式中:σ1为最大主应力,单位为MPa;σ3为最小主应力,单位为MPa;τ为剪应力,单位为MPa;为岩石内摩擦角,单位为弧度;C为岩石内聚力,单位为MPa;f为岩石内摩擦系数;莫尔-库仑准则合理给出了岩石发生剪切破裂所需的应力和剪切破裂方向,直接用于判断裂缝是否产生及裂缝的产状;
(3)、拉张应力状态下采用改进的格里菲斯准则,当(σ1+3σ3)>0时,破裂准则为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>=</mo>
<mn>24</mn>
<msub>
<mi>&sigma;</mi>
<mi>T</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>cos</mi>
<mn>2</mn>
<mi>&theta;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mn>3</mn>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&sigma;</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
当(σ1+3σ3)≤0时,破裂准则为:
θ=0,σ3=-σT
式中:σT为岩石的抗拉强度,单位为MPa;θ为张性破裂角,单位为度;
(4)、对于裂缝性岩体,采用各向异性拉-剪复合破裂准则:
或τ=C+f·σ
C=K·Cj+(1-K)·Cr
f=K·fj+(1-K)·fr
式中:Ft为岩石发生拉张破裂应力条件,单位为MPa;Fs为岩石无法发生拉张破裂后发生剪切破裂的应力条件,单位为Pa;Cj、fj为裂缝面内聚力和内摩擦系数;Cr、fr为基岩内聚力和内摩擦系数;K为裂缝充填率,其值位于0-1之间,以上参数通过野外和单井裂缝统计结果、岩石力学实验获得;
步骤6:进行岩石变形物理测试获取峰值强度,具体方法如下:
通过岩石物理测试定量获得研究区岩石的变形阶段,基于砂岩单轴压缩全应力应变曲线,设置好加载等级,加载等级为单轴抗压强度的10%、20%、30%、40%、50%、60%、70%、80%、90%、95%,同时采用计算机CT扫描仪器对不同受力阶段的裂缝扩展结果进行扫描成像,统计单元体内裂缝的产状、开度、体密度、长度及应变能密度,同时划分砂岩变形-破裂阶段:岩石压密阶段-岩石体积扩容阶段-新增微裂缝阶段-岩石破坏阶段-微裂缝集结产生宏观缝阶段,确定微观裂缝大量集结连通时的峰值强度σp;
步骤7:建立单轴状态应力-应变与裂缝体密度关系模型,具体方法如下:
(1)、基于岩石物理测试结果及应力-应变曲线计算总输入能量,以弹性胡克定律和能量守恒定律为前提,总输入能量=弹性应变能+裂缝摩擦耗能+残余能,其中裂缝面摩擦能占总能量的1/4-1/5,得到用主应力表示的弹性应变能Ue或应变能密度
(2)、选取平行六面表征单元体计算裂缝摩擦耗能,其边长等分别为L1、L2、L3,单位为m;沿σ1方向单元体边长为L1,沿σ2方向单元体边长为L2,沿σ3方向单元体边长为L3,单轴压缩情况下,裂纹尖端表面的剪切滑动位移为:
L=(1-μ)/GKΠ(8r/π)1/2
式中:L为剪切滑移位移,单位为m;G为剪切模量,单位为Gpa;μ为岩石泊松比;KΠ为二类强度因子;因此,岩石从破坏开始产生裂缝并受到摩擦阻力停止前的释放能主要包括裂缝表面能和摩擦耗能:
Ua=σ2acos2θtanθ/(sin2θE)+(1-u2)KΠπhσ2acos2θsin2θ/(4E)
式中:θ为岩石剪裂角,单位为弧度;σ为裂缝面正应力,单位为Mpa;a为裂缝贯穿后的长度,单位为m;
(3)、在各种能量分别表征的基础上,计算单轴压缩状态下裂缝体密度和应变能密度之间的关系:
式中:为新增裂缝表面积的摩擦能密度,单位J/m3;J为产生单位面积裂缝所需要的总能量,单位为J/m2;为要产生裂缝必须克服的弹性应变能密度,单位J/m3;Dvf为裂缝体密度,单位为m2/m3;为系数,直接通过实验数据拟合回归法求得;
步骤8:建立三轴状态应力-应变与裂缝体密度、产状关系模型,具体方法如下:
(1)、基于单轴压缩状态下应力-应变与裂缝密度参数的定量关系表达式,通过区分不同应力状态,计算应变能密度和裂缝体密度之间的定量化关系,关键是确定公式系数,三轴挤压应力状态下,裂缝体密度计算公式:
Dvf=[σ1 2+σ2 2+σ3 2-2μ(σ1+σ2+σ3)-σp 2+2μ(σ2+σ3)σp]/(2EJ)
裂缝线密度计算公式:
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>l</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>D</mi>
<mrow>
<mi>v</mi>
<mi>f</mi>
</mrow>
</msub>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<msub>
<mi>L</mi>
<mn>3</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>-</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>-</mo>
<msub>
<mi>L</mi>
<mn>3</mn>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
</mrow>
<mrow>
<msup>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>+</mo>
<msup>
<msub>
<mi>L</mi>
<mn>3</mn>
</msub>
<mn>2</mn>
</msup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
</mfrac>
</mrow>
(2)、有张应力存在时,裂缝体密度计算公式:
Dvf=[(σ1ε1+σ2ε2+σ3ε3)E-σt 2]/(2EJ);
裂缝线密度计算公式:
当(σ1+3σ3)>0时,裂缝线密度:
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>l</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>D</mi>
<mrow>
<mi>v</mi>
<mi>f</mi>
</mrow>
</msub>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<msub>
<mi>L</mi>
<mn>3</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>-</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>-</mo>
<msub>
<mi>L</mi>
<mn>3</mn>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
</mrow>
<mrow>
<msup>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>+</mo>
<msup>
<msub>
<mi>L</mi>
<mn>3</mn>
</msub>
<mn>2</mn>
</msup>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
</mrow>
</mfrac>
</mrow>
当(σ1+3σ3)≤0时,θ=0,裂缝线密度:
Dlf=Dvf
(3)、基于多级复合破裂准则,获得裂缝的走向和倾角产状数据,在大地三维坐标系中,裂缝倾角是裂缝表面与XY平面的夹角,也就是平面lx+my+nz=0与平面y=0之间的夹角αdip(0°≤αdip≤90°),其计算式为:
步骤9:建立单轴状态应力-应变与裂缝体密度关系模型,具体方法如下:
(1)、将以上得到的裂缝古密度计算公式,基于应力场模拟平台,在步骤1-4获得的古构造地质模型、岩石力学参数、应力状态的基础上,确定模型力学边界,进行力的加载,模拟古构造应力场三维分布,由程序自动判别破裂状态,给出应变能密度,并计算出第一期古裂缝密度、倾角、走向参数;
(2)、以第一期裂缝参数计算结果为基础,继续加载晚期构造作用力,产生新的构造应力场和相应裂缝参数,这时采用裂缝干扰因子,包括裂缝走向夹角、充填程度、应力状态,分两种情况进行叠加计算;
(3)、拉张应力条件下,如果新产生裂缝与原始裂缝夹角超过30°则会发生切割或终止现象,如果夹角小于30°,且裂缝半充填或未充填时,则两者叠加继续延伸,裂缝规模变大,如果夹角小于30°,且裂缝全部充填,则低角度斜交或平行原始裂缝;
(4)、压应力为主的条件下,如果新产生裂缝与原始裂缝夹角超过45°则会发生切割或终止现象,如果夹角大于15°、小于45°,且裂缝半充填或未充填时,则两者叠加继续延伸,裂缝规模变大,如果夹角小于45°,且裂缝全部充填,则低角度斜交或平行原始裂缝;
(5)、晚期裂缝与早期裂缝夹角β≤30°或15°≤β≤45°时,两期叠加后裂缝体密度计算公式为:
晚期裂缝走向与早期裂缝走向夹角β>30°或β>30°时,两期叠加后裂缝体密度计算公式为:
式中:Dvfy为早期裂缝体密度,单位为m2/m3;Dvfx为晚期裂缝体密度,单位为m2/m3;σy1、σy2、σy3分别为早期最大、中间、最小主应力,单位为MPa;σx1、σx2、σx3分别为晚期最大、中间、最小主应力,单位为MPa;εx1、εx2、εx3分别为晚期主应变;
(6)、根据步骤1统计出的裂缝充填规律和平面分布结果,进行多期次裂缝参数叠加,通过前两期构造应力场的模拟计算,运行裂缝计算程序,产生两套单期裂缝参数,接着把后一期的古应力场及裂缝参数数据导出,叠加到前一期应力场和裂缝参数结果中,判断每一节点的充填程度、应力状态及后期最大主应力与早期裂缝走向夹角大小,调用裂缝叠加干扰因子进行判别,最终得到每一节点的叠加裂缝参数值;
步骤10:计算模拟现今条件下裂缝参数,具体方法如下:经过多期裂缝参数叠加之后,再模拟计算现今应力场分布,通过井点应力值的校正,导出每一节点的应力值,叠加到古裂缝参数中,基于“现今应力场不能产生裂缝,仅对裂缝有效性进行改造”的思路,最终获得现今裂缝的密度参数,采用软件三维显示功能,实现裂缝倾角、走向、密度参数的三维空间显示;
步骤11:验证裂缝定量预测结果的可靠性,具体方法如下:针对以上得到的裂缝三维分布预测结果,通过单井岩心统计、成像测井解释和CT扫描裂缝结果进行精度验证,如模拟结果与实际数据吻合度吻合度大于90%,则认为计算模拟结果可靠,否则回到步骤7,重新分析岩石力学及应力测试结果,并重新修正应力、应变与裂缝参数关系模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510752304.8A CN105319603B (zh) | 2015-11-06 | 2015-11-06 | 致密砂岩储层复杂网状裂缝的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510752304.8A CN105319603B (zh) | 2015-11-06 | 2015-11-06 | 致密砂岩储层复杂网状裂缝的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105319603A CN105319603A (zh) | 2016-02-10 |
CN105319603B true CN105319603B (zh) | 2018-04-17 |
Family
ID=55247405
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510752304.8A Active CN105319603B (zh) | 2015-11-06 | 2015-11-06 | 致密砂岩储层复杂网状裂缝的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105319603B (zh) |
Families Citing this family (65)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105842754B (zh) * | 2016-03-24 | 2019-05-07 | 中国石油天然气股份有限公司 | 一种定量确定超深层致密砂岩裂缝有效性的方法 |
CN105866835B (zh) * | 2016-03-28 | 2018-11-02 | 中国石油大学(华东) | 一种基于地应力分布的断层三维封闭性定量评价方法 |
CN107478768B (zh) * | 2016-06-08 | 2019-10-11 | 中国石油化工股份有限公司 | 一种泥页岩古抗剪强度确定方法 |
CN106202737B (zh) * | 2016-07-14 | 2019-08-06 | 中国石油天然气股份有限公司 | 一种致密储层体积压裂适应区筛选方法及装置 |
CN106227957A (zh) * | 2016-07-27 | 2016-12-14 | 长江大学 | 等效裂缝建模的方法 |
CN106290002B (zh) * | 2016-08-03 | 2019-03-12 | 中国矿业大学 | 基于三点弯曲试验的岩石ⅰ型裂纹扩展全过程检测方法 |
CN106404918B (zh) * | 2016-10-09 | 2019-07-02 | 中国矿业大学 | 基于受载煤岩体全应力应变曲线的声发射能量计算方法 |
CN107991188B (zh) * | 2016-10-26 | 2020-04-28 | 中国石油化工股份有限公司 | 一种基于岩心残余应力水平预测水力裂缝复杂性的方法 |
CN106844858B (zh) * | 2016-12-21 | 2020-02-18 | 中国石油天然气股份有限公司 | 地层裂缝发育区带预测方法及装置 |
CN108254262B (zh) * | 2016-12-28 | 2021-03-30 | 中国石油天然气股份有限公司 | 岩石层理裂缝剪切参数预测方法及装置 |
CN107092719B (zh) * | 2017-03-17 | 2020-09-04 | 中国石油天然气股份有限公司 | 注水优势通道识别和微球封堵粒径选择的方法及装置 |
CN106959468A (zh) * | 2017-03-20 | 2017-07-18 | 安徽理工大学 | 一种含破碎带断层模型的建立方法 |
CN107145633B (zh) * | 2017-04-07 | 2018-04-17 | 中国地质大学(武汉) | 一种岩石裂隙网络产状三维统计分布的预测方法 |
CN107091768B (zh) * | 2017-04-24 | 2019-09-03 | 长春师范大学 | 一种测定材料弹性模量的方法 |
CN108957527B (zh) * | 2017-05-27 | 2020-06-12 | 中国石油化工股份有限公司 | 岩层网状裂缝的地震预测方法 |
CN108984804B (zh) * | 2017-05-31 | 2023-05-12 | 中国石油化工股份有限公司 | 一种利用裂缝发育密度评价裂缝性储层质量的方法 |
CN109209356B (zh) * | 2017-07-06 | 2021-08-31 | 中国石油化工股份有限公司 | 一种基于张破裂和剪破裂确定地层可压裂性的方法 |
CN107462190B (zh) * | 2017-07-31 | 2018-06-22 | 中国科学院地质与地球物理研究所 | 一种岩石水力压裂试验裂缝三维形貌高精度成像方法 |
CN109492238B (zh) * | 2017-09-12 | 2022-03-01 | 中国石油天然气股份有限公司 | 一种应用地应力分析断层和裂缝渗流特性的方法 |
CN107957596B (zh) * | 2017-12-20 | 2019-09-06 | 中国石油大学(华东) | 一种基于应力释放率的储层裂缝预测方法 |
US10947841B2 (en) * | 2018-01-30 | 2021-03-16 | Baker Hughes, A Ge Company, Llc | Method to compute density of fractures from image logs |
CN108468538B (zh) * | 2018-03-23 | 2020-04-21 | 西南石油大学 | 一种页岩水压裂缝扩展预测方法 |
CN108590601B (zh) * | 2018-04-08 | 2020-10-23 | 西南石油大学 | 一种注水扩容膨胀施工参数优选的实验方法 |
CN108732092A (zh) * | 2018-05-30 | 2018-11-02 | 山东科技大学 | 一种非透明试件三维裂纹扩展过程的综合反演方法 |
CN108760619B (zh) * | 2018-06-07 | 2019-11-19 | 北京航空航天大学 | 一种软包电池密封寿命的预测方法及系统 |
CN109063250B (zh) * | 2018-06-27 | 2023-04-07 | 青岛理工大学 | 一种混凝土裂缝扩展方向的分析与预测方法 |
CN108843329B (zh) * | 2018-06-27 | 2021-08-20 | 西安智财全技术转移中心有限公司 | 基于围岩分区破裂演化分析的隧道开挖施工方法 |
CN108897066B (zh) * | 2018-08-23 | 2020-01-07 | 中国石油天然气股份有限公司 | 碳酸盐岩裂缝密度定量预测方法及装置 |
CN109459497B (zh) * | 2018-09-14 | 2019-08-16 | 中国石油大学(华东) | 一种含平行裂缝横向各向同性岩石的声波速度计算方法 |
CN109268004B (zh) * | 2018-09-21 | 2020-06-26 | 北京科技大学 | 一种页岩气藏介质结构耦合及缝网形态判识方法 |
CN109469477B (zh) * | 2018-10-18 | 2022-08-02 | 中国海洋石油集团有限公司 | 一种人工裂缝延伸方向的预测方法和装置 |
CN109339759B (zh) * | 2018-10-19 | 2019-12-24 | 中国科学院地质与地球物理研究所 | 高能ct超深层钻井压裂一体化工程试验设备 |
CN109902422B (zh) * | 2019-03-08 | 2021-07-27 | 西南石油大学 | 一种井眼轨迹选取方法及装置 |
CN110031304B (zh) * | 2019-04-25 | 2020-04-21 | 四川大学 | 一种岩石变形破坏模式预测方法 |
CN111859603B (zh) * | 2019-04-30 | 2022-11-01 | 中国石油天然气股份有限公司 | 砂岩裂缝改造方式的评定方法和评定装置 |
CN110286421B (zh) * | 2019-08-09 | 2020-12-18 | 中国石油大学(华东) | 一种致密砂岩储层天然裂缝建模方法 |
CN110850057B (zh) * | 2019-11-15 | 2021-04-09 | 成都理工大学 | 一种基于自相似性理论的储层裂缝建模方法及系统 |
CN110866300B (zh) * | 2019-11-15 | 2022-11-25 | 上海环联生态科技有限公司 | 大型建筑的裂缝预测方法 |
CN111062138B (zh) * | 2019-12-23 | 2021-10-22 | 中国石油大学(北京) | 致密低渗透储层天然裂缝分布样式预测方法及系统 |
CN112145165B (zh) * | 2019-12-26 | 2024-04-23 | 中国海洋石油集团有限公司 | 一种微裂缝-孔隙型储层动静态渗透率转换方法 |
CN111189755B (zh) * | 2020-01-07 | 2022-04-19 | 重庆大学 | 煤矿采动稳定区瓦斯有效储存空间的数值模拟获取方法 |
CN111271055B (zh) * | 2020-02-26 | 2021-10-08 | 中国石油大学(北京) | 页岩的脆性指数确定方法、装置和设备 |
CN111504872B (zh) * | 2020-04-16 | 2021-02-19 | 武汉大学 | 变开度可拆卸的仿真裂隙试验装置及试验方法 |
CN111580166B (zh) * | 2020-05-21 | 2022-03-18 | 中国石油大学(华东) | 一种基于声波远探测和岩石力学的裂缝有效性评价方法 |
US11725483B2 (en) | 2020-06-12 | 2023-08-15 | Saudi Arabian Oil Company | Method and system of fracability measurement based on core fracture density |
CN112065351B (zh) * | 2020-08-25 | 2021-09-17 | 中国石油大学(北京) | 水力裂缝内暂堵体信息的一体化确定方法、装置、设备 |
CN111767323B (zh) * | 2020-09-02 | 2020-12-04 | 西南石油大学 | 一种基于核smo算法的地层裂缝长度实时预测方法 |
CN112329219B (zh) * | 2020-10-26 | 2024-01-26 | 中国科学院力学研究所 | 一种计算巴西劈裂实验中含微孔和微裂缝岩石拉伸损伤区域的方法 |
CN114441301A (zh) * | 2020-11-05 | 2022-05-06 | 中国石油天然气股份有限公司 | 基于构造变形物理模拟的裂缝发育特征分析方法及装置 |
CN112593932B (zh) * | 2020-12-15 | 2022-09-16 | 中国石油大学(华东) | 一种考虑随机裂缝分布的破碎性井眼岩样加工方法 |
CN113281148B (zh) * | 2021-05-19 | 2022-09-09 | 江西理工大学 | 一种岩石蠕变破坏微裂纹贯通阶段起始时间辨识方法 |
CN113341465B (zh) * | 2021-06-11 | 2023-05-09 | 中国石油大学(北京) | 方位各向异性介质的地应力预测方法、装置、介质及设备 |
CN113534291B (zh) * | 2021-07-20 | 2023-02-07 | 中国石油大学(华东) | 岩石力学层约束下的低渗透储层不同尺度裂缝定量预测方法 |
CN113945457B (zh) * | 2021-10-14 | 2023-05-26 | 辽宁科技大学 | 一种分析岩石在复杂卸荷应力条件下破坏机制的方法 |
CN114324002B (zh) * | 2021-11-19 | 2022-06-21 | 中煤科工能源科技发展有限公司 | 一种复杂应力状态下岩石裂隙产状及发育程度的确定方法 |
CN114088516B (zh) * | 2021-11-23 | 2023-07-18 | 贵州大学 | 一种基于岩心损伤各向异性的单孔声发射地应力测量方法 |
CN114235579B (zh) * | 2021-12-17 | 2023-12-19 | 广东石油化工学院 | 天然气水合物岩心抗压、抗剪强度的三轴测试装置及方法 |
CN114778800B (zh) * | 2022-04-28 | 2023-08-08 | 中交第一公路勘察设计研究院有限公司 | 一种基于解析方法的多因素岩爆预测方法 |
CN114622907B (zh) * | 2022-05-12 | 2022-07-19 | 中国科学院地质与地球物理研究所 | 油气储层最小水平主应力剖面模型预测方法、系统和设备 |
CN115079261B (zh) * | 2022-06-06 | 2023-07-11 | 吉林大学 | 一种基于多参数定量解释模板的致密砂岩气储层评价方法 |
CN115267905B (zh) * | 2022-07-28 | 2023-12-19 | 中国石油大学(华东) | 一种复杂构造区钻井工程中裂缝性井漏预测方法 |
CN115856996B (zh) * | 2022-11-08 | 2023-06-30 | 应急管理部国家自然灾害防治研究院 | 一种地震数值预测预报方法和系统 |
CN115906567A (zh) * | 2022-11-21 | 2023-04-04 | 中国矿业大学 | 一种考虑厚度变化的储层构造裂隙模拟预测方法 |
CN115618658B (zh) * | 2022-12-16 | 2023-03-14 | 新疆石油管理局有限公司 | 基于弹塑性本构的地层破裂压力动态预测方法及装置 |
CN115640733B (zh) * | 2022-12-21 | 2023-03-07 | 武汉中旺亿能科技发展有限公司 | 基于全直径岩心ct扫描致密储层压裂缝网数值模拟方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104331745A (zh) * | 2014-10-17 | 2015-02-04 | 成都理工大学 | 油气藏内天然裂缝的分期、分成因预测评价方法 |
CN104500050A (zh) * | 2014-10-31 | 2015-04-08 | 中国石油大学(华东) | 一种裂缝性储层渗透率张量及各向异性定量预测方法 |
CN104730596A (zh) * | 2015-01-25 | 2015-06-24 | 中国石油大学(华东) | 一种基于多尺度因素约束的离散裂缝建模方法 |
-
2015
- 2015-11-06 CN CN201510752304.8A patent/CN105319603B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104331745A (zh) * | 2014-10-17 | 2015-02-04 | 成都理工大学 | 油气藏内天然裂缝的分期、分成因预测评价方法 |
CN104500050A (zh) * | 2014-10-31 | 2015-04-08 | 中国石油大学(华东) | 一种裂缝性储层渗透率张量及各向异性定量预测方法 |
CN104730596A (zh) * | 2015-01-25 | 2015-06-24 | 中国石油大学(华东) | 一种基于多尺度因素约束的离散裂缝建模方法 |
Non-Patent Citations (2)
Title |
---|
储层裂缝预测方法研究;崔健;《重庆科技学院学报(自然科学版)》;20080229;第10卷(第1期);5-8 * |
致密砂岩储层构造裂缝形成机制及定量预测研究进展;徐会永 等;《地质力学学报》;20131231;第19卷(第4期);377-381 * |
Also Published As
Publication number | Publication date |
---|---|
CN105319603A (zh) | 2016-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105319603B (zh) | 致密砂岩储层复杂网状裂缝的预测方法 | |
Wang et al. | Numerical analysis on scale effect of elasticity, strength and failure patterns of jointed rock masses | |
Liu et al. | Experimental investigation of the effects of heterogeneity and geostress difference on the 3D growth and distribution of hydrofracturing cracks in unconventional reservoir rocks | |
CN104361211B (zh) | 岩石统计损伤本构模型的构建和应用方法 | |
CN113534291B (zh) | 岩石力学层约束下的低渗透储层不同尺度裂缝定量预测方法 | |
Liu et al. | Simulation of paleotectonic stress fields and quantitative prediction of multi-period fractures in shale reservoirs: a case study of the Niutitang Formation in the Lower Cambrian in the Cen'gong block, South China | |
CN113820750B (zh) | 基于弹塑性力学定量预测泥岩构造裂缝的方法 | |
Shimizu et al. | A study of the effect of brittleness on hydraulic fracture complexity using a flow-coupled discrete element method | |
Liu et al. | Asymmetric propagation mechanism of hydraulic fracture networks in continental reservoirs | |
CN110864966A (zh) | 一种适用于裂缝型岩石的可压性评价方法 | |
CN104865610A (zh) | 一种深部储层岩心的数值表征方法 | |
Zhao et al. | Numerical study of the effect of natural fractures on shale hydraulic fracturing based on the continuum approach | |
Li et al. | Numerical analysis of the hydrofracturing behaviours and mechanisms of heterogeneous reservoir rock using the continuum-based discrete element method considering pre-existing fractures | |
CN109752262B (zh) | 一种基于原位相对密度确定覆盖层土体动剪模量参数的方法 | |
CN115324557A (zh) | 基于多因素分析预测压裂诱发套管变形风险程度的方法 | |
CN106842316B (zh) | 裂缝确定方法和装置 | |
Fu et al. | Prediction of tectonic fractures in coal reservoirs using geomechanical method | |
CHEN et al. | 3D geological modeling of dual porosity carbonate reservoirs: A case from the Kenkiyak pre-salt oilfield, Kazakhstan | |
Yan et al. | A multiscale analysis of adjacent fault dislocation mechanism induced by tunnel excavation based on continuous-discrete coupling method | |
Feng et al. | Quantitative prediction of 3-D multiple parameters of tectonic fractures in tight sandstone reservoirs based on geomechanical method | |
CN107942381B (zh) | 一种致密油储层层理缝定量预测方法 | |
Chang et al. | Strength anisotropy of jointed rock slope considering mining damage: a numerical analysis method | |
Chen et al. | Numerical modeling of fracture process using a new fracture constitutive model with applications to 2D and 3D engineering cases | |
CN106844858B (zh) | 地层裂缝发育区带预测方法及装置 | |
Pierce et al. | Synthetic rock mass applications in mass mining |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |