• / 8
  • 下载费用:1 下载币  

高精度地震时深转换方法研究及应用_徐立恒_图文

关 键 词:
高精度 地震 转换 方法 研究 应用 徐立恒 图文
资源描述:
 第 44卷第 5期2014年 09月吉 林 大 学 学 报 (地 球 科 学 版 )Journal of Jilin University(Earth Science Edition)Vol.44 No.5Sep.2014 徐 立 恒 ,鲜波 ,薛玉英 ,等 .高精度地震时深转换方法研究及应用 .吉林大学学报 :地球科学版 ,2014,44(5):1712-1719.doi:10.13278/j.cnki.juese.201405306.Xu Liheng,Xian Bo,Xue Yuying,et al.Study and Application on Seismic Time-Depth Conversion with High-Precision.Journal of Jilin University:Earth Science Edition,2014,44(5):1712-1719.doi:10.13278/j.cnki.juese.201405306.高精度地震时深转换方法研究及应用徐立 恒1,鲜波2,薛玉 英3,李彦 斌41.大庆油田勘探开发研究院,黑龙 江 大庆1637122.中油塔里木油田分公司勘探开发研究院,新疆 库尔勒8410003.东方地球物理公司研究院地质研究中心,河 北 涿州0727514.大庆油田第六采油厂,黑龙 江 大庆163000摘 要 :地震体的时深转换是利用地震资料进行构造及储层解释的一个非常关键的环节,而常规时深转换方法精度较低,一定程度上影响了勘探井位布置及开发方案的制定。采用时深一体化网格方法,即顶底界面以多井合成记录标定求取的平面时深关系为准,顶底之间利用时间域和深度域网格一一对应关系进行“物理”搬运的方法,从而减小转换前后的累计误差,使地震体及地震反演体时深转换前后的波形与砂体形态保持一致。利用该方法将喇嘛甸油田北部区块反演结果转换到深度域,精细刻画了断层附近砂体的分布特征,指导了L213等3口水平井的部署,砂体预测准确率达到98%以上,年增油17000t,展现了新的时深转换方法的精度及可靠性。关键词 :高精度;时深转换;地震doi:10.13278/j.cnki.juese.201405306   中图 分类号 :P631.4文献 标志码 :A收 稿 日期 :2014-02-08基金 项目 :国家 “973”计划 项目 (2009CB219307);国家 “十 二 ·五 ”重大专项课题 (2011ZX05010-001)作者 简介 :徐立恒 (1980—),男 ,高级 工程师 ,博士 ,主要从事油气地质与地球物理综合研究 ,E-mail:daqing_daruo@163.com。Study and Application on Seismic Time-DepthConversion with High-PrecisionXu Liheng1,Xian Bo2,Xue Yuying3,Li Yanbin41.Exploration and Development Research Institute,Daqing Oilfield Company Limited,Daqing 163712,Heilongjiang,China2.Research Institute of Exploration and Development,Tarim Oilfield Company,PetroChina,Korla 841000,Xinjiang,China3.Geology Research Center of Geophysical Research Institute,BGP,Zhuozhou 072751,Hebei,China4.Sixth Production Plant of Daqing Oil Field,Daqing 163000,Heilongjiang,ChinaAbstract:Time-depth conversion of seismic bodyis a keypoint in seismic structural interpretationand reservoir prediction.Because of the conventional time-depth conversion is with low-precision,it haspartly influenced the wel-location arrangement and exploration scheme.The time-depth integration gridmethod,i.e.the topand bottom surface horizon time-depth relation was based on multi-wel syntheticsseismogram calibration,is proposed,after this method had carried out,the accumulative conversionerrors can be reduced bythe physical transportation of one to one correspondence on the time domain anddepth domain meshes between the topand bottom,thus the waveforms and sand bodyfeatures of seismicbodyand inversion bodycan be kept consistent.Usingthis method,the inversion results of northernLamadian oilfield were transformed into depth domain,the distribution characters of sand bodynearfault were fine depicted,three horizontal wels,such as L213,were deployed.The results show that theprediction accuracy for the sand body is up to 98%,oil increment is more than 17 000tons per year,andthe accuracyand reliabilityof new time-depth conversion methods are wel showed.Keywords:mesh model;time-depth conversion;seismic时深转换的精度和 速度问题是地震研究非常重要的课题 ,为了提高时深转换精度 ,业界学者们投入了大量的精力 。时深转换是将地震体及地震相关数据 ,如地震反演数据体 ,从时间域信号转变为深度域信号的一个处理步骤 ,是利用地震资料进行构造及储层解释的一个非常关键的环节 。地震体的时深转换对于勘探井位布置 、开发中的储量计算及开发方案的制定具有重大意义 ,不合理的时深转换结果会给油田生产带来极大困难 。时深转换精度主要取决于时深转换速度模型的合理性 ,而速度模型的准确建立与研究地区储层的复杂程度 、所处的勘探开发阶段 、基础资料的多少及采用的速度求取方法息息相关[1-5]。大庆 长垣油田已开发50多年 ,历经4次井网加密 ,平均井距已达到80~200m。目前 ,油 田已进入特高含水开发后期的密井网开采阶段 ,陆相河流-三角洲 沉积微幅构造多 、储层薄 、平面非均质强的特点决定了单纯依靠井点资料进行构造及储层研究存在一定的不足 ,在某种程度上影响了油田开发的效果 。从2007年开 始 ,大庆长垣进行了高密度的三维地震采集 ,利用井震结合手段进行构造及储层研究 ,而地震体的时深转换是其重要的一环 ,成为影响地震在注水开发调整中能否发挥作用的关键因素[6-10]。本次针对大庆长垣 油田密井网区的地质特征及开发实际需求 ,开展地震时深转换方法研究 ,以提高地震体时深转换精度 ,充分发挥地震在储层精细刻画及指导剩余油挖潜中的作用 。1 常规时深转换方法及其局限性1.1 利用 地震叠加速度计算地层平均速度在地震数据处理过程中 ,利用地震勘探的多次覆盖数据求取地下速度信息 ,再利用求取的速度对CMP(common middle point)道 集 进 行NMO(normal moveout)校正 ,使 得CMP道集 的同相轴在各炮检距内都得到校正 ,然后叠加得到叠加剖面 。因此 ,叠加速度既是对CMP道集 进行的最佳叠加 ,又可作为平均速度预测目的层的深度 ,具体求法可参见文献 [11]。一组水平层状介质中某一界面以上介质的平均速度珔v就是地震波垂 直穿过该界面以上各层的总厚度与总传播时间之比 :珔v=∑iHi/∑iti。式中 :Hi为各 层的厚度 ;ti为各层的传播时间 。以平均速度计算的时深转换 公式为∑iHi=珔v∑iti。该方法的优点是可利用地 震的高密度信息求取整个三维空间的速度 ,弥补钻井资料横向信息较少的劣势 ,在勘探初期能够发挥较大的作用 ;缺点是估算的速度值在纵向上的误差较大 ,在应用叠加速度进行时深转换时精度偏低 ,无法满足油田开发中储层精细刻画的需求[12-13]。1.2 VSP速度VSP(垂直 地震剖面 )资料为地表激发 、井中接收的所有波场 ,经特殊处理可将地震剖面和钻井资料结合在一起获得速度信息 。可使用单井VSP资料进 行时深转换 。首先 ,单井VSP资料 与地震资料的时间相位对应 ,通过校正可使二者具有可比性 ;其次 ,VSP资料 与井点合成地震记录的时深关系基本一致 。因此 ,VSP是当 前进行时深转换最可靠的资料 。但VSP测井 的成本高 、测井数量少 、控制点不够 ,而区块平面速度变化又较大 ;因此 ,单纯用少量的VSP资料 进行时深转换达不到开发地震精度要求[14]。1.3 合成记录曲线回归及空间 速度场1.3.1 合成地震记录曲线回归公式该方法一般是利用声波测井曲线合成记录来标定层位 ,从而形成时间和深度的对应关系 ,再通过离散数据拟合出线性转换公式 ,形成整个研究区时深转换标定的速度尺 。该方法的优点是层位界面时深3171 第 5期徐立 恒 ,等 :高精度地震时深转换方法研究及应用转换 前后形态一致 ;但是由于无法保证时深关系在空间上是唯一的固定关系 (拟合公式 ),也就是说 ,所有离散点无法完全落到直线上 ,因此层位转换到深度域后常会出现与实际钻井分层不一致的现象 ,导致时深转换出现误差 。1.3.2 合成地震记录建立空间时深关系与曲线回归公式方法不同 ,该方法是通过多井地震合成记录建立空间的时深关系 ,从而建立空间速度场实现时深转换 ;要求井网密度大 ,从而保证空间速度的稳定性和精度 。该方法适用于密井网区油层组级构造层面的时深转换 ,而不适用于地质沉积单元及油层组内部储层砂体的时深转换 。原因在于 :首先 ,地质沉积单元层面尺度小 ,在地震剖面上无法追踪 ,因此井震标定精度受到影响 ;第二 ,针对油层组内部储层进行时深转换时 ,如微幅构造区域 ,空间速度场精度低 ,时深转换之后地震体及地震反演砂体出现 “变形 ”和 “穿层 ”现象[15]。图 1SⅡ油 层 组顶 (a)、底 (b)界面深 /时值平面分布Fig.1 SⅡoil groupof top(a),bottom(b)interface depth/duration distribution2 一体化网格时深转换方法以上常规方法无法 满足大庆长垣油田开发地震时深转换需求 。因此 ,笔者有针对性地开展了密井网时深转换方法研究 ,形成了高精度的一体化网格时深转换方法 ,从而有利于地震在长垣油田储层精细刻画及挖潜剩余油中更好地发挥作用 。以大庆长垣B区块SⅡ层为例 ,利用实际资料实现时深转换方法研究 。研究区面积20km2,共2 276口井 ,井 网密度达到114口 /km2。2008年完成 高密度三维地震资料采集 ,主频45Hz,频带 宽度8~78Hz,地震 采样频率为1ms×10m×10m。文章所述的一体化网 格时深转换方法包括两部分核心内容 :一是油层组顶底界面利用由多井合成记录标定求取的平面时深关系进行转换 (地震上可连续追踪油组层顶底界面 );二是顶底之间利用时间域和深度域网格一一对应关系进行 “物理 ”搬运 ,使地震波形及地震反演砂体特征在时深转换前后保持一致 。2.1 油层组顶底界面时深转换2.1.1 井震标定及时深关系求取制作全区2 276口井 点的合成记录 ,并与井旁实测地震剖面进行井震标定 ,每口井的合成记录与实测地震剖面的相关性要保证在80%以上 ,在 此基础上求取每口井的时深关系曲线 。利用多井时深关系曲线 ,以地震解释层位为约束 ,通过外部飘移克里金算法空间插值得到深时比值平面分布[16]。图1即为 研究区SⅡ油层组顶底界面SⅡ顶、SⅡ底的时 深关系值 :二者分布相对平滑 ,没有异常值 ,且SⅡ底深 /时值 要大于SⅡ顶 。原因是储层底部的压实程度要大于上部 ,速度较大 。SⅡ顶 、底界面时深关系的准确求取 ,为后续地震体的时深转换打下了好的基础 。2.1.2 构造层面解释及时深转换首先 ,在上述井震精细标定的基础上进行时间域SⅡ顶底界面解释追踪 ,生成时间域SⅡ顶、SⅡ底构造 层面 ;第二 ,利用深 /时平面值将解释的时间域4171   吉 林 大 学 学 报 (地 球 科 学 版 )第 44卷图 2SⅡ油 层 组顶 、底层面时间域 (a)与深度域 (b)连井剖面Fig.2 SⅡoil groupof top,bottom level of time(a)and depth(b)domains connected wels profile层位 进行时深转换 ,生成对应深度域的SⅡ顶、SⅡ底构造 层面 。图2为SⅡ油层组顶 、底层面分别在时间域与深度域的某一连井剖面 。从图2中可以看到 ,时间域和深度域的层面形态一致 ,并且与井点分层符合较好 ,反映了利用井震合成记录进行单一油层组顶 、底层面时深转换的精度 。图 3时 间 域 (a)和深度域 (b)网格模型Fig.3 Mesh model of time domain(a)and depth domain(b)2.2 一体 化网格模型的建立利用上述方法分别生成了时间域和深度域的SⅡ顶、SⅡ底层面 模型 ,在此基础上 ,以SⅡ顶、SⅡ底为界 ,分 别建立一一对应关系的时间域和深度域网格模型 ,从而实现油层组内部地震体及地震反演砂体的转换 。首先 ,利用时间域的SⅡ顶、SⅡ底分别 建立时间域框 架 模 型 ,在 时 间 域 模 型 下 建 立 网 格 ,标 准 为1ms×10m×10m,该值 的大小与实测地震的采样间隔一致 ,并计算框架模型的总网格数N=V(10×10×1)(V为模 型体积 )。其次 ,在深度域SⅡ顶、SⅡ底框架 模型约束下 ,建立与时间域相同网格数N的网 格模型 ,形成与时间域网格在空间上的一一对应关系 。最后 ,利用时间域SⅡ顶、SⅡ底层面 约束控制 ,将地震及地震反演数据体导入时间域网格模型中 ,并根据时深网格一一对应关系进一步导入到深度域网格模型中 ,完成地震体的时深转换处理过程(图3)。从上述的一体化网格时深转换方法流程可知 ,油层组顶 、底界面的生成是利用合成记录标定求取的时深关系进行转换 ,而油层组内部是利用时间域和深度域网格一一对应关系进行 “物理 ”搬运 。其中5171 第 5期徐立 恒 ,等 :高精度地震时深转换方法研究及应用a.时间 域;b.深度域。图 4时深转换前后地震剖面对比分析Fig.4 Seismic section analysis before and after time to depth conversiona.时间 域;b.深度域(空间速度场 );c.深度域(一体化网格 )。图 5空间速度场和一体化网格时深转换方法效果对比剖面Fig.5 Compared profiles of time to depth conversion methods between space velocityand integration grid6171   吉 林 大 学 学 报 (地 球 科 学 版 )第 44卷图 6井震结合储层顶面构造图Fig.6 Structural mapof reservoir topsurface forcombination of wel and seismic涉及 到2个关键环节 :一是时间域网格模型与深度域网格模型 ,二者网格总数必须保证一致 ,这样才能在空间上建立二者一一对应的关系 ;二是时深转换的目的层厚度不宜过大 ,即顶底层位的速度差别不宜过大 。针对大庆长垣来说 ,以单个油层组级别进行时深转换比较合适 ,纵向上网格不会出现拉伸不均 ,从而确保时深转换前后地震及地震相关数据体的一致性 。图 7L213井实 钻轨迹 、测井解释及地震反演剖面图Fig.7 L213wel drilingtrajectory,loggingand seismic inversion profiles3 时深转换精度分析利用一体化网格时深 转换方法 ,将研究区地震体从时间域转换到深度域 ,形成深度域的三维地震体 。图4即为时深转换前后地震对比剖面 ,可见转换前后地震波型特征保持一致 ,确保了深度域进行构造及储层解释的精度 。为了进一步验证一体化网格时深转换方法的精度 ,将该方法与空间速度场模型方法进行对比 。图5即为分别利用上述2种方法将地震反演数据体进行时深转换的效果对比剖面 。在图5中可以看到 ,经一体化网格方法时深转换后 ,地震反演砂体与井匹配程度要高于速度场建模方法 ,时深转化前后反演砂体展布特征保持一致 ,减小了累计误差 ,较好地解决了传统速度场建模方法出现数据失真的问题 。4 开 发 应用实例长垣 油田B区井距达120m,其剩 余油挖潜的重心落实在断层边部 。断层区域砂体井控程度较低 、注采关系不够完善 ,存在剩余油潜力 。首先通过井震结合落实该区域断层位置 、储层顶构造分布 ;其次通过地震反演储层预测方法刻画砂体的分布特征 ,并利用一体化网格时深转换方法将地震反演体转到深度域 ,在深度域内刻画断层附近的砂体发育带 ;最后 ,分析认为断层①和断层②附近为构造高点 ,井 控 程 度 低 ,存 在 剩 余 油 ,因 此 设 计 、部 署 了L213、L211、L221三口 水平井 (图6),以提高油层动用程度 ,实现挖潜断层边部剩余油的目的 。图7列出L213水平 井实钻轨迹及深度域波阻抗地震反演砂体分布剖面 ,同时列出了实测的自然伽马 (GR)、密度 (Den)曲线 。该 井实钻水平段长度292m,钻遇 砂岩厚度277.2m;经与实测曲线对 比 ,显示砂 体 预 测 率 达 到98%。3口 水 平 井 年 产 油17 000t,深度 域反演预测的砂体与实钻轨迹很好地吻合 ,充分体现了地震反演及其时深转换的精度和7171 第 5期徐立 恒 ,等 :高精度地震时深转换方法研究及应用可靠 性 。5 结 论1)地震 叠加速度法 、VSP、空间 速度场等常规时深转换具有局限性 ,无法满足油田开发地震时深转换精度的要求 。一体化网格时深转换方法利用时间域和深度域网格一一对应关系进行 “物理 ”搬运 ,减小了转换前后的累计误差 ,提高了精度 。2)地震体及地震反演体从时间域转到深度域后 ,地震波形特征及反演砂体展布特征时深转换前后保持一致 ,并在油田实际开发中得到应用 ,提高了断层区域水平井的钻遇精度 ,实现了挖潜断层边部剩油的目的 ,展现了良好的应用前景 。参 考 文献 (References):[1]凌云 ,郭 建明 ,郭向宇 ,等 .油藏描述中的井震时深转换技术研究 [J].石油物探 ,2011,50(1):1-13.Ling Yun,Guo Jianming,Guo Xiangyu,et al.Research on Time-Depth Conversion by Wel-to-Seismic in Reservoir Characterization[J].GeophysicalProspecting for Petroleum,2011,50(1):1-13.[2]孟宪 军 ,王延光 ,孙振涛 ,等 .井间地震资料时间域波阻抗反演研究 [J].石油学报 ,2005,26(1):47-54.Meng Xianjun,Wang Yanguang,Sun Zhentao,et al.Research on Impedance Inversion of Cross-HoleSeismic Datas in Time Domain[J].Acta PetroleiSinica,2005,26(1):47-54.[3]尚延 安 .塔里木盆地轮南地区速度场的建立和应用[J].油气地质与采收率,2008,15(4):49-51.Shang Yan’an.The Establishment and Application ofVelocity Field in Lunnan Area of Tarim Basin[J].Oil&Gas Recovery Technology,2008,15(4):49-51.[4]万忠 宏 ,闫玉魁 ,詹世凡 ,等 .属性模型速度建场法[J].石油地球物理勘探,2006,41(3):333-336.Wan Zhonghong,Yan Yukui,Zhan Shifan,et al.Velocity-Field-Built Method of Attributes Models[J].Oil Geophysical Prospecting,2006,41(3):333-336.[5]贾义 蓉 ,贺振华 ,石兰亭 ,等 .变速三维地震速度场的构建与应用 [J].物探化探计算技术 ,2011,33(3):243-247.Jia Yirong,He Zhenhua,Shi Lanting,et al.StructurePapping with Transversely Variable 3D SeismicVelocity Model in Erlian Basin[J].ComputingTechniques  for  Geophysical  and  GeochemicalExploration,2011,33(3):243-247.[6]郝 兰 英 ,郭亚杰 ,李杰 ,等 .地震沉积学在大庆长垣密井网条件下储层精细描述中的初步应用 [J].地学前缘 ,2012,19(3):81-86.Hao Lanying,Guo Yajie,Li Jie,et al.PreliminaryApplication of the Seismic Sedimentology to ReservoirRefined Description Under the Condition of HighDensity Wel-Pattern in DaqingPlacanticline[J].EarthScience Frontiers,2012,19(3):81-86.[7]王琦 ,李 红梅 .开发地震在大庆长垣喇嘛甸油田储层预测中的应用 [J].石油与天然气地质 ,2012,33(3):490-495.Wang Qi,Li Hongmei.Application of DevelopmentSeismology to Reservoir Prediction of LamadianOilfield in Daqing Placanticline[J].Oil & GasGeology,2012,33(3):490-495.[8]史忠 生 ,王天琦 ,王建功 ,等 .地震沉积学在不同地质条件下砂体识别中的应用 :以松辽盆地西斜坡和长垣地区为例 [J].岩性油气藏 ,2012,24(3):5-10.Shi Zhongsheng,Wang Tianqi,Wang Jiangong,et al.Application of Seismic Sedimentology to Sand BodyIdentification Under Different Geological Settings:ACase Study from the Western Slope and Placanticline ofSongliao Basin[J].Lithologic Reservoirs,2012,24(3):5-10.[9]徐立 恒 ,陈显森 ,姜岩 ,等 .不同变差函数变程下随机反演与随机模拟对比分析 [J].物探与化探 ,2012,36(2):224-228.Xu Liheng,Chen Xiansen,Jiang Yan,et al.AComparative Analysis of Stochastic Inversion andStochastic Simulation on Different Ranges ofVariograms[J]. Geophysical and GeochemicalExploration,2012,36(2):224-228.[10]李占 东 ,赵伟 ,李阳 ,等 .开发地震反演可行性研究及应用 :以大庆长垣北部油田为例 [J].石油与天然气地质 ,2011,32(5):797-805.Li Zhandong,Zhao Wei,Li Yang,et al.FeasibilityStudy and Application of Development SeismicInversion:A Case Study from Northern ChangyuanOilfield of Daqing Area[J].Oil & Gas Geology,2011,32(5):797-805.[11]刘军 迎 ,雍学善 ,杨午阳 ,等 .转换波叠加速度分析及动校 正 [J].石 油 地 球 物 理 勘 探 ,2012,47(2):231-235.Liu Junying,Yong Xueshan,Yang Wuyang,et al.Converted-Wave Stacking Velocity Analysis &Moveout  Correction[J]. Oil  GeophysicalProspecting,2012,47(2):231-235.[12]秦宁 ,李 振春 ,杨晓东 ,等 .具共散射点道集与角道8171   吉 林 大 学 学 报 (地 球 科 学 版 )第 44卷集 串 级优化叠前偏移速度分析 [J].吉林大学学报 :地球科学版 ,2013,43(2):623-630.Qin Ning,Li Zhenchun,Yang Xiaodong,et al.Pre-Stack Migration Velocity Analysis Based on SequentOptimization of Common Scattering Point Gathersand Angle Domain Common Imaging Gathers[J].Journal of Jilin University:Earth Science Edition,2013,43(2):623-630.[13]管彦 武 ,管烨 ,高瑞 ,等 .根据地震纵波速度分析青藏高原地壳放射性生热率和地幔热状态 [J].吉林大学学报 :地球科学版 ,2013,42(2):562-568.Guan Yanwu,Guan Ye,Gao Rui,et al.CrustalRadiogenic Heat and Mantle Geothermal Status inTibetan Plateau Based on P-Wave Velocity[J].Journal of Jilin University:Earth Science Edition,2013,42(2):562-568.[14]斯 兴 焱 ,李录明 ,胥良君 ,等 .三维 VSP多波 速度分析方法及应用 [J].石油物探 ,2012,51(1):56-64.Si Xingyan,Li Luming,Xu Liangjun,et al.3D-VSPMultiple Velocity Analysis Method and ItsApplication[J]. Geophysical Prospecting forPetroleum,2012,51(1):56-64.[15]赵 镨 ,杨 艳珍 ,艾劲松 .煤炭三维地震勘探时深转换误差分析及对策 [J].中国煤田地质 ,2006,18(4):49-52.Zhao Pu,Yang Yanzhen,Ai Jinsong.Time-DepthConversion Error Analysis and Countermeasures inCoalfield 3DSeismic Prospecting[J].Coal Geology ofChina,2006,18(4):49-52.[16]刘 文 岭 ,朱庆荣 ,戴晓峰 .具有外部漂移的克里金方法在绘制构造图中的应用 [J].石油物探 ,2004,43(4):404-407.Liu Wenling,Zhu Qingrong,Dai Xiaofeng.Application of Kriging-with-an-External-Drift inStructure Mapping[J].Geophysical Prospecting forPetroleum,2004,43(4):404-407.9171 第 5期徐立 恒 ,等 :高精度地震时深转换方法研究及应用
展开阅读全文
  石油文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

关于本文
本文标题:高精度地震时深转换方法研究及应用_徐立恒_图文
链接地址:http://www.oilwenku.com/p-63651.html
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2016-2020 石油文库网站版权所有
经营许可证编号:川B2-20120048,ICP备案号:蜀ICP备11026253号-10号
收起
展开