CN110969699B - 一种高精度三维地应力建模方法及其应用 - Google Patents
一种高精度三维地应力建模方法及其应用 Download PDFInfo
- Publication number
- CN110969699B CN110969699B CN201811153242.9A CN201811153242A CN110969699B CN 110969699 B CN110969699 B CN 110969699B CN 201811153242 A CN201811153242 A CN 201811153242A CN 110969699 B CN110969699 B CN 110969699B
- Authority
- CN
- China
- Prior art keywords
- stress
- maximum
- rock
- minimum horizontal
- dimensional
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供本发明提供一种高精度三维地应力建模方法及其应用,所述方法以岩石物理建模为基础计算高精度的岩石力学,采用有效应力比值法、应力多边形约束及图像中心点自动分析,提高崩落深度点最大、最小水平主应力的计算精度;结合三维地震速度及三维相控叠前反演,采用弹簧模型计算三维地应力,为三维地应力的进行高精度定量建模,本发明具有算法实现简单、模型精度高和创新性的特点,能方便快捷应用于油气勘探开发中的地应力分析。
Description
技术领域
本发明属于地震信息优化处理与解释技术领域,具体涉及一种高精度三维地应力建模方法及其应用。
背景技术
地应力研究与分析是油气勘探开发中一项基础而且非常重要的研究工作,广泛用于解决诸如油气富集区分布、井壁稳定性、储层改造裂缝分布、水力压裂起裂压力与裂缝延伸压力、套管变形、工程设计优化等问题。在石油工程中,地应力方位往往利用成像测井、偶极快慢横波或者多臂井径资料确定,地应力的方位研究相对简单。
地应力大小包括上覆岩层压力、最大水平主应力、最小水平主应力,最大、最小水平主应力测量方法主要有水压致裂法、声发射Kaiser效应法及差应变法等,此类方法通过岩心或现场试验直接测量离散岩石的地应力大小,为了获得地应力剖面,Zoback提出了有效应力比值法,该方法认为最大或者最小水平主应力有效应力比接近常数,采用有效应力比值法,可以获得基本与岩性无关的连续应力剖面,采用有效应力比值法计算的应力剖面不能有效的指导压裂设计。
学者们更多的考虑利用测井资料和地应力组合弹簧模型,计算最大最小水平主应力,组合弹簧模型的关键是求取构造应力系数,谭强等人利用钻井中的地漏实验(LOT),采用水力压裂井壁诱发拉伸裂缝的条件,简化最大水平主应力公式并求得最大最小水平主应力,以此为基础计算构造应力系数,此种用水力压裂确定最大水平主应力的方法,在绝大多数环境下不能使用,Aadnoy等给出了由崩落的椭圆井眼长短半轴比值,结合最小水平主应力来计算最大水平主应力,张美玲等在Aadnoy基础上,采用综合计算方法计算最大最小水平主应力,并在套损预测中进行了应用,该方法由于采用摩尔库仑准测,适用于脆性较好的地层,王浩等人利用地应力计算结果结合井壁破坏模式合成井壁破坏图像,将合成的图像与电成像测井资料进行对比,当二者反映的井壁破坏情况相符合,此时的水平构造应力系数可作为该区的水平构造应力系数。
对于三维地应力研究,目前多采用ANSYS或者abaqus等有限元分析软件,建立三维地应力场,在实际生产中推广难度较大。
发明内容
为了解决上述问题,本发明提供一种高精度三维地应力建模方法,所述方法以岩石物理建模为基础计算高精度的岩石力学,采用有效应力比值法、应力多边形约束及图像中心点自动分析,提高崩落深度点最大、最小水平主应力的计算精度;结合三维地震速度及三维相控叠前反演,采用弹簧模型计算三维地应力,为三维地应力的进行高精度定量建模;
进一步地,所述方法包括以下步骤;
1):进行井眼扩径校正、曲线标准化、测井曲线Backus平均降频预处理,并根据不同地质环境进行不同方法的岩石物理建模,获得骨架体积模量Ks和剪切模量μs,干岩石体积模量Kdry和剪切模量μdry,及原状地层的体积模量Ksat和剪切模量μsat;
1):根据步骤1)中结果)计算横波速度Vs、杨氏模量E、泊松比v、Biot系数α=1-Kdry/Ks、岩石内摩擦角Φ,采用C0=YM/N*Coeff公式计算单轴抗压强度,其中Coeff为实验室实测C0与测井计算C0的拟合校正系数;
2):根据测井密度计算上覆岩层压力Sv和声波速度,采用Eaton公式、Fillippone公式或者Bowers公式或者岩石物理法,计算地层孔隙压力Pp,其中对Fillippone公式Pp=(Vmax-Vi)/(Vmax-Vmin)*Sv*Coeff,式中Vmax为骨架速度,Vi为流体速度,采用岩石物理建模方法计算Vmax和Vi,Coeff为实测孔隙压力值与预测孔隙压力值的拟合校正系数;
3):根据水力压裂法获得水力压裂点的最小水平主应力σh,采用有效应力比值ESR=(σh-Pp)/(Sv-Pp),确定有效应力比ESR,计算连续的最小水平主应力σh;
4):根据成像测井分析崩落宽度WB0或者钻井过程诱导拉张缝,上覆岩层压力Sv、孔隙压力Pp及根据有效应力比值计算的最小水平主应力σh构建应力多边形;
5):分析重叠图像区域,采用图像分析算法,自动求取多边形图像中心,中心值即为该深度点最终的最大、最小水平主应力,然后对最大、最小水平主应力均采用有效应力比值法,重新计算最大、最小水平主应力连续剖面,采用同样的做法计算其他崩落深度的最大、最小水平主应力,最终获得各崩落深度点的最大最小水平主应力;
7):将各崩落点确定的最大、最小水平主应力带入地应力弹簧模型中,反演构造应变,计算高分辨率最大、最小水平主应力;
式中,σH为最大水平主应力,MPa;σh为最小水平主应力,MPa;Sv为上覆地层压力,MPa;Pp为地层孔隙压力,Mpa;α为Biot系数;E为岩石静态杨氏模量,MPa;v为岩石静态泊松比;εH为最大水平应力方向构造应变;εh为最小水平应力方向构造应变;
8):根据三维地震速度和相控地震叠前反演,计算三维上覆岩层压力Sv、孔隙压力Pp、杨氏模量YM和泊松比v,结合反演的构造应变,带入应力弹簧模型中,得到高精度的三维地应力数据体;
进一步地,所述步骤1)中根据不同地质环境进行的岩石物理建模的方法包括Xu-White及相关修正模型;
进一步地,所述步骤2)中单轴抗压强度C0的计算方法是从杨氏模量和泊松比的物理意义出发,采用C0=YM/N*Coeff公式计算单轴抗压强度,其中Coeff为实验室实测C0与测井计算C0的拟合校正系数,采用岩石物理建模方法,获得诸如Biot系数、单轴抗压强度等岩石力学参数,同时提高岩石力学参数分析精度;
进一步地,所述步骤6)中根据岩石破裂准则对单轴抗压强度进行调整,崩落区域与拉张区域重叠,根据重叠区域进一步约束了岩石破裂准则、单轴抗压强度及最大、最小水平主应力;
进一步地,所述岩石破裂准则为修正Lade准则;
进一步地,一种高精度三维地应力建模方法的应用,所述方法结合三维地震速度及三维地震相控叠前反演,采用弹簧模型计算三维地应力,方便快捷地应用于石油天然气勘探开发中的地应力分析;
本发明的有益效果如下:
1):采用有效应力比值法、应力多边形约束及图像中心点自动分析等,大幅度提高了崩落深度点最大、最小水平主应力的计算精度;
2):结合三维地震速度及三维地震相控叠前反演,采用弹簧模型计算三维地应力,给三维地应力研究开辟了一套新的思路;
3):本发明具有算法实现简单、模型精度高和创新性的特点,能方便快捷应用于油气勘探开发中的地应力分析。
附图说明
图1为本发明所述方法中成像测井与应力多边形约束图(错误的参数);
图2为本发明所述方法中应力多边形约束最大最小水平主应力图(正确的参数);
图3为本发明所述方法中应用弹簧模型计算的高分辨率最大、最小水平主应力;
图4为本发明所述方法中三维最小水平主应力图;
图5为本发明所述方法中三维最大水平主应力图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细描述。应当理解,此处所描述的具体实施例仅仅用于解释本发明,并不用于限定本发明。相反,本发明涵盖任何由权利要求定义的在本发明的精髓和范围上做的替代、修改、等效方法以及方案。进一步,为了使公众对本发明有更好的了解,在下文对本发明的细节描述中,详尽描述了一些特定的细节部分。对本领域技术人员来说没有这些细节部分的描述也可以完全理解本发明。
下面结合附图和具体实施例对本发明作进一步说明,但不作为对本发明的限定。下面为本发明的举出最佳实施例:
如图1-图5所示,本发明提供一种高精度三维地应力建模方法及其应用,本发明所述方法以岩石物理建模为基础,计算高精度的岩石力学参数;采用有效应力比值法、应力多边形约束及图像中心点自动分析等,提高了崩落深度点最大、最小水平主应力的计算精度;结合三维地震速度及三维相控叠前反演,采用弹簧模型计算三维地应力,为三维地应力的高精度定量研究,提供了一种新的思路。
具体方法步骤如下:
1)对测井曲线预处理,包括进行井眼扩径校正、曲线标准化、测井曲线Backus平均降频等预处理,提高测井数据质量,在此基础上开展岩石物理建模,针对不同的地质环境,选用不同的建模方法,比如Xu-White及相关修正模型,岩石物理建模可以获得骨架体积模量Ks和剪切模量μs,干岩石体积模量Kdry和剪切模量μdry,及原状地层的体积模量Ksat和剪切模量μsat,进一步计算横波速度Vs、杨氏模量E、泊松比v、Biot系数α=1-Kdry/Ks、岩石内摩擦角φ等岩石力学参数,对于单轴抗压强度C0,传统做法是分岩性统计C0与孔隙度、声波速度等关系,本发明从杨氏模量和泊松比的物理意义出发,采用C0=YM/N*Coeff公式计算单轴抗压强度,其中Coeff为实验室实测C0与测井计算C0的拟合校正系数,采用岩石物理建模方法,可以非常容易获得诸如Biot系数、单轴抗压强度等岩石力学参数,同时提高岩石力学参数分析精度;
2)根据测井密度计算上覆岩层压力Sv;根据声波速度,采用Eaton公式、Fillippone公式或者Bowers公式或者岩石物理法等孔隙压力计算公式,计算地层孔隙压力Pp,其中对Fillippone公式Pp=(Vmax-Vi)/(Vmax-Vmin)*Sv*Coeff,式中Vmax为骨架速度,Vi为流体速度,传统做法让Vmax和Vi为固定值,或者建立Vmax和Vi与速度趋势的关系,本发明采用岩石物理建模方法计算Vmax和Vi,Coeff为实测孔隙压力值与预测孔隙压力值的拟合校正系数;
3)根据水力压裂法获得水力压裂点的最小水平主应力σh,采用有效应力比值ESR=(σh-Pp)/(Sv-Pp),确定有效应力比ESR,计算连续的最小水平主应力σh;根据成像测井分析崩落宽度WB0或者钻井过程诱导拉张缝,上覆岩层压力Sv、孔隙压力Pp及根据有效应力比值计算的最小水平主应力σh等构建应力多边形,若地层同时发育崩落和拉张缝,可以用应力多边形同时约束岩石破裂准则、单轴抗压强度及最大、最小水平主应力及最大水平主应力σH;为了阐述本方案,图1为某深度段同时发育崩落与拉张缝,如果单轴抗压强度不对,或者岩石破裂准则不对,在应力多边形中,崩落的区域与拉张的区域将不重合,这与事实不符合,图2对单轴抗压强度及岩石破裂准则(修正Lade准则)进行调整,崩落区域与拉张区域重叠,根据重叠区域进一步约束了岩石破裂准则、单轴抗压强度及最大、最小水平主应力。
4)分析图2中的重叠图像区域,采用图像分析算法,自动求取多边形图像中心,中心值即为该深度点最终的最大、最小水平主应力,然后对最大、最小水平主应力均采用有效应力比值法,重新计算最大、最小水平主应力连续剖面,采用同样的做法计算其他崩落深度的最大、最小水平主应力,最终获得各崩落深度点的最大最小水平主应力;
5)将各崩落点确定的最大、最小水平主应力带入地应力弹簧模型中,反演构造应变,计算高分辨率最大、最小水平主应力;
式中,σH为最大水平主应力,MPa;σh为最小水平主应力,MPa;Sv为上覆地层压力,MPa;Pp为地层孔隙压力,Mpa;α为Biot系数;E为岩石静态杨氏模量,MPa;v为岩石静态泊松比;εH为最大水平应力方向构造应变;εh为最小水平应力方向构造应变。
图3为根据前述步骤,计算的各崩落深度点或者拉张裂缝深度点,反演弹簧模型中的应变量,得到的高分辨率的最大、最小水平主应力,分别如图4和图5所示。
6)根据三维地震速度和相控地震叠前反演,计算三维上覆岩层压力Sv、孔隙压力Pp、杨氏模量YM和泊松比v,结合反演的构造应变,带入应力弹簧模型中,得到高精度的三维地应力数据体。
本发明的关键技术包括:
1、对测井曲线预处理提高测井数据质量,在此基础上开展岩石物理建模,提高横波速度Vs、杨氏模量E、泊松比v、Biot系数α、单轴抗压强度C0及岩石内摩擦角φ等岩石力学参数分析精度;
2、根据测井密度计算上覆岩层压力Sv,根据声波速度Vp,采用Eaton公式、Fillippone公式、Bowers公式或者岩石物理公式等计算地层孔隙压力Pp,并用实测点进行标定;
3、根据水力压裂法获得水力压裂点的最小水平主应力σh,采用有效应力比值ESR=(σh-Pp)/(Sv-Pp),确定有效应力比ESR,计算连续的初始最小水平主应力σh;根据成像测井分析崩落宽度WB0或者钻井过程诱导拉张缝,上覆岩层压力Sv、孔隙压力Pp及根据有效应力比值计算的最小水平主应力σh等构建应力多边形,若地层同时发育崩落和拉张缝,在应力多边形里面可以同时约束岩石破裂准则、单轴抗压强度及最大、最小水平主应力,约束最大水平主应力σH;
4、采用图像分析算法,自动求取最大、最小水平主应力多边形图像中心,中心值为该崩落深度点最终的最大、最小水平主应力,然后对最大、最小水平主应力均采用有效应力比值法,重新计算最大最小水平主应力连续剖面,采用同样的做法计算其他崩落深度的最大最小水平主应力,最终获得各崩落深度点的最大、最小水平主应力;
5、将各崩落点计算的最大、最小水平主应力带入地应力弹簧模型中,反演构造应变,最终获得高分辨率最大、最小水平主应力连续剖面,前述步骤为单井高分辨率地应力建模过程;
6、根据前述的模型参数,结合三维地震速度和三维地震相控叠前反演,计算三维上覆岩层压力Sv、孔隙压力Pp、杨氏模量E和泊松比v,结合反演的构造应变,得到高精度的三维地应力数据。
本发明以精细岩石物理建模为基础,采用有效应力比值法、应力多边形约束及图像中心点自动分析等,大幅度提高了崩落深度点最大、最小水平主应力的计算精度,结合三维地震速度及三维地震相控叠前反演,采用弹簧模型计算三维地应力,给三维地应力研究开辟了一套新的思路,本发明具有算法实现简单、模型精度高和创新性的特点,能方便快捷应用于油气勘探开发中的地应力分析
以上所述的实施例,只是本发明较优选的具体实施方式的一种,本领域的技术人员在本发明技术方案范围内进行的通常变化和替换都应包含在本发明的保护范围内。
Claims (2)
1.一种高精度三维地应力建模方法,其特征在于,所述方法以岩石物理建模为基础计算高精度的岩石力学,采用有效应力比值法、应力多边形约束及图像中心点自动分析,提高崩落深度点最大、最小水平主应力的计算精度;结合三维地震速度及三维相控叠前反演,采用弹簧模型计算三维地应力,为三维地应力的进行高精度定量建模;
所述方法包括以下步骤:
1):进行井眼扩径校正、曲线标准化、测井曲线Backus平均降频预处理,并根据不同地质环境进行不同方法的岩石物理建模,获得骨架体积模量Ks和剪切模量μs,干岩石体积模量Kdry和剪切模量μdry,及原状地层的体积模量Ksat和剪切模量μsat;
2):根据步骤1)中结果计算横波速度Vs、杨氏模量E、泊松比v、Biot系数α=1-Kdry/Ks、岩石内摩擦角Φ,采用CO=YM/N*Coeff公式计算单轴抗压强度,其中Coeff为实验室实测CO与测井计算CO的拟合校正系数;
3):根据测井密度计算上覆岩层压力Sv和声波速度,采用Eaton公式、Fillippone公式或者Bowers公式或者岩石物理法,计算地层孔隙压力Pp,其中对Fillippone公式Pp=(Vmax-Vi)/(Vmax-Vmin)*Sv*Coeff,式中Vmax为骨架速度,Vi为流体速度,采用岩石物理建模方法计算Vmax和Vi,Coeff为实测孔隙压力值与预测孔隙压力值的拟合校正系数;
4):根据水力压裂法获得水力压裂点的最小水平主应力σh,采用有效应力比值ESR=(σh-Pp)/(Sv-Pp),确定有效应力比ESR,计算连续的最小水平主应力σh;
5):根据成像测井分析崩落宽度WB0或者钻井过程诱导拉张缝,上覆岩层压力Sv、孔隙压力Pp及根据有效应力比值计算的最小水平主应力σh构建应力多边形;
6):分析重叠图像区域,采用图像分析算法,自动求取多边形图像中心,中心值即为该深度点最终的最大、最小水平主应力,然后对最大、最小水平主应力均采用有效应力比值法,重新计算最大、最小水平主应力连续剖面,采用同样的做法计算其他崩落深度的最大、最小水平主应力,最终获得各崩落深度点的最大最小水平主应力;
7):将各崩落点确定的最大、最小水平主应力带入地应力弹簧模型中,反演构造应变,计算高分辨率最大、最小水平主应力;
式中,σH为最大水平主应力,MPa;σh为最小水平主应力,MPa;Sv为上覆地层压力,MPa;
Pp为地层孔隙压力,Mpa;α为Biot系数;E为岩石静态杨氏模量,MPa;v为岩石静态泊松比;εH为最大水平应力方向构造应变;εh为最小水平应力方向构造应变;
8):根据三维地震速度和相控地震叠前反演,计算三维上覆岩层压力Sv、孔隙压力Pp、杨氏模量YM和泊松比v,结合反演的构造应变,带入应力弹簧模型中,得到高精度的三维地应力数据体;
所述步骤1)中根据不同地质环境进行的岩石物理建模的方法包括Xu-White及相关修正模型;
所述步骤2)中单轴抗压强度CO的计算方法是从杨氏模量和泊松比的物理意义出发,采用CO=YM/N*Coeff公式计算单轴抗压强度,其中Coeff为实验室实测CO与测井计算CO的拟合校正系数,采用岩石物理建模方法,获得诸如Biot系数、单轴抗压强度等岩石力学参数,同时提高岩石力学参数分析精度;
所述步骤6)中根据岩石破裂准则对单轴抗压强度进行调整,崩落区域与拉张区域重叠,根据重叠区域进一步约束了岩石破裂准则、单轴抗压强度及最大、最小水平主应力;
所述岩石破裂准则为修正Lade准则。
2.一种高精度三维地应力建模方法的应用,所述应用基于权利要求1所述的方法实现,其特征在于,所述方法结合三维地震速度及三维地震相控叠前反演,采用弹簧模型计算三维地应力,方便快捷地应用于石油天然气勘探开发中的地应力分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811153242.9A CN110969699B (zh) | 2018-09-29 | 2018-09-29 | 一种高精度三维地应力建模方法及其应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811153242.9A CN110969699B (zh) | 2018-09-29 | 2018-09-29 | 一种高精度三维地应力建模方法及其应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110969699A CN110969699A (zh) | 2020-04-07 |
CN110969699B true CN110969699B (zh) | 2023-08-29 |
Family
ID=70028202
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811153242.9A Active CN110969699B (zh) | 2018-09-29 | 2018-09-29 | 一种高精度三维地应力建模方法及其应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110969699B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111520125B (zh) * | 2020-05-07 | 2021-11-12 | 中国矿业大学 | 一种褶皱发育区现今地应力方位转向的定量分析方法 |
CN112034517B (zh) * | 2020-08-25 | 2021-07-13 | 北京福瑞科润能源科技有限公司 | 一种地应力数据预测方法、装置、储存介质及电子设备 |
CN111980697B (zh) * | 2020-09-23 | 2021-02-19 | 西南石油大学 | 天然裂缝页岩地层水力压裂水平井井筒套变量的计算方法 |
CN112528543B (zh) * | 2020-12-14 | 2022-11-01 | 中国葛洲坝集团第三工程有限公司 | 一种地应力的分析方法、装置以及处理设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102852516A (zh) * | 2012-04-19 | 2013-01-02 | 北京大学 | 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置 |
CN106033127A (zh) * | 2016-06-29 | 2016-10-19 | 中国石油化工股份有限公司 | 基于横波速度变化率的地应力方位地震预测方法 |
CN106546367A (zh) * | 2016-10-20 | 2017-03-29 | 中国石油天然气股份有限公司 | 一种地应力的确定方法和装置 |
CN108009705A (zh) * | 2017-11-07 | 2018-05-08 | 中国石油大学(华东) | 一种基于支持向量机技术的页岩储层可压性评价方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8768671B2 (en) * | 2010-04-26 | 2014-07-01 | Schlumberger Technology Corporation | System for optimizing a drilling operation and method for using same |
-
2018
- 2018-09-29 CN CN201811153242.9A patent/CN110969699B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102852516A (zh) * | 2012-04-19 | 2013-01-02 | 北京大学 | 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置 |
CN106033127A (zh) * | 2016-06-29 | 2016-10-19 | 中国石油化工股份有限公司 | 基于横波速度变化率的地应力方位地震预测方法 |
CN106546367A (zh) * | 2016-10-20 | 2017-03-29 | 中国石油天然气股份有限公司 | 一种地应力的确定方法和装置 |
CN108009705A (zh) * | 2017-11-07 | 2018-05-08 | 中国石油大学(华东) | 一种基于支持向量机技术的页岩储层可压性评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110969699A (zh) | 2020-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110969699B (zh) | 一种高精度三维地应力建模方法及其应用 | |
CN106951660B (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
CN104977618B (zh) | 一种评价页岩气储层及寻找甜点区的方法 | |
CN105386756B (zh) | 一种应用应变量计算脆性地层孔隙度的方法 | |
US20100191511A1 (en) | Method For Multi-Scale Geomechanical Model Analysis By Computer Simulation | |
Liu et al. | Asymmetric propagation mechanism of hydraulic fracture networks in continental reservoirs | |
MX2011001001A (es) | Metodos y sistemas para estimar la tension utilizando datos sismicos. | |
CN113820750A (zh) | 基于弹塑性力学定量预测泥岩构造裂缝的方法 | |
CN106014399A (zh) | 一种非均质地层高精度三维地应力模型建立方法 | |
CN112204225A (zh) | 用于预测岩层的剪切失效的系统和方法 | |
Shi et al. | Characterization of hydraulic fracture configuration based on complex in situ stress field of a tight oil reservoir in Junggar Basin, Northwest China | |
CN111257932B (zh) | 一种更新地层地震速度的方法 | |
CN111257946B (zh) | 一种地球物理指导钻井方法及更新地层地震速度的方法 | |
CN110967742B (zh) | 一种孔隙度反演方法及系统 | |
CN110320562B (zh) | 页岩气储层水平井中声波时差的校正方法 | |
Noufal et al. | Abu Dhabi stress map, UAE | |
CN109339776B (zh) | 一种各向异性地层地应力方位测量方法 | |
Liu et al. | Quantitative prediction of the 3D permeability tensor for a fractured reservoir in the Dingbian oilfield, Ordos basin, Central China | |
Abdideh et al. | Analysis of deep stress field using well log and wellbore breakout data: a case study in Cretaceous oil reservoir, southwest Iran | |
Martemyanov et al. | Analytic modelling for wellbore stability analysis | |
Tao et al. | Application of multipole array sonic logging to acid hydralic fracturing | |
CN109061737A (zh) | 一种基于合成地震记录的储层预测方法及装置 | |
Van Der Zee et al. | Improving Sub Salt Wellbore Stability Predictions Using 3D Geomechanical Simulations | |
Yang et al. | Unreliable determination of in situ stress orientation by cross multipole array acoustic logging in fractured shale reservoirs | |
CN110068862B (zh) | 盐下超深断背斜油气藏优质储层的预测方法及装置 |
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 |