一种面向目的层的vsp观测系统设计方法
【专利摘要】本发明提供了一种面向目的层的VSP观测系统设计方法，属于地震数据采集领域。本发明方法包括：(1)根据勘探任务和目的层深度确定检波器的下放深度和间距；(2)收集工区的现有资料，建立数学模型；(3)根据勘探任务对CRP面元的属性的要求，进行反向射线追踪计算获得理想炮点位置密度分布图；所述CRP面元的属性包括叠加次数、方位角均匀性和偏移距；(4)根据所述理想炮点位置密度分布图确定炮点排列方式和炮点密度；(5)根据步骤(4)确定的炮点排列方式和炮点密度和步骤(1)确定的检波器的下放深度和间距进行射线追踪计算，求取基于CRP的叠加次数统计和方位角统计。
【专利说明】—种面向目的层的VSP观测系统设计方法
【技术领域】
[0001]本发明属于地震数据采集领域，具体涉及一种面向目的层的VSP观测系统设计方法。
【背景技术】
[0002]观测系统设计是VSP野外地震数据采集的施工依据。但长期以来，由于技术的局限，常规的观测系统设计都是假设水平地表和水平介质。随着勘探程度的不断深入，地表平坦和目的层为水平的勘探工区越来越少，大幅度的起伏地表和目的层复杂构造越来越多。在观测系统设计领域一直沿用的手段仍然是水平地表和水平介质的设计方法。近年来，有的软件增加了对起伏地表和复杂构造观测系统设计的评价，但仍然没有提出解决的办法。所以到目前还没有人提出解决起伏地表和复杂构造引起的目的层面元属性不均匀的问题。
【发明内容】
[0003]本发明的目的在于解决上述现有技术中存在的难题，提供一种面向目的层的VSP观测系统设计方法，针对常规VSP观测系统设计中，由于地表起伏和目的层的复杂性引起的实际目的层叠加次数的不均匀性，通过有目的的补炮实现目的层叠加次数的均匀。而这种补炮的位置和密度依据的是通过速度模型的射线反向模拟来获得的，具有更好的针对性和准确性。
[0004]本发明是通过以下技术方案实现的:
[0005]一种面向目的层的VSP观测系统设计方法，所述方法包括以下步骤:
[0006](I)根据勘探任务和目的层深度确定检波器的下放深度和间距；
[0007](2)收集工区的现有资料，建立数学模型；
[0008](3)根据勘探任务对CRP面元的属性的要求，进行反向射线追踪计算获得理想炮点位置密度分布图；所述CRP面元的属性包括叠加次数、方位角均匀性和偏移距；
[0009](4)根据所述理想炮点位置密度分布图确定炮点排列方式和炮点密度；
[0010](5)根据步骤(4)确定的炮点排列方式和炮点密度和步骤(I)确定的检波器的下放深度和间距进行射线追踪计算，求取基于CRP的叠加次数统计和方位角统计；
[0011](6).判断不达标的面元数量是否小于总面元数量的5%，如果是，则转入步骤(8)，如果否，则转入步骤(7);不达标的面元是指叠加次数小于设计次数的面元或方位角分布范围小于设计限值的面元；
[0012](7).对不达标的面元，重新进行反向射线追踪计算获得新的理想炮点位置密度分布图，并以此进行加密补炮设计，然后进行射线追踪计算，求取基于CRP的叠加次数统计和方位角统计，再返回步骤(6)；
[0013](8).结束。
[0014]所述步骤(3)中的所述理想炮点位置密度分布图是这样生成的:
[0015]对于所有的接收点和所有目的层面元都进行反向射线追踪计算求出理想炮点位置，然后把所有的理想炮点位置投射到平面上得到理想炮点位置密度分布图。
[0016]所述步骤(4)是这样实现的:
[0017]在所述理想炮点位置密度分布图中，在密度较大的区域布设更多的炮点，在密度较小的区域布设相对较少的炮点。
[0018]所述步骤(5)中的求取基于CRP的叠加次数统计和方位角统计是这样实现的:统计落入CRP面元的反射点数得到叠加次数统计，根据目的层反射点射线的投影方向确定方位角得到方位角统计。
[0019]所述步骤(7)中进行加密补炮设计是这样实现的:对炮点密度大的区域设计加密炮点，密度越大的区域，设计越多的炮点；加密炮点总数应小于等于N, N =满叠加次数-叠加次数最少的不达标面元的叠加次数；所述满叠加次数是指为完成勘探任务，根据工区具体情况和VSP施工经验确定的有效覆盖目标区的叠加次数，所述叠加次数最少的不达标面元的叠加次数是从步骤(5)中的叠加次数统计中得到的。
[0020]与现有技术相比，本发明的有益效果是:本发明所设计的VSP观测系统，由于是针对目的层的CRP面元属性做了针对性的优化，使得采集得到的VSP单炮记录更真实的反应了目的层的特性。这为后续提高VSP资料处理(包括波场分离、道集叠加和数据成像)质量创造了条件。
【专利附图】
【附图说明】
[0021]图1是复杂构造对VSP传统观测系统的影响示意图。
[0022]图2是本发明观测系统的示意图。
[0023]图3是本发明实施例中的炮点分布密度平面图。
[0024]图4是本发明实施例中的炮线布设示意图。
[0025]图5 (a)是使用传统观测系统接收的VSP资料的成像结果。
[0026]图5(b)是使用本发明设计的观测系统接收的VSP资料的成像结果。
[0027]图6是本发明面向目的层的VSP观测系统设计方法的步骤框图。
【具体实施方式】
[0028]下面结合附图对本发明作进一步详细描述:
[0029]本发明针对VSP采集中接收排列相对固定的特点，用反向射线追踪算法，求取通过目的层CRP面元中心的理想炮点位置，在炮点分布密度平面图的基础上设计炮点排列，从而优化了目的层CRP面元的属性，增强了 CRP面元属性的均匀性。该方法可在普通微机上安装使用，占用计算机资源少、计算量小，方法简单实用。
[0030]复杂构造对VSP观测带来的影响如图1所示，R1、R2为井中检波器接收点，把目的层视为水平层时，经过CRP点的2个叠加道集的射线路径为SI，-CRP-Rl和S2’ -CRP-R2，然而由于目的层的构造起伏，实际的射线轨迹变为SI’ -Pl-Rl和S2’ -P2-R2，显然Pl和P2的位置已经远离CRP点的位置，它们已经不是反映同一个面元的信息。如果把这2道数据进行叠加，就得不到好的结果。图1中SI和S2发出的射线(S1-CRP-R1和S2-CRP-R2)，经过CRP点反射后正好到达Rl和R2点，这2条射线的数据叠加才能增强CRP点的有效信息。
[0031]目的层的构造越复杂，P1、P2点位置偏离CRP点就越远，投射到地面SI’ -SI和S2’ -S2的距离也越远。VSP观测由于检波器放在井中接近目的层的位置，这种偏离越加敏感。而VSP的目的就是查清井旁目的层的精细构造和岩性，对这种偏离是不能容忍的。
[0032]本发明面向目的层的VSP观测系统设计方法的步骤如图6所示，包括:
[0033](I).用常规方法确定检波器的下放深度和间距:使用常规观测系统设计软件，根据勘探任务和目的层深度确定下放深度和间距；
[0034](2).收集工区的现有资料(包括钻进、声波、密度测井资料，已有的地震、地质资料)，建立数学模型；
[0035](3).根据勘探任务对CRP面元(又称为反射面元)的属性(叠加次数、方位角均匀性和偏移距)的要求，进行反向射线追踪计算获得理想炮点位置密度分布图；
[0036](4).根据理想炮点位置密度分布图确定炮点排列方式和炮点密度；
[0037](5).根据步骤(4)确定的炮点排列方式和炮点密度和步骤(1)确定的检波器的下放深度和间距进行射线追踪计算，求取基于CRP的叠加次数统计(即统计落入CRP面元的反射点数)和方位角统计(根据目的层反射点射线的投影方向确定)；
[0038](6).判断不达标的面元数量是否小于总面元数量的5%，如果是，则转入步骤
(8)，如果否，则转入步骤(7);不达标的面元是指叠加次数小于设计次数的面元和方位角分布范围小于设计限值的面元；
[0039](7).对叠加次数和方位角不达标的面元，重新进行反向射线追踪计算获得新的理想炮点位置密度分布图，并以此进行加密补炮设计，然后进行射线追踪计算，求取基于CRP的叠加次数统计和方位角统计，再 返回步骤(6)；
[0040](8).结束。
[0041]所述步骤(7)中进行加密补炮实际是这样实现的:对炮点密度大的区域设计加密炮点，密度越大的区域，设计越多的炮点；炮点排列形式可以采用矩阵、圆形、扇型等规则排列，也可以用不规则排列；加密炮点总数应小于等于N，N =满叠加次数-叠加次数最少的不达标面元的叠加次数。具体实施时，补炮设计的循环次数一般在3~5次，在第5次以后的补炮时，不再考虑方位角和偏移距，而是仅仅考虑叠加次数的均匀性，这样满覆盖面元的比例很快就会达标。最后，当满覆盖面元达到目标范围内面元总量的95%时，设计过程完成。
[0042]所述步骤(3)中的理想炮点位置密度分布图的生成是这样实现的:
[0043]如图2所示，在目的层不为水平界面的情况下，经过面元中心CRP1、CRP2、CRP3到达接收点 Rl、R2 的射线路径分别为:S1-CRP1-R1、S2-CRP1-R2、S3-CRP2-R1、S4-CRP2-R2、S5-CRP3-R1、S6-CRP3-R2。其中 S1、S2、S3、S4、S5、S6 是在已知 RU R2 和 CRP1、CRP2、CRP3的情况下由反向射线追踪计算获得的理想炮点位置。
[0044]对于所有的接收点Rl……R2和所有目的层面元都可以求出其理想炮点位置，把它们投射到平面上就可以获得炮点分布密度平面图，如图3所示。
[0045]所述步骤⑷中的炮点布设排列的确定如下:
[0046]根据上述炮点密度平面分布图，在密度较大的区域布设更多的炮点，在密度较小的区域布设相对较少的炮点。炮点布设的密度与“炮点密度平面分布图”的密度成线性关系。炮点布设的总数根据勘探任务、检波器节数、出射角范围、目的层深度等，用常规的观测系统设计软件(如绿山观测系统设计软件)计算确定。炮点排列形式可以采用矩阵、圆形、扇型等规则排列，也可以用不规则排列。最后可以确定出相应的炮线布设如图4所示:[0047]这样的炮点布设使得基于目的层的CRP面元属性(叠加次数、方位角和偏移距)均匀性大幅度提高.[0048]步骤(5)使用的是正向(常规)的射线追踪计算方法，技术是公开的。步骤(3)和步骤(7)中使用的是反向射线追踪计算方法，两种方法的公式形式都一样，但反向法在过程进行顺序和参数运用设置上是不一样的。
[0049]具体来说，所述步骤(3)和步骤(7)中所用的反向射线追踪计算是这样实现的:
[0050]射线追踪的理论基础是在高频近似的条件下，地震波场主能量沿射线轨迹附近传播。射线追踪方法主要采用费马原理(Fermat’ s principle)、惠更斯原理(Huygens’ principle)、斯捏尔定律(Snell，s law)和程函方程(Eikonal equation)。
[0051]从Fermat原理出发，利用一阶Taylor不完全展开，可得到如下形式的迭代射线追踪算法(关于所有中间路径点的矩阵方程):
【权利要求】
1.一种面向目的层的VSP观测系统设计方法，其特征在于:所述方法包括以下步骤: (1)根据勘探任务和目的层深度确定检波器的下放深度和间距； (2)收集工区的现有资料，建立数学模型； (3)根据勘探任务对CRP面元的属性的要求，进行反向射线追踪计算获得理想炮点位置密度分布图；所述CRP面元的属性包括叠加次数、方位角均匀性和偏移距； (4)根据所述理想炮点位置密度分布图确定炮点排列方式和炮点密度； (5)根据步骤(4)确定的炮点排列方式和炮点密度和步骤(I)确定的检波器的下放深度和间距进行射线追踪计算，求取基于CRP的叠加次数统计和方位角统计； (6).判断不达标的面元数量是否小于总面元数量的5%，如果是，则转入步骤(8)，如果否，则转入步骤(7);不达标的面元是指叠加次数小于设计次数的面元或方位角分布范围小于设计限值的面元； (7).对不达标的面元，重新进行反向射线追踪计算获得新的理想炮点位置密度分布图，并以此进行加密补炮设计，然后进行射线追踪计算，求取基于CRP的叠加次数统计和方位角统计，再返回步骤(6)； (8).结束。
2.根据权利要求1所述的面向目的层的VSP观测系统设计方法，其特征在于:所述步骤(3)中的所述理想炮点位置密度分布图是这样生成的: 对于所有的接收点和所有目的层面元都进行反向射线追踪计算求出理想炮点位置，然后把所有的理想炮点位置投射到平面上得到理想炮点位置密度分布图。
3.根据权利要求2所述的面向目的层的VSP观测系统设计方法，其特征在于:所述:所述步骤(4)是这样实现的: 在所述理想炮点位置密度分布图中，在密度较大的区域布设更多的炮点，在密度较小的区域布设相对较少的炮点。
4.根据权利要求3所述的面向目的层的VSP观测系统设计方法，其特征在于:所述:所述步骤(5)中的求取基于CRP的叠加次数统计和方位角统计是这样实现的:统计落入CRP面元的反射点数得到叠加次数统计，根据目的层反射点射线的投影方向确定方位角得到方位角统计。
5.根据权利要求4所述的面向目的层的VSP观测系统设计方法，其特征在于:所述:所述步骤(7)中进行加密补炮设计是这样实现的:对炮点密度大的区域设计加密炮点，密度越大的区域，设计越多的炮点；加密炮点总数应小于等于N，N =满叠加次数-叠加次数最少的不达标面元的叠加次数；所述满叠加次数是指为完成勘探任务，根据工区具体情况和VSP施工经验确定的有效覆盖目标区的叠加次数，所述叠加次数最少的不达标面元的叠加次数是从步骤(5)中的叠加次数统计中得到的。
【文档编号】G01V1/00GK103777229SQ201210407728
【公开日】2014年5月7日 申请日期:2012年10月23日 优先权日:2012年10月23日
【发明者】王龙泉 申请人:中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院