CN116341174A - 一种裂缝性砂岩储层临界出砂生产压差预测方法及系统 - Google Patents
一种裂缝性砂岩储层临界出砂生产压差预测方法及系统 Download PDFInfo
- 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
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 59
- 239000004576 sand Substances 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 27
- 239000011435 rock Substances 0.000 claims abstract description 75
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 13
- 239000011148 porous material Substances 0.000 claims description 10
- 238000002474 experimental method Methods 0.000 claims description 9
- 230000002093 peripheral effect Effects 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 abstract description 5
- 238000011161 development Methods 0.000 description 6
- 239000007789 gas Substances 0.000 description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 6
- 239000003345 natural gas Substances 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000005483 Hooke's law Effects 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling 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和孔隙压力确定岩石最小和最大水平主应力。
优选的,所述最小和最大水平主应力的计算方法如下:
式中:σh为最小水平主应力,σH为最大水平主应力,α为比奥特系数,Pfo为孔隙压力。
优选的,步骤2中将上覆岩层应力、最小水平主应力、最大水平主应力转化为井周应力场。
优选的,所述井周应力场包括径向主应力σr、周向主应力σθ和垂向主应力σz;
σr=pw
其中,pw为井底流压,θ为方位角。
优选的,步骤3所述周应力场的最大和最小主应力的表达式如下:
σmax=max{σr、σθ、σz}
σmin=min{σr、σθ、σz}
优选的,步骤4中所述临界出砂生产压差的计算方法如下:
式中,σmax为最大主应力,Φ为内摩擦角,σmin为最小主应力。
一种裂缝性砂岩储层临界出砂生产压差预测方法的系统,包括,
模型构建模块,用于根据测井及岩石力学数据建立全井剖面地质力学模型;
井周应力场模块,用于根据全井剖面地质力学模型确定井周应力场;
主应力模块,用于确定井周应力场的最大和最小主应力;
压差计算模块,用于根据井周应力场的最大和最小主应力并结合岩石破坏准则,得到随地层压力衰竭,临界出砂生产压差。
与现有技术相比,本发明具有以下有益的技术效果:
本发明提供的一种裂缝性砂岩储层临界出砂生产压差预测方法,以测井及岩石力学测试实验数据为基础,采用Techlog软件建立全井剖面地质力学模型,并对全井剖面地质力学模型开展出砂预测,计算得到随地层压力衰竭界出砂生产压差变化情况,明确出砂高风险井段,为单井合理配产提供了依据。
附图说明
图1为本发明地质力学剖面图;
图2为本发明临界出砂生产压差剖面图;
图3为本发明高风险井段临界出砂生产压差随地层压力变化图。
具体实施方式
下面结合附图对本发明做进一步的详细说明,所述是对本发明的解释而不是限定。
参阅图1和2,一种裂缝性砂岩储层临界出砂生产压差预测方法,包括以下步骤:
步骤1、根据测井及岩石石力学数据建立全井剖面地质力学模型;
测井及岩石力学数据包括岩石单轴抗压强度UCS、岩石泊松比γ、岩石比奥特系数α、岩石杨氏模量E、岩石的内摩擦系数μ、岩石最小水平主应力方向应变εx、最大水平主应力方向应变εy、上覆岩层压力σv、岩石孔隙压力、岩石最小和最大水平主应力。
S1、获取岩石单轴抗压强度UCS,具体如下:
首先,开展岩心单轴抗压强度实验获取实验测量的岩石单轴抗压强度UCS,修正岩石单轴抗压强度UCS值并计算三维井壁稳定性,直到与成像获得的井壁崩落图像吻合,该时候修正后UCS值则为需要的岩石单轴抗压强度UCS。
S2、获取岩石泊松比γ。
岩石泊松比是岩石在单向受拉或受压时,横向正应变与轴向正应变的绝对值的比值,是一个无量纲量,可以通过下式或岩石压缩实验获得。
式中:vs为横波波速,vp为纵波波速。
S3、获取岩石比奥特系数α,计算公式如下:
式中:kb为岩体体积模量,ks为岩石颗粒体积模量。
S4、获取岩石杨氏模量E。
岩石杨氏模量是描述岩石抵抗形变能力的物理量。当一条长度为L、截面积为S的金属丝在力F作用下伸长ΔL时,F/S叫应力,其物理意义是金属丝单位截面积所受到的力。杨氏模量衡量反映材料的拉伸弹性,定义为在胡克定律适用的范围内,单轴应力和单轴形变之间的比,可以通过下式或进行岩石强度测试实验获得。
式中:vs为横波波速,vp为纵波波速,ρb为岩石密度。
S5、岩石强度测试实验获取岩石的内摩擦系数μ。
S6、获取岩石最小水平主应力方向应变εx、最大水平主应力方向应变εy。
可以通过调研该地区地质力学资料获得。
S7、计算上覆岩层压力。
式中:σv为上覆岩层应力,ρz为密度测井值,g为重力加速度,Z为上覆岩层埋深。
S8、通过井下压力计实测获得孔隙压力。
S9、根据岩石最小水平主应力方向应变εx、最大水平主应力方向应变εy和孔隙压力计算最小和最大水平主应力。
式中:σh为最小水平主应力,σH为最大水平主应力,Pfo为孔隙压力,α为比奥特系数。
步骤2、根据全井剖面地质力学模型确定井周应力场;
将上覆岩层应力、最小水平主应力、最大水平主应力转化为井周应力场。
地层岩石破坏、出砂首先从井周附近开始,因此需要将3个主应力场(上覆岩层应力、最小水平主应力、最大水平主应力)转化为井周应力场(以直井为例):
σr=pw
式中,σr为径向主应力,pw为井底流压,σθ为周向主应力,θ为方位角,σz为垂向主应力。
步骤3、开展井周应力场的最大、最小主应力判别。
σmax=max{σr、σθ、σz}
σmin=min{σr、σθ、σz}
步骤4、将最大和最小主应力代入岩石破坏准则中计算,得到随地层压力衰竭,临界出砂生产压差变化情况,明确出砂高风险井段(图2)。
式中σ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和孔隙压力确定岩石最小和最大水平主应力。
6.根据权利要求2所述的一种裂缝性砂岩储层临界出砂生产压差预测方法,其特征在于,步骤2中将上覆岩层应力、最小水平主应力、最大水平主应力转化为井周应力场。
8.根据权利要求7所述的一种裂缝性砂岩储层临界出砂生产压差预测方法,其特征在于,步骤3所述周应力场的最大和最小主应力的表达式如下:
σmax=max{σr、σθ、σz}
σmin=min{σr、σθ、σz}
10.一种权利要求1-9任一项所述的一种裂缝性砂岩储层临界出砂生产压差预测方法的系统,其特征在于,包括,
模型构建模块,用于根据测井及岩石力学数据建立全井剖面地质力学模型;
井周应力场模块,用于根据全井剖面地质力学模型确定井周应力场;
主应力模块,用于确定井周应力场的最大和最小主应力;
压差计算模块,用于根据井周应力场的最大和最小主应力并结合岩石破坏准则,得到随地层压力衰竭,临界出砂生产压差。
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116894572A (zh) * | 2023-09-11 | 2023-10-17 | 西南石油大学 | 一种超深井考虑岩崩后出砂的合理配产方法 |
-
2021
- 2021-12-24 CN CN202111605136.1A patent/CN116341174A/zh active Pending
Cited By (2)
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 |