CN103592680A - 一种基于正反演的测井数据和深度域地震剖面合成方法 - Google Patents
一种基于正反演的测井数据和深度域地震剖面合成方法 Download PDFInfo
- Publication number
- CN103592680A CN103592680A CN201210286553.9A CN201210286553A CN103592680A CN 103592680 A CN103592680 A CN 103592680A CN 201210286553 A CN201210286553 A CN 201210286553A CN 103592680 A CN103592680 A CN 103592680A
- Authority
- CN
- China
- Prior art keywords
- depth domain
- seismic
- depth
- obtains
- logging speed
- 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
Images
Abstract
本发明是地震勘探技术的一种地震数据处理和地质解释中的基于正反演的测井数据和深度域地震剖面合成方法。利用野外实测的声波测井数据以及岩性速度数据,通过波动方程正反演方法制作深度域合成记录,深度域合成记录用于叠前深度偏移地震剖面的地震地质层位标定,确定地震剖面中反射同相轴的地质层位。本发明适用于深度域复杂构造地震地质解释,是精细构造解释的关键性方法,可提高深度域构造解释的正确性。
Description
技术领域
本发明涉及地震勘探技术,是一种地震数据处理和地质解释中的基于正反演的测井数据和深度域地震剖面合成方法。
背景技术
反射波地震勘探通过在地面激发地震波,利用检波器接受来自于地下地质体的反射信息。通过对已得到的反射信息进行地震资料数据处理得到可用于解释的地震数据。随着地震数据处理技术的发展,地震数据处理结果已经由时间域数据转向深度域数据。时间域地震数据的地震地质层位标定和精细构造解释已经是比较成熟的技术了,但国内外还没有被大家普遍认可的利用测井数据制作深度域合成记录的方法来确定深度域地震数据的地震地质层位。以往深度域的层位标定工作是将钻井分层直接投到深度域地震剖面上,在时间域完成,合成记录过程分以下几个步骤:
第一步:将深度域的声波测井和密度测井数据转换到时间域,并计算反射系数序列:r(tj),j=1,nt。
第二步:从时间域地震剖面中提取地震子波或提供理论子波:w(tk),k=1,nw。
第三步:将第二步得到的子波和第一步得到的反射系数序列进行褶积(式1),得到时间域合成记录序列:s(tj),j=1,nt。
褶积公式:
第四步:把合成记录和井旁地震道进行对比,调整测井速度重新计算反射系数序列
第五步:重复第三步和第四步,直到合成记录和井旁地震道具有较高的一致性。
由于测井和地震速度差异,进行深度域的地震地质综合解释,可能导致钻井分层与地震数据上的同相轴不是同一套地层,如果层位标定不准确,就不能保障解释精度,加大了开发的风险。
发明内容
本发明目的在于提供一种用于深度域层位标定的,基于正反演的合成记录制作方法。测井数据和深度域地震剖面合成方法
本发明通过以下技术步骤实现:
1)钻井取不同岩性的岩石样本,实测的岩石速度,在深度域岩性分界面约束下,离散为深度域低频背景速度;
2)测井并提取高频测井速度和平均密度;
(1)实测得到测井速度,在深度域岩性分界面约束下,计算测井速度在层内的平均速度并离散;
(2)在深度域岩性分界面约束下,计算测井速度在离散采样间隔内的平均速度;
(3)用采样间隔内的平均速度减去测井速度在层内的平均速度,得到高频测井速度;
(4)实测密度测井数据,计算离散采样间隔内的平均密度;
3)利用高频测井速度、低频背景速度和平均密度,用下式计算得到深度域反射系数r(zi),i=1,nz:
式中:vl为步骤1)得到的低频背景速度,vh为步骤2)得到的高频测井速度,ρ为步骤2)得到的平均密度;zi为不同深度的采样点,i是不同深度采样点的编号,nz是采样点的总个数,r(zi),vh(zi),vl(zi),ρ(zi)分别为深度域反射系数、高频测井速度、低频背景速度和平均密度的第i个样点的值;
4)从时间域地震剖面中提取地震子波;
5)用步骤1)得到的低频背景速度,步骤3)得到的反射系数;步骤4)得到的地震子波,用时间二阶差分、空间十阶差分来求解一维初值问题声波方程,得到时间域合成记录;
6)用步骤5)得到的时间域合成记录和步骤1)得到的低频背景速度,利用与步骤5)相同的方式,求解一维初边值问题声波方程,得到深度域合成记录;
7)将步骤6)得到的深度域合成记录和深度偏移剖面的井旁地震道进行对比,对步骤1)得到的低频背景速度进行修正;
8)重复步骤6)和步骤7),直到深度域合成记录和井旁地震道一致。
本发明使得测井数据和深度域地震剖面结合,得到的深度域合成记录包含了界面位置、波数变化和能量衰减信息,和深度偏移剖面特征更为接近,提升了深度域构造解释的可靠性,可直接在深度域进行地震地质综合解释,提高储层预测的准确性。
附图说明
图1深度域低频背景速度;
图2测井速度;
图3测井分层平均速度;
图4测井速度的高频部分;
图5反射系数;
图6本发明深度域合成记录。
具体实施方式
本发明利用修正的反射系数计算公式求取反射系数,通过波动方程正反演来制作深度域合成记录,经过多次迭代,得到与深度偏移剖面中井旁地震道一致的深度域地震剖面合成方法。
本发明的具体实施方式为:
1)在实测岩性分界面深度的约束下,将实测岩性速度离散为深度域低频背景速度(vl(zi),i=1,nz)。
对钻井得到的不同岩性的岩石样本,实测的岩石速度,在深度域岩性分界面约束下,离散为深度域低频背景速度(vl(zi),i=1,nz),见图1。
2)在实测岩性分界面深度约束下,计算测井速度在层内的平均速度并离散(va(i),i=1,nz);计算测井速度在离散网格内的平均速度(vw(zi),i=1,nz)。
从测井数据提取高频测井速度和平均密度:
(1)实测得到测井速度[图2],在深度域岩性分界面约束下,计算测井速度在层内的平均速度并离散(va(i),i=1,nz),见图3;
(2)在深度域岩性分界面约束下,计算测井速度在离散采样间隔内的平均速度(vw(zi),i=1,nz);
(3)用采样间隔内的平均速度减去测井速度在层内的平均速度,得到高频测井速度(vh(zi)=vw(zi)-va(zi),i=1,nz),见图4;
(4)实测密度测井数据,计算离散采样间隔内的平均密度
(ρ(zi),i=1,nz)。
3)利用测井速度的高频分量(vh(zi)=vw(zi)-va(zi),i=1,nz)和步骤1得到的低频背景速度以及密度测井数据(ρ(zi),i=1,nz),利用反射系数计算公式(2),计算深度域反射系数(r(zi),i=1,nz)。
利用高频测井速度、低频背景速度和平均密度,用下式计算得到深度域反射系数(r(zi),i=1,nz),见图5:
4)从时间域地震剖面中提取地震子波或提供理论子波(w(tj),j=1,nt)。
从时间域地震剖面中提取地震子波(w(tj),j=1,nt);
5)利用步骤1)得到的低频背景速度,步骤3)得到的反射系数;步骤4)得到的子波,子波主频(ω)以及品质因子(Q),利用时间二阶差分、空间十阶差分来求解一维初值问题粘滞声波方程(公式3),得到时间域合成记录(trace(tj)=p(z1,tj),j=1,nt)。
p(zi,t-1)=0
p(zi,t0)=0
(3)
式中:Δt,Δx分别为时间采样间隔和深度采样间隔, 为差分系数。
6)利用步骤5)得到的合成记录和步骤1)得到的低频背景速度,利用与步骤5)相同的方式,求解一维初边值问题粘滞声波方程(公式4),得到深度域合成记录(p(zi,t1),i=1,nz),见图6。
p(zi,tnt+2)=0
p(zi,tnt+1)=0
p(z1,tj)=trace(tj)
(4)
7)将步骤6)得到的合成记录和深度偏移剖面的井旁地震道进行对比,对步骤1)得到的低频背景速度进行修正。
8)重复步骤6)和步骤7),直到合成记录和井旁地震道具有较高的一致性。
本发明使得测井数据和深度域地震剖面结合,得到的深度域合成记录包含了界面位置、波数变化和能量衰减信息,和深度偏移剖面特征更为接近,提升了深度域构造解释的可靠性。可直接在深度域进行地震地质综合解释,提高储层预测的准确性,从而提高钻井的成功率,为滚动勘探开发奠定了坚实的基础。
Claims (2)
1.一种基于正反演的测井数据和深度域地震剖面合成方法,特点是通过以下步骤实现:
1)钻井取不同岩性的岩石样本,实测的岩石速度,在深度域岩性分界面约束下,离散为深度域低频背景速度;
2)测井并提取高频测井速度和平均密度;
3)利用高频测井速度、低频背景速度和平均密度,用下式计算得到深度域反射系数r(zi),i=1,nz:
式中:vl为步骤1)得到的低频背景速度,vh为步骤2)得到的高频测井速度,ρ为步骤2)得到的平均密度;zi为不同深度的采样点,i是不同深度采样点的编号,nz是采样点的总个数,r(zi),vh(zi),vl(zi),ρ(zi)分别为深度域反射系数、高频测井速度、低频背景速度和平均密度的第i个样点的值;
4)从时间域地震剖面中提取地震子波;
5)用步骤1)得到的低频背景速度,步骤3)得到的反射系数;步骤4)得到的地震子波,用时间二阶差分、空间十阶差分来求解一维初值问题声波方程,得到时间域合成记录;
6)用步骤5)得到的时间域合成记录和步骤1)得到的低频背景速度,利用与步骤5)相同的方式,求解一维初边值问题声波方程,得到深度域合成记录;
7)将步骤6)得到的深度域合成记录和深度偏移剖面的井旁地震道进行对比,对步骤1)得到的低频背景速度进行修正;
8)重复步骤6)和步骤7),直到深度域合成记录和井旁地震道一致。
2.根据权利要求1的方法,特点是步骤2)所述的高频测井速度和平均密度采用如下步骤得到;
(1)实测得到测井速度,在深度域岩性分界面约束下,计算测井速度在层内的平均速度并离散;
(2)在深度域岩性分界面约束下,计算测井速度在离散采样间隔内的平均速度;
(3)用采样间隔内的平均速度减去测井速度在层内的平均速度,得到高频测井速度;
(4)实测密度测井数据,计算离散采样间隔内的平均密度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210286553.9A CN103592680B (zh) | 2012-08-13 | 2012-08-13 | 一种基于正反演的测井数据和深度域地震剖面合成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210286553.9A CN103592680B (zh) | 2012-08-13 | 2012-08-13 | 一种基于正反演的测井数据和深度域地震剖面合成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103592680A true CN103592680A (zh) | 2014-02-19 |
CN103592680B CN103592680B (zh) | 2016-09-07 |
Family
ID=50082889
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210286553.9A Active CN103592680B (zh) | 2012-08-13 | 2012-08-13 | 一种基于正反演的测井数据和深度域地震剖面合成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103592680B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105353405A (zh) * | 2014-08-21 | 2016-02-24 | 中国石油化工股份有限公司 | 一种全波形反演方法和系统 |
CN106094019A (zh) * | 2016-04-12 | 2016-11-09 | 中国石油化工股份有限公司 | 基于地质信息映射的深度域地层结构反演方法 |
CN106199694A (zh) * | 2016-06-22 | 2016-12-07 | 中国石油化工股份有限公司 | 基于深变子波的合成记录制作方法 |
CN107193040A (zh) * | 2017-06-27 | 2017-09-22 | 中国石油天然气股份有限公司 | 深度域合成地震记录的确定方法和装置 |
CN107229075A (zh) * | 2017-05-02 | 2017-10-03 | 中国石油天然气股份有限公司 | 深度域地震子波的确定方法和装置 |
CN110088647A (zh) * | 2016-09-30 | 2019-08-02 | 艾奎诺能源公司 | 改进的结构建模 |
CN111208564A (zh) * | 2020-02-28 | 2020-05-29 | 中国煤炭地质总局地球物理勘探研究院 | 一种深度域层位标定方法及装置 |
US11333783B2 (en) * | 2018-03-07 | 2022-05-17 | Chengdu University Of Technology | Integrated method for estimation of seismic wavelets and synthesis of seismic records in depth domain |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100142321A1 (en) * | 2008-12-10 | 2010-06-10 | Toshi Thomas Chang | Hybrid modeling in the tau-p domain |
CN102565856A (zh) * | 2010-12-29 | 2012-07-11 | 中国石油天然气集团公司 | 一种基于波动方程正演的近地表噪音压制方法 |
-
2012
- 2012-08-13 CN CN201210286553.9A patent/CN103592680B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100142321A1 (en) * | 2008-12-10 | 2010-06-10 | Toshi Thomas Chang | Hybrid modeling in the tau-p domain |
CN102565856A (zh) * | 2010-12-29 | 2012-07-11 | 中国石油天然气集团公司 | 一种基于波动方程正演的近地表噪音压制方法 |
Non-Patent Citations (3)
Title |
---|
NOZOMU TAKEUCHI EL.: ""Optimally accurate second order time-domain finite difference scheme for computing synthetic seismograms in 2-D and 3-D media"", 《PHYSICS OF THE EARTH AND PLANETARY INTERIORS》, vol. 119, no. 1, 30 April 2000 (2000-04-30), pages 100 - 105 * |
刘财 等: ""用Stokes方程的差分方法制作合成地震记录"", 《石油地球物理勘探》, vol. 37, no. 3, 30 June 2002 (2002-06-30), pages 232 - 236 * |
张雪建 等: ""深度域合成地震记录的制作方法研究"", 《石油地球物理勘探》, vol. 35, no. 3, 30 June 2000 (2000-06-30), pages 377 - 383 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105353405A (zh) * | 2014-08-21 | 2016-02-24 | 中国石油化工股份有限公司 | 一种全波形反演方法和系统 |
CN106094019A (zh) * | 2016-04-12 | 2016-11-09 | 中国石油化工股份有限公司 | 基于地质信息映射的深度域地层结构反演方法 |
CN106199694A (zh) * | 2016-06-22 | 2016-12-07 | 中国石油化工股份有限公司 | 基于深变子波的合成记录制作方法 |
CN110088647A (zh) * | 2016-09-30 | 2019-08-02 | 艾奎诺能源公司 | 改进的结构建模 |
CN107229075A (zh) * | 2017-05-02 | 2017-10-03 | 中国石油天然气股份有限公司 | 深度域地震子波的确定方法和装置 |
CN107229075B (zh) * | 2017-05-02 | 2019-06-11 | 中国石油天然气股份有限公司 | 深度域地震子波的确定方法和装置 |
CN107193040A (zh) * | 2017-06-27 | 2017-09-22 | 中国石油天然气股份有限公司 | 深度域合成地震记录的确定方法和装置 |
US11333783B2 (en) * | 2018-03-07 | 2022-05-17 | Chengdu University Of Technology | Integrated method for estimation of seismic wavelets and synthesis of seismic records in depth domain |
CN111208564A (zh) * | 2020-02-28 | 2020-05-29 | 中国煤炭地质总局地球物理勘探研究院 | 一种深度域层位标定方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103592680B (zh) | 2016-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103592680A (zh) | 一种基于正反演的测井数据和深度域地震剖面合成方法 | |
CN103454685B (zh) | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 | |
US10386531B2 (en) | Geological model analysis incorporating cross-well electromagnetic measurements | |
CN101329407B (zh) | 一种快速p-sv转换波直接模拟以确定地层岩性、岩相变化的方法 | |
CN104267429A (zh) | 确定地层压力的方法及装置 | |
CN104656142A (zh) | 一种利用垂直地震剖面与测井联合的地震层位标定方法 | |
CA2818790C (en) | Seismic trace attribute | |
WO2017035104A1 (en) | Velocity model seismic static correction | |
US8694262B2 (en) | System and method for subsurface characterization including uncertainty estimation | |
CN103135135A (zh) | 一种基于疏松砂岩模型进行烃类定量预测的方法和装置 | |
US20140336940A1 (en) | Estimation of q-factor in time domain | |
CN105527653A (zh) | 一种基于地质信息的虚拟井构建方法 | |
CN103149588B (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
CN103576200A (zh) | 低信噪比区浅层波阻抗界面静校正方法 | |
CN102943669B (zh) | 适用于水平井轨迹设计的高精度深度域随机模拟反演方法 | |
CN104122582B (zh) | 利用叠加速度求取高精度地震波速度的方法 | |
Yusof et al. | Uncertainty analysis of hydrocarbon in place calculation using 3D seismic and well data during appraisal stage–Case study of Goldie Field, offshore Sarawak | |
CN104237946A (zh) | 基于井控的单层反射纵波和反射转换横波的振幅匹配方法 | |
CN111077578B (zh) | 岩层分布预测方法和装置 | |
CN107238859A (zh) | 计算虚拟井合成地震记录的方法及系统 | |
CN104280767A (zh) | 一种基于柯西分布的稀疏脉冲反演方法 | |
CN105445795A (zh) | 倾斜界面的速度分析方法及其系统 | |
CN104345337B (zh) | 一种用于地震反演的时控储层参数建模方法 | |
CN103576189B (zh) | 一种随钻地震直达波组搜索方法 | |
Gu et al. | Investigation of fractures using seismic computerized crosshole tomography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |