基于地质模型的全方位观测系统参数论证方法
【技术领域】
[0001] 本发明涉及地震勘探观测系统参数论证，特别是涉及到一种基于地质模型的全方 位观测系统参数论证方法。
【背景技术】
[0002] 以往观测系统参数论证只是在全工区选择2-3个参数论证点进行论证，所使用的 方法也比较简单，都是基于各向同性介质的假设下计算的，而在西部非常复杂的山前带区 域，由于地质构造已经不满足各向同性介质假设的条件，因此这种方法已经不能满足技术 设计的需要，首先2-3个论证点不能满足需要，所计算出来的结果具有偏颇性，另外计算方 法不合适，计算结果不满足实际生产需要。为此我们发明了一种新的基于地质模型的全方 位观测系统参数论证方法，解决了以上技术问题。
【发明内容】
[0003] 本发明的目的是提供一种以地质模型为基础，在模型的基础上计算论证每个物理 点采集参数需求的基于地质模型的全方位观测系统参数论证方法。
[0004] 本发明的目的可通过如下技术措施来实现：基于地质模型的全方位观测系统参数 论证方法，该基于地质模型的全方位观测系统参数论证方法包括：步骤1，采集多种地质数 据，建立工区地质模型；步骤2,选择地质模型中的一个层位，采用公式计算层位上每个物 理点的面元的需求；步骤3,将每个物理点的面元需求显示在模型中对应的层面上；步骤4， 选择地质模型中的一个层位，采用公式计算层位上每个物理点的最大炮检距的需求；步骤 5,将每个物理点的最大炮检距需求显示在模型中对应的层面上；以及步骤6,确定满足生 产需要的观测系统参数。
[0005] 本发明的目的还可通过如下技术措施来实现：
[0006] 在步骤1中，采集的地质数据包括层位数据、速度资料、密度资料、深度资料、断层 资料、频率资料这些资料。
[0007] 在步骤1中，层位数据为txt格式，包含各个层位的每个物理点的X横坐标、Y纵坐 标、Z深度信息，各个层位中每个离散物理点的密度小于该区地震勘探所需成像面元大小， 速度资料、频率资料、密度资料、断层资料与层位数据是对应的。
[0008] 在步骤1中，采用层位数据建立工区地质模型，并将各层位对应的速度资料、密度 资料、断层资料导入，每个物理点赋予其地质属性。
[0009]，在步骤2中，在采用公式计算层位上每个物理点的面元的需求时，采用两种公式 分别进行计算，分别是满足最高无混叠频率法则和满足横向分辨率法则公式；
[0010] 满足最高无混叠频率法则的公式为：b=Vint/(4XFmaxXSin0 )
[0011] 式中：b-面兀大小，Vint-上一层层速度，Fmax-最商无混叠频率，Q地层 倾角；
[0012] 满足横向分辨率要求，依据保证良好横向分辨率面元大小的经验公式为：
[0013] b=Vint/(2XFdom)
[0014] 式中：b-面元大小，Vint-上一层层速度，Fdom-反射层视主频，采用这两种公 式，分别进行计算。
[0015] 在步骤3中，原有地质模型不变,将每个物理点的面元需求叠合显示在层面上,采 用颜色代表面元需求的大小，根据步骤2所采用的两种公式计算出的两种面元大小，分别 显示这两种面元大小。
[0016] 在步骤4中，在采用公式计算层位上每个物理点的最大炮检距的需求时，采用两 种公式分别进行计算，分别是满足动校正拉伸需要法则和满足速度分析精度法则公式；
[0017] 满足动校正拉伸需要法则的公式为：
[0019]式中：
[0020] Vrms--均方根速度；
[0021]Xnax-最大炮检距；
[0022] t〇--双程旅行时；
[0023]D--拉伸系数；
[0024] 满足速度分析精度法则的公式为：
[0026] 式中：Av-待测速度增量，
[0027]Xmax-最大炮检距
[0028]At--速度分析中可以检测出的最小正常时差值
[0029]t〇--双程旅行时；
[0030] fp反射波主频；
[0031]vRMS--均方根速度。
[0032] 在步骤5中，原有地质模型不变，将每个物理点的最大炮检距需求叠合显示在层 面上，采用颜色代表最大炮检距需求的大小，根据步骤4所采用的两种公式计算出的两种 最大炮检距大小，分别显示这两种最大炮检距大小。
[0033] 在步骤6中，通过对比分析各个层位的参数需求大小，并结合施工成本需求，合理 选择满足全区生产需要的观测系统参数。
[0034] 本发明中的基于地质模型的全方位观测系统参数论证方法，通过收集工区的以往 地质资料，建立精细的地质模型，在该模型的基础上采用公式计算每个物理点观测系统参 数需求，然后将参数需求叠合显示在地质模型上，这种显示方式比较直观的显示出全区的 参数需求，因此更加容易确定出了适合于该层位地震采集需求的参数，可以较为直观、全面 的分析全区的参数需求，更加科学合理。该发明采用地质模型来计算每个物理点的参数需 求，克服了以往参数论证方法中，只选择几个（一般为3-4个）物理点进行论证，无法满足 全区所有物理点参数需求的矛盾，论证结果更加全面，可以弥补以往方法中论证点少的不 足，和弥补以往参数论证方法无法了解全区参数需求的缺点，该方法适合于所有二维、三维 地震勘探的观测系统参数论证。
【附图说明】
[0035] 图1为本发明的基于地质模型的全方位观测系统参数论证方法的一具体实施例 的流程图；
[0036] 图2是本发明的一具体实施例中某一层位地质模型；
[0037] 图3是本发明的一具体实施例中某一层位中考虑动校拉伸的最大炮检距分析展 示图；
[0038] 图4是本发明的一具体实施例中某一层位中考虑速度分析精度的最大炮检距分 析展示图。
【具体实施方式】
[0039] 为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施 例，并配合所附图式，作详细说明如下。
[0040] 如图1所示，图1为本发明的基于地质模型的全方位观测系统参数论证方法的流 程图。
[0041] 在步骤101，采集多种地质数据，建立精确的工区地质模型。采集的地质数据包括 层位数据、速度、密度、深度、断层、频率等资料，采用层位数据建立工区地质模型，并将各 层位对应的速度、密度、断层信息导入，每个物理点赋予其地质属性。如图2所示，图2是本 发明的一具体实施例中某一层位地质模型。模型是计算的基础，模型要求精确才能保证参 数计算的正确性。所述的收集工区资料中，层位数据的特征是txt格式，包含各个层位的每 个物理点的X(横坐标）、Y(纵坐标）、Z(深度）信息，各个层位中每个离散物理点的密度应 该小于该区地震勘探所需成像面元大小，速度资料、频率资料、密度资料、断层资料应该与 层位数据是对应的。所述的建立模型中，各个层位中的每个物理点需要赋予对应的速度、密 度，在后续的公式计算中需要调用这些参数进行计算，为使模型建立的精细，需要导入断层 信息。流程进入到步骤102。
[0042] 在步骤102,选择其中的一个层位，采用公式计算层位上每个物理点的面元的需 求，采用两种公式分别进行计算，分别是满足最高无混叠频率法则和满足横向分辨率法则 公式。所述的计算每个物理点的参数需求中，面元的选择主要是依据两种法则，
[0043] 第一是满足最高无混叠频率需求，满足偏移成像时不产生偏移噪音，即满足最高 无混叠频率法则，依据公式：b=Vint/(4XFmaxXSin0 )
[0044] 式中：b-面兀大小，Vint-上一层层速度，Fmax-最商无混叠频率，9 地层 倾角。
[0045] 根据以上公式，计算各层位上每个物理点的面元需求。
[0046] 第二是满足横向分辨率要求，依据保证良好横向分辨率面元大小的经验公式：
[0047]b=Vint/