CN103790569A - 一种利用vsp时深关系对声波测井资料进行校正的方法 - Google Patents
一种利用vsp时深关系对声波测井资料进行校正的方法 Download PDFInfo
- Publication number
- CN103790569A CN103790569A CN201210418575.6A CN201210418575A CN103790569A CN 103790569 A CN103790569 A CN 103790569A CN 201210418575 A CN201210418575 A CN 201210418575A CN 103790569 A CN103790569 A CN 103790569A
- Authority
- CN
- China
- Prior art keywords
- depth
- vsp
- data
- degree
- correct
- 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
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明提供了一种利用VSP时深关系对声波测井资料进行校正的方法,属于地震勘探领域。本方法包括:(1)输入校正前的声波数据和该井的VSP时深数据;(2)将所述声波数据转换成时深数据;(3)按VSP分层解释深度或设定的固定深度分别计算出声波测井和VSP的时差ΔT0ac和ΔT0vsp,同时计算出T0ac在深度1和深度2内的声波个数n1;(4)比较声波测井与VSP的时差ΔT0ac和ΔT0vsp,得到声波测井和VSP在同一段深度中的差值:ΔT=ΔT0ac-ΔT0vsp,然后把它均匀地校正到该深度段内每个点上;(5)在完成步骤(4)的校正后,用VSP时间对校正好的声波时间进行标定;(6)按VSP分层解释深度或设定的固定深度类推下去,对每对新的深度1和深度2,重复步骤(3)到(5)。
Description
技术领域
本发明属于地震勘探领域,具体涉及一种利用VSP时深关系对声波测井资料进行校正的方法。
背景技术
地震资料解释前的地层标定,人们常用声波测井资料来计算平均速度,做合成地震记录,然而平均速度的准确对层位的标定结果,合成地震记录与实际地震记录是否符合,这些在很大程度上都与测井资料有关。声波测井资料由于受地层、井径、井孔扩大、泥浆渗入地层等诸多因素的影响,会有不同程度的失真。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种利用VSP时深关系对声波测井资料进行校正的方法,利用VSP资料获得的T0时间准确这一成果对声波测井进行校正,提高了声波测井资料在处理解释应用中的可信度。
本发明是通过以下技术方案实现的:
一种利用VSP时深关系对声波测井资料进行校正的方法,包括以下步骤:
(1)输入校正前的声波数据和该井的VSP时深数据;
(2)将所述声波数据转换成时深数据;
(3)按VSP分层解释深度或设定的固定深度分别计算出声波测井和VSP的时差ΔT0ac和ΔT0vsp,同时计算出T0ac在深度1和深度2内的声波个数n1;
(4)比较声波测井与VSP的时差ΔT0ac和ΔT0vsp,得到声波测井和VSP在同一段深度中的差值:ΔT=ΔT0ac-ΔT0vp,然后把它均匀地校正到该深度段内每个点上;
(5)在完成步骤(4)的校正后,用VSP时间对校正好的声波时间进行标定;
(6)按VSP分层解释深度或设定的固定深度类推下去,对每对新的深度1和深度2,重复步骤(3)到(5),直到完成所有的数据校正后,再把校正好的声波时间转换成声波数据,即得到校正后的声波数据。
所述步骤(2)是利用下式将声波数据转换成时深数据的:
当声波测井采用英制单位,即为微秒/英尺时,k=0.3048,
当声波测井采用米制单位,即为微秒/米时,k=1,
其中H为深度,单位是米,AC为声波测井值,n为声波测井的点数,T0ac(i)为第i点处的时间,单位是毫秒。
所述步骤(3)具体如下:
ΔT0ac=T0ac(深度1)-T0ac(深度2),ΔT0vsp=T0vsp(深度1)-T0vsp(深度2),T0vsp指第(1)步输入的VSP时深数据。
所述步骤(5)具体如下:
Δtt=T0vsp(深度1)-TT0ac(深度1),TT0ac(ii)=TT0ac(ii)+Δtt。
所述步骤(6)中的所述把校正好的声波时间转换成声波数据是采用下式完成的:
其中,acnew(i)为校正后的声波数据,单位为微秒/米。
与现有技术相比,本发明的有益效果是:本发明是VSP资料处理成果应用的又一拓展。垂直地震剖面(VSP)技术是将检波器置于井中接收来自地面激发的地震信号的一种地震勘探技术。由于偏移距小,近似于垂直入射,因此它提供的时深关系既直接又准确。利用VSP资料获得的T0时间准确这一成果对声波测井进行校正,提高了声波测井资料在处理解释应用中的可信度,且本发明方法无需迭代和求解大型方程组,占用计算机资源少、计算量小。
附图说明
图1是本发明利用VSP时深关系对声波测井资料进行校正的方法的步骤框图。
图2是本发明实施例中声波校正前后的测井曲线。
图3-1是本发明实施例中的地面资料。
图3-2是本发明实施例中的校正后的声波合成记录。
图3-3是本发明实施例中的校正前的声波合成记录。
具体实施方式
下面结合附图对本发明作进一步详细描述:
VSP技术是将检波器置于井中接收来自地面激发的地震信号的一种地震勘探技术。由于偏移距小,近似于垂直入射,因此它提供的时深关系既直接又准确。随着VSP(垂直地震剖面)处理和解释技术的不断进步,用VSP资料获得的时深关系的准确度一直为大家所认可。本发明是利用VSP时深关系对声波测井资料进行校正,即利用VSP资料获得的垂直时间T0准确这一成果对声波测井进行校正,提高了声波测井资料在处理解释应用中的可信度。
如图1所示,本发明利用VSP时深关系对声波测井资料进行校正的方法包括以下步骤:
(1)首先输入校正前的声波数据和该井的VSP时深数据(这是VSP资料处理成果之一);
(2)将所述声波数据转换成时间数据:
当声波测井采用英制单位,即为微秒/英尺时,k=0.3048,
当声波测井采用米制单位,即为微秒/米时,k=1,
其中H为深度(单位米),AC为声波测井值,n为声波测井的点数,T0ac(i)为第i点处的时间(单位毫秒)。
(3)利用VSP垂直时间T0准确的优势,按VSP分层解释深度或某个固定深度分别计算出声波测井和VSP的时差ΔT0ac、ΔT0vsp,即:ΔT0ac=T0ac(深度1)-T0ac(深度2),ΔT0vsp=T0vsp(深度1)-T0vsp(深度2),T0vsp指第(1)步输入的VSP时深数据,同时计算出T0ac在深度1和深度2内的声波个数n1(通过编程,条件判断语句和计数器来完成)。
(4)比较声波测井与VSP的时差ΔT0ac、ΔT0vsp,得到声波测井和VSP在同一段深度中的差值:ΔT=ΔT0ac-ΔT0vsp,可以认为声波测井和VSP在这段深度段内时差ΔT,是多种原因的局部系统误差,因此把它均匀地校正到该深度段内每个点上,校正量为: ii从1到n1,TT0ac(ii)为声波在某个点处校正后的时间;
(5)在做好上面的校正后,用VSP时间去标定上面校正好的声波时间,即Δtt=T0vsp(深度1)-TT0ac(深度1),TT0ac(ii)=TT0ac(ii)+Δtt,保证声波时间和VSP时间在控制点上一致,消除累计误差。
(6)按VSP分层解释深度或某个固定深度类推下去(如每隔200米),对新的深度1和深度2,重复步骤(3)到(5),完成所有的数据校正后,再把校正好的声波时间转换成声波数据: acnew(i)为校正后的声波数据,单位:微秒/米。
下面以一口井的声波测井资料为例来说明本发明方法的效果。图2是该井校正前、后的声波测井资料,图2中的实线ac为原始声波曲线,虚线ac-new为校正后的声波曲线,随着井深加深,校正量增大,在这口近5000米深的井中,声波时间的最大校正量为13ms。图3-1至图3-2是用校正前、后声波测井资料做的合成记录和地面地震资料的对比图,其中图3-1是地面资料,图3-3是校正前的声波合成记录,图3-2是校正后的声波合成记录,从图中可以看出1.6秒后,校正后的合成记录与地面资料的吻合度明显好于校正前的。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (6)
1.一种利用VSP时深关系对声波测井资料进行校正的方法,其特征在于:所述方法包括以下步骤;
(1)输入校正前的声波数据和该井的VSP时深数据;
(2)将所述声波数据转换成时深数据;
(3)按VSP分层解释深度或设定的固定深度分别计算出声波测井和VSP的时差ΔT0ac和ΔT0vsp,同时计算出T0ac在深度1和深度2内的声波个数n1;
(4)比较声波测井与VSP的时差ΔT0ac和ΔT0vsp,得到声波测井和VSP在同一段深度中的差值:ΔT=ΔT0ac-ΔT0vsp,然后把它均匀地校正到该深度段内每个点上;
(5)在完成步骤(4)的校正后,用VSP时间对校正好的声波时间进行标定;
(6)按VSP分层解释深度或设定的固定深度类推下去,对每对新的深度1和深度2,重复步骤(3)到(5),直到完成所有的数据校正后,再把校正好的声波时间转换成声波数据,即得到校正后的声波数据。
2.根据权利要求1所述的利用VSP时深关系对声波测井资料进行校正的方法,其特征在于:所述步骤(2)是利用下式将声波数据转换成时深数据的:
当声波测井采用英制单位,即为微秒/英尺时,k=0.3048,
当声波测井采用米制单位,即为微秒/米时,k=1,
其中H为深度,单位是米,AC为声波测井值,n为声波测井的点数,T0ac(i)为第i点处的时间,单位是毫秒。
3.根据权利要求2所述的利用VSP时深关系对声波测井资料进行校正的方法,其特征在于:所述步骤(3)具体如下:
ΔT0ac=T0ac(深度1)-T0ac(深度2),ΔT0vsp=T0vsp(深度1)-T0vsp(深度2),T0vsp指第(1)步输入的VSP时深数据。
5.根据权利要求4所述的利用VSP时深关系对声波测井资料进行校正的方法,其特征在于:所述步骤(5)具体如下:
Δtt=T0vsp(深度1)-TT0ac(深度1),TT0ac(ii)=TT0ac(ii)+Δtt。
6.根据权利要求5所述的利用VSP时深关系对声波测井资料进行校正的方法,其特征在于:所述步骤(6)中的所述把校正好的声波时间转换成声波数据是采用下式完成的:
其中,acnew(i)为校正后的声波数据,单位为微秒/米。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210418575.6A CN103790569B (zh) | 2012-10-26 | 2012-10-26 | 一种利用vsp时深关系对声波测井资料进行校正的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210418575.6A CN103790569B (zh) | 2012-10-26 | 2012-10-26 | 一种利用vsp时深关系对声波测井资料进行校正的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103790569A true CN103790569A (zh) | 2014-05-14 |
CN103790569B CN103790569B (zh) | 2016-12-21 |
Family
ID=50666581
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210418575.6A Active CN103790569B (zh) | 2012-10-26 | 2012-10-26 | 一种利用vsp时深关系对声波测井资料进行校正的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103790569B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104133249A (zh) * | 2014-07-31 | 2014-11-05 | 中国石油天然气集团公司 | 一种微测井资料与声波测井资料联合解释的方法及装置 |
CN105182408A (zh) * | 2015-08-28 | 2015-12-23 | 中国石油天然气集团公司 | 一种合成地震记录的制作方法和装置 |
CN105445802A (zh) * | 2014-09-04 | 2016-03-30 | 中国石油化工股份有限公司 | 一种基于时频谱约束的高精度时深校正方法 |
CN106154329A (zh) * | 2015-04-16 | 2016-11-23 | 中国石油化工股份有限公司 | 利用vsp速度校正声波均方根速度的方法及其系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1882852A (zh) * | 2003-10-24 | 2006-12-20 | 国际壳牌研究有限公司 | 储层区的延时地震勘测 |
CN101231346A (zh) * | 2008-02-26 | 2008-07-30 | 中国矿业大学(北京) | 应用地震波速度预测煤、岩体物理力学参数方法 |
CN101625418A (zh) * | 2009-07-17 | 2010-01-13 | 中国石化集团胜利石油管理局 | 井间地震反射波优选角度成像方法及装置 |
WO2010118624A1 (zh) * | 2009-04-13 | 2010-10-21 | 中国石油集团川庆钻探工程有限公司 | 地震波速度场构建中井约束横向可变h-v曲线构建方法 |
CN102053267A (zh) * | 2010-10-22 | 2011-05-11 | 中国石油化工股份有限公司 | 一种地震剖面资料处理中基于参数反演的vsp波场分离方法 |
-
2012
- 2012-10-26 CN CN201210418575.6A patent/CN103790569B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1882852A (zh) * | 2003-10-24 | 2006-12-20 | 国际壳牌研究有限公司 | 储层区的延时地震勘测 |
CN101231346A (zh) * | 2008-02-26 | 2008-07-30 | 中国矿业大学(北京) | 应用地震波速度预测煤、岩体物理力学参数方法 |
WO2010118624A1 (zh) * | 2009-04-13 | 2010-10-21 | 中国石油集团川庆钻探工程有限公司 | 地震波速度场构建中井约束横向可变h-v曲线构建方法 |
CN101625418A (zh) * | 2009-07-17 | 2010-01-13 | 中国石化集团胜利石油管理局 | 井间地震反射波优选角度成像方法及装置 |
CN102053267A (zh) * | 2010-10-22 | 2011-05-11 | 中国石油化工股份有限公司 | 一种地震剖面资料处理中基于参数反演的vsp波场分离方法 |
Non-Patent Citations (2)
Title |
---|
高志凌等: "用VSP时深关系对声波测井资料进行校正", 《中国地球物理学会第二十三届年会论文集》, 1 October 2007 (2007-10-01) * |
高荣涛: "井地资料多尺度匹配方法研究", 《中国石油大学硕士研究生学位论文》, 28 December 2007 (2007-12-28) * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104133249A (zh) * | 2014-07-31 | 2014-11-05 | 中国石油天然气集团公司 | 一种微测井资料与声波测井资料联合解释的方法及装置 |
CN104133249B (zh) * | 2014-07-31 | 2016-08-31 | 中国石油天然气集团公司 | 一种微测井资料与声波测井资料联合解释的方法及装置 |
CN105445802A (zh) * | 2014-09-04 | 2016-03-30 | 中国石油化工股份有限公司 | 一种基于时频谱约束的高精度时深校正方法 |
CN106154329A (zh) * | 2015-04-16 | 2016-11-23 | 中国石油化工股份有限公司 | 利用vsp速度校正声波均方根速度的方法及其系统 |
CN105182408A (zh) * | 2015-08-28 | 2015-12-23 | 中国石油天然气集团公司 | 一种合成地震记录的制作方法和装置 |
CN105182408B (zh) * | 2015-08-28 | 2017-09-01 | 中国石油天然气集团公司 | 一种合成地震记录的制作方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103790569B (zh) | 2016-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102967882B (zh) | 地层的层速度模型的建模方法 | |
CN106842313B (zh) | 基于方位叠前地震数据的各向异性参数反演方法 | |
CN101630016B (zh) | 一种提高垂直地震剖面成像质量的方法 | |
CN104730579B (zh) | 一种基于表层横波速度反演的纵横波联合静校正方法 | |
CN104656142A (zh) | 一种利用垂直地震剖面与测井联合的地震层位标定方法 | |
US20190345816A1 (en) | Methods and apparatus to measure formation features | |
CN104570079A (zh) | 一种纵波、转换横波地震资料的时间匹配方法 | |
CN105182408A (zh) | 一种合成地震记录的制作方法和装置 | |
CN104570102B (zh) | 一种近地表速度模型与中深层速度模型的融合方法 | |
CN102879819B (zh) | 保持地震波场运动学特征的地震数据处理方法及装置 | |
CN105589100A (zh) | 一种微地震震源位置和速度模型同时反演方法 | |
CN104375172B (zh) | 火山岩下伏地层的构造形态正确成像的方法 | |
CN103592680B (zh) | 一种基于正反演的测井数据和深度域地震剖面合成方法 | |
CN104459782A (zh) | 使用薄层的沿层层析反演地震速度建模方法以及建模单元 | |
CN103149588B (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
CN103576200A (zh) | 低信噪比区浅层波阻抗界面静校正方法 | |
CN103790569A (zh) | 一种利用vsp时深关系对声波测井资料进行校正的方法 | |
CN102707318B (zh) | 虚共炮点射线追踪动态静校正方法 | |
CN102901985A (zh) | 一种适用于起伏地表的深度域层速度修正方法 | |
CN105093319A (zh) | 基于三维地震数据的地面微地震静校正方法 | |
CN103576197A (zh) | 一种转换波角道集抽取方法 | |
CN101201409A (zh) | 一种地震数据变相位校正方法 | |
CN105093326A (zh) | 一种地震勘探资料中波长剩余静校正处理方法 | |
CN105626056A (zh) | 一种确定目标区三维地层孔隙压力的方法和装置 | |
CN104570119A (zh) | 一种三维垂直地震剖面反射波拉伸校正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |