CN112241025B - 一种井震联合地层压力确定方法及系统 - Google Patents

一种井震联合地层压力确定方法及系统 Download PDF

Info

Publication number
CN112241025B
CN112241025B CN201910648700.4A CN201910648700A CN112241025B CN 112241025 B CN112241025 B CN 112241025B CN 201910648700 A CN201910648700 A CN 201910648700A CN 112241025 B CN112241025 B CN 112241025B
Authority
CN
China
Prior art keywords
seismic
data
obtaining
velocity
longitudinal wave
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
Application number
CN201910648700.4A
Other languages
English (en)
Other versions
CN112241025A (zh
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 CN201910648700.4A priority Critical patent/CN112241025B/zh
Publication of CN112241025A publication Critical patent/CN112241025A/zh
Application granted granted Critical
Publication of CN112241025B publication Critical patent/CN112241025B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters
    • 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
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种井震联合地层压力确定方法及系统，所述方法包括：通过地震资料和地质数据得到地震速度分布；根据所述地震速度分布得到实测纵波速度；通过拟合趋势线得到正常压实下的纵波速度，并根据所述实测纵波速度得到地层压力，本发明引入对地层压力具有敏感性的地震多属性参数，以提高地层压力预测精度，降低压力预测多解性。

Description

一种井震联合地层压力确定方法及系统
技术领域
本发明涉及石油勘探技术领域，尤其涉及一种井震联合地层压力确定方法及系统。
背景技术
一些最可靠的定量超压检测和评价技术，例如地层压力确定，是以地球物理电缆测井为基础的，但这些技术多数是已经钻开井孔获得了测井和实际测试资料，是一种“事后”技术，而在实际的勘探开发过程中，往往需要在钻开井孔之前了解地下地层压力状况，以便降低钻井工程风险。在海域勘探早期，由于勘探技术、高昂成本等因素的限制，区块内通常无钻井或仅有少量钻井，可供研究的井筒资料较匮乏，前期研究主要以地震资料为主，钻前压力的预测和计算主要是利用地震相关资料。
当前钻前压力预测通常采用的方法有经验关系法、等效深度法和有效应力法等，主要是利用到钻井声波、地震层速度等相关资料，其理论基础是，速度通常会随深度增加而增加，但是当出现超压带时,常常伴随有速度反转特征，表现为低速异常响应。Eaton(1972，1976)依据有效应力理论提出了Eaton公式法，该方法通过建立正常压实趋势线，利用实测数据与正常趋势线间的偏离关系计算地层压力，取得了一定的成效。然而地震速度的变化是受地下多种地质因素的影响，如果把低速异常作为超压层的唯一地震响应，会导致地震地层压力结果具有多解性和近似性。且随着钻井技术的不断发展,其局限性也逐渐暴露出来：随着压实程度的提高，声波速度增加的程度将逐渐变慢，达到一定的深度后声波速度将几乎不再随深度变化，而不是沿传统的正常趋势线无限延伸下去，同时对正常压实趋势线的求取，往往带有主观性，实践经验的多少和正确与否对检测地层压力结果影响很大。同时，Hunt(1990)通过大量实践数据研究，认为正常压实趋势线的拟合往往不是线性或半对数线性，有时需要建立多条趋势线，有时在某些层段甚至还无法建立趋势线。
发明内容
本发明的一个目的在于提供一种井震联合地层压力确定方法，引入对地层压力具有敏感性的地震多属性参数，以提高地层压力预测精度，降低压力预测多解性。本发明的另一个目的在于提供一种井震联合地层压力确定系统。本发明的再一个目的在于提供一种计算机设备。本发明的还一个目的在于提供一种可读介质。
为了达到以上目的，本发明一方面公开了一种井震联合地层压力确定方法，包括：
通过地震资料和地质数据得到地震速度分布；
根据所述地震速度分布得到实测纵波速度；
通过拟合趋势线得到正常压实下的纵波速度，并根据所述实测纵波速度得到地层压力。
优选地，所述通过地震资料和地质数据得到地震速度分布具体包括：
根据地震资料和地质数据建立目标地层的速度场模型，得到层速度数据体；
根据所述层速度数据体提取伪井数据；
通过对目标地层的已钻探井或所述伪井数据进行合成记录标定得到反演的子波；
根据已钻探井或所述伪井数据建立初始波阻抗模型；
根据所述子波和所述初始波阻抗模型反演得到所述地震速度分布。
优选地，所述通过拟合趋势线得到正常压实下的纵波速度，并根据所述实测纵波速度得到地层压力具体包括：
通过幂函数或二次多项式拟合趋势线；
根据所述趋势线得到对应实测纵波速度反转段的正常压实下的纵波速度；
通过所述实测纵波速度和所述正常压实下的纵波速度的偏离关系得到地层压力。
优选地，所述地震资料包括砂岩测井资料和泥岩油井资料。
优选地，所述方法进一步包括在通过地震资料和地质数据得到地震速度分布之前：
通过DST数据对所述地质数据进行修正。
优选地，所述方法进一步包括：
根据地层压力与所述地震频率的对应关系，通过目标地层对应的地震瞬时频率去除所述地层压力的异常值。
本发明还公开了一种井震联合地层压力确定系统，包括：
速度分布确定单元，用于通过地震资料和地质数据得到地震速度分布；
实测速度确定单元，用于根据所述地震速度分布得到实测纵波速度；
地层压力确定单元，手于通过拟合趋势线得到正常压实下的纵波速度，并根据所述实测纵波速度得到地层压力。
优选地，所述速度分布确定单元用于根据地震资料和地质数据建立目标地层的速度场模型，得到层速度数据体，根据所述层速度数据体提取伪井数据，通过对目标地层的已钻探井或所述伪井数据进行合成记录标定得到反演的子波，根据已钻探井或所述伪井数据建立初始波阻抗模型，根据所述子波和所述初始波阻抗模型反演得到所述地震速度分布。
优选地，所述地层压力确定单元用于通过幂函数或二次多项式拟合趋势线，根据所述趋势线得到对应实测纵波速度反转段的正常压实下的纵波速度，通过所述实测纵波速度和所述正常压实下的纵波速度的偏离关系得到地层压力。
优选地，所述地震资料包括砂岩测井资料和泥岩测井资料。
优选地，所述系统进一步包括参数修正单元，用于通过DST数据对所述地质数据进行修正。
优选地，所述系统进一步包括异常排除单元，用于根据地层压力与所述地震频率的对应关系，通过目标地层对应的地震瞬时频率去除所述地层压力的异常值。
本发明还公开了一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上所述方法。
本发明还公开了一种计算机可读介质，其上存储有计算机程序，该程序被处理器执行时实现如上所述方法。
本发明通过震资料和地质数据得到地震速度分布，得到实测纵波速度，并进一步通过拟合趋势线得到正常压实下的纵波速度，并根据所述实测纵波速度得到地层压力。本发明通过改进Eaton公式等手段联合，提高地震资料计算地层压力的精度和可靠性，以期用最大限度排除异常压力预测过程中的多解性，更好的服务于油田勘探生产。进一步还综合分岩性压力预测和地震多属性约束等方法，以进一步提高地层压力确定的准确度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出本发明一种井震联合地层压力确定方法一个具体实施例的流程图之一；
图2示出地震均方根速度转换的层速度剖面的示意图；
图3示出本发明一种井震联合地层压力确定方法一个具体实施例波阻抗反演得到的层速度剖面图；
图4示出本发明一种井震联合地层压力确定方法一个具体实施例的流程图之二；
图5示出本发明一种井震联合地层压力确定方法一个具体实施例井筒声波速度与井旁地震反演层速度对比图；
图6示出本发明一种井震联合地层压力确定方法一个具体实施例的流程图之三；
图7示出本发明一种井震联合地层压力确定方法一个具体实施例井筒声波测井数据及正常趋势线拟合图；
图8示出本发明一种井震联合地层压力确定方法一个具体实施例井点实测泥浆比重、孔隙压力与破裂压力的示意图；
图9示出本发明一种井震联合地层压力确定方法一个具体实施例的流程图之四；
图10示出本发明一种井震联合地层压力确定方法一个具体实施例测井速度(分砂、泥岩)、地震层速度分别计算地层压力与DST实测压力系数对比图；
图11示出本发明一种井震联合地层压力确定方法一个具体实施例的流程图之五；
图12示出本发明一种井震联合地层压力确定方法一个具体实施例计算的压力系数剖面图；
图13a示出本发明一种井震联合地层压力确定方法一个具体实施例瞬时频率与有效压力关系图一；
图13b示出本发明一种井震联合地层压力确定方法一个具体实施例瞬时频率与有效压力关系图二；
图14a示出本发明一种井震联合地层压力确定方法一个具体实施例随有效压力增加瞬时频率和纵波速度增长量、变化率关系图之一；
图14b示出本发明一种井震联合地层压力确定方法一个具体实施例随有效压力增加瞬时频率和纵波速度增长量、变化率关系图之二；
图15示出本发明一种井震联合地层压力确定方法一个具体实施例井点测线提取瞬时频率剖面图；
图16示出本发明一种井震联合地层压力确定方法一个具体实施例井点测线提取瞬时频率剖面图的局部放大图；
图17示出本发明一种井震联合地层压力确定系统一个具体实施例的结构图之一；
图18示出本发明一种井震联合地层压力确定系统一个具体实施例的结构图之二；
图19示出本发明一种井震联合地层压力确定系统一个具体实施例的结构图之三；
图20示出适于用来实现本发明实施例的计算机设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
根据本发明的一个方面，本实施例公开了一种井震联合地层压力确定方法。如图1所示，所述方法包括：
S100：通过地震资料和地质数据得到地震速度分布；
S200：根据所述地震速度分布得到实测纵波速度；
S300：通过拟合趋势线得到正常压实下的纵波速度，并根据所述实测纵波速度得到地层压力。
本发明通过震资料和地质数据得到地震速度分布，得到实测纵波速度，并进一步通过拟合趋势线得到正常压实下的纵波速度，并根据所述实测纵波速度得到地层压力。本发明通过改进Eaton公式等手段联合，提高地震资料计算地层压力的精度和可靠性，以期用最大限度排除异常压力预测过程中的多解性，更好的服务于油田勘探生产。
在可选的实施方式中，地震速度一般可以从VSP(垂直地震剖面)资料、叠加速度谱或地震反演层速度获取。其中，VSP资料费用高、相对稀少，一般很难获取，而且是在钻开井孔后获得的。地震叠加速度谱来源于地震资料处理，主要是为了满足叠加偏移需要，会将一些比较关键的变化切除或平滑掉，且速度拾取横纵向间隔范围大、采样稀，速度在横向上的变化也会受到数据的平均效应影响而减弱，而且在构造相对复杂、纵横向速度变化较大地区，在转为地层层速度时会存在较大误差，从而影响到压力预测的垂向精度，仅用地震叠加速度预测地层压力会存在诸多缺陷。地震反演层速度方法，可以把高分辨率井资料、地震速度谱的速度变化趋势、地震数据的横纵向密度高等信息结合起来，利用少量井和地震速度谱约束建立地质模型，通过地震反演获得相对较高精度的地震层速度。例如，图2示出了地震均方根速度转换的层速度剖面。井点处中上中新统～下中新统存在低速异常，但由于地震速度谱纵向和横向采样点稀疏，直接用于预测压力精度低，难以满足压力预测需求。图3示出了波阻抗反演得到的层速度剖面。由于结合了高分辨率井资料、地震速度谱的速度变化趋势、地震数据的横纵向密度高等信息，对应的层速度剖面分辨率有明显的提高，剖面中层速度的成层性更明显，井点处的中上中新统和下中新统地层中可见两个清晰的低速带。
在优选的实施方式中，可采用地震反演层速度方法得到地震速度分布。如图4的所示，所述S100具体可包括：
S110：根据地震资料和地质数据建立目标地层的速度场模型，得到层速度数据体。具体的，通过速度直方图分析，对三维速度谱资料中的奇异值进行剔除，通过Dix公式转换为层速度，进行三维空间插值，建立本区的速度场模型，得到层速度体，如图3所示。
S120：根据所述层速度数据体提取伪井数据。具体的，依据计算的层速度体来提取伪井，伪井井点坐标选择的原则为：井点所处的线道为过原始的采样点，尽量避免井位处提取的层速度数据不是通过原始速度谱点插值出来的数据值，提高数据的准确性，提取多口伪井时，井点分布在研究区的不同构造区带上，能反映不同构造带上的曲线特征作为后续参与反演的井曲线。
S130：通过对目标地层的已钻探井或所述伪井数据进行合成记录标定得到反演的子波。具体的，对已钻探井(通常较少或没有)或伪井进行合成记录标定，提取子波，同时从多口井提取的子波中，挑选波形较好相对较平稳的子波取平均值作为反演的子波。
S140：根据已钻探井或所述伪井数据建立初始波阻抗模型。具体的，对实钻井和伪井数据沿层内插外推，建立地质框架模型，产生一个平滑的初始波阻抗模型，初始波阻抗模型是补偿低频成分。
S150：根据所述子波和所述初始波阻抗模型反演得到所述地震速度分布。其中，优选的，可采用稀疏脉冲反演得到地震速度分布。约束稀疏反演的关键参数λ反映波阻抗值和子波褶积产生的合成地震道与实际地震道匹配程度的好坏，λ值太大，过分强调地震残差最小，一味地使合成记录与原始地震道吻合，会使一些噪音也加到了反演剖面中，同时也忽略了反射系数的稀疏，即忽略了波阻抗变化的低频成分，通过反演参数尝试确定λ的取值，使反演结果中能补充井上和叠加速度的低频背景及地震有限带宽部分，有效提高反演层速度的分辨率和精度。例如图5井筒声波速度与井旁地震反演层速度对比图。由浅至深地震层速度总体变化趋势与测井速度是基本吻合，在约800米处二者均出现速度反转。
在优选的实施方式中，如图6的所示，所述S300具体可包括：
S310：通过幂函数或二次多项式拟合趋势线。
对于地层正常压实时速度Vn的求取，传统方式是以正常压力段的声波测井数据和对应深度以线性或半对数线性进行拟合，通过拟合趋势线向下外推，从而得到对应实测纵波速度反转段的正常压实下的纵波速度，通过实测数据与正常趋势线间的偏离关系求得地层压力。
随着深度的增加，声波速度增加的程度是逐渐变慢的，达到一定的深度后声波速度将几乎不再随深度变化，以线性或半对数线性求取Vn值，是有理论缺陷的，当地层埋藏较深时，该方法计算的地层压力值会明显偏大于地层实际压力，误差较大。本发明通过研究后发明正常趋势线更偏向于为幂函数(Vn＝A×ZB)或二次多项式关系(Vn＝A×Z2+B×Z+C)。本实施例中，以幂函数方式拟合出正常趋势线，计算对应的Vn值，其中，A、B、C为拟合系数，常量。例如图7示出井筒声波测井数据及正常趋势线拟合图。井数据为实钻测井声波数据，正常趋势线为以上部正常压实段声波数据按幂函数关系拟合。
S320：根据所述趋势线得到对应实测纵波速度反转段的正常压实下的纵波速度。即通过拟合趋势线向下外推，从而得到对应实测纵波速度反转段的正常压实下的纵波速度，通过实测数据与正常趋势线间的偏离关系求得地层压力。
S330：通过所述实测纵波速度和所述正常压实下的纵波速度的偏离关系得到地层压力。
具体的，Eaton(1976)提出利用声波速度计算地层压力公式：
Vn＝Vo+kZ或log(Vn)＝Vo+k×Z
其中：Pp——预测的孔隙压力，MPa；Pr——静岩压力，MPa；Ph——正常的静水压力，MPa；Vn——地层正常压实时的速度，m/s；V——实测地层速度，m/s；N——实验系数，通常为1-3；Vo——海底对应速度，通常为1524m/s；k——拟合系数，常量；Z——从海底开始的测量深度，m。从该方程可以看出，Eaton认为有效应力与纵波速度的多次方成正比，通过上述公式可得到实时的地层压力。
在优选的实施方式中，所述地震资料包括砂岩测井资料和泥岩油井资料。为了消除岩性对速度的影响，对井上声波测井速度按砂岩和泥岩进行过滤分选，然后分别运用砂、泥岩测井速度采用Eaton改进方法计算地层压力，计算结果与DST测试结果对比，测井速度分岩性计算的地层压力与实测结果吻合。例如图8示出井点实测泥浆比重、孔隙压力与破裂压力。钻井深为3000米，900米以上泥浆比重基本维持在1.1s.g.左右；自900米开始至约1050米，泥浆开始发生气侵，泥浆比重增大到1.22s.g.；在约1050米处，泥浆比重增大到1.30s.g.，后泥浆比重增大到1.6s.g.，最后泥浆比重被调配到1.67s.g.，后经测试此处真实的地层压力梯度为1.80MPa/100m。到达完钻深度时，泥浆比重增至2.03s.g.。该井约900米以上为正常压力地层，900米以下地层发育异常高压。
在优选的实施方式中，如图9的所示，所述方法进一步包括在通过地震资料和地质数据得到地震速度分布之前：
S000：通过DST(有线随钻测斜仪)数据对所述地质数据进行修正。以实测DST数据进行刻度，对实验参数值进行微调，以此参数用于井旁地震反演的层速度资料中进行压力计算，井震数据所得结果较吻合，但是基于声波测井的预测结果更精细。例如图10示出测井速度(分砂、泥岩)、地震层速度分别计算地层压力与DST实测压力系数对比。以测井速度分岩性计算的地层压力与实测结果吻合，地震速度与测井泥岩速度对应的计算结果吻合较好。分析认为砂岩中的异常高压一般都是源自附近的高压泥岩，高压泥岩中的异常高压会通过各种形式向附近地层释放，导致砂岩压力高于周围泥岩压力。
在优选的实施方式中，如图11的所示，所述方法进一步包括：
S400：根据地层压力与所述地震频率的对应关系，通过目标地层对应的地震瞬时频率去除所述地层压力的异常值。基于地层压力变化与岩石物理参数的实验分析，认为地层超压与地震频率存在较好的对应关系，总体上为超压对应地震低频响应，提取该井附近地震瞬时频率属性作约束参数进行辅助表征，排除其多解性并验证压力计算结果的可信度。例如图12中，通过Eaton改进法计算压力系数剖面。井点附近从下上新统开始出现地层超压，到中上中新统，压力系数主要在1.4至1.6之间，但其下部有一个很明显的强超压仓，仓内压力系数超过1.6；下中新统底部也有一个近似的强超压仓。而在井点左侧约3000公里以外，在中上上新统和下上新统中发育强超压，至下部地层反而压力系数减小；井点右侧约1500公里以外，在下上新统和中上中新统中发育几个独立的强超压仓，可能计算结果异常。图13a和图13b示出瞬时频率与有效压力关系图。从图13a和图13b可以看出，瞬时频率随着有效压力的降低而减小。图14a和图14b示出随有效压力增加瞬时频率和纵波速度增长量、变化率关系图。随着有效压力的增加，瞬时频率增长量与纵波速度增长量呈线性关系，瞬时频率的变化率相对大于纵波速度的变化率。
图15和图16示出了井点测线提取瞬时频率剖面图。井点两侧的中上上新统和下上新统地层中基本上为中高频响应，未见明显大范围的低频特征，但在井点到断层发育范围附近，这两套地层处的瞬时频率属性相对两侧是明显变低的，频率的降低可能是由于断层附近地层较破碎及深部异常压力以断层为通道泄压所致。井点瞬时频率剖面的局部放大图上，可见分岩性计算地层压力的结果与地震频率匹配较好，较高压力的砂岩对应的地震瞬时频率较低，较低压力的泥岩对应的地震瞬时频率较高。
本实施例中通过地震反演层速度方法得到目标地层的速度分布，能有效提高地震速度分辨率和精度。并对Eaton经验公式法的改进，完善了参数求取方法，可应用于地震层速度中进行压力计算，井震计算结果匹配较好。进一步地，本实施例还通过地震多属性约束表征，能有效排除其多解性并验证压力计算结果的可信度。
地层压力的预测在油气勘探生产中意义重大，地层压力的研究可适用到石油勘探至开发的整个过程：在地质评价阶段，地层压力可作为研究油气成藏流体动力学的基础、分析圈闭和盆地几何形状的依据、评价储集层资源潜力的重要参数；在钻井工程阶段，地层压力的预测可准确确定套管程序、合理选择泥浆密度、有效降低储层伤害、避免钻井重大事故；在开发生产阶段，地层压力可以确定油藏压力、保障科学开采、指导合理注水等。
现有技术特征是直接运用前人已发明的经验公式法，在相关应用上存在理论缺陷。应用效果上存在多解性或相似性，地层较浅时预测精度较高，地层较深时会存在较大误差。本发明运用地震反演层速度、压力计算数学模型改进、分岩性压力预测、地震多属性约束的组合方法。在相关的发明技术中，通过地震反演层速度方法有效提高了地震速度分辨率和精度；通过对Eaton经验公式法的改进，完善了参数求取方法，应用于地震层速度中进行压力计算，井震计算结果匹配好；通过地震多属性约束表征，有效排除其多解性并验证压力计算结果的可信度。
基于相同原理，本实施例还公开了一种井震联合地层压力确定系统。如图17所示，所述系统包括速度分布确定单元11、实测速度确定单元12和地层压力确定单元13。
速度分布确定单元11用于通过地震资料和地质数据得到地震速度分布。
实测速度确定单元12用于根据所述地震速度分布得到实测纵波速度。
地层压力确定单元13用于通过拟合趋势线得到正常压实下的纵波速度，并根据所述实测纵波速度得到地层压力。
在优选的实施方式中，所述速度分布确定单元11用于根据地震资料和地质数据建立目标地层的速度场模型，得到层速度数据体，根据所述层速度数据体提取伪井数据，通过对目标地层的已钻探井或所述伪井数据进行合成记录标定得到反演的子波，根据已钻探井或所述伪井数据建立初始波阻抗模型，根据所述子波和所述初始波阻抗模型反演得到所述地震速度分布。
在优选的实施方式中，所述地层压力确定单元13用于通过幂函数或二次多项式拟合趋势线，根据所述趋势线得到对应实测纵波速度反转段的正常压实下的纵波速度，通过所述实测纵波速度和所述正常压实下的纵波速度的偏离关系得到地层压力。
在优选的实施方式中，所述地震资料包括砂岩测井资料和泥岩油井资料。
在优选的实施方式中，如图18所示，所述系统进一步包括参数修正单元14，用于通过DST数据对所述地质数据进行修正。
在优选的实施方式中，如图19所示，所述系统进一步包括异常排除单元15，用于根据地层压力与所述地震频率的对应关系，通过目标地层对应的地震瞬时频率去除所述地层压力的异常值。
本实施例中系统的实现原理与上述方法相同，系统的实施可参照地层压力确定方法的实施，在此不再赘述。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机设备，具体的，计算机设备例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
在一个典型的实例中计算机设备具体包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上所述的方法。
下面参考图20，其示出了适于用来实现本申请实施例的计算机设备600的结构示意图。
如图20所示，计算机设备600包括中央处理单元(CPU)601，其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM))603中的程序而执行各种适当的工作和处理。在RAM603中，还存储有系统600操作所需的各种程序和数据。CPU601、ROM602、以及RAM603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。
以下部件连接至I/O接口605：包括键盘、鼠标等的输入部分606；包括诸如阴极射线管(CRT)、液晶反馈器(LCD)等以及扬声器等的输出部分607；包括硬盘等的存储部分608；以及包括诸如LAN卡，调制解调器等的网格接口卡的通信部分609。通信部分609经由诸如因特网的网格执行通信处理。驱动器610也根据需要连接至I/O接口606。可拆卸介质611，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器610上，以便于从其上读出的计算机程序根据需要被安装如存储部分608。
特别地，根据本发明的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本发明的实施例包括一种计算机程序产品，其包括有形地包含在机器可读介质上的计算机程序，所述计算机程序包括用于执行流程图所示的方法的程序代码。在这样的实施例中，该计算机程序可以通过通信部分609从网格上被下载和安装，和/或从可拆卸介质611被安装。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网格而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (10)

1.一种井震联合地层压力确定方法，其特征在于，包括：
通过地震资料和地质数据得到地震速度分布；
根据所述地震速度分布得到实测纵波速度；
通过拟合趋势线得到正常压实下的纵波速度，并根据所述实测纵波速度得到地层压力；
所述通过地震资料和地质数据得到地震速度分布具体包括：
根据地震资料和地质数据建立目标地层的速度场模型，得到层速度数据体；
根据所述层速度数据体提取伪井数据；
通过对目标地层的已钻探井或所述伪井数据进行合成记录标定得到反演的子波；
根据已钻探井或所述伪井数据建立初始波阻抗模型；
根据所述子波和所述初始波阻抗模型反演得到所述地震速度分布；
所述通过拟合趋势线得到正常压实下的纵波速度，并根据所述实测纵波速度得到地层压力具体包括：
通过幂函数或二次多项式拟合趋势线；
根据所述趋势线得到对应实测纵波速度反转段的正常压实下的纵波速度；
通过所述实测纵波速度和所述正常压实下的纵波速度的偏离关系得到地层压力。
2.根据权利要求1所述的井震联合地层压力确定方法，其特征在于，所述地震资料包括砂岩测井资料和泥岩测井资料。
3.根据权利要求1所述的井震联合地层压力确定方法，其特征在于，所述方法进一步包括在通过地震资料和地质数据得到地震速度分布之前：
通过有线随钻测斜仪数据对所述地质数据进行修正。
4.根据权利要求1所述的井震联合地层压力确定方法，其特征在于，所述方法进一步包括：
根据地层压力与地震频率的对应关系，通过目标地层对应的地震瞬时频率去除所述地层压力的异常值。
5.一种井震联合地层压力确定系统，其特征在于，包括：
速度分布确定单元，用于通过地震资料和地质数据得到地震速度分布；
实测速度确定单元，用于根据所述地震速度分布得到实测纵波速度；
地层压力确定单元，用于通过拟合趋势线得到正常压实下的纵波速度，并根据所述实测纵波速度得到地层压力；
所述速度分布确定单元用于根据地震资料和地质数据建立目标地层的速度场模型，得到层速度数据体，根据所述层速度数据体提取伪井数据，通过对目标地层的已钻探井或所述伪井数据进行合成记录标定得到反演的子波，根据已钻探井或所述伪井数据建立初始波阻抗模型，根据所述子波和所述初始波阻抗模型反演得到所述地震速度分布；
所述地层压力确定单元用于通过幂函数或二次多项式拟合趋势线，根据所述趋势线得到对应实测纵波速度反转段的正常压实下的纵波速度，通过所述实测纵波速度和所述正常压实下的纵波速度的偏离关系得到地层压力。
6.根据权利要求5所述的井震联合地层压力确定系统，其特征在于，所述地震资料包括砂岩测井资料和泥岩测井资料。
7.根据权利要求5所述的井震联合地层压力确定系统，其特征在于，所述系统进一步包括参数修正单元，用于通过有线随钻测斜仪数据对所述地质数据进行修正。
8.根据权利要求5所述的井震联合地层压力确定系统，其特征在于，所述系统进一步包括异常排除单元，用于根据地层压力与地震频率的对应关系，通过目标地层对应的地震瞬时频率去除所述地层压力的异常值。
9.一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，
所述处理器执行所述程序时实现如权利要求1-4任一项所述方法。
10.一种计算机可读介质，其上存储有计算机程序，其特征在于，
该程序被处理器执行时实现如权利要求1-4任一项所述方法。
CN201910648700.4A 2019-07-18 2019-07-18 一种井震联合地层压力确定方法及系统 Active CN112241025B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910648700.4A CN112241025B (zh) 2019-07-18 2019-07-18 一种井震联合地层压力确定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910648700.4A CN112241025B (zh) 2019-07-18 2019-07-18 一种井震联合地层压力确定方法及系统

Publications (2)

Publication Number Publication Date
CN112241025A CN112241025A (zh) 2021-01-19
CN112241025B true CN112241025B (zh) 2023-11-28

Family

ID=74167446

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910648700.4A Active CN112241025B (zh) 2019-07-18 2019-07-18 一种井震联合地层压力确定方法及系统

Country Status (1)

Country Link
CN (1) CN112241025B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113514890B (zh) * 2021-03-18 2022-12-09 中国石油大学(华东) 一种使用地震数据预测地层压力系数的方法、装置及设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106094027A (zh) * 2016-06-01 2016-11-09 中国海洋石油总公司 一种垂直地震剖面vsp钻前压力预测方法和系统
CN107843927A (zh) * 2016-09-20 2018-03-27 中国石油化工股份有限公司 基于井震联合速度的页岩地层压力预测方法及装置
CN109143373A (zh) * 2017-06-28 2019-01-04 中国石油化工股份有限公司 页岩气储层地层压力计算方法及计算机可读存储介质
CN109667573A (zh) * 2018-12-12 2019-04-23 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 三维页岩储层孔隙压力预测方法、装置和电子设备

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106094027A (zh) * 2016-06-01 2016-11-09 中国海洋石油总公司 一种垂直地震剖面vsp钻前压力预测方法和系统
CN107843927A (zh) * 2016-09-20 2018-03-27 中国石油化工股份有限公司 基于井震联合速度的页岩地层压力预测方法及装置
CN109143373A (zh) * 2017-06-28 2019-01-04 中国石油化工股份有限公司 页岩气储层地层压力计算方法及计算机可读存储介质
CN109667573A (zh) * 2018-12-12 2019-04-23 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 三维页岩储层孔隙压力预测方法、装置和电子设备

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
井约束地震反演预测地层压力的方法――以渤海湾盆地某凹陷为例;管红;《天然气地球科学》;19(2);第276-279页 *
地层压力预测方法及其在莺歌海盆地的应用;刘天佑 等;《中国海上油气(地质)》;9(5);第333-339页 *
声波测井资料地质应用的某些实例及探索;楚泽涵 等;《地球物理学报》;第32卷(第6期);第700-710页 *
用地震层速度预测地层压力的方法;郭燕珩 等;《内蒙古石油化工》(5);第66-67页 *

Also Published As

Publication number Publication date
CN112241025A (zh) 2021-01-19

Similar Documents

Publication Publication Date Title
US11591900B2 (en) Method to predict overpressure uncertainty from normal compaction trendline uncertainty
US11098565B2 (en) Method for estimating permeability of fractured rock formations from induced slow fluid pressure waves
AU2013384285B2 (en) System, method and computer program product for predicting well production
US8788208B2 (en) Method to estimate pore pressure uncertainty from trendline variations
US20130176822A1 (en) Method for estimating elastic parameters by inverting 4d seismic measurements
CN109138975B (zh) 一种基于时移测井数据的求解相渗特征曲线的新方法
Johnson et al. Statistical comparison of methods for estimating sediment thickness from horizontal-to-vertical spectral ratio (HVSR) seismic methods: An example from Tylerville, Connecticut, USA
US20180164463A1 (en) Verifying Measurements of Elastic Anisotropy Parameters In An Anisotropic Wellbore Environment
CN110095811B (zh) 膏岩层速度模型构建处理方法及装置
CN110826157B (zh) 一种多压力成因超压贡献比例的计算方法
CN104090303A (zh) 一种地震反演方法和装置
CN112241025B (zh) 一种井震联合地层压力确定方法及系统
CN113279746A (zh) 一种套管变形风险区域的确定方法及应用
Anbazhagan et al. Establishing empirical correlation between sediment thickness and resonant frequency using HVSR for the Indo-Gangetic Plain
Gherabati et al. New engineering and geologic parameters to predict infill well performance in the Wolfcamp of the Delaware Basin
CN112147700A (zh) 速度异常区的低频模型构建方法及系统
El Faidouzi Contribution of Tidal Analysis to Reservoir Monitoring: A Field Case Study in a Fractured Reservoir Offshore Abu Dhabi, UAE
US8141259B2 (en) Method of determining the dip of a formation
CN110133721B (zh) 一种水力压裂过程监测方法及系统
US11035966B2 (en) Method for estimating reservoir formation quality using seismic velocities
CN111502647A (zh) 一种钻井地质环境因素确定方法、装置及存储介质
CN107589457A (zh) 测井资料尺度粗化的方法和装置
US11702924B2 (en) Determination of stimulated reservoir volume and estimated ultimate recovery of hydrocarbons for unconventional reservoirs
CN114563816B (zh) 油气藏评价阶段建立地震解释速度模型的方法及装置
Chen et al. Application of prediction techniques in carbonate karst reservoir in tarim basin

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