确定三维地震观测系统面元属性均匀程度及均匀化的方法
【技术领域】
[0001]本发明涉及三维地震勘探技术领域，特别涉及一种确定三维地震观测系统面元属性均匀程度及均匀化的方法。
【背景技术】
[0002]在石油工业中，三维地震勘探所起的作用日益明显，有关于三维地震资料采集观测系统方面的研究也逐渐增多，其中，对三维地震观测系统的属性:覆盖次数、炮检距和方位角的分析和对该观测系统的评价得到了越来越多的关注。由于一个面元内通常只有一个覆盖次数，但通常存在多个炮检距、方位角，依照目前存在的方法，能够定量地统计一个面元内炮检距和方位角，但是不能够对一个面元内炮检距和方位角的分布进行定量直观的评价，这样就不能合理的确定三维地震观测系统面元属性均匀程度。目前，确定三维地震观测系统面元属性均匀程度的方法，一般是通过设计人员的观察和实际经验来确定，还没有一个能够较合理地确定三维地震观测系统面元属性均匀程度的方法。
[0003]进一步地，随着岩性油气藏和隐蔽油气藏勘探对地震成像分辨率要求的提高，规则的三维地震观测系统、规则的空间采样同样会产生“采集脚印”，也会影响地质目标的高精度、高质量地震成像。由于地震采集过程中，三维地震观测系统中炮线距、接收线距以及排列片的滚动距的不同，使得覆盖次数、炮检距、方位角分布不均匀，并且呈周期性变化，从而导致CDP面元或CRP面元叠加、偏移的振幅、相位出现不均匀性，在地震叠加或偏移数据体上出现周期性的振幅变化，这样的“采集脚印”对储层预测有一定的影响，对油藏评价和储层预测产生影响。目前，也没有对三维地震观测系统面元属性均匀化的方法。
[0004]综上所述，目前，三维地震观测系统评价和应用方面存在的主要问题是:1、没有一个合理的确定三维地震观测系统面元属性均匀程度的方法；2、没有很好的对三维地震观测系统面元属性均匀化，来压制三维地震观测系统引起“采集脚印”问题的方法。
【发明内容】
[0005]本发明实施例提供了一种确定三维地震观测系统面元属性均匀程度的方法，用以确定三维地震观测系统面元属性的均匀程度，该方法包括:
[0006]根据三维地震观测系统的面元大小，将三维地震观测系统的目标区域均匀分成多个面元；
[0007]计算目标区域内每个面元的覆盖次数、炮检距范围大小和方位角范围大小；
[0008]根据目标区域内每个面元的覆盖次数，计算目标区域内所有面元的覆盖次数的均值、方差和方差与均值的比值；根据目标区域内每个面元的炮检距范围大小，计算目标区域内所有面元的炮检距范围大小的均值、方差和方差与均值的比值；根据目标区域内每个面元的方位角范围大小，计算目标区域内所有面元的方位角范围大小的均值、方差和方差与均值的比值；
[0009]根据目标区域内所有面元的覆盖次数的均值、方差和方差与均值的比值，炮检距范围大小的均值、方差和方差与均值的比值，方位角范围大小的均值、方差和方差与均值的比值，确定三维地震观测系统面元属性均匀程度。
[0010]在一个实施例中，计算目标区域内每个面元的炮检距范围大小，包括:
[0011]根据目标区域内每个面元的覆盖次数，将每个面元的炮检距按Cl1= Xmax/N等间距分成N等份炮检距段；所述Xmax为最大炮检距，(I1为每段炮检距段的大小，N为每个面元的覆盖次数，其为大于等于I的正整数；
[0012]根据每个面元对应的炮检距在N等份炮检距段内的分布情况，确定每个面元的炮检距范围大小。
[0013]在一个实施例中，计算目标区域内每个面元的方位角范围大小，包括:
[0014]根据目标区域内每个面元的覆盖次数，将每个面元的方位角按d2= 360° /N等间距分成N等份方位角段；所述d2为每段方位角段的大小，N为每个面元的覆盖次数，其为大于等于I的正整数；
[0015]根据每个面元对应的方位角在N等份方位角段内的分布情况，确定每个面元的方位角范围大小。
[0016]基于同一发明构思，本发明实施例还提供了一种三维地震观测系统面元属性均匀化的方法，用以减小面元之间的属性差别，使得面元属性分布均匀，该方法包括:
[0017]利用上述确定三维地震观测系统面元属性均匀程度的方法，确定多种三维地震观测系统面元属性均匀程度；
[0018]对面元属性均匀程度最优的三维地震观测系统，根据三维地震观测系统的目标区域内每个面元的覆盖次数、炮检距范围大小和方位角范围大小，对三维地震观测系统的面元属性进行均匀化。
[0019]在一个实施例中，对面元属性均匀程度最优的三维地震观测系统，根据三维地震观测系统的目标区域内每个面元的覆盖次数、炮检距范围大小和方位角范围大小，对三维地震观测系统的面元属性进行均匀化，包括:
[0020]确定目标区域内所有面元的覆盖次数的限定范围为&?N2，其中，N1为目标区域内面元的最小限定覆盖次数，N2为目标区域内面元的最大限定覆盖次数，NJPN2S大于等于I的正整数；
[0021]对目标区域内面元覆盖次数大于队的面元按炮检距均匀分布的原则进行抽道集处理。
[0022]在一个实施例中，对目标区域内面元覆盖次数大于N2的面元按炮检距均匀分布的原则进行抽道集处理，包括:
[0023]将覆盖次数大于N2的面元的炮检距按d3= Xmax/N2等间距分成N2等份炮检距段；所述Xmax为最大炮检距，(13为每段炮检距段的大小；
[0024]若一个面元对应道集的炮检距分布在仏个炮检距段内，所述M1为大于等于I的正整数；当1=1时，每个炮检距段内抽一道，并确保所抽道集的炮检距与所在炮检距段的中心值最近。
[0025]在一个实施例中，对目标区域内面元覆盖次数大于N2的面元按炮检距均匀分布的原则进行抽道集处理，还包括:
[0026]当，确定每个炮检距段对应的道集的炮检距分布情况；
[0027]根据每个炮检距段对应的道集的炮检距分布情况，从道集分布密集的炮检距段内抽取数量多的道集；
[0028]从道集分布稀少的炮检距段内抽取数量少的道集，直至总抽取道数为N2道。
[0029]在一个实施例中，对面元属性均匀程度最优的三维地震观测系统，根据三维地震观测系统的目标区域内每个面元的覆盖次数、炮检距范围大小和方位角范围大小，对三维地震观测系统的面元属性进行均匀化，包括:
[0030]确定目标区域内所有面元的覆盖次数的限定范围为&?N2，其中，N1为目标区域内面元的最小限定覆盖次数，N2为目标区域内面元的最大限定覆盖次数，NJPN2S大于等于I的正整数；
[0031]对目标区域内面元覆盖次数大于队的面元按方位角均匀分布的原则进行抽道集处理。
[0032]在一个实施例中，对目标区域内面元覆盖次数大于N2的面元按方位角均匀分布的原则进行抽道集处理，包括:
[0033]将覆盖次数大于N2的面元的方位角按d 4= 360° /N 2等间距分成N 2等份方位角段；所述山为每段方位角段的大小；
[0034]若一个面元对应道集的方位角分布在％个方位角段内，所述M2为大于等于I的正整数；当％=1时，每个方位角段内抽一道，并确保所抽道集的方位角与所在方位角段的中心值最近。
[0035]在一个实施例中，对目标区域内面元覆盖次数大于N2的面元按方位角均匀分布的原则进行抽道集处理，还包括:
[0036]当112〈队时，确定每个方位角段对应的道集的方位角分布情况；
[0037]根据每个方位角段对应的道集的方位角分布情况，从道集分布密集的方位角段内抽取数量多的道集；
[0038]从道集分布稀少的方位角段内抽取数量少的道集，直至总抽取道数为队道。
[0039]与现有技术中仅通过设计人员的观察和实际经验来确定三维地震观测系统面元属性均匀程度的方法相比较，本发明实施例，首先通过对三维地震观测系统的属性进行量化，即计算目标区域内每个面元的覆盖次数、炮检距范围大小和方位角范围大小，这样将面元的覆盖次数、炮检距范围和方位角范围分布情况量化，易于后续对面元的属性进行评价；其次，根据目标区域内每个面元的覆盖次数、炮检距范围大小和方位角范围大小，分别计算它们的均值、方差和方差与均值的比值；最后，分别根据目标区域内所有面元的覆盖次数、炮检距范围大小和方位角范围大小的均值、方差和方差与均值，来确定三维地震观测系统面元属性均匀程度，从而确定合适的三维地震观测系统，这样利用均值、方差值、方差与均值的比值，对面元属性进行统计和分析，能够全面地对三维观测系统进行整体评价，因此，利用本发明提供的技术方案确定的三维地震观测系统，来进行地震勘探采集数据等，利于后续油藏评价、储层预测和地质目标的高精度、高质量地震成像。
[0040]进一步地，本发明实施例