CN107038313A - 基于数值岩芯的分层地应力精细描述方法 - Google Patents

基于数值岩芯的分层地应力精细描述方法 Download PDF

Info

Publication number
CN107038313A
CN107038313A CN201710282601.XA CN201710282601A CN107038313A CN 107038313 A CN107038313 A CN 107038313A CN 201710282601 A CN201710282601 A CN 201710282601A CN 107038313 A CN107038313 A CN 107038313A
Authority
CN
China
Prior art keywords
stress
core
dimension
static
rock mechanics
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.)
Granted
Application number
CN201710282601.XA
Other languages
English (en)
Other versions
CN107038313B (zh
Inventor
李连崇
朱万成
李明
杨天鸿
于永军
刘洪磊
于庆磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeastern University China
Original Assignee
Northeastern University China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northeastern University China filed Critical Northeastern University China
Priority to CN201710282601.XA priority Critical patent/CN107038313B/zh
Publication of CN107038313A publication Critical patent/CN107038313A/zh
Application granted granted Critical
Publication of CN107038313B publication Critical patent/CN107038313B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Mining & Mineral Resources (AREA)
  • Health & Medical Sciences (AREA)
  • Geometry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Agronomy & Crop Science (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明涉及低渗多薄层油气储层精细分层压裂相关的技术领域,具体为一种基于数值岩芯的分层地应力精细描述方法,主要基于连续测井数据、岩芯岩石力学参数物理测试数据、全井段岩芯岩石力学参数数值模拟测试数据。本发明提出了运用数模测试进行岩芯静态岩石力学参数模拟分析,校正得到近似真实的岩石力学参数,并结合测井资料测算得到的储层动态岩石力学参数,建立起动、静态力学参数的关系,并结合压裂施工数据来确定分层地应力,从而获得精细地应力剖面的新方法。通过精细的岩性识别和精细的岩石力学参数计算,可提高储层纵向上的有效改造程度。

Description

基于数值岩芯的分层地应力精细描述方法
技术领域:
本发明涉及低渗多薄层油气储层精细分层压裂相关的技术领域,具体为一种基于数值岩芯的分层地应力精细描述方法,主要基于连续测井数据、岩芯岩石力学参数物理测试数据、全井段岩芯岩石力学参数数值模拟测试数据。
背景技术:
薄互层低渗油气藏通常需要压裂改造后才能经济动用,在薄互层发育的低渗透油气储层开发时,储层多且薄、隔层效果差、非均质性强等原因,最终导致地层滤失严重、缝高难以控制、施工困难等问题。其中储存纵向上分层地应力分布状态对水力压裂裂缝的产生及其形态、方位和高度等都有重要影响。地应力剖面可以反映地应力场在纵向上的变化规律,准确获取分层地应力参数可以为钻井工程、油气藏工程和采油气工程等各个环节的决策和设计提供基础参数。因此,获取精细的地应力剖面对于油气田的开发具有重要意义。
现在通常采用的测井动态参数,不能完全反映地层真实情况,导致复杂岩性储层岩石力学解释存在局限性,需要动静态岩石力学参数转化。但是在进行静态岩石力学参数室内物模实验时,受控于复杂的岩层地质情况、施工设施以及成本问题,实际钻取的岩芯数量有限、取芯也是非连续的,不可能涵盖工程设计所需要的各个层位,而且物模实验数据的离散性较大;此外,受限于当前的真三轴岩芯物理实验设备的局限性,当前的岩石力学实验大多是在单轴或假三轴应力状态下进行的,其得到的岩芯力学参数也较难反映实际岩芯所处的真三维应力状态。因此,可在有限的物模实验基础之上,利用数值模拟方法,开展岩芯的数值化,进行岩芯岩石力学参数的数值模拟分析,获得连续的、精细的岩芯静态岩石力学参数。
发明内容:
本发明的目的在于提供一种基于数值岩芯的分层地应力精细描述方法,该方法通过精细的岩性识别以及静态岩石力学参数计算,可获得连续的、精细的应力分层剖面,有助于提高储层纵向上的有效改造程度。
本发明的技术方案是:
一种基于数值岩芯的分层地应力精细描述方法,包括下列步骤:
a、利用测井资料获取连续的动态岩石力学参数;
b、获取非连续静态岩石力学参数;
c、将步骤b中测试得到的岩芯岩石力学静态参数与步骤a中的连续测井动态岩石力学参数相结合,回归形成动-静态岩石泊松比、弹性模量关系式,见公式(1)、(2),获得精细的、连续的、纵向上不同储层深度的静态岩石力学参数;
岩石静-动态泊松比的关系:
υs=f(υd) (1)
岩石静-动态弹性模量的关系:
Es=f(Ed) (2)
式中:υs,υd分别是岩石静态、动态泊松比,无量纲;Es,Ed分别是岩石静态、动态弹性模量,量纲:MPa;
d、根据水力压裂施工压力曲线计算施工点的最大、最小水平主应力:按照水力压裂施工压力曲线,得到原地应力与施工压力的表征关系,见公式(3)、(4):
σh=p闭合 (3)
σH=3p闭合-p重张-αp重张 (4)
σh是最小水平主应力,量纲:MPa;σH是最大水平主应力,量纲:MPa;
e、计算连续的地应力剖面数据:首先根据步骤c得到的静态弹性模量Es、静态泊松比υs,以及步骤d得到的最大主应力σH、最小主应力σh,以及经典的地应力组合弹簧计算模型,见公式(5)、(6)、(7),计算得到最大、最小水平构造应力系数KH和Kh,构造应力系数KH和Kh不随井深和计算地点发生变化;最后根据公式(5)、(6)、(7)计算出该区块不同储层深度、连续的、精细的地应力分布;
式中:是地层i点深度的上覆压力,量纲:Pa;ρ0是无密度测井数据段平均密度,量纲:kg/m3;H0是无密度测井数据段长度,量纲:m;ρi是测井密度数据,量纲:kg/m3;dhi是对应ρi的测井层段厚度,量纲:m;g是重力加速度,9.8m/s2;σv是垂向主应力,量纲:MPa;H是储层深度,量纲:m;α是有效应力系数,无量纲;pp是孔隙压力,量纲:MPa;Kh是最小水平主应力方向的构造系数,在同一断块内为常数,量纲:m-1;KH为最大水平主应力方向的构造系数,在同一断块内为常数,量纲:m-1
所述的基于数值岩芯的分层地应力精细描述方法,步骤a中,声波测井资料作为油气田的常规资料,由连续的声波测井资料实现对储存纵向地应力剖面的连续计算,即:获得区块储层纵向上的、连续的动态弹性模量Ed和动态泊松比υd的分布数据。
所述的基于数值岩芯的分层地应力精细描述方法,声波测井资料采用声波时差。
所述的基于数值岩芯的分层地应力精细描述方法,步骤b中,首先按照国家标准GB/T50266-99、国际岩石力学学会的岩石力学实验建议方法,进行有限数量岩芯的单轴和三轴岩石力学物模实验,获得岩芯静态弹性模量Es和静态泊松比υs
所述的基于数值岩芯的分层地应力精细描述方法,对于典型的需要进一步深入校核的储层、以及尚未取到实际岩芯的储层,参考物模流程,建立数值岩芯,进行任意层位储层、任意孔隙度、任意围压的单轴和三轴岩石力学模拟测试,得到岩芯的全程应力应变曲线,进而纵向上储层的精细静态岩石力学参数。
本发明的优点及有益效果是:
本发明提出了运用数模测试进行岩芯静态岩石力学参数模拟分析,校正得到近似真实的岩石力学参数,并结合测井资料测算得到的储层动态岩石力学参数,建立起动、静态力学参数的关系,并结合压裂施工数据来确定分层地应力,从而获得精细地应力剖面的新方法。通过精细的岩性识别和精细的岩石力学参数计算,可提高储层纵向上的有效改造程度。
附图说明:
图1是基于数值模拟的岩芯单轴/三轴岩石力学参数获取示意图。
图2是岩芯的应力应变图。
图3是经校正后的精细连续储层静态岩石力学参数示意图。
图4是水力压裂施工压力曲线示意图。
图5是精细连续分层地应力曲线。
具体实施方式:
在具体实施过程中,本发明基于数值岩芯的分层地应力精细描述方法,实施过程包括下列步骤:
a、利用测井资料获取连续的动态岩石力学参数:声波测井资料(如:声波时差)作为油气田的常规资料,丰富并且便于获取,由连续的声波测井资料实现对储存纵向地应力剖面的连续计算,即:获得区块储层纵向上的、连续的动态弹性模量Ed(量纲:MPa)和动态泊松比υd(无量纲)的分布数据。
b、获取非连续静态岩石力学参数:首先按照国家标准GB/T50266-99(工程岩体实验方法标准)、国际岩石力学学会(ISRM)的岩石力学实验建议方法,进行有限数量岩芯的单轴和三轴岩石力学物模实验,获得岩芯静态弹性模量Es(量纲:MPa)和静态泊松比υs(无量纲)。对于典型的需要进一步深入校核的储层、以及尚未取到实际岩芯的储层,参考物模流程,建立数值岩芯,如图1,进行任意层位储层、任意孔隙度、任意围压的单轴和三轴岩石力学模拟测试,得到岩芯的全程应力应变曲线,一般的岩石岩芯应力-应变曲线见图2,进而纵向上储层的精细静态岩石力学参数。
c、将步骤b中测试得到的岩芯岩石力学静态参数与步骤a中的连续测井动态岩石力学参数相结合,回归形成动-静态岩石泊松比、弹性模量关系式,如:公式(1)、(2),这样就获得了精细的、连续的、纵向上不同储层深度的静态岩石力学参数,如图3所示。
岩石静-动态泊松比的关系:
υs=f(υd) (1)
岩石静-动态弹性模量的关系:
Es=f(Ed) (2)
式中:υs,υd分别是岩石静态、动态泊松比,无量纲;Es,Ed分别是岩石静态、动态弹性模量,量纲:MPa;
d、根据水力压裂施工压力曲线计算施工点(任意储层深度)的最大、最小水平主应力:按照一般的水力压裂施工压力曲线,见图4,得到原地应力与施工压力的表征关系,见公式(3)、(4):
σh=p闭合 (3)
σH=3p闭合-p重张-αp重张 (4)
σh是最小水平主应力,量纲:MPa;σH是最大水平主应力,量纲:MPa;
e、计算连续的地应力剖面数据:首先根据步骤c得到的静态弹性模量Es、静态泊松比υs,以及步骤d得到的最大主应力σH、最小主应力σh,以及经典的地应力组合弹簧计算模型(见公式(5)、(6)、(7)),计算得到最大、最小水平构造应力系数KH和Kh,构造应力系数KH和Kh不随井深和计算地点发生变化;最后根据公式(5)、(6)、(7)计算出该区块不同储层深度、连续的、精细的地应力分布,如图5所示。
式中:是地层i点深度的上覆压力,量纲:Pa;ρ0是无密度测井数据段平均密度,量纲:kg/m3;H0是无密度测井数据段长度,量纲:m;ρi是测井密度数据,量纲:kg/m3;dhi是对应ρi的测井层段厚度,量纲:m;g是重力加速度,9.8m/s2;σv是垂向主应力,量纲:MPa;H是储层深度,量纲:m;α是有效应力系数,无量纲;pp是孔隙压力,量纲:MPa;Kh是最小水平主应力方向的构造系数,在同一断块内为常数,量纲:m-1;KH为最大水平主应力方向的构造系数,在同一断块内为常数,量纲:m-1
如图1所示,本发明基于数值模拟的岩芯单轴/三轴岩石力学参数获取过程主要包括:基于实际物理岩芯的几何尺寸、物理力学参数,进行网格剖分,建立起数值岩芯,并进行压缩、拉伸测试,测试结果的合理性首先通过破裂面模式进行校核,数值岩芯的破裂面模式要符合压缩或拉伸破坏模式,并且与物理岩芯的破裂面模式基本一致。
如图2所示,在岩芯破裂面模式基本一致的基础之上,进一步校核岩芯的应力应变曲线,特别关注弹性变形阶段、峰值强度,保证数模测试结果与物理模拟结果的一致性,例如实施例的弹性变形阶段和峰值强度,数模测试结果与物理模拟结果的一致性较好。在数模测试结果与物理模拟结果的对比标定基础之上,进行任意层位岩芯岩石力学参数数模测试,特别是对于没有取得实际物理岩芯的层位,可进行多次测试校正。
如图3所示,是经校正后的实施例的精细连续储层静态岩石力学参数,这些参数是在岩石力学参数数模测试数据、有限的物理岩芯岩石力学参数测试、连续测井动态岩石力学参数数据基础之上,并利用公式(1)、公式(2)计算求得的,弹性模量和泊松比将为下一步计算分层地应力剖面奠定基础。
如图4所示,是水力压裂施工压力曲线,该曲线是压裂施工中实际测得的,利用该压力曲线中的关键点:P闭合压力、P重张压力,并结合公式(3)、公式(4)计算出任意储层深度的最大、最小水平主应力。
如图5所示,是精细连续分层地应力曲线,该曲线有效表征了实施例不同储层深度的垂直地应力、最大水平主应力和最小水平主应力的大小,该曲线上的每个数据点是根据连续的储层静态岩石力学参数和施工压力曲线的特征点压力值,并结合公式(5)、公式(6)和公式(7)计算求得的。
结果表明,本发明提出的基于数值岩芯的分层地应力精细描述方法,可以获得精细的连续地应力剖面,对于压裂改造工程设计和提高储层纵向上的有效改造程度,具有较高的实用意义和经济价值。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (5)

1.一种基于数值岩芯的分层地应力精细描述方法,其特征在于,包括下列步骤:
a、利用测井资料获取连续的动态岩石力学参数;
b、获取非连续静态岩石力学参数;
c、将步骤b中测试得到的岩芯岩石力学静态参数与步骤a中的连续测井动态岩石力学参数相结合,回归形成动-静态岩石泊松比、弹性模量关系式,见公式(1)、(2),获得精细的、连续的、纵向上不同储层深度的静态岩石力学参数;
岩石静-动态泊松比的关系:
υs=f(υd) (1)
岩石静-动态弹性模量的关系:
Es=f(Ed) (2)
式中:υs,υd分别是岩石静态、动态泊松比,无量纲;Es,Ed分别是岩石静态、动态弹性模量,量纲:MPa;
d、根据水力压裂施工压力曲线计算施工点的最大、最小水平主应力:按照水力压裂施工压力曲线,得到原地应力与施工压力的表征关系,见公式(3)、(4):
σh=p闭合 (3)
σH=3p闭合-p重张-αp重张 (4)
σh是最小水平主应力,量纲:MPa;σH是最大水平主应力,量纲:MPa;
e、计算连续的地应力剖面数据:首先根据步骤c得到的静态弹性模量Es、静态泊松比υs,以及步骤d得到的最大主应力σH、最小主应力σh,以及经典的地应力组合弹簧计算模型,见公式(5)、(6)、(7),计算得到最大、最小水平构造应力系数KH和Kh,构造应力系数KH和Kh不随井深和计算地点发生变化;最后根据公式(5)、(6)、(7)计算出该区块不同储层深度、连续的、精细的地应力分布;
式中:是地层i点深度的上覆压力,量纲:Pa;ρ0是无密度测井数据段平均密度,量纲:kg/m3;H0是无密度测井数据段长度,量纲:m;ρi是测井密度数据,量纲:kg/m3;dhi是对应ρi的测井层段厚度,量纲:m;g是重力加速度,9.8m/s2;σv是垂向主应力,量纲:MPa;H是储层深度,量纲:m;α是有效应力系数,无量纲;pp是孔隙压力,量纲:MPa;Kh是最小水平主应力方向的构造系数,在同一断块内为常数,量纲:m-1;KH为最大水平主应力方向的构造系数,在同一断块内为常数,量纲:m-1
2.按照权利要求1所述的基于数值岩芯的分层地应力精细描述方法,其特征在于,步骤a中,声波测井资料作为油气田的常规资料,由连续的声波测井资料实现对储存纵向地应力剖面的连续计算,即:获得区块储层纵向上的、连续的动态弹性模量Ed和动态泊松比υd的分布数据。
3.按照权利要求2所述的基于数值岩芯的分层地应力精细描述方法,其特征在于,声波测井资料采用声波时差。
4.按照权利要求1所述的基于数值岩芯的分层地应力精细描述方法,其特征在于,步骤b中,首先按照国家标准GB/T50266-99、国际岩石力学学会的岩石力学实验建议方法,进行有限数量岩芯的单轴和三轴岩石力学物模实验,获得岩芯静态弹性模量Es和静态泊松比υs
5.按照权利要求4所述的基于数值岩芯的分层地应力精细描述方法,其特征在于,对于典型的需要进一步深入校核的储层、以及尚未取到实际岩芯的储层,参考物模流程,建立数值岩芯,进行任意层位储层、任意孔隙度、任意围压的单轴和三轴岩石力学模拟测试,得到岩芯的全程应力应变曲线,进而纵向上储层的精细静态岩石力学参数。
CN201710282601.XA 2017-04-26 2017-04-26 基于数值岩芯的分层地应力精细描述方法 Expired - Fee Related CN107038313B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710282601.XA CN107038313B (zh) 2017-04-26 2017-04-26 基于数值岩芯的分层地应力精细描述方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710282601.XA CN107038313B (zh) 2017-04-26 2017-04-26 基于数值岩芯的分层地应力精细描述方法

Publications (2)

Publication Number Publication Date
CN107038313A true CN107038313A (zh) 2017-08-11
CN107038313B CN107038313B (zh) 2020-03-20

Family

ID=59535230

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710282601.XA Expired - Fee Related CN107038313B (zh) 2017-04-26 2017-04-26 基于数值岩芯的分层地应力精细描述方法

Country Status (1)

Country Link
CN (1) CN107038313B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109113742A (zh) * 2018-08-02 2019-01-01 中国矿业大学 一种煤储层现今地应力预测方法
CN110501758A (zh) * 2019-07-30 2019-11-26 东北大学 一种砂砾岩储层纵向连续脆性指数预方法
CN110954419A (zh) * 2019-11-15 2020-04-03 山东大学 一种预置裂隙水力疲劳压裂设计方法
CN112394397A (zh) * 2019-08-13 2021-02-23 中国石油化工股份有限公司 一种页岩气藏三维岩石力学参数场建模方法
CN112412434A (zh) * 2020-06-02 2021-02-26 中国海洋石油集团有限公司 一种改进的疏松砂岩地应力计算方法
CN112557185A (zh) * 2019-09-26 2021-03-26 中国石油天然气股份有限公司 储层地应力状态的测定方法及装置
CN113868923A (zh) * 2021-10-13 2021-12-31 西南石油大学 油气储层压前三维地质评价方法
CN114112651A (zh) * 2020-08-27 2022-03-01 中国石油化工股份有限公司 一种用于人造岩心的岩石动静态力学参数转换方法及系统
CN114184764A (zh) * 2020-09-15 2022-03-15 中国石油化工股份有限公司 致密碳酸盐岩储层岩石力学层划分方法及系统
CN115586569A (zh) * 2022-09-27 2023-01-10 西南石油大学 基于理论模型约束下数据驱动的地层水平地应力计算方法
CN117633998A (zh) * 2023-12-13 2024-03-01 石家庄铁道大学 一种隧道支护体系数值试验模拟实现方法
CN117634228A (zh) * 2024-01-26 2024-03-01 中国科学院地质与地球物理研究所 储层可控冲击波致裂的非连续变形分析模拟计算方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100319909A1 (en) * 2006-10-13 2010-12-23 Symington William A Enhanced Shale Oil Production By In Situ Heating Using Hydraulically Fractured Producing Wells
CN103206206A (zh) * 2012-01-16 2013-07-17 中国石油化工股份有限公司 用于计算超深井地层破裂压力的方法
CN103267678A (zh) * 2013-04-28 2013-08-28 北京大学 一种动静态岩石力学参数同步测量方法及装置
CN103278389A (zh) * 2013-04-28 2013-09-04 北京大学 一种岩石动、静态弹性参数同步测量的方法
KR101415196B1 (ko) * 2013-10-15 2014-07-04 한국지질자원연구원 전기비저항 및 밀도 검층 데이터를 이용한 셰일가스전 스윗스팟 수평시추구간의 파쇄구간 선정방법
CN105114047A (zh) * 2015-07-22 2015-12-02 西南石油大学 一种利用测井资料优选体积压裂射孔层位的方法
CN105866835A (zh) * 2016-03-28 2016-08-17 中国石油大学(华东) 一种基于地应力分布的断层三维封闭性定量评价方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100319909A1 (en) * 2006-10-13 2010-12-23 Symington William A Enhanced Shale Oil Production By In Situ Heating Using Hydraulically Fractured Producing Wells
CN103206206A (zh) * 2012-01-16 2013-07-17 中国石油化工股份有限公司 用于计算超深井地层破裂压力的方法
CN103267678A (zh) * 2013-04-28 2013-08-28 北京大学 一种动静态岩石力学参数同步测量方法及装置
CN103278389A (zh) * 2013-04-28 2013-09-04 北京大学 一种岩石动、静态弹性参数同步测量的方法
KR101415196B1 (ko) * 2013-10-15 2014-07-04 한국지질자원연구원 전기비저항 및 밀도 검층 데이터를 이용한 셰일가스전 스윗스팟 수평시추구간의 파쇄구간 선정방법
CN105114047A (zh) * 2015-07-22 2015-12-02 西南石油大学 一种利用测井资料优选体积压裂射孔层位的方法
CN105866835A (zh) * 2016-03-28 2016-08-17 中国石油大学(华东) 一种基于地应力分布的断层三维封闭性定量评价方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王丽忱: "页岩气藏岩石力学性质的测井评价方法与应用", 《中国优秀硕士学位论文全文数据库》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109113742B (zh) * 2018-08-02 2019-12-31 中国矿业大学 一种煤储层现今地应力预测方法
CN109113742A (zh) * 2018-08-02 2019-01-01 中国矿业大学 一种煤储层现今地应力预测方法
CN110501758A (zh) * 2019-07-30 2019-11-26 东北大学 一种砂砾岩储层纵向连续脆性指数预方法
CN110501758B (zh) * 2019-07-30 2021-06-18 东北大学 一种砂砾岩储层纵向连续脆性指数预测方法
CN112394397A (zh) * 2019-08-13 2021-02-23 中国石油化工股份有限公司 一种页岩气藏三维岩石力学参数场建模方法
CN112557185A (zh) * 2019-09-26 2021-03-26 中国石油天然气股份有限公司 储层地应力状态的测定方法及装置
CN112557185B (zh) * 2019-09-26 2023-01-10 中国石油天然气股份有限公司 储层地应力状态的测定方法及装置
CN110954419B (zh) * 2019-11-15 2020-11-10 山东大学 一种预置裂隙水力疲劳压裂设计方法
CN110954419A (zh) * 2019-11-15 2020-04-03 山东大学 一种预置裂隙水力疲劳压裂设计方法
CN112412434A (zh) * 2020-06-02 2021-02-26 中国海洋石油集团有限公司 一种改进的疏松砂岩地应力计算方法
CN112412434B (zh) * 2020-06-02 2024-03-26 中国海洋石油集团有限公司 一种改进的疏松砂岩地应力计算方法
CN114112651A (zh) * 2020-08-27 2022-03-01 中国石油化工股份有限公司 一种用于人造岩心的岩石动静态力学参数转换方法及系统
CN114184764A (zh) * 2020-09-15 2022-03-15 中国石油化工股份有限公司 致密碳酸盐岩储层岩石力学层划分方法及系统
CN114184764B (zh) * 2020-09-15 2024-06-04 中国石油化工股份有限公司 致密碳酸盐岩储层岩石力学层划分方法及系统
CN113868923A (zh) * 2021-10-13 2021-12-31 西南石油大学 油气储层压前三维地质评价方法
CN115586569A (zh) * 2022-09-27 2023-01-10 西南石油大学 基于理论模型约束下数据驱动的地层水平地应力计算方法
CN117633998A (zh) * 2023-12-13 2024-03-01 石家庄铁道大学 一种隧道支护体系数值试验模拟实现方法
CN117634228A (zh) * 2024-01-26 2024-03-01 中国科学院地质与地球物理研究所 储层可控冲击波致裂的非连续变形分析模拟计算方法
CN117634228B (zh) * 2024-01-26 2024-04-02 中国科学院地质与地球物理研究所 储层可控冲击波致裂的非连续变形分析模拟计算方法

Also Published As

Publication number Publication date
CN107038313B (zh) 2020-03-20

Similar Documents

Publication Publication Date Title
CN107038313A (zh) 基于数值岩芯的分层地应力精细描述方法
CN105866835B (zh) 一种基于地应力分布的断层三维封闭性定量评价方法
Lu et al. Model development and analysis of the evolution of coal permeability under different boundary conditions
CN112502701B (zh) 一种低渗透储层综合地质-工程的分类评价方法
CN106501856B (zh) 亚地震断层定量预测方法
CN106295095B (zh) 基于常规测井资料预测低渗透砂岩储层产能的方法
CN106154351A (zh) 一种低孔渗储层渗透率的估算方法
CN110130882A (zh) 一种基于测井测试资料的油藏区域评价方法
CN103761448B (zh) 煤层底板突水变权脆弱性评价法中调权参数的确定方法
CN103913772B (zh) 基于储层地质力学参数的微地震事件正演模拟方法
CN105527652B (zh) 一种岩石脆性的测井方法和装置
CN104500050A (zh) 一种裂缝性储层渗透率张量及各向异性定量预测方法
CN110147561A (zh) 一种含天然裂缝致密油气储层体积压裂缝网预测方法
CN103713327A (zh) 基于测井及地震资料的页岩气储层最小闭合压力评价方法
CN109583113A (zh) 一种岩石地层压实系数和有效孔隙体积压缩系数计算方法
Zhang et al. Numerical modeling of land subsidence resulting from oil production
TORRES et al. Characterization of complex hydraulic fractures in Eagle Ford shale oil development through embedded discrete fracture modeling
CN110501758A (zh) 一种砂砾岩储层纵向连续脆性指数预方法
CN105804732B (zh) 一种基于相对电阻率随机模拟反演的油气分布预测方法
Zhang et al. Experimental study on natural gas migration and accumulation mechanism in sweet spots of tight sandstones
CN108459357B (zh) 地层欠压实与流体膨胀超压的评价方法
Ren et al. The influence of argillaceous content in carbonate rocks on the 3D modeling and characterization of tectonic fracture parameters—example from the carboniferous and ordovician formations in the hetianhe gas field, Tarim Basin, NW China
Maxwell et al. Calibrated microseismic geomechanical modeling of a Horn River basin hydraulic fracture
Jambayev Discrete fracture network modeling for a carbonate reservoir
Zhang et al. Geomechanical evaluation enabled successful stimulation of a HPHT tight gas reservoir in western China

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200320