【摘要】 本发明是利用垂直地震剖面数据对各向异性介质成像的方法，先建立初 始的深度-速度模型，再利用Walkaway VSP或三维VSP初至时间反演模型 的层速度和各向异性系数建立初始速度模型，将三维空间内任意法线方向的 平面转换为法线平行于Z轴的平面，加载观测系统信息，确定从炮点到接收 点射线所穿过的反射层，求取近似射线参数，确定反射点在反射界面上的坐 标并记录相对应的接收点序号，确定反射点到接收点的旅行时间及射线在反 射点的出射角，输出成像道集数据。本发明计算简化，运算效率高，提高了 VSP成像结果的精度与分辨率，可方便地对Walkaway VSP及三维VSP纵波、 横波或转换波数据进行偏移成像。 【专利类型】发明申请 【申请人】中国石油天然气集团公司; 中国石油集团东方地球物理勘探有限责任公司 【申请人类型】企业 【申请人地址】100011北京市东城区安德路16号洲际大厦519室 【申请人地区】中国 【申请人城市】北京市 【申请人区县】东城区 【申请号】CN200810116707.3 【申请日】2008-07-16 【申请年份】2008 【公开公告号】CN101630014A 【公开公告日】2010-01-20 【公开公告年份】2010 【授权公告号】CN101630014B 【授权公告日】2011-10-05 【授权公告年份】2011.0 【发明人】张文波; 侯爱源; 耿伟峰; 王文闯; 吴恒 【主权项内容】1、一种利用垂直地震剖面数据对各向异性介质成像的方法，其特征在于通过以下具体步骤实现： 1)采集VSP叠前地震记录； 2)波场分离，将VSP叠前地震记录分离成不同类型的波场； 3)用零井源距VSP初至时间求得VSP层速度，建立初始的速度-深度模型，再利用Walkaway?VSP或三维VSP的多条初至时间-深度时距曲线反演垂直对称轴各向异性(VTI)模型的层速度和各向异性系数； 4)用每个反射界面的深度、每一层的纵波速度、纵波各向异性系数、横波速度、横波各向异性系数，反射界平面的与Z轴的夹角 反射界面的与X轴的夹角θ建立初始速度模型；根据坐标变换公式将三维空间内任意法线方向的平面转换为法线平行于Z轴的平面； 5)加载要偏移的地震数据并利用地震数据道头加载观测系统信息； 6)对于每一个炮点、每一个接收点和模型的每一个反射层，根据炮点、接收点和当前反射层在模型内的位置，确定从炮点到接收点射线所穿过的反射层； 7)用射线参数P(θ)的函数 利用二分法求取近似射线参数P(θ)， 上式中，k为模型层序号，Hk为第k层的厚度，Vk为第k层的垂向层速度，Ak为第k层的速度各向异性系数，L为炮点到接收井的水平距离； 8)做当前炮点到各个接收点的射线追踪，利用公式 确定反射点在第k个反射界面上的坐标，并记录相应的接收点序号，上式中，k为模型层序号，Hk为第k层的厚度，Vk为第k层的垂向层速度，Ak为第k层 的速度各向异性系数，P(θ1)为第1层的射线参数； 9)根据给定的成像面元大小，把第k个反射界面剖分成一系列等间隔的点，将这些点当作界面上的反射点，计算出每一个反射点的坐标，查找邻近的通过射线追踪计算出的反射点，并记录相对应的接收点序号； 10)做从当前炮点到第k个反射界上等间隔的反射点的射线追踪； 11)做从第k个反射界上等间隔的反射点到与之相对应的接收点的射线追踪； 12)根据步骤10)和步骤11)计算得到炮点到反射点的旅行时间与反射点到接收点的旅行时间之和，将加载的地震数据相应接收点上相应时间上的幅值乘以与射线传播方向相关的系数置于相应的成像点上； 13)做从第k个反射界上等间隔的反射点到与井径上各个接收点的射线追踪，采用与步骤11)相同的方法确定反射点到接收点的旅行时间及射线在反射点的出射角，如果步骤13)计算得到的射线出射角与步骤11)计算得到的射线出射角之差小于给定的偏移孔径角度，则按照步骤12)抽取地震数据进行成像； 14)根据给定的成像网格大小，将成像范围内的数据选排成道集，为各道成像数据置道头信息，输出成像道集数据。 【当前权利人】中国石油天然气集团公司; 中国石油集团东方地球物理勘探有限责任公司 【当前专利权人地址】北京市东城区安德路16号洲际大厦519室; 河北省保定市涿州市范阳西路189号 【专利权人类型】国有企业; 有限责任公司(非自然人投资或控股的法人独资) 【统一社会信用代码】91110000100010433L; 91130000710926983Y 【被引证次数】36 【被自引次数】15.0 【家族被引证次数】40
