CN116068616B - 一种储层异常地层压力研究方法及系统 - Google Patents
一种储层异常地层压力研究方法及系统 Download PDFInfo
- Publication number
- CN116068616B CN116068616B CN202310203688.2A CN202310203688A CN116068616B CN 116068616 B CN116068616 B CN 116068616B CN 202310203688 A CN202310203688 A CN 202310203688A CN 116068616 B CN116068616 B CN 116068616B
- Authority
- CN
- China
- Prior art keywords
- pressure
- point
- formation
- rock
- depth
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Physics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- Geology (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Geophysics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种储层异常地层压力研究方法及系统,方法包括:使用超声波设备测量地层中H深度范围内声波传播时间差△t;通过静水压力模型获取地层中H深度范围内各深度点的静岩压力值;使用计算机绘制H与Ln△t关系图,其中:关系图中斜线表示静岩压力值,关系图中点状曲线表示深度H与Ln声波传播时间差△t变化曲线;在计算机根据斜线和点状曲线在关系图中确定各个点异常地层压力。系统基于上述方法。本发明通过泥岩声波时差确定孔隙压力,根据静水压力计算公式确定静水压力,将孔隙压力和静水压力绘制在关系图,根据关系图确定各点异常地层压力,能够准确方便的预测出异常地层压力。
Description
技术领域
本发明涉及油气勘探技术领域,尤其涉及一种储层异常地层压力研究方法及系统。
背景技术
地层压力与地层温度是开发油气田的能量,也是油气田开发中重要的基础参数。油气藏地层压力和温度的高低,不仅决定着油气流体的性质,还决定着油气田开发的方式、油气开采的技术特点与经济成本,以及最终的采收率。地层压力即岩石孔隙空间中的流体压力,对于石油勘探与开发具有十分重要的意义。
地层压力之一是异常地层压力,异常地层压力研究不仅是石油地质学的一个重要理论问题,而且也是长期困扰钻井工程界的重大技术难题。异常地层压力在世界范围内普遍存在,由于其形成机理和影响因素很多,国内外的研究一直未找到一种理想的方法,能够准确地对异常地层压力进行预测。
有鉴于此,有必要研究出一种储层异常地层压力研究方法及系统,以解决异常地层压力预测问题。
发明内容
本发明提供了一种储层异常地层压力研究方法及系统,通过泥岩声波时差确定孔隙压力,根据静水压力计算公式确定静水压力,将孔隙压力和静水压力绘制在关系图,根据关系图确定各点异常地层压力,能够准确方便的预测出异常地层压力。
实现本发明目的的技术方案如下:
一方面,提供一种储层异常地层压力研究方法,包括:
使用超声波设备测量地层中H深度范围内声波传播时间差△t;
通过静水压力模型获取地层中H深度范围内各深度点的静岩压力值;
使用计算机绘制H与Ln△t关系图,其中:关系图中斜线表示静岩压力值,关系图中点状曲线表示深度H与Ln声波传播时间差△t变化曲线;
在计算机根据斜线和点状曲线在关系图中确定各个点异常地层压力。
基于一方面,在一种可能的实现方式中,所述测量地层中H深度范围内声波传播时间差△t,包括:
使用超声波设备测量声波传播时间差△t是利用声波测井原理测量的;
声波传播时间差△t匹配有地层孔隙度Φ;
根据声波传播时间差△t和地层孔隙度Φ确定孔隙压力Ln△t;
孔隙压力Ln△t随H深度呈指数关系变化。
基于一方面,在一种可能的实现方式中,所述使用超声波设备测量声波传播时间差△t是利用声波测井原理测量的,包括:
泥岩孔隙度Φ随深度H的增加而减小,其减小程度呈指数关系变化,数学表达式为:
dΦ/Φ = -kdH (1)
Φ =Φ0e –k H (2)
Φ =(△t-△tma)/(△tf -△tma) (3)
式中(1)-公式(3)中: Φ-地层孔隙度;Φ0-上覆地层孔隙度;△t-地层声波时差,单位μs/m;△tf-孔隙流体声波时差,单位μs/m;△tma-岩石颗粒声波时差,单位μs/m;k-斜线斜率;H-目的层深度,单位m。
基于一方面,在一种可能的实现方式中,所述通过静水压力模型获取地层中H深度范围内各深度点的静岩压力值,包括:
静水压力模型输出连通岩石表面及地表的开放体系下的水柱压力;
参考水柱压力和静岩压力计算公式确定地层中H深度范围内各深度点的静岩压力值;
静岩压力值随着埋深增加而增加。
优选静水压力为连通岩石表面及地表的开放体系下的水柱压力;静岩压力由某一深度的上覆岩石和岩石中的流体所形成的压力;静水压力的压力值随深度值而增加;静水压力的压力变化率仅依赖于水的密度变化;静水压力代表压力增加最大速率的方向向量总是垂直的。
基于一方面,在一种可能的实现方式中,所述静水压力模型的计算公式为:
PH = ρw•H (4)
公式(4)中,PH-静水压力;ρw -水的平均密度; H-水柱高度;
所述静岩压力计算公式为:
Pob = Pe + P p (5)
公式(5)中,Pob-静岩压力,静岩压力的值随深度增加而增加;Pe-粒间接触所施加应力,粒间接触所施加应力等于上覆岩柱的重量;Pp-孔隙流体压力,孔隙流体压力取决于岩石骨架密度、孔隙流体密度及岩石孔隙度。
基于一方面,在一种可能的实现方式中,所述使用计算机绘制H与Ln△t关系图,包括:
H与Ln△t关系图基于压实理论,在地层岩石中传播时间相等的区域基岩压实程度相同、骨架应力相等;
利用计算机绘图绘制的H与Ln△t关系图上呈线性变化关系;
通过该直线在两坐标轴上截距可求出k和Ln△t0的值。
基于一方面,在一种可能的实现方式中,所述在计算机根据斜线和点状曲线在关系图中确定各个点异常地层压力,包括:
由A点作纵轴的平行线,交压实趋势线于点B;
△tA=△tB,ΦA=ΦB,PA=PB;
其中,△tA为A点声波传播时间差,△tB为B点声波传播时间差,ΦA为A点泥岩孔隙度,ΦB为B点泥岩孔隙度,PA为A点地层压力,PB为B点地层压力;
A点较B点所增加的上覆地层压力即A点与B点之间的覆盖层压力就加在了A点的孔隙流体之上,引起流储层异常地层压力。
另一方面,提供一种储层异常地层压力研究系统,包括:
测量单元,使用超声波设备测量地层中H深度范围内声波传播时间差△t;
处理单元,通过静水压力模型获取地层中H深度范围内各深度点的静岩压力值;
绘图单元,使用计算机绘制H与Ln△t关系图,其中:关系图中斜线表示静岩压力值,关系图中点状曲线表示深度H与Ln声波传播时间差△t变化曲线;
计算单元,在计算机根据斜线和点状曲线在关系图中确定各个点异常地层压力。
基于另一方面,在一种可能的实现方式中,所述测量单元是利用声波测井原理测量的,包括:
泥岩孔隙度Φ随深度H的增加而减小,其减小程度呈指数关系变化,数学表达式为:
dΦ/Φ = -kdH (1)
Φ =Φ0e –k H (2)
Φ =(△t-△tma)/(△tf -△tma) (3)
式中(1)-(3)中: Φ-地层孔隙度;Φ0-上覆地层孔隙度;△t-地层声波时差,单位μs/m;△tf-孔隙流体声波时差,单位μs/m;△tma-岩石颗粒声波时差,单位μs/m;k-斜线斜率;H-目的层深度,单位m。
基于另一方面,在一种可能的实现方式中,所述静水压力模型的计算公式为:
PH = ρw•H (4)
公式(4)中,PH-静水压力;ρw -水的平均密度; H-水柱高度;
所述静岩压力计算公式为:
Pob = Pe + P p (5)
公式(5)中,Pob-静岩压力,静岩压力的值随深度增加而增加;Pe-粒间接触所施加应力,粒间接触所施加应力等于上覆岩柱的重量;Pp-孔隙流体压力,孔隙流体压力取决于岩石骨架密度、孔隙流体密度及岩石孔隙度。
与现有技术相比,本发明的有益效果是:
本发明通过泥岩声波时差确定孔隙压力,根据静水压力计算公式确定静水压力,将孔隙压力和静水压力绘制在关系图,根据关系图确定各点异常地层压力,能够准确方便的预测出异常地层压力。
附图说明
图1为本发明提供的一种储层异常地层压力研究方法流程图;
图2为本发明提供的H与Ln△t关系图;
图3为本发明提供的静水压力值与静岩压力值对比图示。
具体实施方式
下面结合附图所示的各实施方式对本发明进行详细说明,但应当说明的是,这些实施方式并非对本发明的限制,本领域普通技术人员根据这些实施方式所作的功能、方法、或者结构上的等效变换或替代,均属于本发明的保护范围之内。
现有的异常压力研究方法有以下三种:(1)直接利用等效泥浆比重、反应岩石可钻性的dc指数、试油测压资料、钻井过程中的中途测试资料、重复地层测试资料等,这类资料不必经过计算处理就可以获得地层压力参数,其缺点是钻井分布范围有限,难以获得纵向上和横向上的全局面貌。此外,还存在钻井深度较深,测试结果不准确等缺点。(2)利用地球物理测井资料,主要是采用声波时差和电阻率曲线进行地层流体压力计算,该方法仅限于钻井中的泥岩层段。(3)利用地震资料,即利用地震速度谱、波阻抗以及振幅来预测异常地层压力。
在描述本发明实施例的方案之前,请参阅图3,孔隙压力与静岩压力之差称为异常地层压力。正值为异常高压(超压),负值为异常低压(负压)。静岩压力与孔隙压力之差称为有效应力,表示由岩石颗粒所支撑的那部分静岩压力。在孔隙压力为常压的情况下,有效应力随深度的增加而增大。异常高压介于静水压力与静岩压力之间。孔隙压力与同深度静水压力的比值称为压力系数。当地层压力达到某一临界值时地层将破裂,将这一临界值称为破裂压力。由于超强压力可以引起地层破裂。
孔隙压力即岩石孔隙空间中的流体压力,静水压力指的是与岩石表面及地表连通的开放体系下的水柱压力。通常具有以下性质:①压力值随深度而增加;②压力变化率只依赖于水的密度变化;③代表压力增加最大速率的方向向量总是垂直的;④压力-深度的关系与流体容器的形状完全无关。
本发明实施例不同于上述三种现有的异常压力研究方法,本发明实施例具体为:
请参阅图1和图2,本发明实施例提供一种储层异常地层压力研究方法,包括:
使用超声波设备测量地层中H深度范围内声波传播时间差△t;
通过静水压力模型获取地层中H深度范围内各深度点的静岩压力值;
使用计算机绘制H与Ln△t关系图,其中:关系图中斜线表示静岩压力值,关系图中点状曲线表示深度H与Ln声波传播时间差△t变化曲线;
计算机根据斜线和点状曲线在关系图中确定各个点异常地层压力。
本发明实施例的声波测井是测量在一定距离地层中声波传播时间,记录声波传播速度之倒数,单位是μs/m,其大小取决于岩性、压实程度、孔隙度及孔隙空间流体含量。在岩性、地层水性质变化不大时,声波时差主要反映地层孔隙度大小,而上覆地层压力为孔隙度的单值函数,因此,可将声波时差转换为地层压力。
请参阅图2,本发明实施例的储层异常地层压力研究方法中,储层异常地层压力分类主要依据地层压力梯度(Pa/m)或压力系数,压力系数(Pc)是指地层压力与静水压力的比值。正常地层压力与异常压力之间的压力递变带称为压力过渡带。
请参阅图3,本发明实施例的静水压力模型获取静水压力,静水压力的计算公式为:
PH = ρw·H
式中,PH-静水压力;ρw -水的平均密度; H-水柱高度。
静岩压力:由某一深度的上覆岩石和岩石中的流体所形成的压力。即:
Pob = Pe + P p
式中Pob代表静岩压力,该值随深度的增加而增加;Pe是粒间接触所施加的应力,等于上覆岩柱的重量;Pp是孔隙流体压力。该压力主要取决于岩石骨架密度、孔隙流体密度及岩石孔隙度。对于整个沉积层段来说,沉积岩的密度不是固定不变的,而是随着埋深的增加,由低到高,直至大于2.6 g/cm3。而孔隙流体的密度与温度、压力、流体成分(包括溶解气和固体)及流体相态(气体或液体)有关,所以其变化范围也较大。
孔隙压力与静岩压力之差称为异常地层压力。正值为异常高压(超压),负值为异常低压(负压)。静岩压力与孔隙压力之差称为有效应力,表示由岩石颗粒所支撑的那部分静岩压力。在孔隙压力为常压的情况下,有效应力随深度的增加而增大。异常高压介于静水压力与静岩压力之间。孔隙压力与同深度静水压力的比值称为压力系数。当地层压力达到某一临界值时地层将破裂,将这一临界值称为破裂压力。由于超强压力可以引起地层破裂,因此,多数情况下地层压力低于静岩压力,即上覆地层压力。
基于上述实施例,在上述方案中,测量地层中H深度范围内声波传播时间差△t,包括:
使用超声波设备测量声波传播时间差△t是利用声波测井原理测量的;
声波传播时间差△t匹配有地层孔隙度Φ;
根据声波传播时间差△t和地层孔隙度Φ确定孔隙压力Ln△t;
孔隙压力Ln△t随H深度呈指数关系变化。
基于上述实施例,在上述方案中,使用超声波设备测量声波传播时间差△t是利用声波测井原理测量的,包括:
泥岩孔隙度Φ随深度H的增加而减小,其减小程度呈指数关系变化,数学表达式为:
dΦ/Φ = -kdH (1)
Φ =Φ0e –k H (2)
Φ =(△t-△tma)/(△tf -△tma) (3)
式中(1)-公式(3)中: Φ-地层孔隙度;Φ0-上覆地层孔隙度;△t-地层声波时差,单位μs/m;△tf-孔隙流体声波时差,单位μs/m;△tma-岩石颗粒声波时差,单位μs/m;k-斜线斜率;H-目的层深度,单位m。
基于上述实施例,在上述方案中,通过静水压力模型获取地层中H深度范围内各深度点的静岩压力值,包括:
静水压力模型输出连通岩石表面及地表的开放体系下的水柱压力;
参考水柱压力和静岩压力计算公式确定地层中H深度范围内各深度点的静岩压力值;
静岩压力值随着埋深增加而增加。
优选静水压力为连通岩石表面及地表的开放体系下的水柱压力;静岩压力由某一深度的上覆岩石和岩石中的流体所形成的压力;静水压力的压力值随深度值而增加;静水压力的压力变化率仅依赖于水的密度变化;静水压力代表压力增加最大速率的方向向量总是垂直的。
基于上述实施例,在上述方案中,静水压力模型的计算公式为:
PH = ρw•H (4)
公式(4)中,PH-静水压力;ρw -水的平均密度; H-水柱高度;
静岩压力计算公式为:
Pob = Pe + P p (5)
公式(5)中,Pob-静岩压力,静岩压力的值随深度增加而增加;Pe-粒间接触所施加应力,粒间接触所施加应力等于上覆岩柱的重量;Pp-孔隙流体压力,孔隙流体压力取决于岩石骨架密度、孔隙流体密度及岩石孔隙度。
基于上述实施例,在上述方案中,使用计算机绘制H与Ln△t关系图,包括:
H与Ln△t关系图基于压实理论,在地层岩石中传播时间相等的区域基岩压实程度相同、骨架应力相等;
利用计算机绘图绘制的H与Ln△t关系图上呈线性变化关系;
通过该直线在两坐标轴上截距可求出k和Ln△t0的值。
需指出,由于本发明实施例需要采用计算机绘制图2所示的H与Ln△t关系图,故本发明实施例研究的储层异常地层压力研究方法基于以下的假设条件:①不考虑地层水热作用对欠压实地层中异常压力的贡献;②不考虑地层矿物组分,尤其是粘土矿物在地质埋藏过程中的脱水作用对欠压实地层中异常压力的贡献;③假设沉积地层的性质基本相同;④假设不同深度上岩石的孔隙度相同,则其岩石颗粒间的骨架应力也相同。
基于上述实施例,在上述方案中,在计算机根据斜线和点状曲线在关系图中确定各个点异常地层压力,包括:
由A点作纵轴的平行线,交压实趋势线于点B;
△tA=△tB,ΦA=ΦB,PA=PB;
其中,△tA为A点声波传播时间差,△tB为B点声波传播时间差,ΦA为A点泥岩孔隙度,ΦB为B点泥岩孔隙度,PA为A点地层压力,PB为B点地层压力;
A点较B点所增加的上覆地层压力即A点与B点之间的覆盖层压力就加在了A点的孔隙流体之上,引起流储层异常地层压力。
需要说明的是:在正常压实情况下,泥岩孔隙度Φ随深度H的增加而减小,其减小程度呈指数关系变化,数学表达式为:
dΦ/Φ = -kdH (a)
对上式进行积分可得:
Φ =Φ0e –k H (b)
又Φ =(△t-△tma)/(△tf -△tma) (c)
式中: Φ,Φ0-地层、上覆地层孔隙度;△t-地层声波时差,μs/m;△tf-孔隙流体声波时差,μs/m;△tma-岩石颗粒声波时差,μs/m;k-正常压实趋势线斜率;H-目的层深度,m。
特殊情况下,△t和△tf对于同一区域可视为常数,孔隙度Φ与声波时差△t存在着线性关系。(c)式可变为:
Φ=△t/(△tf -△tma)-△tma/(△tf -△tma)= A·△t + B(d)
由(b)、(d)可得:
△t =△t0e–k H (e)
式中:△t0-深度H为零时泥岩的声波时差,μs/m;k-正常压实趋势线斜率。
对上式两边取对数则有:
H = -1/k(Ln△t-Ln△t0) (f)
由此可以看出,在正常压实情况下,泥岩的声波时差随着深度的增加而呈指数性减小,请参阅图2,利用计算机绘图绘制的H与Ln△t关系图上呈线性变化关系,通过该直线在两坐标轴上截距可求出k和Ln△t0的值。
压实理论,认为在地层岩石中传播时间相等的区域,其基岩压实程度也相同,骨架应力相等。具体作法是由A点作纵轴的平行线,交压实趋势线于点B,由于△tA=△tB,所以ΦA=ΦB,PA=PB,换言之,A点较B点所增加的上覆地层压力,即A点与B点之间的覆盖层压力就加在了A点的孔隙流体之上,从而引起了流体的异常压力。请参阅图2,整个压实曲线均包括正常压实段和异常压实段,计算方法如下:
在A点,PA=ρwgHB+ρsg (HA-HB) (g)
式(g)中:PA -A点处地层压力;HA、HB-AB两点的深度;ρw-地层水的密度;ρs-上覆沉积岩的平均密度;g-重力加速度。取上覆岩柱的平均密度为2.31g/㎝3,地层水密度为1.02g/㎝3。由此可以计算出目标层位的地层流体压力。
A点剩余压力:PpA=PA-ρwgHB (h)。
本发明的另一实施例提供一种储层异常地层压力研究系统,包括:
测量单元,使用超声波设备测量地层中H深度范围内声波传播时间差△t;
处理单元,通过静水压力模型获取地层中H深度范围内各深度点的静岩压力值;
绘图单元,使用计算机绘制H与Ln△t关系图,其中:关系图中斜线表示静岩压力值,关系图中点状曲线表示深度H与Ln声波传播时间差△t变化曲线;
计算单元,在计算机根据斜线和点状曲线在关系图中确定各个点异常地层压力。
基于另一实施例,在上述方案中,测量单元是利用声波测井原理测量的,包括:
泥岩孔隙度Φ随深度H的增加而减小,其减小程度呈指数关系变化,数学表达式为:
dΦ/Φ = -kdH (1)
Φ =Φ0e –k H (2)
Φ =(△t-△tma)/(△tf -△tma) (3)
式中(1)-(3)中: Φ-地层孔隙度;Φ0-上覆地层孔隙度;△t-地层声波时差,单位μs/m;△tf-孔隙流体声波时差,单位μs/m;△tma-岩石颗粒声波时差,单位μs/m;k-斜线斜率;H-目的层深度,单位m。
基于另一实施例,在上述方案中,静水压力模型的计算公式为:
PH = ρw•H (4)
公式(4)中,PH-静水压力;ρw -水的平均密度; H-水柱高度;
静岩压力计算公式为:
Pob = Pe + P p (5)
公式(5)中,Pob-静岩压力,静岩压力的值随深度增加而增加;Pe-粒间接触所施加应力,粒间接触所施加应力等于上覆岩柱的重量;Pp-孔隙流体压力,孔隙流体压力取决于岩石骨架密度、孔隙流体密度及岩石孔隙度。
需要说明的是,本发明实施例的测量单元在读取测井数据时,要综合各类电测曲线并结合录井资料,尽可能选择厚层纯泥岩段,如自然伽马处于基线、低电阻、高电导率、无严重扩径现象,单井泥岩厚度要大于2m时读取声波时差数据。然后,建立各单井正常压实曲线。在正常压实带中,声波时差传播时间与深度呈指数关系。数据的准确读取是获得正确结果的前提。
本发明实施例通过泥岩声波时差确定孔隙压力,根据静水压力计算公式确定静水压力,将孔隙压力和静水压力绘制在关系图,根据关系图确定各点异常地层压力,能够准确方便的预测出异常地层压力。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (8)
1.一种储层异常地层压力研究方法,其特征在于,包括:
使用超声波设备测量地层中H深度范围内声波传播时间差△t;
通过静水压力模型获取地层中H深度范围内各深度点的静岩压力值;
使用计算机绘制H与ln△t关系图,其中:关系图中纵坐标为H,横坐标为ln△t,斜线表示静岩压力值,关系图中点状曲线表示深度H与ln△t变化曲线;
采用计算机根据斜线和点状曲线在关系图中确定各个点异常地层压力;
所述使用计算机绘制H与ln△t关系图,包括:
H与ln△t关系图基于压实理论,在地层岩石中传播时间相等的区域基岩压实程度相同、骨架应力相等;
利用计算机绘图绘制的H与ln△t关系图上呈线性变化关系;
通过该直线在两坐标轴上截距可求出k和ln△t0的值;
所述采用计算机根据斜线和点状曲线在关系图中确定各个点异常地层压力,包括:
基于压实理论,在关系图确定正常压实带和异常压实带,在异常压实带从点状曲线取A点;
由A点作纵轴的平行线,交正常压实带的点状曲线于点B;
△tA=△tB,ΦA=ΦB,PA=PB;
其中,△tA为A点声波传播时间差,△tB为B点声波传播时间差,ΦA为A点地层孔隙度,ΦB为B点地层孔隙度,PA为A点地层压力,PB为B点地层压力;
A点较B点所增加的上覆地层压力,即A点与B点之间的覆盖层压力,就加在了A点的孔隙流体之上,引起储层异常地层压力。
2.根据权利要求1所述的一种储层异常地层压力研究方法,其特征在于,所述测量地层中H深度范围内声波传播时间差△t,包括:
使用超声波设备测量声波传播时间差△t是利用声波测井原理测量的;
声波传播时间差△t匹配有地层孔隙度Φ;
根据声波传播时间差△t和地层孔隙度Φ确定孔隙压力;
孔隙压力随H深度呈指数关系变化。
3.根据权利要求2所述的一种储层异常地层压力研究方法,其特征在于,所述使用超声波设备测量声波传播时间差△t是利用声波测井原理测量的,包括:
地层孔隙度Φ随深度h的增加而减小,其减小程度呈指数关系变化,数学表达式为:
dΦ/Φ = -kdh (1)
Φ =Φ0e –kh (2)
Φ =(△t-△tma)/(△tf -△tma) (3)
式中(1)-(3)中: Φ-地层孔隙度;Φ0-上覆地层孔隙度;△t-地层声波时差,单位μs/m;△tf-孔隙流体声波时差,单位μs/m;△tma-岩石颗粒声波时差,单位μs/m;k-斜线斜率;h-目的层深度,单位m。
4.根据权利要求1所述的一种储层异常地层压力研究方法,其特征在于,所述通过静水压力模型获取地层中H深度范围内各深度点的静岩压力值,包括:
静水压力模型输出连通岩石表面及地表的开放体系下的水柱压力;
参考水柱压力和静岩压力计算公式确定地层中H深度范围内各深度点的静岩压力值;
静岩压力值随着埋深增加而增加。
5.根据权利要求4所述的一种储层异常地层压力研究方法,其特征在于,所述静水压力模型的计算公式为:
PH = ρw•h' (4)
公式(4)中,PH-静水压力;ρw -水的平均密度; h'-水柱高度;
所述静岩压力计算公式为:
Pob = Pe + P p (5)
公式(5)中,Pob-静岩压力,静岩压力的值随深度增加而增加;Pe-粒间接触所施加应力,粒间接触所施加应力等于上覆岩柱的重量;Pp-孔隙流体压力,孔隙流体压力取决于岩石骨架密度、孔隙流体密度及岩石孔隙度。
6.一种储层异常地层压力研究系统,其特征在于,包括:
测量单元,使用超声波设备测量地层中H深度范围内声波传播时间差△t;
处理单元,通过静水压力模型获取地层中H深度范围内各深度点的静岩压力值;
绘图单元,使用计算机绘制H与ln△t关系图,其中:关系图中纵坐标为H,横坐标为ln△t,关系图中斜线表示静岩压力值,关系图中点状曲线表示深度H与ln△t变化曲线;
计算单元,在计算机根据斜线和点状曲线在关系图中确定各个点异常地层压力;
所述使用计算机绘制H与ln△t关系图,包括:
H与ln△t关系图基于压实理论,在地层岩石中传播时间相等的区域基岩压实程度相同、骨架应力相等;
利用计算机绘图绘制的H与ln△t关系图上呈线性变化关系;
通过该直线在两坐标轴上截距可求出k和ln△t0的值;
所述采用计算机根据斜线和点状曲线在关系图中确定各个点异常地层压力,包括:
基于压实理论,在关系图确定正常压实带和异常压实带,在异常压实带从点状曲线取A点;
由A点作纵轴的平行线,交正常压实带的点状曲线于点B;
△tA=△tB,ΦA=ΦB,PA=PB;
其中,△tA为A点声波传播时间差,△tB为B点声波传播时间差,ΦA为A点地层孔隙度,ΦB为B点地层孔隙度,PA为A点地层压力,PB为B点地层压力;
A点较B点所增加的上覆地层压力,即A点与B点之间的覆盖层压力,就加在了A点的孔隙流体之上,引起储层异常地层压力。
7.根据权利要求6所述的一种储层异常地层压力研究系统,其特征在于,所述测量单元是利用声波测井原理测量的,包括:
地层孔隙度Φ随深度h的增加而减小,其减小程度呈指数关系变化,数学表达式为:
dΦ/Φ = -kdh (1)
Φ =Φ0e –kh (2)
Φ =(△t-△tma)/(△tf -△tma) (3)
式中(1)-(3)中: Φ-地层孔隙度;Φ0-上覆地层孔隙度;△t-地层声波时差,单位μs/m;△tf-孔隙流体声波时差,单位μs/m;△tma-岩石颗粒声波时差,单位μs/m;k-斜线斜率;h-目的层深度,单位m。
8.根据权利要求6所述的一种储层异常地层压力研究系统,其特征在于,所述静水压力模型的计算公式为:
PH = ρw•h' (4)
公式(4)中,PH-静水压力;ρw -水的平均密度; h'-水柱高度;
所述静岩压力计算公式为:
Pob = Pe + P p (5)
公式(5)中,Pob-静岩压力,静岩压力的值随深度增加而增加;Pe-粒间接触所施加应力,粒间接触所施加应力等于上覆岩柱的重量;Pp-孔隙流体压力,孔隙流体压力取决于岩石骨架密度、孔隙流体密度及岩石孔隙度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310203688.2A CN116068616B (zh) | 2023-03-06 | 2023-03-06 | 一种储层异常地层压力研究方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310203688.2A CN116068616B (zh) | 2023-03-06 | 2023-03-06 | 一种储层异常地层压力研究方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116068616A CN116068616A (zh) | 2023-05-05 |
CN116068616B true CN116068616B (zh) | 2023-06-20 |
Family
ID=86180318
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310203688.2A Active CN116068616B (zh) | 2023-03-06 | 2023-03-06 | 一种储层异常地层压力研究方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116068616B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109917489A (zh) * | 2019-03-22 | 2019-06-21 | 西北大学 | 一种地下承压水位确定的新方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6751558B2 (en) * | 2001-03-13 | 2004-06-15 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
FR2893421B1 (fr) * | 2005-11-14 | 2007-12-21 | Inst Francais Du Petrole | Methode d'evaluation quantitative des pressions de fluides et de detection des surpressions d'un milieu souterrain. |
US9279897B2 (en) * | 2005-12-12 | 2016-03-08 | Hannes Zuercher | Locating oil or gas passively by observing a porous oil and gas saturated system giving off its characteristic resonance response to ambient background noise, including optional differentiation of oil, locatinggas and water |
US10385678B2 (en) * | 2014-03-21 | 2019-08-20 | Conocophillips Company | Method for analysing pore pressure in shale formations |
CN104977618B (zh) * | 2014-04-09 | 2017-09-22 | 中国石油集团东方地球物理勘探有限责任公司 | 一种评价页岩气储层及寻找甜点区的方法 |
CN104853822A (zh) * | 2014-09-19 | 2015-08-19 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
CN106285622B (zh) * | 2015-06-01 | 2020-05-15 | 中国石油化工股份有限公司 | 用于校正压实曲线的方法 |
CN112651098A (zh) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | 基于胡克定律的地层压力计算方法及系统 |
CN115220095A (zh) * | 2021-04-14 | 2022-10-21 | 中国石油化工股份有限公司 | 含水圈闭静动态相结合的断层垂向封闭性评价方法 |
CN115680638A (zh) * | 2021-07-26 | 2023-02-03 | 中国石油化工股份有限公司 | 一种利用压力衰减梯度识别超压顶封层的方法 |
CN113790044B (zh) * | 2021-09-08 | 2023-06-30 | 西安石油大学 | 一种多机制复合构成超压地层孔隙压力的评价方法 |
-
2023
- 2023-03-06 CN CN202310203688.2A patent/CN116068616B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109917489A (zh) * | 2019-03-22 | 2019-06-21 | 西北大学 | 一种地下承压水位确定的新方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116068616A (zh) | 2023-05-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mikhailov et al. | Using borehole electroseismic measurements to detect and characterize fractured (permeable) zones | |
Vernik et al. | Velocity anisotropy in shales: A petrophysical study | |
Geertsma et al. | Some aspects of elastic wave propagation in fluid-saturated porous solids | |
RU2602409C2 (ru) | Инверсия анизотропии многоскважинной системы | |
Wang et al. | Poroelasticity: Efficient modeling of strongly coupled, slow deformation processes in a multilayered half-space | |
Nur et al. | Critical porosity: A key to relating physical properties to porosity in rocks | |
Kuster et al. | Velocity and attenuation of seismic waves in two-phase media: Part I. Theoretical formulations | |
Wang et al. | Pore pressure prediction using geophysical methods in carbonate reservoirs: Current status, challenges and way ahead | |
Hu et al. | Acoustic response of gas hydrate formation in sediments from South China Sea | |
CN111460602A (zh) | 一种基于岩石物理建模的横观各向同性地层地应力预测方法 | |
CN106295119B (zh) | 一种页岩气地层地应力计算方法 | |
CN103713327B (zh) | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 | |
Jeppson et al. | San Andreas fault zone velocity structure at SAFOD at core, log, and seismic scales | |
CN101285381B (zh) | 一种泄漏模式波反演软地层横波速度的方法 | |
AU3396199A (en) | Method of imaging the permeability and fluid content structure within sediment | |
King | Recent developments in seismic rock physics | |
CN107703543B (zh) | 一种基于广义泊松比的地层孔隙压力预测方法及系统 | |
ZHAO et al. | Seismic rock physics characteristics and modeling of ultra-deep carbonate reservoirs | |
CN116068616B (zh) | 一种储层异常地层压力研究方法及系统 | |
Han et al. | Understanding how overpressure affects the physical properties of sandstones | |
King | Rock-physics developments in seismic exploration: A personal 50-year perspective | |
RU2761935C1 (ru) | Способ локализации перспективных зон в нефтематеринских толщах | |
CN113589385B (zh) | 一种基于地震散射波场分析的储层特征反演方法 | |
CN112965106B (zh) | 利用纵横波速度的奇异性识别天然气水合物赋存类型方法 | |
YE et al. | Experimental research on the relationship between gas hydrate saturation and acoustic parameters |
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 |