CN107065011B - 一种应用于陆相盆地储层反演的曲线频率融合方法 - Google Patents
一种应用于陆相盆地储层反演的曲线频率融合方法 Download PDFInfo
- Publication number
- CN107065011B CN107065011B CN201710264472.1A CN201710264472A CN107065011B CN 107065011 B CN107065011 B CN 107065011B CN 201710264472 A CN201710264472 A CN 201710264472A CN 107065011 B CN107065011 B CN 107065011B
- Authority
- CN
- China
- Prior art keywords
- curve
- frequency
- amplitude
- seismic
- log
- 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
- 238000007500 overflow downdraw method Methods 0.000 title claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 12
- 230000004927 fusion Effects 0.000 claims abstract description 10
- 239000000203 mixture Substances 0.000 claims abstract description 6
- 230000008569 process Effects 0.000 claims abstract description 5
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 4
- 230000002269 spontaneous effect Effects 0.000 claims description 39
- 238000012937 correction Methods 0.000 claims description 38
- 238000001228 spectrum Methods 0.000 claims description 30
- 238000000265 homogenisation Methods 0.000 claims description 17
- 239000004575 stone Substances 0.000 claims description 13
- 239000002131 composite material Substances 0.000 claims description 12
- 238000001914 filtration Methods 0.000 claims description 9
- 238000011161 development Methods 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 230000006641 stabilisation Effects 0.000 claims description 3
- 238000011105 stabilization Methods 0.000 claims description 3
- 230000008021 deposition Effects 0.000 abstract description 4
- 239000003921 oil Substances 0.000 description 12
- 239000010410 layer Substances 0.000 description 11
- 244000018633 Prunus armeniaca Species 0.000 description 10
- 235000009827 Prunus armeniaca Nutrition 0.000 description 10
- 238000009826 distribution Methods 0.000 description 6
- 239000004576 sand Substances 0.000 description 6
- 239000011435 rock Substances 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 235000009754 Vitis X bourquina Nutrition 0.000 description 1
- 235000012333 Vitis X labruscana Nutrition 0.000 description 1
- 240000006365 Vitis vinifera Species 0.000 description 1
- 235000014787 Vitis vinifera Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000155 melt Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000009466 transformation Effects 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
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
本发明公开一种应用于陆相盆地储层反演的曲线频率融合方法，包括如下步骤：（a）基础曲线预处理；（b）地震合成记录，建立时深关系曲线；（c）曲线频谱分析及滤波处理；（d）多频率曲线融合；（e）基于频率融合曲线的地震反演。本发明提出了一种基于频率融合的曲线融合方法，在陆相盆地储层储层反演中根据测井曲线的不同频率尺度反映不同的地质与沉积特征的特点，利用频率融合的方法建立拟波阻抗测井曲线。将曲线频率融合成果应用于地震反演，解决了地震反演精细储层刻画的难题，有利于提高地震反演储层预测的精度和效率。
Description
技术领域
本发明内容涉及陆相河流储层精细研究范畴，具体涉及一种利用地震联合测井资料的综合频率融合方法，主要利用自然电位低频分量及电阻率曲线高频信息综合应用到储层预测，以提高储层刻画精度。
背景技术
目前地震反演方法广泛应用于储层预测之中，其主要根据地震储层可识别的波阻抗差异，实现反演成果能够有效体现储集层与围岩边界。单一测井曲线对陆相碎屑岩储层与围岩差异响应具有局限性，难以有效体现两者地球物理性质差异。采用测井曲线多频段融合方法构建对储层特征敏感的曲线，突出储集层的地球物理响应，以达到精细储层刻画目的。
综合测井、地震、地质等多种资料开展地震反演工作。最终实现地球物理多信息融合，突出各信息的优势分量。各种地球物理资料中，声波阻抗与声波时差测井曲线起到了关键作用。但波阻抗和声波速度曲线受各种因素限制往往与真实地层条件匹配不好，单一曲线的精度局限性影响到地震反演储集层的精确分布。自然电位、电阻率等非速度类曲线在实际研究中对岩性识别发挥明显优势，虽然与地震资料没有与地震反射没有直接对应关系，但可探求一种非速度类测井曲线与速度类测井曲线融合曲线方法，综合体现地层背景与岩性背景信息，相当于在地震反演过程中加入了丰富的地质先验信息和岩石物性信息，使得反演结果更加准确。
发明内容
本发明的目的在于解决目前技术中存在的不足，提供一种应用于陆相盆地储层反演的曲线频率融合方法，提高地震反演预测的精度。
本发明的技术方案是这样实现的，一种应用于陆相盆地储层反演的曲线频率融合方法，包括如下步骤：
（a）基础曲线预处理；
（b）地震合成记录，建立时深关系曲线；
（c）曲线频谱分析及滤波处理；
（d）多频率曲线融合；
（e）基于频率融合曲线的地震反演。
(a1)对声波时差曲线开展标准化处理，选取全区稳定发育的泥岩段作为 标准层，计算标准井标准层段的测井曲线幅度平均值同时计算每口井 目的层段测井曲线平均值求出校正差，从测井曲线中减去校正量，即完成 校正。校正公式为：
其中，AC标为标准化后测井曲线值；AC为标准化前测井曲线值；为待标 准化井标准层测井曲线平均值；为标准井标准层测井曲线平均值。
(a2)对自然电位曲线开展基线校正，选取全区稳定发育、测井响应具有 一致性的稳定的泥岩段作为标准层，读取标准层自然电位曲线平均值作为自然 电位曲线基线漂移的基值，选取目的层内各个泥岩层深度及其对自然电位幅值， 按照自然电位曲线采样率进行垂向上线性插值，生成校正基线曲线SP基，对自 然电位曲线进行基线校正，利用公式：
其中，为第k个采样点处基线校正后的自然电位幅度值；Spk为第k 个采样点处基线校正前的自然电位幅度值；为第k个采样点处校正基准曲 线幅度值；SP标准为标准泥岩段自然电位平均幅值。
(a3)将步骤(a2)校正结果SPNEW曲线归一化到波阻抗曲线量纲范围内， 根据读取波阻抗曲线开展频率直方图统计，得出波阻抗曲线的最大幅值最 小幅值由于自然电位幅值极性与波阻抗曲线幅值极性相反，对步骤(a2) 校正结果SPNEW曲线开展反极性求取SP-NEW＝SPNEW，并对SP-NEW曲线进行频率直方 图统计，统计SP-NEW曲线最大值和最小值将SPNEW曲线量纲到波 阻抗曲线的范围之内，均一化公式如下：
其中，为第k个采样点处均一化后的自然电位幅度值；为第k 个采样点处均一化前的自然电位幅度值；和为波阻抗曲线的最大幅 值和最小幅值；和为井待校正自然电位曲线的最大幅值和最小幅 值。
(a4)引用步骤(a3)所读取的波阻抗量曲线最大幅值和最小幅值 对步骤(a2)校正结果SPNEW曲线进行频率直方图统计，得出电阻率曲 线最大值和最小值将电阻率曲线量纲到波阻抗曲线的范围 之内，均一化公式如下：
其中，为第k个采样点处均一化后的电阻率幅度值；为第k 个采样点处均一化前的电阻率幅度值；和为波阻抗曲线的最大幅值 和最小幅值；和为井待校正电阻率曲线的最大幅值和最小幅 值。
通过步骤(a3)和步骤(a4)将完成自然电位曲线和电阻率曲线的归一化处 理，将自然电位曲线和电阻率曲线归一到波阻抗曲线量纲范围，以此达到各井 曲线的一致性和可对比性。
上述一种应用于陆相盆地储层反演的曲线频率融合方法，在步骤(b)中：
(b1)制作地震合成记录，建立时深关系曲线，利用步骤(a1)标准化的 声波曲线和密度曲线，根据公式：
其中，Rk为第k个采样点处换算反射系数值；ρk为第k个采样点处密度曲 线幅度值；vk为第k个采样点处速度曲线幅值，其大小为
换算出反射系数。并选用合适频率的雷克子波δ，将反射系数曲线与子波褶 积得到合成记录：
其中，Xi为第i个采样点合成记录值；δk为第k个采样点处子波值；Ri为第 i个采样点处反射系数；M1、M2分别为子波的起始样点和末尾样点序号。
以合成记录为基础，建立地震数据的时间域与测井深度域之间的对应关系。
（b2）应用步骤（b1）所建立地震数据的时间域与测井深度域之间的对应关系，将步骤（a3）中所得到SPlog曲线转换为地震反射时间与自然电位幅值相对应的SPtime曲线。
（b3）应用步骤（b1）所地震数据的时间域与测井深度域之间的对应关系，将步骤（a3）中所得到RLLDlog曲线转换为地震反射时间与电阻率幅值相对应的RLLDtime曲线。
上述一种应用于陆相盆地储层反演的曲线频率融合方法，在步骤(c)中：
(c1)对步骤(b1)中的得到SPtime曲线进行快速傅里叶变换(FFT)，得到 SPtime曲线频率谱H(f)SP，利用低通滤波器公式：
H(f)SPlow＝C(f)×H(f)SP
其中，C(f)为低通滤波函数；f为频率；fL为低通滤波器的截止频率；H(f)SP为自然电位曲线频谱函数；H(f)SPlow为自然电位曲线低频部分。
通过低通滤波处理得到低频分量频谱曲线H(f)SPlow。
(c2)对步骤(b2)中的得到RLLDlog曲线进行快速傅里叶变换(FFT)得 到RLLDtime曲线频谱H(f)RLLD，根据高通滤波器公式：
H(f)RLLDhigh＝C(f)×H(f)RLLD
其中，C(f)为高通滤波函数；f为频率；fH为高通滤波器的截止频率； H(f)RLLD为自然电位曲线频谱函数；H(f)RLLDhigh为电阻率曲线高频部分。
通过高通滤波处理得到高频分量频谱曲线H(f)RLLDhigh。
截止频率fL和fH的选取为地震反演所应用到的地震数据的有效频带的最高 值。
上述一种应用于陆相盆地储层反演的曲线频率融合方法，在步骤(d)中：
通过融合技术将步骤(c2)得到的自然电位低频分量频谱曲线H(f)SPlow和 步骤(c3)得到的电阻率高频分量频谱曲线H(f)RLLDhigh进行融合生成频谱曲 线H(f)Ip。针对频谱曲线H(f)Ip开展反傅里叶变换，即得到基于频率融合的 拟波阻抗曲线Ip拟合。
上述一种应用于陆相盆地储层反演的曲线频率融合方法，在步骤(e)中： 将步骤(d)所得频率融合拟声波阻抗曲线Ip拟合应用于地震反演运算过程。
本发明的有益效果是：
本发明利用“低频补偿高频恢复”的思路开展曲线频率融合，能够更好的识别薄厚不同的砂岩，突破单一因素的局限性。根据测井曲线的不同频率尺度反映不同的地质与沉积特征的特点，采用自然电位曲线的低频分量与电阻率曲线高频分量相融合，融合出拟波阻抗曲线。融合后的拟波阻抗曲线对砂岩和泥岩的区分较之前有很大改进，对薄厚不同尺度砂岩响应非常敏感，可提高地震反演储层横向预测精度。
附图说明
图1a为X4-3-B602井标准井标准层选取图。
图1b为杏北东部过渡带地区声波时差曲线校正前直方图。
图1c为杏北东部过渡带地区声波时差曲线校正后直方图。
图2a为X4-4-B431井基线校正前曲线自然电位曲线图。
图2b为X4-4-B431井基线校正前曲线自然电位曲线图。
图3为杏北东部过渡带地区波阻抗曲线直方图。
图4a为X4-4-B431井均一化校正后自然电位曲线图。
图4b为X4-4-B431井均一化校正后电阻率曲线图。
图5为X4-4-B431井合成记录图。
图6为依据X4-4-B431井合成记录所求时深转换关系公式图。
图7a为时深转化后X4-4-B431井时间域自然电位曲线图。
图7b为时深转化后X4-4-B431井时间域电阻率曲线图。
图8a为X4-4-B431井自然电位曲线低频分量曲线图。
图8b为X4-4-B431井电阻率曲线高频分量曲线图。
图8c为X4-4-B431井曲线融合成果图。
图9a为基于频率融合曲线的地震反演过井剖面。
图9b为基于模型的地震波阻抗反演过井剖面。
图10为本发明所述方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合附图对本发明做进一步的说明：
本实施例一种应用于陆相盆地储层反演的曲线频率融合方法在杏北东部过渡带地区萨尔图油层进行试验：
一、试验背景介绍：
大庆油田杏北开发区东部过渡带位于杏树岗背斜构造北部的东翼，是指葡Ⅰ内油水边界至萨Ⅱ外油水边界地区，包括四个条带，总面积31.6km2，原始地质储量3175.6×104t。共发育萨尔图、葡萄花两套油层，5个油层组、15个砂岩组、36个小层、51个沉积单元。油层埋藏深度900～1200m左右，砂泥质交互分布，非均质严重。东部过渡带萨尔图油层位于坳陷期姚家组沉积时期，该时期以三角洲外前缘沉积为主，主要是席状砂与泥岩的薄互层，层位相对稳定。岩性以细砂岩、粉砂岩和泥质粉砂岩为主，渗透率较低，岩石颗粒较细。分流河道砂相对较发育，储层中表内薄层和表外储层主要以局部变差和镶边搭桥形式充填于大面积分布的表内厚层之间，其中表内层钻遇率较高，在80%以上。区域沉积的砂体分选好、砂质纯、密度低，砂岩以薄互层发育为主，因此将杏北东部过渡带地区确定为目标区。
二、试验步骤
(a)基础曲线预处理
(a1)如图1所示对杏北东部过渡带地区萨尔图油层声波时差曲线开展标 准化处理，选取全区稳定发育的泥岩段S1-S2夹层作为标准层(图1a)，计算标 准井X4-3-B602井标准层的声波时差曲线平均幅值为同 时计算区内各井在标准层段声波时差曲线平均值求出校正差，从测井曲线 中减去校正量，即完成校正。校正公式为：
其中，AC标为标准化后测井曲线值；AC为标准化前测井曲线值；为待 标准化井标准层测井曲线平均值；为标准井标准层测井曲线平均值。
通过上述公式完成杏北东部过渡带地区萨尔图油层声波时差曲线标准化处理，通过对比声波时差曲线校正前后直方图即图1b、图1c，可知各井标准化处理后曲线正态分布，则说明各井声波时差曲线已经校正到同一刻度水平上。
(a2)如图2所示对X4-4-B431井自然电位曲线开展基线校正，选取全区 稳定发育、测井响应具有一致性的稳定的S1-S2夹层作为标准层，如图2a读取 标准层段自然电位曲线平均值SP基＝22.2716mv作为自然电位曲线基线漂移的基 值，选取萨尔图油层内各个泥岩层深度及其对应自然电位幅值，按照自然电位 曲线采样率进行垂向上线性插值，生成校正基线曲线SP基图2a中红色曲线，对 自然电位曲线进行基线校正，利用公式：
其中，为第k个采样点处基线校正后的自然电位幅度值；Spk为第k 个采样点处基线校正前的自然电位幅度值；为第k个采样点处校正基准曲 线幅度值；SP标准为标准泥岩段自然电位平均幅值。
完成了如图2b所示X4-4-B431井自然电位曲线基线校正，将萨尔图油层内各个泥岩层位置校正到相同基线位置，消除了自然电位曲线由于水淹等影响产生的基线偏移现象，增强了单井自然电位曲线对砂泥岩的区分能力。
(a3)如图3杏北东部过渡带地区波阻抗曲线直方图所示，将步骤(a2)校 正所得SPNEW曲线归一化到波阻抗曲线量纲范围内，根据读取波阻抗曲线开展 频率直方图统计，得出波阻抗曲线的最大幅值最 小幅值由于自然电位幅值极性与波阻抗曲线幅值 极性相反，对步骤(a2)校正结果SPNEW曲线反极性求取SP-NEW＝SPNEW，并对SP-NEW曲线进行频率直方图统计，统计SP-NEW曲线最大值和最小值将 SPNEW曲线量纲到波阻抗曲线的范围之内，均一化公式如下：
其中，为第k个采样点处均一化后的自然电位幅度值；为第k 个采样点处均一化前的自然电位幅度值；和为波阻抗曲线的最大幅值 和最小幅值；和为井待校正自然电位曲线的最大幅值和最小幅值。
完成了如图4a所示X4-4-B431井量纲化校正后自然电位曲线成果。
(a4)如图4b所示将X4-4-B431井电阻率曲线量纲到波阻抗范围内，对电 阻率曲线进行频率直方图统计，得出电阻率曲线最大值和最小值 利用步骤(a3)所读取的波阻抗量曲线最大幅值和最小幅值将电阻率曲线量纲到波阻抗曲线的范围之内，均一化公式如下：
其中，为第k个采样点处均一化后的电阻率幅度值；为第k 个采样点处均一化前的电阻率幅度值；和为波阻抗曲线的最大幅值 和最小幅值；和为待校正电阻率曲线的最大幅值和最小幅 值。
通过步骤（a3）和步骤（a4）将完成自然电位曲线和电阻率曲线的均一化处理，将自然电位和电阻率归一到波阻抗曲线量纲范围，增加了各井曲线的一致性和可对比性。
(b)地震合成记录，建立时深关系曲线
(b1)如图5X4-4-B431井合成记录图所示，利用步骤(a1)标准化的声波 曲线和密度曲线利用公式：
其中，Rk为第k个采样点处换算反射系数值；ρk为第k个采样点处密度曲 线幅度值；vk为第k个采样点处速度曲线幅值，其大小为
换算出反射系数。并选用合适频率的雷克子波δ，将反射系数曲线与子波褶 积得到合成记录：
其中，Xi为第i个采样点合成记录值；δk为第k个采样点处子波值；Ri为第 i个采样点处反射系数；M1、M2分别为子波的起始样点和末尾样点序号。
以合成记录为基础，建立图6依据X4-4-B431井合成记录所求时深转换关系公式图，确定地震数据的时间域与测井测量深度域之间的对应关系。
（b2）应用图6所建立时深关系，将步骤（a3）中所得到SPlog曲线转换为如图7a所示地震反射时间与自然电位幅值相对应的SPtime曲线。
（b3）应用图6所建立时深关系，将步骤（a4）中所得到RLLDlog曲线转换为如图7b所示地震反射时间与电阻率幅值相对应的RLLDtime曲线。
(c)曲线频谱分析及滤波处理
(c1)对步骤(b1)中的得到SPtime曲线进行快速傅里叶变换(FFT)得到 SPtime曲线频率谱H(f)SP，分析地震反演所应用到的地震数据有效频带分布，选 取fL＝80Hz作为低通滤波器截止频率；利用低通滤波器公式：
H(f)SPlow＝C(f)×H(f)SP
其中，C(f)为低通滤波函数；f为频率；fL为低通滤波器的截止频率；H(f)SP为自然电位曲线频谱函数；H(f)SPlow为自然电位曲线低频部分。
通过低通滤波处理得到低频分量频谱曲线H(f)SPlow，对H(f)SPlow进行反 傅里叶变换得到图8a所示X4-4-B431井自然电位曲线低频分量曲线SPlog。
(c2)对步骤(b2)中的得到RLLDlog曲线进行快速傅里叶变换(FFT)得 到RLLDtime曲线频谱H(f)RLLD，分析地震反演所应用到的地震数据有效频带分 布，选取fH＝80Hz作为高通滤波器截止频率，根据高通滤波器公式：
H(f)RLLDhigh＝C(f)×H(f)RLLD
其中，C(f)为高通滤波函数；f为频率；fH为高通滤波器的截止频率； H(f)RLLD为自然电位曲线频谱函数；H(f)RLLDhigh为电阻率曲线高频部分。
通过高通滤波处理得到高频分量频谱曲线H(f)RLLDhigh，对H(f)RLLDhigh进行反 傅里叶变换得到图8b所示X4-4-B431井自然电位曲线低频分量曲线RLLDlog。
(d)多频率曲线融合
通过融合技术将步骤(c2)得到的自然电位低频分量频谱曲线H(f)SPlow和 步骤(c3)得到的电阻率高频分量频谱曲线H(f)RLLDhigh进行融合生成频谱曲线 H(f)Ip。针对频谱曲线H(f)Ip开展反傅里叶变换，即得到如图8c所示既具有地 层背景的低频信息又能反映地层岩性变化的频率融合拟波阻抗曲线Ip拟合。
(e)基于频率融合曲线的地震反演
应用步骤(d)所得频率融合拟波阻抗曲线Ip拟合约束下应用稀疏脉冲反演 方法开展杏北东部过渡带地区萨尔图油层储层预测。
三、实验效果分析
从预测效果来看，应用图8a基于频率融合的反演方法预测储层结果，较图8b基于模型的地震波阻抗反演预测相比，储层与测井钻遇储层一致性更好，薄差储层识别精度更高，储层与非储层之间差异明显，有效地提高薄储层预测精度，增加反演预测可信度。表明了一种应用于陆相盆地储层反演的曲线频率融合方法的技术先进性和实用性。
上述实施例仅仅是为清楚地说明本发明创造所作的举例，而并非对本发明创造具体实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所引伸出的任何显而易见的变化或变动仍处于本发明创造权利要求的保护范围之中。
Claims (1)
1.一种应用于陆相盆地储层反演的曲线频率融合方法，包括如下步骤：
(a)基础曲线预处理；
(b)地震合成记录，建立时深关系曲线；
(c)曲线频谱分析及滤波处理；
(d)多频率曲线融合；
(e)基于频率融合曲线的地震反演；
其中，所述步骤(a)包括如下内容：
(a1)对声波时差曲线开展标准化处理，选取全区稳定发育的泥岩段作为标准层，计算标准井标准层段的测井曲线平均值同时计算每口井目的层段测井曲线平均值求出校正差，从测井曲线中减去校正量，即完成校正，校正公式为：
其中，AC标为标准化后测井曲线值；AC为标准化前测井曲线值；为目的层段测井曲线平均值；为标准井标准层段的测井曲线平均值；
(a2)对自然电位曲线开展基线校正，选取全区稳定发育、测井响应具有一致性的稳定的泥岩段作为标准层，读取标准层自然电位曲线平均值作为自然电位曲线基线漂移的基值，选取目的层内各个泥岩层深度及其对自然电位幅值，按照自然电位曲线采样率进行垂向上线性插值，生成校正基线曲线SP基，对自然电位曲线进行基线校正，利用公式：
其中，为第k个采样点处基线校正后的自然电位幅度值；Spk为第k个采样点处基线校正前的自然电位幅度值；为第k个采样点处校正基准曲线幅度值；SP标准为标准泥岩段自然电位平均幅值；
(a3)将步骤(a2)校正结果SPNEW曲线归一化到波阻抗曲线量纲范围内，根据读取波阻抗曲线开展频率直方图统计，得出波阻抗曲线的最大幅值最小幅值由于自然电位幅值极性与波阻抗曲线幅值极性相反，对步骤(a2)校正结果SPNEW曲线开展反极性求取SP-NEW＝SPNEW，并对SP-NEW曲线进行频率直方图统计，统计SP-NEW曲线最大值和最小值将SPNEW曲线量纲到波阻抗曲线的范围之内，均一化公式如下：
其中，为第k个采样点处均一化后的自然电位幅度值；为第k个采样点处均一化前的自然电位幅度值；和为波阻抗曲线的最大幅值和最小幅值；和为井待校正自然电位曲线的最大幅值和最小幅值；
(a4)引用步骤(a3)所读取的波阻抗量曲线最大幅值和最小幅值对步骤(a2)校正结果SPNEW曲线进行频率直方图统计，得出电阻率曲线最大值和最小值将电阻率曲线量纲到波阻抗曲线的范围之内，均一化公式如下：
其中，为第k个采样点处均一化后的电阻率幅度值；为第k个采样点处均一化前的电阻率幅度值；和为波阻抗曲线的最大幅值和最小幅值；和为井待校正电阻率曲线的最大幅值和最小幅值；
通过步骤(a3)和步骤(a4)将完成自然电位曲线和电阻率曲线的归一化处理，将自然电位曲线和电阻率曲线归一到波阻抗曲线量纲范围，以此达到各井曲线的一致性和可对比性；
所述步骤(b)包括如下内容：
(b1)制作地震合成记录，建立时深关系曲线，利用步骤(a1)标准化的声波时差曲线和密度曲线，根据公式：
其中，Rk为第k个采样点处换算反射系数值；ρk为第k个采样点处密度曲线幅度值；vk为第k个采样点处速度曲线幅值，其大小为
换算出反射系数后，并选用合适频率的雷克子波δ，将反射系数曲线与子波褶积得到合成记录：
其中，Xi为第i个采样点合成记录值；δk为第k个采样点处子波值；Ri为第i个采样点处反射系数；M1、M2分别为子波的起始样点和末尾样点序号；
以合成记录为基础，建立地震数据的时间域与测井深度域之间的对应关系；
(b2)应用步骤(b1)所建立地震数据的时间域与测井深度域之间的对应关系，将步骤(a3)中所得到SPlog曲线转换为地震反射时间与自然电位幅值相对应的SPtime曲线；
(b3)应用步骤(b1)所建立地震数据的时间域与测井深度域之间的对应关系，将步骤(a3)中所得到RLLDlog曲线转换为地震反射时间与电阻率幅值相对应的RLLDtime曲线；
所述步骤(c)包括如下内容：
(c1)对步骤(b1)中的得到SPtime曲线进行快速傅里叶变换(FFT)，得到SPtime曲线频率谱H(f)SP，利用低通滤波器公式：
H(f)SPlow＝C(f)×H(f)SP
其中，C(f)为低通滤波函数；f为频率；fL为低通滤波器的截止频率；H(f)SP为自然电位曲线频谱函数；H(f)SPlow为自然电位曲线低频分量频谱函数；
通过低通滤波处理得到自然电位曲线低频分量频谱函数H(f)SPlow；
(c2)对步骤(b2)中的得到RLLDlog曲线进行快速傅里叶变换(FFT)得到RLLDtime曲线频谱H(f)RLLD，根据高通滤波器公式：
H(f)RLLDhigh＝C(f)×H(f)RLLD
其中，C(f)为高通滤波函数；f为频率；fH为高通滤波器的截止频率；H(f)RLLD为自然电位曲线频谱函数；H(f)RLLDhigh为电阻率曲线高频分量频谱函数；
通过高通滤波处理得到电阻率曲线高频分量频谱函数H(f)RLLDhigh；
截止频率fL和fH的选取为地震反演所应用到的地震数据的有效频带的最高值；
所述步骤(d)包括如下内容：
通过融合技术将步骤(c1)得到的自然电位曲线低频分量频谱函数H(f)SPlow和步骤(c2)得到的电阻率曲线高频分量频谱函数H(f)RLLDhigh进行融合生成频谱曲线H(f)Ip；针对频谱曲线H(f)Ip开展反傅里叶变换，即得到基于频率融合拟波阻抗曲线Ip拟合；
所述步骤(e)包括：将步骤(d)所得频率融合拟波阻抗曲线Ip拟合应用于地震反演运算过程。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710264472.1A CN107065011B (zh) | 2017-06-22 | 2017-06-22 | 一种应用于陆相盆地储层反演的曲线频率融合方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710264472.1A CN107065011B (zh) | 2017-06-22 | 2017-06-22 | 一种应用于陆相盆地储层反演的曲线频率融合方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107065011A CN107065011A (zh) | 2017-08-18 |
| CN107065011B true CN107065011B (zh) | 2018-11-02 |
Family
ID=59600092
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710264472.1A Active CN107065011B (zh) | 2017-06-22 | 2017-06-22 | 一种应用于陆相盆地储层反演的曲线频率融合方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107065011B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109541685B (zh) * | 2017-09-22 | 2020-06-09 | 中国石油天然气股份有限公司 | 一种河道砂体识别方法 |
| CN108756867B (zh) * | 2018-05-11 | 2021-11-19 | 中国地质调查局油气资源调查中心 | 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 |
| CN109001815B (zh) * | 2018-08-27 | 2019-12-31 | 中国石油化工股份有限公司 | 基于地震和井径约束的声波曲线校正方法 |
| CN109001803B (zh) * | 2018-09-26 | 2021-02-23 | 中国石油化工股份有限公司 | 采用拟合微电位曲线提高地震分辨率的方法 |
| CN112824940A (zh) * | 2019-11-20 | 2021-05-21 | 中国石油天然气集团有限公司 | 反射纵波与反射转换横波的振幅匹配方法及装置 |
| CN111827966B (zh) * | 2020-03-25 | 2022-04-15 | 大庆油田有限责任公司 | 一种多井声波测井曲线一致性处理方法、装置及存储介质 |
| CN113759419A (zh) * | 2020-06-04 | 2021-12-07 | 中国石油化工股份有限公司 | 一种储层预测方法、装置、存储介质及电子设备 |
| CN115903026B (zh) * | 2023-01-09 | 2023-05-09 | 东北石油大学三亚海洋油气研究院 | 复合砂体构型解析方法、设备及介质 |
| CN117153289B (zh) * | 2023-09-14 | 2024-04-05 | 大庆油田有限责任公司 | 一种适用于窄薄砂体的油层水淹程度预测方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101206264A (zh) * | 2007-11-08 | 2008-06-25 | 符力耘 | 高分辨率非线性地震波阻抗反演方法 |
| GB2459887A (en) * | 2008-05-09 | 2009-11-11 | Total Sa | Combinating seismic data sets with overlapping bandwidths |
| CN101881153A (zh) * | 2010-06-04 | 2010-11-10 | 中国石油天然气股份有限公司 | 一种常规测井信息融合可视化方法及系统 |
| CN103510947A (zh) * | 2012-06-21 | 2014-01-15 | 中国石油化工股份有限公司 | 滩坝砂岩微相识别模式的建立及其在微相识别中的应用 |
| CN103711482A (zh) * | 2014-01-06 | 2014-04-09 | 山东科技大学 | 一种基于两种测井曲线重叠图快速识别碎屑岩储层的方法 |
| CN103993871A (zh) * | 2014-05-08 | 2014-08-20 | 中国石油化工股份有限公司 | 针对薄互层地层的测井资料标准化处理方法及装置 |
| CN105223618A (zh) * | 2015-09-22 | 2016-01-06 | 中国石油天然气集团公司 | 薄层沉积微相反演方法和装置 |
| CN105988137A (zh) * | 2015-02-11 | 2016-10-05 | 中国石油化工股份有限公司 | 砂砾岩体基于岩心刻度测井的测井特征曲线重构方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105786484B (zh) * | 2016-01-29 | 2019-05-10 | 华为技术有限公司 | 一种应用融合方法及装置 |
| CN106223942A (zh) * | 2016-08-26 | 2016-12-14 | 中国石油新疆油田分公司勘探开发研究院 | 一种基于测井曲线重构的砾岩油藏泥质含量计算方法 |
-
2017
- 2017-06-22 CN CN201710264472.1A patent/CN107065011B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101206264A (zh) * | 2007-11-08 | 2008-06-25 | 符力耘 | 高分辨率非线性地震波阻抗反演方法 |
| GB2459887A (en) * | 2008-05-09 | 2009-11-11 | Total Sa | Combinating seismic data sets with overlapping bandwidths |
| CN101881153A (zh) * | 2010-06-04 | 2010-11-10 | 中国石油天然气股份有限公司 | 一种常规测井信息融合可视化方法及系统 |
| CN103510947A (zh) * | 2012-06-21 | 2014-01-15 | 中国石油化工股份有限公司 | 滩坝砂岩微相识别模式的建立及其在微相识别中的应用 |
| CN103711482A (zh) * | 2014-01-06 | 2014-04-09 | 山东科技大学 | 一种基于两种测井曲线重叠图快速识别碎屑岩储层的方法 |
| CN103993871A (zh) * | 2014-05-08 | 2014-08-20 | 中国石油化工股份有限公司 | 针对薄互层地层的测井资料标准化处理方法及装置 |
| CN105988137A (zh) * | 2015-02-11 | 2016-10-05 | 中国石油化工股份有限公司 | 砂砾岩体基于岩心刻度测井的测井特征曲线重构方法 |
| CN105223618A (zh) * | 2015-09-22 | 2016-01-06 | 中国石油天然气集团公司 | 薄层沉积微相反演方法和装置 |
Non-Patent Citations (2)
| Title |
|---|
| "井震联合多学科油藏技术及在大庆长垣油田中的应用";李占东 等;《地球物理学进展》;20151231;第30卷(第1期);第242-248页 * |
| "地震属性融合技术及其在储层描述中的应用";于正军 等;《特种油气藏》;20131231;第20卷(第6期);第6-9页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107065011A (zh) | 2017-08-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107065011B (zh) | 一种应用于陆相盆地储层反演的曲线频率融合方法 | |
| CN109061765B (zh) | 非均质薄砂岩互层油藏的圈闭评价方法 | |
| CN104047597B (zh) | 含油气泥页岩地层测井曲线标准化方法 | |
| CN103454685B (zh) | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 | |
| CN103135135B (zh) | 一种基于疏松砂岩模型进行烃类定量预测的方法和装置 | |
| CN106405642B (zh) | 基于脱压实声波速度的地震反演储层预测方法 | |
| CN103954995B (zh) | 一种砂岩型铀矿勘查中的砂体识别方法 | |
| CN104898161B (zh) | 一种基于测井响应模拟体的有效砂岩预测方法 | |
| CN105759312A (zh) | 一种储层特征法井震标定方法 | |
| CN105317435B (zh) | 一种水平井裂缝识别方法 | |
| CN101738637B (zh) | 一种基于速度随频率变化信息的油气检测方法 | |
| CN101556337A (zh) | 一种确定地下深层特殊岩性体的方法 | |
| CN101907729A (zh) | 射线弹性参数反演方法 | |
| CN106021793A (zh) | 基于存储系数和渗流系数的低渗储层甜点评价方法 | |
| CN106526671B (zh) | 一种碳酸盐岩储层地震波响应特征分析方法 | |
| CN108072903A (zh) | 一种测井曲线重构方法 | |
| CN106405639A (zh) | 一种叠前地震储层岩性参数的反演方法 | |
| Cyz et al. | Brittleness prediction for the Lower Paleozoic shales in northern Poland | |
| CN107762483A (zh) | 一种基于测井曲线的相关系数及包络面积的流体识别方法 | |
| CN116027416A (zh) | 一种砂砾岩体岩性圈闭识别方法 | |
| CN105158803A (zh) | 一种选取优势储层的方法 | |
| CN106353813A (zh) | 基于阵列声波测井的流体性质识别方法 | |
| Ekine et al. | Delineation of hydrocarbon bearing reservoirs from surface seismic and well log data (Nembe Creek) in Niger Delta oil field | |
| CN109492938A (zh) | 一种基于甜点指示因子的深层碳酸盐岩储层品质评价方法 | |
| Afife et al. | Petrophysical reservoir-rock properties and source-rock characterization of Abu Roash Formation in Wadi El-Rayan oil field, Western Desert, 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 |