CN107605469B - 预测地层孔隙压力的方法 - Google Patents
预测地层孔隙压力的方法 Download PDFInfo
- Publication number
- CN107605469B CN107605469B CN201610548005.7A CN201610548005A CN107605469B CN 107605469 B CN107605469 B CN 107605469B CN 201610548005 A CN201610548005 A CN 201610548005A CN 107605469 B CN107605469 B CN 107605469B
- Authority
- CN
- China
- Prior art keywords
- work area
- interval velocity
- velocity
- seismic interval
- functions
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及地震资料解释领域,具体公开了一种预测地层孔隙压力的方法。该方法包括:得到工区内的地震层速度v并基于测井资料得到多个位置处的伊顿指数ni;基于使E2最小的原则确定工区内伊顿指数n与地震层速度v间的关系n=p(v);在工区内,基于地震层速度v和确定的函数关系n=p(v)得到工区范围内的伊顿指数n;进而预测工区内的地层孔隙压力Pp。应用本发明,能够大大提高预测地层孔隙压力的精度。
Description
技术领域
本发明涉及地震资料解释领域,更具体地,涉及一种预测地层孔隙压力的方法。
背景技术
地层孔隙压力是指地层孔隙或裂缝中流体所具有的压力。它是合理确定套管程序的基础,也是合理选择泥浆密度,实现安全高效钻井的关键。同时地层孔隙压力是油气成藏与分布的主控因素之一,是油气成藏流体动力学研究的依据。
地层孔隙压力的研究已有40多年的历史,但是严格讲来至今并未完全得到解决。进入本世纪以来,地层孔隙压力预测在西方国家地质、钻井、测井、物探等领域再次成为研究的热点。具体的预测方法有利用地震层速度预测、测井资料解释(检测)、钻井资料解释(检测)、实测等。
这些传统方法是基于正常压实趋势线的经验半经验方法,总的来讲存在以下共同的缺限:仅适用于“不平衡压实过程导致的地层欠压实”高压的情况;绝大部分方法仅限于在纯泥页岩中使用;都需要建立正常压实趋势线,且假定半对数坐标系中为直线;因建立经验图版的压力来源于渗透性地层,反过来预测泥岩地层,结果往往偏低;在定量化方面是经验和半经验的方法,缺乏理论基础。
地层孔隙压力预测的理论基础是有效应力定理。依据该定理,地层孔隙压力等于上覆岩层压力与垂直有效应力之差。上覆岩层压力可以通过已钻井的视密度测井资料等多种途径求得,因此只要设法求出垂直有效应力即可以确定地层孔隙压力。饱和多孔介质的有效应力是物理学上无法直接测量的参数,只有通过测量其产生的结果(如孔隙度的变化)来间接计算其值。在检索到的国内外文献中,基于有效应力定理利用地震速度资料来求取孔隙压力的方法比较有代表性的有Foster&Whalen法、Eaton法、Fan简易法、基于岩石物理方法等。Foster&Whalen法是一种地区半经验方法,Fan简易法更多适用于欠压实地层,基于岩石物理的方法主要针对砂泥岩剖面且计算繁琐。因此优选Eaton法来进行地层孔隙压力预测。Eaton公式中的重要参量Eaton指数n通常是根据经验确定的,由于不同工区的情况不同,因此,根据经验确定的Eaton指数n往往不能很好地适应工区的实际情况,从而严重影响孔隙压力预测精度。
发明内容
本发明的目的是得到一种能够克服上述缺陷的利用Eaton法准确预测孔隙压力的方法。
为了实现上述目的,本发明提供一种预测地层孔隙压力的方法,该方法包括:得到工区内的地震层速度v并基于测井资料得到多个位置处的伊顿指数ni,下标i表示各个位置的编号,i=1,2,...,N,N表示所述多个位置中包含的所述位置的数量;基于使下式中E2最小的原则确定工区内伊顿指数n与地震层速度v间的关系n=p(v):其中,vi表示相应位置处的地震层速度,n=p(v)从预设的函数集合中选择;在工区内,基于地震层速度v和确定的函数关系n=p(v)得到工区范围内的伊顿指数n;基于下式预测工区内的地层孔隙压力Pp:Pp=Sv-(Sv-Ph)(v/vnorm),其中,Sv为上覆地层压力,Ph为静水压力,vnorm为正常压实情况下的层速度。
本发明中,首先优选Eaton公式法作为基于有效应力理论进行地层孔隙压力预测的有效方法;然后针对页岩气等非常规储层的实际地层复杂性,充分考虑各向异性地层对孔隙压力预测公式中关键参数的影响,从而大大提高了预测精度。
本发明的其它特征和优点将在随后具体实施方式部分予以详细说明。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施方式中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的一个实施例的预测地层孔隙压力的方法的流程图。
图2示出了根据本发明的孔隙压力预测前地震速度精细处理、反演技术的示意图。
图3示出了根据本发明的应用Gardner公式求取速度进而得到上覆地层压力的示意图。
图4示出了根据本发明确定的Eaton指数的误差率示意图。
具体实施方式
下面将参照附图更详细地描述本发明的优选实施方式。虽然附图中显示了本发明的优选实施方式,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
图1示出了一种预测地层孔隙压力的方法,该方法包括:
步骤101,得到工区内的地震层速度v并基于测井资料得到多个位置处的伊顿指数ni,下标i表示各个位置的编号,i=1,2,...,N,N表示所述多个位置中包含的所述位置的数量;
步骤102,基于使下式中E2最小的原则确定工区内伊顿指数n与地震层速度v间的关系n=p(v):
步骤103,在工区内,基于地震层速度v和确定的函数关系n=p(v)得到工区内的伊顿指数n;
步骤104,基于下式预测工区内的地层孔隙压力Pp:
Pp=Sv-(Sv-Ph)(v/vnorm),其中,Sv为上覆地层压力,Ph为静水压力,vnorm为正常压实情况下的层速度。
本实施例中,借助Eaton公式预测地层孔隙压力,并充分考虑各向异性地层对孔隙压力预测公式中关键参数(例如Eaton指数)的影响,大大提高了预测精度。
在一种实施方式中,所述预设的函数集合中可以包括下列中的部分或全部:线性拟合函数、二次函数、指数函数、幂函数、对数函数。本领域技术人员可根据经验和实际情况设置函数集合中所包含的函数形式。
在一种实施方式中,可以采用最小二乘法确定工区内伊顿指数n与地震层速度v间的关系n=p(v)。
在一种实施方式中,所述地震层速度v可以是通过下列方法得到的:进行速度谱处理及调整解释;进行剩余速度分析;进行约束速度反演;进行层速度矫正,得到所述地震层速度v。
在一种实施方式中,所述上覆地层压力Sv可以是通过下列方法得到的:基于测井资料拟合得到Gardner系数;基于所述地震层速度并利用Gardner公式得到工区内的密度;基于所述密度得到所述上覆地层压力Sv。
应用实例
图2示出了根据本发明的孔隙压力预测前地震速度精细处理、反演技术的示意图。图3示出了根据本发明的应用Gardner公式求取速度进而得到上覆地层压力的示意图。图4示出了根据本发明确定的Eaton指数的误差率示意图。
本发明可以是***、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质,其上载有用于使处理器实现本发明的各个方面的计算机可读程序指令。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (4)
1.一种预测地层孔隙压力的方法,该方法包括:
得到工区内的地震层速度v并基于测井资料得到多个位置处的伊顿指数ni,下标i表示各个位置的编号,i=1,2,...,N,N表示所述多个位置中包含的所述位置的数量;
基于使下式中E2最小的原则确定工区内伊顿指数n与地震层速度v间的关系n=p(v):
在工区内,基于地震层速度v和确定的函数关系n=p(v)得到工区范围内的伊顿指数n;
基于下式预测工区内的地层孔隙压力Pp:
Pp=Sv-(Sv-Ph)(v/vnorm)n,其中,Sv为上覆地层压力,Ph为静水压力,vnorm为正常压实情况下的层速度,n为伊顿指数;
采用最小二乘法确定工区内伊顿指数n与地震层速度v间的关系n=p(v)。
2.根据权利要求1所述的方法,其中,所述预设的函数集合中包括下列中的部分或全部:线性拟合函数、二次函数、指数函数、幂函数、对数函数。
3.根据权利要求1所述的方法,其中,所述地震层速度v是通过下列方法得到的:
进行速度谱处理及调整解释;
进行剩余速度分析;
进行约束速度反演;
进行层速度矫正,得到所述地震层速度v。
4.根据权利要求1所述的方法,其中,所述上覆地层压力Sv是通过下列方法得到的:
基于测井资料拟合得到Gardner系数;
基于所述地震层速度v并利用Gardner公式得到工区内的密度v;
基于所述密度得到所述上覆地层压力Sv。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610548005.7A CN107605469B (zh) | 2016-07-12 | 2016-07-12 | 预测地层孔隙压力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610548005.7A CN107605469B (zh) | 2016-07-12 | 2016-07-12 | 预测地层孔隙压力的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107605469A CN107605469A (zh) | 2018-01-19 |
CN107605469B true CN107605469B (zh) | 2021-03-12 |
Family
ID=61055141
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610548005.7A Active CN107605469B (zh) | 2016-07-12 | 2016-07-12 | 预测地层孔隙压力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107605469B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108562938B (zh) * | 2018-03-23 | 2019-09-06 | 中国石油天然气股份有限公司 | 一种消除频散效应的方法、装置及*** |
CN108594305A (zh) * | 2018-08-13 | 2018-09-28 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 页岩储层含气性分析方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5615115A (en) * | 1994-12-15 | 1997-03-25 | Atlantic Richfield Company | Method of determining pore pressure and fracture gradient profiles using seismic transit times |
CN105134189A (zh) * | 2015-08-24 | 2015-12-09 | 西南石油大学 | 一种LogGMIR测井地应力识别优质储层的方法 |
CN105445791A (zh) * | 2015-11-25 | 2016-03-30 | 成都理工大学 | 一种基于多种地震属性的地层孔隙压力预测方法 |
CN105626056A (zh) * | 2015-12-22 | 2016-06-01 | 中国石油天然气集团公司 | 一种确定目标区三维地层孔隙压力的方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10385678B2 (en) * | 2014-03-21 | 2019-08-20 | Conocophillips Company | Method for analysing pore pressure in shale formations |
-
2016
- 2016-07-12 CN CN201610548005.7A patent/CN107605469B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5615115A (en) * | 1994-12-15 | 1997-03-25 | Atlantic Richfield Company | Method of determining pore pressure and fracture gradient profiles using seismic transit times |
CN105134189A (zh) * | 2015-08-24 | 2015-12-09 | 西南石油大学 | 一种LogGMIR测井地应力识别优质储层的方法 |
CN105445791A (zh) * | 2015-11-25 | 2016-03-30 | 成都理工大学 | 一种基于多种地震属性的地层孔隙压力预测方法 |
CN105626056A (zh) * | 2015-12-22 | 2016-06-01 | 中国石油天然气集团公司 | 一种确定目标区三维地层孔隙压力的方法和装置 |
Non-Patent Citations (4)
Title |
---|
Eaton法预测M油田地层孔隙压力;杨振平等;《石油天然气学报》;20120930;第34卷(第9期);第181-182页 * |
利用Eaton法计算地层孔隙压力的不确定性分析;臧艳彬等;《西南石油大学学报(自然科学版)》;20120704(第4期);第55-61页 * |
利用已钻井资料构建区域地层压力剖面的方法;管志川等;《中国石油大学学报(自然科学版)》;20131031;第37卷(第5期);第71-73页 * |
孔隙压力地震预测技术综述;张荣忠等;《勘探地球物理进展》;20050430;第28卷(第2期);第90-96页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107605469A (zh) | 2018-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11591900B2 (en) | Method to predict overpressure uncertainty from normal compaction trendline uncertainty | |
CN108713089B (zh) | 基于钻孔流体和钻探录井估计地层性质 | |
US10400572B2 (en) | Apparatus and methods using drillability exponents | |
CN108301825B (zh) | 一种高温高压储层孔隙压力分析方法 | |
CN106401574B (zh) | 一种钻前高温地热井地层压力的预测方法 | |
CN104863577B (zh) | 利用地震纵波传播时间预测地层孔隙压力的方法 | |
US20150370934A1 (en) | Completion design based on logging while drilling (lwd) data | |
CN107503739B (zh) | 一种用于水平井判别来水方向的压力监测方法 | |
WO2018212781A1 (en) | Predictive lithology and formation type for downhole drilling | |
NO20140166A1 (no) | Estimering av usikkerhet ved undergrunnens poretrykk ut ifra trendlinjevariasjoner | |
US11841476B2 (en) | Methods and systems of determining parameters characterizing porous media from data gathered by a plurality of different tools | |
CN108304959B (zh) | 提高地层流体压力预测精度的方法 | |
CN107605469B (zh) | 预测地层孔隙压力的方法 | |
CN103790579B (zh) | 随钻地质导向中确定钻头与地层界面距离的方法及装置 | |
CN105626056A (zh) | 一种确定目标区三维地层孔隙压力的方法和装置 | |
US10401808B2 (en) | Methods and computing systems for processing and transforming collected data to improve drilling productivity | |
WO2014077833A1 (en) | Triple phase evaluation of formation fluids | |
Almalikee et al. | Overburden stress and pore pressure prediction for the North Rumaila oilfield, Iraq | |
Feng et al. | Evaluation of reservoir permeability using array induction logging | |
US10802177B2 (en) | Evaluating hydrocarbon reserves using tool response models | |
Alvarez et al. | Evaluation of a Fractured Tight Reservoir in Real-Time: The importance of Detecting Open Fractures While Drilling with Accurate Mud Flow Measurement | |
US20230349286A1 (en) | Geologic formation characterization | |
Vandenbohede et al. | Double forced gradient tracer test: Performance and interpretation of a field test using a new solute transport model | |
WO2018132103A1 (en) | Detecting a flood front in a formation | |
Noah | Pore pressure evaluation from well logging and drilling exponent at Amal field, Gulf of Suez area, Egypt |
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 |