CN106646603A - 一种实际时移地震资料处理差异的可靠性判断方法 - Google Patents
一种实际时移地震资料处理差异的可靠性判断方法 Download PDFInfo
- Publication number
- CN106646603A CN106646603A CN201710004041.1A CN201710004041A CN106646603A CN 106646603 A CN106646603 A CN 106646603A CN 201710004041 A CN201710004041 A CN 201710004041A CN 106646603 A CN106646603 A CN 106646603A
- Authority
- CN
- China
- Prior art keywords
- velocity
- seismic
- data
- difference
- actual
- 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
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000012545 processing Methods 0.000 title claims abstract description 34
- 238000004088 simulation Methods 0.000 claims description 44
- 239000004576 sand Substances 0.000 claims description 28
- 230000004044 response Effects 0.000 claims description 23
- 230000008859 change Effects 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 15
- 239000004575 stone Substances 0.000 claims description 12
- 230000015572 biosynthetic process Effects 0.000 claims description 11
- 239000012530 fluid Substances 0.000 claims description 9
- 238000005516 engineering process Methods 0.000 claims description 8
- 108010046685 Rho Factor Proteins 0.000 claims description 4
- HJUFTIJOISQSKQ-UHFFFAOYSA-N fenoxycarb Chemical compound C1=CC(OCCNC(=O)OCC)=CC=C1OC1=CC=CC=C1 HJUFTIJOISQSKQ-UHFFFAOYSA-N 0.000 claims description 4
- 238000002637 fluid replacement therapy Methods 0.000 claims description 4
- 239000011435 rock Substances 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 5
- 238000012937 correction Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/612—Previously recorded data, e.g. time-lapse or 4D
-
- 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/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- 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/622—Velocity, density or impedance
- G01V2210/6224—Density
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (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
本发明涉及一种实际时移地震资料处理差异的可靠性判断方法,包括以下步骤:1)获得开发前后油藏模型正演模拟的地震记录振幅能量差异;2)获得开发前后实际地震资料处理结果的振幅能量差异;3)计算实际地震处理结果与油藏模型的差异地震数据的相关性;4)判断实际地震资料处理的结果对时移地震差异的解释是否可靠。
Description
技术领域
本发明涉及一种实际时移地震资料处理差异的可靠性判断方法,属于石油天然气地震勘探开发领域。
背景技术
时移地震的差异主要是油藏饱和度和压力变化引起的地震响应特征的变化。饱和度变化的解释可用于确定流体边界等目的,压力变化解释可用来描述连通性、断层封闭等目的。通常,两种变化会在同一个储层内发生,因此时移地震解释需要联系生产数据和油藏模型地震正演模拟,来解决油藏饱和度和压力变化共同产生的储层变化的不确定性问题。
基于地下地质结构和油藏开发前后相应的地球物理物性参数建立油藏地质模型,通过正演模拟合成油藏开发前后的地震记录,由此可以得到由于油藏参数变化引起的地震响应的变化。但两次实际地震资料采集的不一致性、各种噪声的干扰、以及处理技术的影响等,使得实际地震资料的处理结果有时不能正确反映油藏开采引起的地震响应的变化。现有技术中尚不存在判断实际地震资料的处理是否能够正确反映油藏开采引起的地震响应的方法。
发明内容
针对上述问题,本发明的目的是提供一种能够判断实际地震资料的处理是否准确反映油藏参数变化的基于模型的实际时移地震资料处理差异的可靠性判断方法。
为实现上述目的,本发明采用以下技术方案:一种实际时移地震资料处理差异的可靠性判断方法,包括以下步骤:
1)获得开发前后油藏模型正演模拟的地震记录振幅能量差异;
2)获得开发前后实际地震资料处理结果的振幅能量差异;
3)计算实际地震资料处理结果与油藏模型正演模拟的地震记录的振幅能量差异的相关性;
4)判断实际地震资料处理的结果对时移地震差异的解释是否可靠。
所述步骤1)中,开发前后油藏模型正演模拟的地震记录振幅能量差异获取过程如下:①通过油藏数模结果得到开发前后两个时间点的孔隙度和饱和度模型;②将开发前后两个时间点的孔隙度和饱和度这两个三维油藏数模网格属性转换为三维油藏弹性属性,三维油藏弹性属性包括纵波速度、横波速度和密度;③利用正演模拟方法,获得开发前后三维油藏模型的地震响应记录;④利用开发前后三维油藏模型的模拟地震响应记录直接相减得到模拟数据的差异;⑤计算开发前后三维油藏模型的模拟数据差异的均方根振幅作为其振幅能量。
所述步骤①中,开发前后两个时间点的孔隙度和饱和度模型的获取方法如下:
利用ECLIPSE数值模拟软件,针对目标地质模型,经过数值计算,得以获取开发前后两个时间点的孔隙度和饱和度两个三维油藏网格属性参数。
所述步骤②中,通过开放前后两个时间点的孔隙度和饱和度这两个三维油藏数模网格属性获取纵波速度、横波速度和密度的三维油藏数模网格属性方法如下:
根据测井数据区分砂层和泥层,分别统计出砂层和泥层的泥质含量和孔隙度的关系;
对于砂层,根据其流体情况,应用流体替换技术,把砂岩的纵波速度和密度的数据转化为饱含水砂岩的纵波速度和密度的数据,统计出饱含水砂岩的纵波速度和密度与孔隙度的关系;对于泥层,直接统计出其纵波速度和密度与孔隙度之间的关系;
对于砂层,应用泥岩公式得到饱含水砂岩的横波速度与纵波速度的关系;对于泥层,应用泥岩公式得到岩石的横波速度和纵波速度的关系;
对于砂岩,将统计出的饱含水砂岩的纵波速度、横波速度、密度和泥质含量的公式输入到三维油藏数模网格属性数据中,应用流体替换计算得到实际地层的纵波速度、横波速度和密度;对于泥岩,统计出的纵波速度、横波速度和密度即为实际地层数据。
所述步骤③中,利用正演模拟方法获得开发前后三维油藏模型的地震响应记录方法如下:
根据式(1)计算开发前后三维油藏模型中各层不同位置的反射系数:
其中,ti代表第i网格点的时间,R(ti)是第i网格点的反射系数,ρ(ti+1)是第i网格点界面上部的密度,v(ti+1)是第i网格点界面上部的纵波速度,ρ(ti)是第i网格点界面下部的密度,v(ti)是第i网格点界面下部的纵波速度,t代表时间;
选取与实际地震资料主频一致的子波;
根据式(2),通过子波与反射系数褶积得到开发前后三维油藏模型的模拟地震响应记录:
其中,S(t)是模拟地震响应记录,b(t)代表地震子波。
在所述步骤2)中,获得开发前后实际地震资料处理结果的振幅能量差异的过程如下:
①采集开发前后实际地震资料;
②利用开发前后实际地震资料进行数据处理得到开发前后实际地震的地震响应记录;
③计算开发前后实际地震资料处理结果的差异;
④计算实际地震资料差异数据的均方根振幅作为其振幅能量差异。
所述步骤3)中,计算实际地震资料处理结果与油藏模型正演模拟的地震记录的振幅能量差异的相关性的方法如下:
比较开发前后油藏模型正演模拟与实际地震资料处理结果沿层振幅能量差异的一致性,如果两者振幅能量分布、强弱关系等一致性较高,利用下述相关性值公式(3)计算对应井点位置实际资料处理和油藏模型模拟的差异地震数据在目的层段的相关性,
其中,f1是井点位置实际资料处理的目的层段差异数据,f2为井点位置油藏模型模拟的目的层段差异数据,t1、t2为目的层段的起止时间。
所述步骤4)中,判断实际地震资料处理的结果对时移地震差异的解释是否可靠的方法如下:
如果步骤3)中计算结果相关性值若大于等于90%,则判定实际地震资料处理的结果能够反映油藏参数的变化;如果计算结果表明相关性低于90%,则判定实际地震资料处理的结果不能反映油藏参数的变化。
本发明由于采取以上技术方案,其具有以下优点:1、本发明可为后续时移地震差异数据解释和反演提供可靠的数据基础。2、本发明可提高油田剩余油预测精度。3、本发明根据实际地震资料的时移处理结果能够反映油藏参数的变化,满足后续时移地震解释的需求。4、本发明可用来解决一些由于油藏饱和度和压力变化共同产生的储层变化的不确定性问题。
附图说明
图1是本发明三维油藏模型,(A)图是2003年的孔隙度模型图,(B)图是2013年的孔隙度模型图;
图2是本发明三维油藏模型,(A)图是2003年的饱和度模型图,(B)图是2013年的饱和度模型图;
图3是本发明三维油藏模型中H层的属性模型,(A)图是2003年的孔隙度属性模型图,(B)是2013年的孔隙度属性模型图;
图4是本发明三维油藏模型中H层的属性模型,(A)图是2003年的饱和度属性模型图,B是2013年的饱和度属性模型图;
图5是本发明三维油藏模型,(A)图是2003年的纵波速度模型图,(B)图是2003年的纵波速度模型图;
图6是本发明三维油藏模型,(A)图是2003年横波速度模型图,(B)图是2013年横波速度模型图;
图7是本发明三维油藏模型,(A)图是2003年密度的模型图,(B)图是2013年密度的模型图;
图8是本发明三维油藏模型中H层的属性模型,(A)图是2003年的纵波速度模型图,(B)图是2003年的纵波速度模型图;
图9是本发明三维油藏模型中H层的属性模型,(A)图是2003年的横波速度模型图,(B)图是2003年的横波速度模型图;
图10是本发明三维油藏模型中H层的属性模型,(A)图是2003年的密度的模型图,(B)图是2013年的密度的模型图;
图11是本发明正演模拟得到的地震记录,(A)图是基于2003年三维油藏弹性模型得到的地震记录图,(B)图是基于2013年的三维油藏弹性模型得到的地震记录图,图(C)是图(A)的记录数值减去图(B)的记录数值的差值绘制成的记录图;
图12是本发明时移地震资料处理结果及差异,(A)图为2003年数据图;(B)图为2013年数据图;图(C)为图A的记录数值减去图B的记录数值的差值绘制成的记录图;
图13是本发明H层模拟数据振幅能量差异与实际地震资料振幅能量差异的对比,(A)图是模拟数据振幅能量差异图,(B)图是实际地震资料振幅能量差异图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明提出了一种实际时移地震资料处理差异的可靠性判断方法,包括以下步骤:
1)获得开发前后油藏模型正演模拟的地震记录振幅能量差异,具体过程如下:
①通过油藏数模结果得到开发前后两个时间点的孔隙度和饱和度模型。
利用ECLIPSE数值模拟软件,针对目标地质模型,经过数值计算,可以输出任意开发阶段的三维油藏网格属性参数(孔隙度、渗透率、饱和度、压力、油藏深度等),获取开发前后两个时间点的孔隙度和饱和度两个三维油藏网格属性参数。
②将开发前后两个时间点的孔隙度和饱和度这两个三维油藏数模网格属性转换为三维油藏弹性属性,三维油藏弹性属性包括纵波速度、横波速度和密度。
根据测井数据区分砂层和泥层,分别统计出砂层和泥层的泥质含量和孔隙度的关系。
对于砂层,根据其流体情况,应用流体替换技术,把砂岩的纵波速度和密度的数据转化为饱含水砂岩的纵波速度和密度的数据,统计出饱含水砂岩的纵波速度和密度与孔隙度的关系;对于泥层,直接统计出其纵波速度和密度与孔隙度之间的关系。
对于砂层,应用泥岩公式得到饱含水砂岩的横波速度与纵波速度的关系;对于泥层,应用泥岩公式得到岩石的横波速度和纵波速度的关系。
对于砂岩,将统计出的饱含水砂岩的纵波速度、横波速度、密度和泥质含量的公式输入到三维油藏数模网格属性数据中,应用流体替换计算得到实际地层的纵波速度、横波速度和密度;对于泥岩,统计出的纵波速度、横波速度和密度即为实际地层数据。
③利用正演模拟方法,获得开发前后三维油藏模型的地震响应记录;
根据式(1)计算开发前后三维油藏模型中各层不同位置的反射系数:
其中,ti代表第i网格点的时间,R(ti)是第i网格点的反射系数,ρ(ti+1)是第i网格点界面上部的密度,v(ti+1)是第i网格点界面上部的纵波速度,ρ(ti)是第i网格点界面下部的密度,v(ti)是第i网格点界面下部的纵波速度,t代表时间。
选取与实际地震资料主频一致的子波;
根据式(2),通过子波与反射系数褶积得到开发前后三维油藏模型的模拟地震响应记录:
其中,S(t)是模拟地震响应记录,b(t)代表地震子波。
④利用开发前后三维油藏模型的模拟地震响应记录直接相减得到模拟数据的差异;
⑤计算开发前后三维油藏模型的模拟数据差异的均方根振幅作为其振幅能量。
2)获得开发前后实际地震资料处理结果的振幅能量差异。
①采集开发前后实际地震资料;
②利用下述现有的时移地震一致性处理流程对采集到的资料进行数据处理得到开发前后实际地震资料的地震响应记录,时移地震一致性处理流程如下;
a、地震数据解编、导航数据合并;
b、多次波等噪音衰减;
c、潮汐校正;
d、船速校正;
e、振幅补偿;
f、子波整形;
g、预测反褶积;
h、一致性抽取;
i、面元均化、速度分析、叠加;
j、拉冬变换压制远道多次波;
k、针对相位的Q补偿;
l、叠前偏移、剩余多次波衰减;
m、全区叠加。
③利用下述现有的时移地震匹配处理流程得到开发前后实际地震资料处理结果的差异,时移地震匹配处理流程如下;
a、输入开发前和开发后的叠后数据;
b、整体大尺度能量匹配;
c、频率特征匹配;
d、开发后数据体基于标准层时间相位匹配;
e、开发后数据体基于标准层动态拉伸校正;
f、开发后数据体纵向能量小尺度匹配;
g、开发后数据体基于标志层的匹配滤波。
④计算实际地震资料差异数据的均方根振幅作为其振幅能量差异。
3)计算实际地震资料处理结果与油藏模型正演模拟的地震记录的振幅能量差异的相关性。
比较开发前后油藏模型正演模拟与实际地震资料处理结果沿层振幅能量差异的一致性,如果两者振幅能量分布、强弱关系等一致性较高,利用下述相关性值公式(3)计算对应井点位置实际资料处理和油藏模型模拟的差异地震数据在目的层段的相关性,
其中,f1是井点位置实际资料处理的目的层段差异数据,f2为井点位置油藏模型模拟的目的层段差异数据,t1、t2为目的层段的起止时间。
4)判断实际地震资料处理的结果对时移地震差异的解释是否可靠。
如果步骤3)中计算结果相关性值若大于等于90%,则说明实际地震资料处理的结果能够反映油藏参数的变化,如果计算结果表明相关性低于90%,则说明实际地震资料处理的结果不能反映油藏参数的变化,需要进一步改进。
下面通过一个具体的实施例,用以说明本发明的效果。根据某地区2003年和2013年的开发油藏模型模拟的差异和实际地震资料处理的差异两者比较对本发明所提方法进行验证。
1)获得开发前后油藏模型正演模拟的地震记录振幅能量差异。
①通过油藏数模结果得到开发前后两个时间点的孔隙度和饱和度模型;
利用ECLIPSE数值模拟软件,针对目标地质模型,经过数值计算,可以输出任意开发阶段的三维油藏网格属性参数,其中开发前后两个时间点的孔隙度和饱和度的图像模型如图1~图4所示。
②将开发前后两个时间点的孔隙度和饱和度这两个三维油藏数模网格属性转换为三维油藏弹性属性,三维油藏弹性属性包括纵波速度、横波速度和密度,2003年和2013年三维油藏模型中的纵波速度、横波速度、密度的模型图如图5~图7所示;2003年和2013年三维油藏模型中H层的纵波速度、横波速度、密度的模型图如图8~图10所示。
根据测井数据区分砂层和泥层,分别分析砂层和泥层的泥质含量和孔隙度的关系。
对于砂层,根据其流体情况,应用流体替换技术,把砂岩的纵波速度和数据转化为饱含水砂岩的纵波速度和密度,统计出其纵波速度和密度与孔隙度的关系;对于泥层,直接统计出其纵波速度和密度与孔隙度之间的关系。
对于砂层,应用泥岩公式得到饱含水砂岩的横波数据与纵波速度的关系;对于泥层,应用泥岩公式得到岩石的横波数据和纵波速度的关系。
对于砂岩,将统计出的饱含水砂岩的纵波速度、横波速度、密度和泥质含量的公式输入到三维油藏数模网格属性数据中,应用流体替换计算得到实际地层的纵波速度、横波速度和密度的数据;对于泥岩,统计出的纵波速度、横波速度和密度的数据即为实际地层数据。
重复上述步骤,可得到开发后的实际地层数据。
③利用正演模拟方法,获得开发前后三维油藏模型的地震响应记录,如图11所示;
根据式(1)计算开发前后三维油藏模型中各层不同位置的反射系数:
其中,ti代表第i网格点的时间,R(ti)是第i网格点的反射系数,ρ(ti+1)是第i网格点界面上部的密度,v(ti+1)是第i网格点界面上部的纵波速度,ρ(ti)是第i网格点界面下部的密度,v(ti)是第i网格点界面下部的纵波速度,t代表时间。
选取与实际地震资料主频一致的子波;
根据式(2),通过子波与反射系数褶积得到开发前后三维油藏模型的模拟地震响应记录:
其中,S(t)是模拟地震响应记录,b(t)代表地震子波。
④利用开发前后三维油藏模型的模拟地震响应记录直接相减得到模拟数据的差异,其差异如图11(C)所示;
⑤计算开发前后三维油藏模型的模拟数据差异的均方根振幅作为其振幅能量,其数据差异的振幅能量如图13(A)所示。
2)获得两次实际地震资料处理结果的振幅能量差异。
①采集两次实际地震资料;
②利用时移地震一致性处理流程对采集到的资料进行数据处理得到开发前后实际地震资料的地震响应记录,其记录如图12(A)、12(B)所示;
③利用时移地震匹配处理流程得到开发前后实际地震资料处理结果的差异如图12(C)所示;
④计算实际地震资料差异数据的均方根振幅作为其振幅能量差异,其差异如图13(B)所示。
3)计算实际地震资料处理结果与油藏模型正演模拟的地震记录的振幅能量差异的相关性。
比较开发前后油藏模型正演模拟与实际地震资料处理结果沿层振幅能量差异的一致性,如果两者振幅能量分布、强弱关系等一致性较高,利用下述相关性值公式(3)计算对应井点位置实际资料处理和油藏模型模拟的差异地震数据在目的层段的相关性,
其中,f1是井点位置实际资料处理的目的层段差异数据,f2为井点位置油藏模型模拟的目的层段差异数据,t1、t2为目的层段的起止时间。
4)判断实际地震资料处理的结果对时移地震差异的解释是否可靠。
根据前序步骤中提供的数据在步骤3)中进行计算得出其结果相关性值大于90%,则说明实际地震资料处理的结果能够反映油藏参数的变化。
上述各实施例仅用于对本发明的目的、技术方案和有益效果进行了进一步详细说明,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种实际时移地震资料处理差异的可靠性判断方法,包括以下步骤:
1)获得开发前后油藏模型正演模拟的地震记录振幅能量差异;
2)获得开发前后实际地震资料处理结果的振幅能量差异;
3)计算实际地震资料处理结果与油藏模型正演模拟的地震记录的振幅能量差异的相关性;
4)判断实际地震资料处理的结果对时移地震差异的解释是否可靠。
2.如权利要求1所述的一种实际时移地震资料处理差异的可靠性判断方法,其特征在于:所述步骤1)中,开发前后油藏模型正演模拟的地震记录振幅能量差异获取过程如下:①通过油藏数模结果得到开发前后两个时间点的孔隙度和饱和度模型;②将开发前后两个时间点的孔隙度和饱和度这两个三维油藏数模网格属性转换为三维油藏弹性属性,三维油藏弹性属性包括纵波速度、横波速度和密度;③利用正演模拟方法,获得开发前后三维油藏模型的地震响应记录;④利用开发前后三维油藏模型的模拟地震响应记录直接相减得到模拟数据的差异;⑤计算开发前后三维油藏模型的模拟数据差异的均方根振幅作为其振幅能量。
3.如权利要求2所述的一种实际时移地震资料处理差异的可靠性判断方法,其特征在于:所述步骤①中,开发前后两个时间点的孔隙度和饱和度模型的获取方法如下:
利用ECLIPSE数值模拟软件,针对目标地质模型,经过数值计算,得以获取开发前后两个时间点的孔隙度和饱和度两个三维油藏网格属性参数。
4.如权利要求2所述的一种实际时移地震资料处理差异的可靠性判断方法,其特征在于:所述步骤②中,通过开放前后两个时间点的孔隙度和饱和度这两个三维油藏数模网格属性获取纵波速度、横波速度和密度的三维油藏数模网格属性方法如下:
根据测井数据区分砂层和泥层,分别统计出砂层和泥层的泥质含量和孔隙度的关系;
对于砂层,根据其流体情况,应用流体替换技术,把砂岩的纵波速度和密度的数据转化为饱含水砂岩的纵波速度和密度的数据,统计出饱含水砂岩的纵波速度和密度与孔隙度的关系;对于泥层,直接统计出其纵波速度和密度与孔隙度之间的关系;
对于砂层,应用泥岩公式得到饱含水砂岩的横波速度与纵波速度的关系;对于泥层,应用泥岩公式得到岩石的横波速度和纵波速度的关系;
对于砂岩,将统计出的饱含水砂岩的纵波速度、横波速度、密度和泥质含量的公式输入到三维油藏数模网格属性数据中,应用流体替换计算得到实际地层的纵波速度、横波速度和密度;对于泥岩,统计出的纵波速度、横波速度和密度即为实际地层数据。
5.如权利要求2所述的一种实际时移地震资料处理差异的可靠性判断方法,其特征在于:所述步骤③中,利用正演模拟方法获得开发前后三维油藏模型的地震响应记录方法如下:
根据式(1)计算开发前后三维油藏模型中各层不同位置的反射系数:
其中,ti代表第i网格点的时间,R(ti)是第i网格点的反射系数,ρ(ti+1)是第i网格点界面上部的密度,v(ti+1)是第i网格点界面上部的纵波速度,ρ(ti)是第i网格点界面下部的密度,v(ti)是第i网格点界面下部的纵波速度,t代表时间;
选取与实际地震资料主频一致的子波;
根据式(2),通过子波与反射系数褶积得到开发前后三维油藏模型的模拟地震响应记录:
其中,S(t)是模拟地震响应记录,b(t)代表地震子波。
6.如权利要求1所述的一种实际时移地震资料处理差异的可靠性判断方法,其特征在于:在所述步骤2)中,获得开发前后实际地震资料处理结果的振幅能量差异的过程如下:
①采集开发前后实际地震资料;
②利用开发前后实际地震资料进行数据处理得到开发前后实际地震的地震响应记录;
③计算开发前后实际地震资料处理结果的差异;
④计算实际地震资料差异数据的均方根振幅作为其振幅能量差异。
7.如权利要求1所述的一种实际时移地震资料处理差异的可靠性判断方法,其特征在于:所述步骤3)中,计算实际地震资料处理结果与油藏模型正演模拟的地震记录的振幅能量差异的相关性的方法如下:
比较开发前后油藏模型正演模拟与实际地震资料处理结果沿层振幅能量差异的一致性,如果两者振幅能量分布、强弱关系等一致性较高,利用下述相关性值公式(3)计算对应井点位置实际资料处理和油藏模型模拟的差异地震数据在目的层段的相关性,
其中,f1是井点位置实际资料处理的目的层段差异数据,f2为井点位置油藏模型模拟的目的层段差异数据,t1、t2为目的层段的起止时间。
8.如权利要求1所述的一种实际时移地震资料处理差异的可靠性判断方法,其特征在于:所述步骤4)中,判断实际地震资料处理的结果对时移地震差异的解释是否可靠的方法如下:
如果步骤3)中计算结果相关性值若大于等于90%,则判定实际地震资料处理的结果能够反映油藏参数的变化;如果计算结果表明相关性低于90%,则判定实际地震资料处理的结果不能反映油藏参数的变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710004041.1A CN106646603B (zh) | 2017-01-04 | 2017-01-04 | 一种实际时移地震资料处理差异的可靠性判断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710004041.1A CN106646603B (zh) | 2017-01-04 | 2017-01-04 | 一种实际时移地震资料处理差异的可靠性判断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106646603A true CN106646603A (zh) | 2017-05-10 |
CN106646603B CN106646603B (zh) | 2018-11-27 |
Family
ID=58843125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710004041.1A Active CN106646603B (zh) | 2017-01-04 | 2017-01-04 | 一种实际时移地震资料处理差异的可靠性判断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106646603B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108802816A (zh) * | 2018-04-11 | 2018-11-13 | 中石化石油工程技术服务有限公司 | 城市地下空间勘探方法及系统 |
CN111239805A (zh) * | 2020-02-13 | 2020-06-05 | 中国石油大学(北京) | 基于反射率法的块约束时移地震差异反演方法及系统 |
CN112014881A (zh) * | 2020-08-27 | 2020-12-01 | 中海石油(中国)有限公司 | 基于时移地震的水驱速度预测方法 |
CN112649844A (zh) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | 基于地震振幅差的反演方法及系统 |
CN113484914A (zh) * | 2021-07-13 | 2021-10-08 | 中海石油(中国)有限公司 | 海上风浪一致性影响量板制作方法、系统、介质及设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101937101A (zh) * | 2010-08-04 | 2011-01-05 | 中国海洋石油总公司 | 一种鉴定能否实施时移地震的方法 |
CN102508294A (zh) * | 2011-10-20 | 2012-06-20 | 西北大学 | 一种利用时移地震勘探资料进行差异avo分析的方法 |
US20140343859A1 (en) * | 2013-05-15 | 2014-11-20 | Conocophillips Company | Time-lapse 4d scattering for imaging hydraulically induced fractures |
CN104536046A (zh) * | 2014-12-26 | 2015-04-22 | 吉林大学 | 基于地震记录的震源激发信号一致性评价方法 |
US20160097870A1 (en) * | 2014-10-03 | 2016-04-07 | Partha S. Routh | Seismic Survey Design Using Full Wavefield Inversion |
-
2017
- 2017-01-04 CN CN201710004041.1A patent/CN106646603B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101937101A (zh) * | 2010-08-04 | 2011-01-05 | 中国海洋石油总公司 | 一种鉴定能否实施时移地震的方法 |
CN102508294A (zh) * | 2011-10-20 | 2012-06-20 | 西北大学 | 一种利用时移地震勘探资料进行差异avo分析的方法 |
US20140343859A1 (en) * | 2013-05-15 | 2014-11-20 | Conocophillips Company | Time-lapse 4d scattering for imaging hydraulically induced fractures |
US20160097870A1 (en) * | 2014-10-03 | 2016-04-07 | Partha S. Routh | Seismic Survey Design Using Full Wavefield Inversion |
CN104536046A (zh) * | 2014-12-26 | 2015-04-22 | 吉林大学 | 基于地震记录的震源激发信号一致性评价方法 |
Non-Patent Citations (1)
Title |
---|
李景叶 等: "时移地震油藏监测可行性分析评价技术", 《石油物探》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108802816A (zh) * | 2018-04-11 | 2018-11-13 | 中石化石油工程技术服务有限公司 | 城市地下空间勘探方法及系统 |
CN112649844A (zh) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | 基于地震振幅差的反演方法及系统 |
CN112649844B (zh) * | 2019-10-10 | 2023-08-15 | 中国石油化工股份有限公司 | 基于地震振幅差的反演方法及系统 |
CN111239805A (zh) * | 2020-02-13 | 2020-06-05 | 中国石油大学(北京) | 基于反射率法的块约束时移地震差异反演方法及系统 |
CN111239805B (zh) * | 2020-02-13 | 2021-02-05 | 中国石油大学(北京) | 基于反射率法的块约束时移地震差异反演方法及系统 |
CN112014881A (zh) * | 2020-08-27 | 2020-12-01 | 中海石油(中国)有限公司 | 基于时移地震的水驱速度预测方法 |
CN112014881B (zh) * | 2020-08-27 | 2023-10-27 | 中海石油(中国)有限公司 | 基于时移地震的水驱速度预测方法 |
CN113484914A (zh) * | 2021-07-13 | 2021-10-08 | 中海石油(中国)有限公司 | 海上风浪一致性影响量板制作方法、系统、介质及设备 |
CN113484914B (zh) * | 2021-07-13 | 2023-09-12 | 中海石油(中国)有限公司 | 海上风浪一致性影响量板制作方法、系统、介质及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN106646603B (zh) | 2018-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106646603A (zh) | 一种实际时移地震资料处理差异的可靠性判断方法 | |
US5475589A (en) | System for evaluating seismic sequence lithology and property, and for evaluating risk associated with predicting potential hydrocarbon reservoir, seal, trap or source | |
US5586082A (en) | Method for identifying subsurface fluid migration and drainage pathways in and among oil and gas reservoirs using 3-D and 4-D seismic imaging | |
CN109425896A (zh) | 白云岩油气储层分布预测方法及装置 | |
Eberli et al. | Seismic Imaging of Carbonate Reservoirs and Systems: AAPG Memoir 81 | |
CA2246151C (en) | A method for determining, displaying and utilizing the relative changes of elastic moduli and that of density of geologic formations | |
WO2006031383A2 (en) | Geologic models of subsurface sedimentray volumes | |
CN104678438B (zh) | 一种co2地质封存中四维地震资料co2分布预测的方法 | |
CN105319585B (zh) | 一种利用薄层干涉振幅恢复识别油气藏的方法 | |
Li | Geophysical exploration technology: Applications in lithological and stratigraphic reservoirs | |
Kianoush et al. | Compressional and shear interval velocity modeling to determine formation pressures in an oilfield of SW Iran | |
CN103293553B (zh) | 一种复杂海底上下缆地震采集数据边界元延拓校正方法 | |
CN107462924A (zh) | 一种不依赖于测井资料的绝对波阻抗反演方法 | |
CN106842326B (zh) | 无横波速度测井时砂泥互层co2地质封存时移地震正演模拟方法 | |
CN107765308A (zh) | 基于褶积思想与精确震源的重构低频数据频域全波形反演方法 | |
CN104330826A (zh) | 一种去除复杂地表条件下多种噪音的方法 | |
CN103543478A (zh) | 地质形态插值的km方法 | |
CN101937101B (zh) | 一种鉴定能否实施时移地震的方法 | |
CN111983683A (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
CN113219531A (zh) | 致密砂岩气水分布的识别方法及装置 | |
Reda et al. | Hydrocarbon Reservoir Characterization in the Challenging Structural Setting of Southern Gulf of Suez: Synergistic Approach of Well Log Analyses and 2D Seismic Data Interpretation | |
Denham | Seismic interpretation | |
Pennington et al. | Calibration of seismic attributes for reservoir characterization | |
CN104316959B (zh) | 一种基于等效流体声波阻抗的流体识别方法 | |
CN109613615A (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 |