一种地质构造成像的方法及装置的制造方法
【技术领域】
[0001] 本发明涉及石油地质勘探技术领域，具体涉及一种地质构造成像的方法及装置。
【背景技术】
[0002] 通常的地震勘探资料处理和成像是基于均匀或水平层状的地震波传播介质得到 的，并要求地表面为水平面，在此基础上引出了速度分析、叠加以及偏移。这样，在资料处理 开始就要定义基准面（地震剖面的起始零线），地震剖面的时间都要以基准面做参考。根据 地形的具体情况，结合物探的技术要求，目前基准面的选择主要有两种方法：水平基准面法 和地表圆滑面法。
[0003] 由于静校正一般都只做垂直方向的时移校正，不考虑波的传播方向。当静校正量 较小时，这种做法的问题不大，但当静校正量较大时，这种做法就会影响静校正的精度和反 射波时距曲线的性质。然而，真正的基准面校正需要按射线路径把炮点、检波点校正到基准 面上，这样才能确保基准面校正后的反射波时距曲线的性质。因此，通常的基准面静校正量 就存在着静校正量的误差。
[0004] 针对上述减小静校正量误差的问题，目前尚未提出有效的解决方案。
【发明内容】
[0005] 本发明目的在于提供一种地质构造成像的方法及装置，利用本发明地质构造成像 的方法及装置可以减小静校正量，进而减小静校正量的误差，提高地质构造成像的精度。
[0006] 本发明的上述目的可采用下列技术方案来实现：
[0007] 本发明提供一种地质构造成像的方法，包括：
[0008] 拾取工区的地震资料和大地测量数据；
[0009] 基于所述大地测量数据，对一次曲面方程进行地表拟合，得到所述一次曲面方程 中的方程系数，将得到的方程系数代入所述一次曲面方程中，得到拟合面；
[0010] 对所述方程系数中影响拟合面高程的方程系数进行调整，使得拟合面位于所述工 区的最大高程以上，将调整后的拟合面作为基准面；
[0011] 将炮点或检波点的坐标数据代入所述基准面中，得到炮点或检波点的基准面高 程，根据所述基准面高程得到静校正量；
[0012] 基于所述静校正量将所述炮点或检波点校正至所述基准面，得到带有炮点或检波 点信息的基准面；
[0013] 将所述带有炮点或检波点信息的基准面转换成水平基准面；
[0014] 根据所述地震资料在所述水平基准面上形成地质构造图。
[0015] 在一个实施方式中，按照以下公式计算得到所述静校正量：
[0017] 其中，T为炮点或检波点的静校正量，单位ms也为炮点或检波点的第i层低降速 层的厚度，单位m;Vi为炮点或检波点的第i层低降速层的速度，单位m/sJ为炮点的井口 时间，单位ms此为炮点或检波点的基准面高程，单位m;Hg为炮点或检波点的地表高程，单 位m;VS为替换速度，单位m/s。
[0018] 在一个实施方式中，将所述带有炮点或检波点信息的基准面转换成水平基准面， 包括：
[0019] 保持所述带有炮点或检波点信息的基准面的位置不变，旋转所述带有炮点或检波 点信息的基准面所在的第一坐标系至第二坐标系，其中，所述带有炮点或检波点信息的基 准面在所述第二坐标系中为水平状态。
[0020] 在一个实施方式中，根据所述地震资料在所述水平基准面上形成地质构造图，包 括：
[0021] 根据所述地震资料，在所述第二坐标系中形成第一地质构造图；
[0022] 将所述第一地质构造图所在的第二坐标系反旋至第一坐标系，得到第二地质构造 图。
[0023] 在一个实施方式中，按照以下公式表示所述一次曲面方程：
[0024] ax+by+cz+d= 0
[0025] 其中，a、b、c、d表示方程系数；X、y、z表示坐标轴。
[0026] 本发明还提供了一种地质构造成像的装置，包括：
[0027] 拾取模块，被配置为拾取工区的地震资料和大地测量数据；
[0028] 拟合模块，被配置为基于所述大地测量数据，对一次曲面方程进行地表拟合，得到 所述一次曲面方程中的方程系数，将得到的方程系数代入所述一次曲面方程中，得到拟合 面；
[0029] 调整模块，被配置为对所述方程系数中影响拟合面高程的方程系数进行调整，使 得拟合面位于所述工区的最大高程以上，将调整后的拟合面作为基准面；
[0030] 代入模块，被配置为将炮点或检波点的数据代入所述基准面中，得到炮点或检波 点的基准面高程，根据所述基准面高程得到静校正量；
[0031] 校正模块，被配置为基于所述静校正量将所述炮点或检波点校正至所述基准面， 得到带有炮点或检波点信息的基准面；
[0032] 转换模块，被配置为将所述带有炮点或检波点信息的基准面转换成水平基准面；
[0033] 成像模块，被配置为根据所述地震资料在所述水平基准面上形成地质构造图。
[0034] 在一个实施方式中，所述代入模块具体被配置为按照以下公式计算得到所述静校 正量：
[0036] 其中，T为炮点或检波点的静校正量，单位ms也为炮点或检波点的第i层低降速 层的厚度，单位m ;Vi为炮点或检波点的第i层低降速层的速度，单位m/sJ为炮点的井口 时间，单位ms此为炮点或检波点的基准面高程，单位m;Hg为炮点或检波点的地表高程，单 位m;VS为替换速度，单位m/s。
[0037] 在一个实施方式中，所述转换模块包括：
[0038] 坐标系转换单元，被配置为保持所述带有炮点或检波点信息的基准面的位置不 变，旋转所述带有炮点或检波点信息的基准面所在的第一坐标系至第二坐标系，其中，所述 带有炮点或检波点信息的基准面在所述第二坐标系中为水平状态。
[0039] 在一个实施方式中，所述成像模块包括：
[0040] 第一成像单元，被配置为根据所述地震资料，在所述第二坐标系中形成第一地质 构造图；
[0041] 第二成像单元，被配置为将所述第一地质构造图所在的第二坐标系反旋至第一坐 标系，得到第二地质构造图。
[0042] 在一个实施方式中，按照以下公式表示所述一次曲面方程：
[0043] ax+by+cz+d = 0
[0044] 其中，a、b、c、d表示方程系数；X、y、z表示坐标轴。
[0045] 综上所述，本发明通过大地测量数据与一次曲面方程进行地表拟合后，调整方程 系数使拟合后的拟合面转换为基准面后，再利用基准面高程得到较小的静校正量，基于该 较小的静校正量将炮点或检波点校正到基准面，得到带有炮点或检波点信息的基准面；将 带有炮点或检波点信息的基准面转换成水平基准面后成像，可有效提高地质勘探成像精 度。
【附图说明】
[0046] 在此描述的附图仅用于解释目的，而不意图以任何方式来限制本发明公开的范 围。另外，图中的各部件的形状和比例尺寸等仅为示意性的，用于帮助对本发明的理解，并 不是具体限定本发明各部件的形状和比例尺寸。本领域的技术人员在本发明的教导下，可 以根据具体情况选择各种可能的形状和比例尺寸来实施本发明。
[0047] 图1示出了基准面静校正量示意图；
[0048] 图2示出了基准面静校正量的误差和地表与基准面的高差以及炮检距之间的关 系曲线图；
[0049] 图3示出了本发明的地质构造成像的方法的流程图；
[0050] 图4示出了本发明在第一坐标系下工区的地形图；
[0051] 图5不出了本发明在第一坐标系下工区地形图的拟合面图；
[0052] 图6示出了本发明在第一坐标系下基准面的静校正图；
[0053] 图7示出了本发明在第二坐标系下基准面的静校正图；
[0054] 图8示出了本发明在第二坐标系下的工区地形图；
[0055] 图9示出了本发明在第一坐标系下的基准面叠后偏移剖面；
[0056] 图10示出了本发明在第二坐标系下水平基准面的叠后偏移剖面图；
[0057]图11示出了本发明地质构造成像的装置的模块图。
【具体实施方式】
[0058] 结合附图和本发明【具体实施方式】的描述，能够更加清楚地了解本发明的细节。但 是，在此描述的本发明的【具体实施方式】，仅用于解释本发明的目的，而不能以任何方式理解 成是对本发明的限制。在本发明的教导下，技术人员可以构想基于本发明的任意可能的变 形，这些都应被视为属于本发明的范围。
[0059] 如图1所示，地表和基准面都是平面，静校正相当于把在地表激发和接收的共中 心点（CMP)各道延迟一个时间Atl，所以基准面静校正量是：
[0061] 其中，Z表示地表和基准面的高差，V表示地表和基准面间的速度；
[0062] 而实际的基准面静校正量是：
[0064] 其中，Z表示地表和基准面的高差，X为炮检距，Η表示地下某一界面的深度，v是 Η深度界面的速度；
[0065] 这样静校正量的误差为：Δt=Δt「Δt2。
[0066] 图2表示基准面静校正量的误差Δt和地表与基准面的高差Z以及炮检距x之间 的关系曲线。其中，横坐标表示炮检距X(也叫做偏移距），纵坐标表示静校正量的误差At。
[0067] 从图中可以看出，随着炮检距x的增加，静校正量的误差Δt也增加；随着地表与 基准面