CN110069751A - 一种非欠压实地层孔隙压力预测方法 - Google Patents
一种非欠压实地层孔隙压力预测方法 Download PDFInfo
- Publication number
- CN110069751A CN110069751A CN201910373439.1A CN201910373439A CN110069751A CN 110069751 A CN110069751 A CN 110069751A CN 201910373439 A CN201910373439 A CN 201910373439A CN 110069751 A CN110069751 A CN 110069751A
- Authority
- CN
- China
- Prior art keywords
- formation
- pore pressure
- mud stone
- pressure
- density
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Computer Hardware Design (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种非欠压实地层孔隙压力预测方法。利用现场钻井来建立纵波速度与密度交汇图,根据现场钻井的不同深度的泥岩,在纵波速度与密度交汇图中建立这些泥岩的速度与密度交汇点,然后将偏离加载趋势线的速度与密度交汇点和经典压力卸载模型进行对比,通过对比来判断非欠压地层异常压力成因,利用地层速度计算不同深度的泥岩孔隙压力,统计计算的泥岩孔隙压力和实测的泥岩孔隙压力之间的压力系数误差,利用地层深度与压力系数误差来建立拟合函数,利用拟合函数得到的泥岩孔隙压力与得到的泥岩孔隙压力来实现非欠压实地层孔隙压力预测。本发明的有益效果:实现非欠压实地层孔隙压力预测,实现非欠压地层异常压力成因的定性判断。
Description
技术领域
本发明涉及一种非欠压实地层孔隙压力预测方法。
背景技术
国内外地层孔隙压力预测技术发展,包括20世纪70年代至80年代,研究者Fillippone提出不依赖余正常压实速度趋势,而利用地层层速度预测压力的Fillione公式;研究者Bowers在1995年指出:由于不同的增压作用和压力演化过程会使地层孔隙度经常不随着地层压力而变化;研究者Djevanshir和研究者Akhnerdiev在1998年利用已钻井的实际测量压力来建立经验关系,并对等效深度法的估算结果进行修正,以及根据各自研究地区的情况提出改进的公式等;在1990年,国内研究者通过对辽东湾辽西凹陷压力测试数据的分析中发现,地层压力和速度之间的关系不是简单的线性内插关系,提出了修正后的Fillippone公式,进一步提高了解释精度。
现阶段,国内外已经提出了许多预测地层孔隙压力的方法,但是基本都是基于传统的泥岩欠压实模式,对于非欠压实地层的压力成因只能定性分析而无法定量计算,随着钻井勘探领域和层系的扩大,欠压实地层孔隙压力预测技术已经无法实现高精度的预测地层孔隙压力,极大制约了油气勘探的进程。
发明内容
本发明的目的在于提供一种非欠压实地层孔隙压力预测方法,能够在定性判断非欠压地层异常压力成因的基础上实现非欠压地层孔隙压力的定量预测,在高温超压的地层钻前和随钻过程中实现高精度的地层孔隙压力预测。
为了实现上述目的,本发明采用如下技术方案:
一种非欠压实地层孔隙压力预测方法,包括以下步骤:
步骤一:利用现场钻井来建立纵波速度与密度交汇图,根据现场钻井的不同深度的泥岩,在纵波速度与密度交汇图中建立这些泥岩的速度与密度交汇点,然后将偏离加载趋势线的速度与密度交汇点和经典压力卸载模型进行对比,通过两者的对比来判断非欠压地层异常压力成因;
步骤二:利用地层速度计算不同深度的泥岩孔隙压力;
步骤三:统计在所述的步骤二中得到的泥岩孔隙压力和实际测量的泥岩孔隙压力之间的压力系数误差;
步骤四:利用地层深度与在所述的步骤三中得到的压力系数误差来建立拟合函数;
步骤五:利用拟合函数得到的泥岩孔隙压力与在所述的步骤二中得到的泥岩孔隙压力来实现非欠压实地层孔隙压力预测。
作为本发明的进一步方案,在所述的步骤一中,不同深度的泥岩的速度与密度交汇点明显偏离加载趋势线,不同深度的泥岩的速度与密度交汇点和经典压力卸载模型线接近。
作为本发明的进一步方案,在所述的步骤一中,非欠压实地层异常压力成因是水热增压或者生烃增压。
本发明的有益效果如下:算法简单,能够在地层钻前和随钻过程中实现地层孔隙压力预测;操作方便,能够利用自主研发的软件实现所需的计算;应用前景良好,能够广泛应用于高温超压的非欠压实地层,实现非欠压实地层孔隙压力预测;能够实现非欠压地层异常压力成因的定性判断,对地层的适应性强;拟合函数能通过新钻井数据更新完善。
附图说明
图1是本发明实施例中纵波速度与密度交汇图。
图2是本发明实施例中速度谱速度计算压力误差函数图。
图3是本发明实施例中声波速度计算压力误差函数图。
具体实施方式
下面将结合的附图,对本发明实施例中的技术方案进行描述。
一种非欠压实地层孔隙压力预测方法,包括以下步骤:
步骤一:利用现场钻井来建立纵波速度与密度交汇图,根据现场钻井的不同深度的泥岩,在纵波速度与密度交汇图中建立这些泥岩的速度与密度交汇点,然后将偏离加载趋势线的速度与密度交汇点和经典压力卸载模型进行对比,通过两者的对比来判断非欠压地层异常压力成因;
步骤二:利用地层速度计算不同深度的泥岩孔隙压力;
步骤三:统计在所述的步骤二中得到的泥岩孔隙压力和实际测量的泥岩孔隙压力之间的压力系数误差;
步骤四:利用地层深度与在所述的步骤三中得到的压力系数误差来建立拟合函数;
步骤五:利用拟合函数得到的泥岩孔隙压力与在所述的步骤二中得到的泥岩孔隙压力来实现非欠压实地层孔隙压力预测。
实施例
图1为利用现场钻井建立的纵波速度与密度交汇图,在图1中,曲线为现场钻井目的地层的加载趋势线,大部分在加载趋势线附近的泥岩速度与密度交汇点均为欠压实成因,明显偏离加载趋势线且与经典压力卸载模型线接近的泥岩速度与密度交汇点均为非欠压实异常压力成因,该非欠压实异常压力成因为水热增压或者生烃增压。
图2和图3都为速度与压力的误差交汇图,利用地层速度计算不同深度的泥岩孔隙压力,统计不同深度的泥岩孔隙压力和实际测量的泥岩孔隙压力之间的压力系数误差,通过压力系数误差和地层深度来构建拟合函数,最后利用拟合函数预测非欠压实异常压力和欠压实异常压力。
Claims (3)
1.一种非欠压实地层孔隙压力预测方法,其特征是包括以下步骤:
步骤一:利用现场钻井来建立纵波速度与密度交汇图,根据现场钻井的不同深度的泥岩,在纵波速度与密度交汇图中建立这些泥岩的速度与密度交汇点,然后将偏离加载趋势线的速度与密度交汇点和经典压力卸载模型进行对比,通过两者的对比来判断非欠压地层异常压力成因;
步骤二:利用地层速度计算不同深度的泥岩孔隙压力;
步骤三:统计在所述的步骤二中得到的泥岩孔隙压力和实际测量的泥岩孔隙压力之间的压力系数误差;
步骤四:利用地层深度与在所述的步骤三中得到的压力系数误差来建立拟合函数;
步骤五:利用拟合函数得到的泥岩孔隙压力与在所述的步骤二中得到的泥岩孔隙压力来实现非欠压实地层孔隙压力预测。
2.根据权利要求1所述的非欠压实地层孔隙压力预测方法,其特征是在所述的步骤一中,不同深度的泥岩的速度与密度交汇点明显偏离加载趋势线,不同深度的泥岩的速度与密度交汇点和经典压力卸载模型线接近。
3.根据权利要求2所述的非欠压实地层孔隙压力预测方法,其特征是在所述的步骤一中,非欠压实地层异常压力成因是水热增压或者生烃增压。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910373439.1A CN110069751B (zh) | 2019-05-07 | 2019-05-07 | 一种非欠压实地层孔隙压力预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910373439.1A CN110069751B (zh) | 2019-05-07 | 2019-05-07 | 一种非欠压实地层孔隙压力预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110069751A true CN110069751A (zh) | 2019-07-30 |
CN110069751B CN110069751B (zh) | 2022-05-03 |
Family
ID=67370004
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910373439.1A Active CN110069751B (zh) | 2019-05-07 | 2019-05-07 | 一种非欠压实地层孔隙压力预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110069751B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110529104A (zh) * | 2019-08-29 | 2019-12-03 | 中国海洋石油集团有限公司 | 一种改进的异常高压泥岩速度拟合方法 |
CN112034521A (zh) * | 2020-09-07 | 2020-12-04 | 中国海洋石油集团有限公司 | 一种欠压实与生烃混合成因地层超压的预测方法 |
CN112925018A (zh) * | 2019-12-06 | 2021-06-08 | 中国石油天然气股份有限公司 | 一种异常压力成因判别方法及装置 |
CN113187463A (zh) * | 2021-04-14 | 2021-07-30 | 中海石油(中国)有限公司海南分公司 | 一种基于地层超压单成因贡献率的随钻孔隙压力预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080033704A1 (en) * | 2006-08-07 | 2008-02-07 | Schlumberger Technology Corporation | Method and system for pore pressure prediction |
CN104863577A (zh) * | 2015-04-09 | 2015-08-26 | 中国石油大学(北京) | 利用地震纵波传播时间预测地层孔隙压力的方法 |
CN107703543A (zh) * | 2017-09-11 | 2018-02-16 | 武汉拓盟能源科技有限公司 | 一种基于广义泊松比的地层孔隙压力预测方法及系统 |
CN108150160A (zh) * | 2017-12-21 | 2018-06-12 | 西安石油大学 | 一种地层中欠压实作用超压的求取方法 |
CN108301825A (zh) * | 2018-01-23 | 2018-07-20 | 中国海洋石油集团有限公司 | 一种高温高压储层孔隙压力分析方法 |
-
2019
- 2019-05-07 CN CN201910373439.1A patent/CN110069751B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080033704A1 (en) * | 2006-08-07 | 2008-02-07 | Schlumberger Technology Corporation | Method and system for pore pressure prediction |
CN104863577A (zh) * | 2015-04-09 | 2015-08-26 | 中国石油大学(北京) | 利用地震纵波传播时间预测地层孔隙压力的方法 |
CN107703543A (zh) * | 2017-09-11 | 2018-02-16 | 武汉拓盟能源科技有限公司 | 一种基于广义泊松比的地层孔隙压力预测方法及系统 |
CN108150160A (zh) * | 2017-12-21 | 2018-06-12 | 西安石油大学 | 一种地层中欠压实作用超压的求取方法 |
CN108301825A (zh) * | 2018-01-23 | 2018-07-20 | 中国海洋石油集团有限公司 | 一种高温高压储层孔隙压力分析方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110529104A (zh) * | 2019-08-29 | 2019-12-03 | 中国海洋石油集团有限公司 | 一种改进的异常高压泥岩速度拟合方法 |
CN112925018A (zh) * | 2019-12-06 | 2021-06-08 | 中国石油天然气股份有限公司 | 一种异常压力成因判别方法及装置 |
CN112925018B (zh) * | 2019-12-06 | 2023-09-26 | 中国石油天然气股份有限公司 | 一种异常压力成因判别方法及装置 |
CN112034521A (zh) * | 2020-09-07 | 2020-12-04 | 中国海洋石油集团有限公司 | 一种欠压实与生烃混合成因地层超压的预测方法 |
CN113187463A (zh) * | 2021-04-14 | 2021-07-30 | 中海石油(中国)有限公司海南分公司 | 一种基于地层超压单成因贡献率的随钻孔隙压力预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110069751B (zh) | 2022-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110069751A (zh) | 一种非欠压实地层孔隙压力预测方法 | |
CN110005389B (zh) | 一种基于热流固耦合作用的超深层砂岩缝网改造评价方法 | |
CN108009705B (zh) | 一种基于支持向量机技术的页岩储层可压性评价方法 | |
CN104727798B (zh) | 一种低渗透气藏转向重复压裂工艺方法 | |
Pan et al. | Coalbed methane (CBM) exploration, reservoir characterisation, production, and modelling: A collection of published research (2009–2015) | |
CN104632187B (zh) | 一种确定产水煤层气井动态储量的方法 | |
CN108894777B (zh) | 一种分层压裂多层合采油气藏储层物性及裂缝特性参数的确定方法 | |
CN110735635B (zh) | 一种页岩含气量测试中损失气含量的确定方法 | |
CN107725034B (zh) | 一种用于多级压裂水平井判别来水方向的压力监测方法 | |
CN111577256B (zh) | 一种穿层钻孔水力冲孔增透效果定量评价方法 | |
CN104695950A (zh) | 火山岩油藏产能预测方法 | |
CN104373118B (zh) | 钻孔有效抽采半径的测定方法 | |
CN109033504B (zh) | 一种油水井套管损坏预测方法 | |
CN114427432B (zh) | 一种气藏剩余气开发潜力确定方法 | |
Shi et al. | An accurate method for permeability evaluation of undersaturated coalbed methane reservoirs using early dewatering data | |
CN111927417A (zh) | 一种页岩气分段压裂水平井组储量动用状况评价方法 | |
CN112766630A (zh) | 一种评价低渗气藏气井无阻流量的方法 | |
CN112836442B (zh) | 一种页岩油水平井网老井水力裂缝注液量的确定方法 | |
CN112647930B (zh) | 一种水平井油藏工程找水方法 | |
Ding et al. | Study on the Development Height of Overburden Water‐Flowing Fracture Zone of the Working Face | |
Jang et al. | Effect of fracture design parameters on the well performance in a hydraulically fractured shale gas reservoir | |
CN109899050A (zh) | 形成复杂裂缝网络的煤层气藏压裂方法 | |
CN112035993A (zh) | 一种底部定压的碳酸盐岩储层测试评价方法 | |
Cheng et al. | Production Data Analysis and Practical Applications in the Sulige Tight Gas Reservoir, Ordos Basin, China | |
CN111982567A (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 |