CN105388522B - 一种井震联合的沉积旋回匹配分析方法 - Google Patents
一种井震联合的沉积旋回匹配分析方法 Download PDFInfo
- Publication number
- CN105388522B CN105388522B CN201410449179.9A CN201410449179A CN105388522B CN 105388522 B CN105388522 B CN 105388522B CN 201410449179 A CN201410449179 A CN 201410449179A CN 105388522 B CN105388522 B CN 105388522B
- Authority
- CN
- China
- Prior art keywords
- data
- well
- cycle
- sedimentation
- decomposition
- 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
Abstract
本发明提供了一种井震联合的沉积旋回匹配分析方法,属于地球物理勘探地震解释及综合研究领域。首先通过数据的预处理得到相同采样率的地震和测井数据,在此基础上分别对两者采用经验模态分解的方法进行处理,得到不同频谱范围的分解曲线,对分解曲线进行相关性分析,并通过地震相位旋转调整优化相关性,最后得到一个较为满意的匹配结果。该方法能够解决不同尺度资料的地层联合对比问题,可用于井震联合地层对比等勘探综合研究工作。
Description
技术领域
本发明属于地球物理勘探地震数据解释及综合研究领域,具体涉及一种井震联合的沉积旋回匹配分析方法。
背景技术
旋回是地层学和沉积学中表示具有重视性或周期性沉积记录的术语,包含描述性和解释性的双重含义。地层记录的旋回性是沉积作用、沉积环境或驱动机制的不同,地层旋回通常区分为自旋回和它旋回两类。自旋回主要由环境旋回所致,具有局域性和区域性,如河道、三角洲朵叶体的侧向迁移形成的地层旋回。它旋回通常是由环境旋回因素以外的驱动机制所致,具有大区域、跨地区,甚至全球可对比性,如海平面变换旋回和米兰科维奇旋回等。在沉积旋回划分对比过程中,常常受到很多问题的制约。应用地震资料划分沉积旋回,由于受分辨率的限制。只能实现大尺度旋回的追踪,如典型的不整合面或大的海(湖)泛面,不能实现精细的等时对比,而测井数据是载有地层沉积旋回且分辨率更高的信号,但其横向对比难度大。因此,需要通过井震联合的对比方法将两者结合起来,取长补短,才能更好的进行沉积旋回的解释。井震联合对比面临的主要难题是不同尺度的地震旋回和测井旋回该如何匹配?本发明提出了一种井震联合的沉积旋回匹配方法,能够较好解决这一难题,有助于开展井震联合的沉积旋回对比解释研究工作。
发明内容
本发明的目的是针对现有技术的不足,提供一种井震联合的沉积旋回匹配分析方法,以解决由于地震与测井数据分辨率不同造成的对应沉积旋回对比困难的问题。
本发明采用的技术方案是:
首先通过数据的预处理得到相同采样率的地震和测井数据,在此基础上分别对两者采用经验模态分解的方法进行处理,得到不同频谱范围的分解曲线,对分解曲线进行相关性分析,并通过地震相位旋转调整优化相关性,最后得到一个较为满意的匹配结果。
依据上述技术方案的包括以下步骤:
A:加载测井数据与对应井旁道地震数据,作为沉积旋回匹配分析的初始数据;
B:选择一个时深关系将测井数据转换到时间域,按照井旁道地震数据采样率将测井数据重采样;
C:通过经验模态分解方法分别对测井和地震数据进行处理得到不同频率范围的分析曲线;
D:对各个分解曲线进行两两相关性分析,如果得到的相关性较差,旋转地震相位后再重复步骤C,直到得到一个比较满意的匹配结果;
E:利用得到的分解曲线,进行地层沉积旋回对比。
上述步骤C包括下列步骤:
C1:找出测井地震信号的局部最大最小值,将这些最大、最小值用三次样条插值函数拟合形成信号的上下包络;
C2:计算上下包络的平均值,将测井地震信号与平均值相减,得到一个新的信号;
C3:重复上述步骤C1、C2直到新的信号不存在负的局部极大值和正的局部极小值。
上述步骤D的相关性分析方法包括如下步骤:
D1:分别计算井分解曲线、地震分解曲线的数学期望;
D2:分别计算井分解曲线、地震分解曲线的方差;
D3:计算井分解曲线与地震分解曲线的协方差;
D4:计算井分解曲线与地震分解曲线的相关系数。
上述步骤D的地震相位旋转方法包括如下步骤:
DD1:对数据进行离散傅里叶变换,得到数据的实部和虚部;
DD2:将实部数据的前半部分乘以cos(θ),虚部数据的前半部分乘以sin(θ),其中θ为旋转角度;
DD3:实部和虚部后半部分赋零值;
DD4:将实部和虚部进行离散傅里叶反变换,其中实部数据为相位旋转后的数据。
本发明利用经验模态分解将不同频率的地震数据和测井数据进行处理,得到能够反映不同级别沉积旋回的分解曲线,通过给定的时深关系对井震各个分解曲线进行相关性分析,选取相关性最高的井震分解曲线可用于联合开展沉积旋回的对比划分。该方法能够解决不同尺度资料的地层联合对比问题,可用于井震联合地层对比等勘探综合研究工作。
附图说明
图1井震联合沉积旋回匹配分析流程
图2井震数据分解曲线对比图
图3地震相位旋转前井第4分量与地震第2分量对比图
图4地震相位90°旋转后井第4分量与地震第2分量对比图
图5地震道第2分量曲线地层旋回解释示意图
具体实施方式
本发明的一种井震联合的沉积旋回匹配分析方法较佳实施方式为:
a.加载测井数据与对应井旁道地震数据,作为沉积旋回匹配分析的初始数据;选择一个时深关系将测井数据转换到时间域,按照井旁道地震数据采样率将测井数据重采样。具体步骤如下:首先将地震起始时间为开始点,查找离该点最近的对应深度点,将该深度点对应的曲线值作为井数据的起始点,将起始时间按照地震采样率依次查找,并记录对应深度点的曲线值,直到完成井曲线的重采样。
b.通过经验模态分解方法分别对测井和地震数据进行处理得到不同频率范围的分析曲线。主要输入参数是样点个数、数据数组、采样率和终止条件值(SD)。具体步骤如下:找出数据的局部最大最小值,如果极大极小值个数小于4则终止,将这些最大、最小值用三次样条插值函数拟合形成信号的上下包络;求取上下包络的平均值并用原始数据减去该平均值,利用积分法计算终止条件值(SD),如果不满足终止条件则将减去包络平均值的数据作为原始数据继续上一步操作,直到满足终止条件,得到最终分解曲线。
c.对各个分解曲线进行两两相关性分析,如果得到的相关性较差,旋转地震相位后再重复步骤C。直到得到一个比较满意的匹配结果。相关性分析具体步骤如下:分别计算井分解曲线、地震分解曲线的数学期望和方差;利用期望和方差计算协方差;利用协方差计算相关系数。地震相位旋转具体步骤如下:对数据进行离散傅里叶变换,得到数据的实部和虚部;将实部数据的前半部分乘以cos(θ),虚部数据的前半部分乘以sin(θ),其中θ为旋转角度;将实部和虚部后半部分赋零值;将实部和虚部进行离散傅里叶反变换,其中实部数据为相位旋转后的数据。
d.利用得到的分解曲线,进行地层沉积旋回对比。
下面结合附图对本发明作进一步详细描述:
图1为本方法主要步骤和流程,其具体实施方式如下:
(1)加载需要进行匹配分析的数据,通过选择给定的时深关系将数据 匹配到统一的时间域,并按照地震数据采样率进行重采样。
(2)通过经验模态分解方法分别对测井和地震数据进行处理得到不同频率范围的分解曲线。图2为N107井声波测井与过井地震数据的分解曲线图,主要参数设置为采样率2ms、终止条件SD=0.2。
(3)对分解后的曲线进行两两相关性分析,表1为地震相位旋转前各分量间的相关系数,从表中可以看出相关性较好的有(粗体):地震第2分量与井的第4分量相关性为0.21左右,地震第3分量与井的第4分量相关性为0.21,地震第4分量与井的第5分量相关性为0.4,地震第6分量与井第7分量相关性为0.59,地震第6分量与井第8分量相关系数为0.98。其中地震第2分量的频率较高,与井上的沉积旋回较为匹配,可以考虑将地震第2分量曲线作为旋回对比的辅助曲线。
表1地震相位旋转前各分量间的相关系数
(4)对地震数据进行90°相位旋转后,再进行经验模态分解,与井的分量进行相关性对比。表2为地震相位90°旋转后的相关系数,可以看出主要分量间的相关性有较大提高,其中地震第2分量与井第4分量相关性为0.41。图3和图4分别为相位旋转前与旋转后地震第2分量与第4分量曲线对比图,可以看出相位旋转后的井震分量曲线有较好的匹配关系。
表2地震相位90°旋转后各分量间的相关系数
利用与井沉积旋回较为匹配的地震第2分量,可以再无井区进行地层沉积旋回对比与解释,图5为地震第2分量曲线对比示意图。
Claims (6)
1.一种井震联合的沉积旋回匹配分析方法,其特征在于:该方法通过经验模态分解对井震数据进行处理,得到能够反映不同级别沉积旋回的井震分解曲线,通过给定的时深关系对井震各个分解曲线进行相关性分析,通过地震数据相位旋转增强相关性,最终得到井震分辨尺度一致的分解曲线,用于沉积旋回的对比。
2.如权利要求1所述的井震联合的沉积旋回匹配分析方法,其特征在于,该方法包括下列步骤:
A:加载测井数据与对应井旁道地震数据,作为沉积旋回匹配分析的初始数据;
B:选择一个时深关系将测井数据转换到时间域,按照井旁道地震数据采样率将测井数据重采样;
C:通过经验模态分解方法分别对测井和地震数据进行处理得到不同频率范围的分解曲线;
D:对各个分解曲线进行两两相关性分析,如果得到的相关性较差,旋转地震相位后再重复步骤C,直到得到一个比较满意的匹配结果;
E:利用得到的分解曲线,进行地层沉积旋回对比。
3.根据权利要求2所述的井震联合的沉积旋回匹配分析方法,其特征在于,上述步骤C包括下列步骤:
C1:找出测井地震信号的局部最大最小值,将这些最大、最小值用三次样条插值函数拟合形成信号的上下包络;
C2:计算上下包络的平均值,将测井地震信号与平均值相减,得到一个新的信号;
C3:重复上述步骤C1、C2直到新的信号不存在负的局部极大值和正的局部极小值。
4.根据权利要求2或3所述的井震联合的沉积旋回匹配分析方法,其特征在于,上述步骤D的相关性分析方法包括如下步骤:
D1:分别计算井分解曲线、地震分解曲线的数学期望;
D2:分别计算井分解曲线、地震分解曲线的方差;
D3:计算井分解曲线与地震分解曲线的协方差;
D4:计算井分解曲线与地震分解曲线的相关系数。
5.根据权利要求2或3所述的井震联合的沉积旋回匹配分析方法,其特征在于,上述步骤D的地震相位旋转方法包括如下步骤:
DD1:对数据进行离散傅里叶变换,得到数据的实部和虚部;DD2:将实部数据的前半部分乘以cos(θ),虚部数据的前半部分乘以sin(θ),其中θ为旋转角度;
DD3:实部和虚部后半部分赋零值;
DD4:将实部和虚部进行离散傅里叶反变换,其中实部数据为相位旋转后的数据。
6.根据权利要求4所述的井震联合的沉积旋回匹配分析方法,其特征在于,上述步骤D的地震相位旋转方法包括如下步骤:
DD1:对数据进行离散傅里叶变换,得到数据的实部和虚部;DD2:将实部数据的前半部分乘以cos(θ),虚部数据的前半部分乘以sin(θ),其中θ为旋转角度;
DD3:实部和虚部后半部分赋零值;
DD4:将实部和虚部进行离散傅里叶反变换,其中实部数据为相位旋转后的数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410449179.9A CN105388522B (zh) | 2014-09-04 | 2014-09-04 | 一种井震联合的沉积旋回匹配分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410449179.9A CN105388522B (zh) | 2014-09-04 | 2014-09-04 | 一种井震联合的沉积旋回匹配分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105388522A CN105388522A (zh) | 2016-03-09 |
CN105388522B true CN105388522B (zh) | 2018-10-12 |
Family
ID=55421003
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410449179.9A Active CN105388522B (zh) | 2014-09-04 | 2014-09-04 | 一种井震联合的沉积旋回匹配分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105388522B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106019402B (zh) * | 2016-05-12 | 2018-01-30 | 山东科技大学 | 综合判别与验证巨厚煤层内部古气候旋回的方法 |
CN107884830B (zh) * | 2016-09-29 | 2020-03-24 | 中国石油化工股份有限公司 | 沉积地层旋回划分方法及装置 |
CN111830559B (zh) * | 2019-04-19 | 2023-01-10 | 中国石油天然气股份有限公司 | 一种地震资料拓频处理的方法及装置 |
CN112180465B (zh) * | 2020-11-03 | 2022-06-21 | 中国石油化工股份有限公司 | 一种井震联合的地层展布确定方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6246963B1 (en) * | 1999-01-29 | 2001-06-12 | Timothy A. Cross | Method for predicting stratigraphy |
CN102520448A (zh) * | 2012-01-06 | 2012-06-27 | 中国石油化工股份有限公司胜利油田分公司西部新区研究中心 | 年代地层旋回域井震标定系统 |
CN102536195A (zh) * | 2011-12-19 | 2012-07-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 测井沉积层序自动划分方法 |
CN103698808A (zh) * | 2012-09-28 | 2014-04-02 | 中国石油天然气集团公司 | 一种对地震和测井数据波形极值特征点分离与波形重构方法 |
CN103901478A (zh) * | 2012-12-28 | 2014-07-02 | 中国石油天然气集团公司 | 一种井震信息联合确定储层沉积特征和分布的方法 |
-
2014
- 2014-09-04 CN CN201410449179.9A patent/CN105388522B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6246963B1 (en) * | 1999-01-29 | 2001-06-12 | Timothy A. Cross | Method for predicting stratigraphy |
CN102536195A (zh) * | 2011-12-19 | 2012-07-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 测井沉积层序自动划分方法 |
CN102520448A (zh) * | 2012-01-06 | 2012-06-27 | 中国石油化工股份有限公司胜利油田分公司西部新区研究中心 | 年代地层旋回域井震标定系统 |
CN103698808A (zh) * | 2012-09-28 | 2014-04-02 | 中国石油天然气集团公司 | 一种对地震和测井数据波形极值特征点分离与波形重构方法 |
CN103901478A (zh) * | 2012-12-28 | 2014-07-02 | 中国石油天然气集团公司 | 一种井震信息联合确定储层沉积特征和分布的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105388522A (zh) | 2016-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2019062655A1 (zh) | 薄夹层的确定方法和装置 | |
CN103842853B (zh) | 针对覆盖层和记录效应来校正时移地震数据 | |
CN101598805B (zh) | 一种多分量地震资料层位对比和标定方法 | |
CN105388522B (zh) | 一种井震联合的沉积旋回匹配分析方法 | |
van Hoek et al. | Geometric attributes for seismic stratigraphic interpretation | |
CN105954802A (zh) | 一种岩性数据体的转换方法及装置 | |
CN105527653A (zh) | 一种基于地质信息的虚拟井构建方法 | |
Aleardi et al. | Characterisation of shallow marine sediments using high‐resolution velocity analysis and genetic‐algorithm‐driven 1D elastic full‐waveform inversion | |
CN110824563A (zh) | 一种基于Xgboost算法的储层岩性预测方法 | |
Lyu et al. | Multispectral coherence: Which decomposition should we use? | |
CN103364826A (zh) | 基于独立分量分析的地震盲源反褶积方法 | |
Zeng et al. | Recent progress in analysis of seismically thin beds | |
CN102721979A (zh) | 一种基于地震资料的薄层自动解释及厚度预测方法和装置 | |
CN106199694A (zh) | 基于深变子波的合成记录制作方法 | |
CN104199098B (zh) | 利用反射转换横波的偏振特性预测地下裂缝的方法 | |
Shakir et al. | Selection of sensitive post-stack and pre-stack seismic inversion attributes for improved characterization of thin gas-bearing sands | |
RU2490677C2 (ru) | Способ комплексной обработки геофизических данных и технологическая система "литоскан" для его осуществления | |
CN107991716B (zh) | 薄层的沉积相图和沉积体厚度的确定方法和装置 | |
CN111983683B (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
CN107884830B (zh) | 沉积地层旋回划分方法及装置 | |
Pereira et al. | Iterative geostatistical seismic inversion incorporating local anisotropies | |
CN106610506B (zh) | 地震勘探薄层识别方法 | |
CN110988996A (zh) | 一种基于Adaboost回归算法的储层孔隙度预测方法 | |
Lubo-Robles | DEVELOPMENT OF INDEPENDENT COMPONENT ANALYSIS FOR RESERVOIR GEOMORPHOLOGY AND UNSUPERVISED SEISMIC FACIES CLASSIFICATION IN THE TARANAKI BASIN, NEW ZEALAND. | |
Li* et al. | VMD based sedimentary cycle division for unconventional facies analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |