CN116341174A - 一种裂缝性砂岩储层临界出砂生产压差预测方法及系统 - Google Patents

一种裂缝性砂岩储层临界出砂生产压差预测方法及系统 Download PDF

Info

Publication number
CN116341174A
CN116341174A CN202111605136.1A CN202111605136A CN116341174A CN 116341174 A CN116341174 A CN 116341174A CN 202111605136 A CN202111605136 A CN 202111605136A CN 116341174 A CN116341174 A CN 116341174A
Authority
CN
China
Prior art keywords
rock
stress
sand production
minimum
well
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.)
Pending
Application number
CN202111605136.1A
Other languages
English (en)
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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202111605136.1A priority Critical patent/CN116341174A/zh
Publication of CN116341174A publication Critical patent/CN116341174A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提供的一种裂缝性砂岩储层临界出砂生产压差预测方法及系统,以测井及岩石力学测试实验数据为基础,采用Techlog软件建立全井剖面地质力学模型,并对全井剖面地质力学模型开展出砂预测,计算得到随地层压力衰竭界出砂生产压差变化情况,明确出砂高风险井段,为单井合理配产提供了依据。

Description

一种裂缝性砂岩储层临界出砂生产压差预测方法及系统
技术领域
本发明涉及石油开采领域,具体为一种裂缝性砂岩储层临界出砂生产压差预测方法及系统。
背景技术
随着石油天然气工业的不断发展和国内天然气需求量的日益攀升,超高压裂缝性致密砂岩气藏因其储量巨大、产能高,已经成为非常规天然气勘探开发的重要领域之一。随着此类气藏的大规模开采,出砂是油气井生产过程中面临的重要问题之一。出砂多发生于岩石孔隙度大、胶结程度低、岩石强小的疏松砂岩油气藏,致密砂岩油气藏通常不易出砂。然而随着致密气开发规模的不断扩张,部分超高压致密裂缝性砂岩气藏出现了大量出砂的情况。
塔里木油田库车山前地下储层致密、岩石强度高、发育裂缝,岩性为砂岩,实际开采过程中出现大面积出砂问题,而传统出砂预测方法认为该地区不出砂或临界出砂生产压差原高于现场实际出砂生产压差,与该地区实际生产情况不符,无法指导该地区合理开发。
因此,建立一套裂缝性砂岩出砂预测方法对于开发生产显得尤为必要。
发明内容
针对现有技术中存在的问题,本发明提供一种裂缝性砂岩储层临界出砂生产压差预测方法及系统,为单井合理配产提供了依据。
本发明是通过以下技术方案来实现:
一种裂缝性砂岩储层临界出砂生产压差预测方法,包括以下步骤:
步骤1、根据测井及岩石力学数据建立全井剖面地质力学模型;
步骤2、根据全井剖面地质力学模型确定井周应力场;
步骤3、确定井周应力场的最大和最小主应力;
步骤4、根据井周应力场的最大和最小主应力并结合岩石破坏准则,得到随地层压力衰竭,临界出砂生产压差。
优选的,步骤1中所述测井及岩石力学数据包括岩石单轴抗压强度UCS、岩石泊松比γ、岩石比奥特系数α、岩石杨氏模量E、岩石的内摩擦系数μ、岩石最小水平主应力方向应变εx、最大水平主应力方向应变εy、上覆岩层压力σv、岩石孔隙压力、岩石最小和最大水平主应力。
优选的,所述岩石单轴抗压强度UCS的确定方法如下:
首先,开展岩心单轴抗压强度实验获取实验测量的岩石单轴抗压强度UCS,对测量的岩石单轴抗压强度UCS修正并计算三维井壁稳定性,直到与成像获得的井壁崩落图像吻合,得到岩石单轴抗压强度UCS。
优选的,根据岩石最小水平主应力方向应变εx、最大水平主应力方向应变εy和孔隙压力确定岩石最小和最大水平主应力。
优选的,所述最小和最大水平主应力的计算方法如下:
Figure BDA0003433446870000021
Figure BDA0003433446870000022
式中:σh为最小水平主应力,σH为最大水平主应力,α为比奥特系数,Pfo为孔隙压力。
优选的,步骤2中将上覆岩层应力、最小水平主应力、最大水平主应力转化为井周应力场。
优选的,所述井周应力场包括径向主应力σr、周向主应力σθ和垂向主应力σz
σr=pw
Figure BDA0003433446870000031
Figure BDA0003433446870000032
其中,pw为井底流压,θ为方位角。
优选的,步骤3所述周应力场的最大和最小主应力的表达式如下:
σmax=max{σr、σθ、σz}
σmin=min{σr、σθ、σz}
优选的,步骤4中所述临界出砂生产压差的计算方法如下:
Figure BDA0003433446870000033
式中,σmax为最大主应力,Φ为内摩擦角,σmin为最小主应力。
一种裂缝性砂岩储层临界出砂生产压差预测方法的系统,包括,
模型构建模块,用于根据测井及岩石力学数据建立全井剖面地质力学模型;
井周应力场模块,用于根据全井剖面地质力学模型确定井周应力场;
主应力模块,用于确定井周应力场的最大和最小主应力;
压差计算模块,用于根据井周应力场的最大和最小主应力并结合岩石破坏准则,得到随地层压力衰竭,临界出砂生产压差。
与现有技术相比,本发明具有以下有益的技术效果:
本发明提供的一种裂缝性砂岩储层临界出砂生产压差预测方法,以测井及岩石力学测试实验数据为基础,采用Techlog软件建立全井剖面地质力学模型,并对全井剖面地质力学模型开展出砂预测,计算得到随地层压力衰竭界出砂生产压差变化情况,明确出砂高风险井段,为单井合理配产提供了依据。
附图说明
图1为本发明地质力学剖面图;
图2为本发明临界出砂生产压差剖面图;
图3为本发明高风险井段临界出砂生产压差随地层压力变化图。
具体实施方式
下面结合附图对本发明做进一步的详细说明,所述是对本发明的解释而不是限定。
参阅图1和2,一种裂缝性砂岩储层临界出砂生产压差预测方法,包括以下步骤:
步骤1、根据测井及岩石石力学数据建立全井剖面地质力学模型;
测井及岩石力学数据包括岩石单轴抗压强度UCS、岩石泊松比γ、岩石比奥特系数α、岩石杨氏模量E、岩石的内摩擦系数μ、岩石最小水平主应力方向应变εx、最大水平主应力方向应变εy、上覆岩层压力σv、岩石孔隙压力、岩石最小和最大水平主应力。
S1、获取岩石单轴抗压强度UCS,具体如下:
首先,开展岩心单轴抗压强度实验获取实验测量的岩石单轴抗压强度UCS,修正岩石单轴抗压强度UCS值并计算三维井壁稳定性,直到与成像获得的井壁崩落图像吻合,该时候修正后UCS值则为需要的岩石单轴抗压强度UCS。
S2、获取岩石泊松比γ。
岩石泊松比是岩石在单向受拉或受压时,横向正应变与轴向正应变的绝对值的比值,是一个无量纲量,可以通过下式或岩石压缩实验获得。
Figure BDA0003433446870000051
式中:vs为横波波速,vp为纵波波速。
S3、获取岩石比奥特系数α,计算公式如下:
Figure BDA0003433446870000052
式中:kb为岩体体积模量,ks为岩石颗粒体积模量。
S4、获取岩石杨氏模量E。
岩石杨氏模量是描述岩石抵抗形变能力的物理量。当一条长度为L、截面积为S的金属丝在力F作用下伸长ΔL时,F/S叫应力,其物理意义是金属丝单位截面积所受到的力。杨氏模量衡量反映材料的拉伸弹性,定义为在胡克定律适用的范围内,单轴应力和单轴形变之间的比,可以通过下式或进行岩石强度测试实验获得。
Figure BDA0003433446870000053
式中:vs为横波波速,vp为纵波波速,ρb为岩石密度。
S5、岩石强度测试实验获取岩石的内摩擦系数μ。
S6、获取岩石最小水平主应力方向应变εx、最大水平主应力方向应变εy。
可以通过调研该地区地质力学资料获得。
S7、计算上覆岩层压力。
Figure BDA0003433446870000061
式中:σv为上覆岩层应力,ρz为密度测井值,g为重力加速度,Z为上覆岩层埋深。
S8、通过井下压力计实测获得孔隙压力。
S9、根据岩石最小水平主应力方向应变εx、最大水平主应力方向应变εy和孔隙压力计算最小和最大水平主应力。
Figure BDA0003433446870000062
Figure BDA0003433446870000063
式中:σh为最小水平主应力,σH为最大水平主应力,Pfo为孔隙压力,α为比奥特系数。
步骤2、根据全井剖面地质力学模型确定井周应力场;
将上覆岩层应力、最小水平主应力、最大水平主应力转化为井周应力场。
地层岩石破坏、出砂首先从井周附近开始,因此需要将3个主应力场(上覆岩层应力、最小水平主应力、最大水平主应力)转化为井周应力场(以直井为例):
σr=pw
Figure BDA0003433446870000064
Figure BDA0003433446870000071
式中,σr为径向主应力,pw为井底流压,σθ为周向主应力,θ为方位角,σz为垂向主应力。
步骤3、开展井周应力场的最大、最小主应力判别。
σmax=max{σr、σθ、σz}
σmin=min{σr、σθ、σz}
步骤4、将最大和最小主应力代入岩石破坏准则中计算,得到随地层压力衰竭,临界出砂生产压差变化情况,明确出砂高风险井段(图2)。
Figure BDA0003433446870000072
式中σmax为最大主应力,Φ为内摩擦角,σmin为最小主应力。
步骤5、根据随地层压力衰竭,临界出砂生产压差变化情况,生成高风险段临界出砂生产压差随地层压力变化图版(图3)。
一种裂缝性砂岩储层临界出砂生产压差预测方法的系统,包括,
模型构建模块,用于根据测井及岩石力学数据建立全井剖面地质力学模型;
井周应力场模块,用于根据全井剖面地质力学模型确定井周应力场;
主应力模块,用于确定井周应力场的最大和最小主应力;
压差计算模块,用于根据井周应力场的最大和最小主应力并结合岩石破坏准则,得到随地层压力衰竭,临界出砂生产压差。
本发明提供的一种裂缝性砂岩储层临界出砂生产压差预测方法,以测井及岩石力学测试实验数据为基础,采用Techlog软件建立全井剖面地质力学模型,并对全井剖面地质力学模型开展出砂预测,计算得到随地层压力衰竭界出砂生产压差变化情况,明确出砂高风险井段,为单井合理配产提供了依据。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。

Claims (10)

1.一种裂缝性砂岩储层临界出砂生产压差预测方法,其特征在于,包括以下步骤:
步骤1、根据测井及岩石力学数据建立全井剖面地质力学模型;
步骤2、根据全井剖面地质力学模型确定井周应力场;
步骤3、确定井周应力场的最大和最小主应力;
步骤4、根据井周应力场的最大和最小主应力并结合岩石破坏准则,得到随地层压力衰竭,临界出砂生产压差。
2.根据权利要求1所述的一种裂缝性砂岩储层临界出砂生产压差预测方法,其特征在于,步骤1中所述测井及岩石力学数据包括岩石单轴抗压强度UCS、岩石泊松比γ、岩石比奥特系数α、岩石杨氏模量E、岩石的内摩擦系数μ、岩石最小水平主应力方向应变εx、最大水平主应力方向应变εy、上覆岩层压力σv、岩石孔隙压力、岩石最小和最大水平主应力。
3.根据权利要求1所述的一种裂缝性砂岩储层临界出砂生产压差预测方法,其特征在于,所述岩石单轴抗压强度UCS的确定方法如下:
首先,开展岩心单轴抗压强度实验获取实验测量的岩石单轴抗压强度UCS,对测量的岩石单轴抗压强度UCS修正并计算三维井壁稳定性,直到与成像获得的井壁崩落图像吻合,得到岩石单轴抗压强度UCS。
4.根据权利要求2所述的一种裂缝性砂岩储层临界出砂生产压差预测方法,其特征在于,根据岩石最小水平主应力方向应变εx、最大水平主应力方向应变εy和孔隙压力确定岩石最小和最大水平主应力。
5.根据权利要求4所述的一种裂缝性砂岩储层临界出砂生产压差预测方法,其特征在于,所述最小和最大水平主应力的计算方法如下:
Figure FDA0003433446860000021
Figure FDA0003433446860000022
式中:σh为最小水平主应力,σH为最大水平主应力,α为比奥特系数,Pfo为孔隙压力。
6.根据权利要求2所述的一种裂缝性砂岩储层临界出砂生产压差预测方法,其特征在于,步骤2中将上覆岩层应力、最小水平主应力、最大水平主应力转化为井周应力场。
7.根据权利要求6所述的一种裂缝性砂岩储层临界出砂生产压差预测方法,其特征在于,所述井周应力场包括径向主应力σr、周向主应力σθ和垂向主应力σz
σr=pw
Figure FDA0003433446860000023
Figure FDA0003433446860000024
其中,pw为井底流压,θ为方位角。
8.根据权利要求7所述的一种裂缝性砂岩储层临界出砂生产压差预测方法,其特征在于,步骤3所述周应力场的最大和最小主应力的表达式如下:
σmax=max{σr、σθ、σz}
σmin=min{σr、σθ、σz}
9.根据权利要求8所述的一种裂缝性砂岩储层临界出砂生产压差预测方法,其特征在于,步骤4中所述临界出砂生产压差的计算方法如下:
Figure FDA0003433446860000031
式中,σmax为最大主应力,Φ为内摩擦角,σmin为最小主应力。
10.一种权利要求1-9任一项所述的一种裂缝性砂岩储层临界出砂生产压差预测方法的系统,其特征在于,包括,
模型构建模块,用于根据测井及岩石力学数据建立全井剖面地质力学模型;
井周应力场模块,用于根据全井剖面地质力学模型确定井周应力场;
主应力模块,用于确定井周应力场的最大和最小主应力;
压差计算模块,用于根据井周应力场的最大和最小主应力并结合岩石破坏准则,得到随地层压力衰竭,临界出砂生产压差。
CN202111605136.1A 2021-12-24 2021-12-24 一种裂缝性砂岩储层临界出砂生产压差预测方法及系统 Pending CN116341174A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111605136.1A CN116341174A (zh) 2021-12-24 2021-12-24 一种裂缝性砂岩储层临界出砂生产压差预测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111605136.1A CN116341174A (zh) 2021-12-24 2021-12-24 一种裂缝性砂岩储层临界出砂生产压差预测方法及系统

Publications (1)

Publication Number Publication Date
CN116341174A true CN116341174A (zh) 2023-06-27

Family

ID=86879491

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111605136.1A Pending CN116341174A (zh) 2021-12-24 2021-12-24 一种裂缝性砂岩储层临界出砂生产压差预测方法及系统

Country Status (1)

Country Link
CN (1) CN116341174A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116894572A (zh) * 2023-09-11 2023-10-17 西南石油大学 一种超深井考虑岩崩后出砂的合理配产方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116894572A (zh) * 2023-09-11 2023-10-17 西南石油大学 一种超深井考虑岩崩后出砂的合理配产方法
CN116894572B (zh) * 2023-09-11 2023-12-15 西南石油大学 一种超深井考虑岩崩后出砂的合理配产方法

Similar Documents

Publication Publication Date Title
CN108868748B (zh) 一种页岩气水平井重复压裂裂缝开启压力的计算方法
CN108952700B (zh) 一种各向异性地层井壁破裂压力确定方法
CN104806233B (zh) 一种预测弱面地层坍塌压力当量密度窗口的方法
CN109356567B (zh) 深水浅部地层井壁稳定性预测方法
CN107038313B (zh) 基于数值岩芯的分层地应力精细描述方法
EP2583214B1 (en) Systems and methods for wellbore optimization
WO2022193655A1 (zh) 一种基于井壁力学失稳垮塌的地应力大小评价方法
CN107169248B (zh) 一种特殊地层安全泥浆密度窗口确定方法
CN109684685B (zh) 一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法
CN110566171A (zh) 一种超高压致密裂缝性砂岩气藏出砂预测方法
CN106289964B (zh) 一种泥页岩古抗压强度确定方法
CN108304959B (zh) 提高地层流体压力预测精度的方法
CN116341174A (zh) 一种裂缝性砂岩储层临界出砂生产压差预测方法及系统
Yu et al. On how asymmetric stimulated rock volume in shales may impact casing integrity
CN114547906A (zh) 一种深部含软弱结构面地层的井壁稳定测井解释方法
CN113109162B (zh) 一种基于热流固耦合的岩石起裂压力计算方法
CN112412434B (zh) 一种改进的疏松砂岩地应力计算方法
CN116205163B (zh) 一种天然气水合物藏数值模拟中的地层稳定性分析方法
CN108708715A (zh) 快速预测页岩气层地层破裂压力梯度的测录井方法
CN114033356B (zh) 一种煤系地层地应力计算方法及装置
CN115951422A (zh) 构建天然裂缝漏失压力模型的方法
CN117313508A (zh) 一种砂砾岩水平井地层可钻性分析方法
CN107016219B (zh) 一种碳酸盐岩储层钻井放空的预警方法及系统
CN112651113A (zh) 适用于东海压力衰竭地层的安全泥浆密度窗口判定方法
CN114647913A (zh) 一种基于Mohr-Coulomb准则的最大水平主地应力获取方法

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