在地震剖面上显示测井信息的方法及装置与流程
本申请涉及地震解释技术领域，尤其是涉及一种在地震剖面上显示测井信息的方法及装置。
背景技术：
目前传统的测井信息绘制是将测井信息垂直显示(如图2所示，其中直线为井轨迹，曲线为测井信息)，忽略井轨迹的影响，以关注单井的数据变化。当在地震剖面上显示测井信息时，是把每个深度上的数据点进行横向偏移，从而绘制出测井信息沿着井轨迹变化的趋势(如图3所示，其中斜线为井轨迹，曲线为测井信息)。
如图1所示，根据井轨迹类型不同测井分为直井、斜井、水平井。而在实现本申请的过程中，本申请的发明人发现：随着大角度斜井和水平井的出现，进入井的大斜率段或者水平段，一般沿井轨迹延伸，而在垂直的深度上几乎没有变化，这样，绘制出的很多数据点可能会挤在一起，从而产生数据挤压变形和显示交叠等情况。此外，本申请的发明人还发现：目前传统的在井轨迹上水平绘制分层的方案是单纯地沿水平方向绘制分层(如图4所示，其中s1、s2和s3为划分出的各个水平层)，其没有结合井和地震反射特征的走向信息，无法准确确定地质分层的走势。
技术实现要素：
本申请实施例的目的在于提供一种在地震剖面上显示测井信息的方法及装置，以提高测井信息在地震剖面上的显示质量。
为达到上述目的，一方面，本申请实施例提供了一种在地震剖面上显示测井信息的方法，包括：
在给定的投影距离下将三维空间中的测井轨迹投射到地震剖面上，获得测井投射轨迹；
确定所述测井投射轨迹上每个测量点的法向，并按照每个测量点的法向绘制对应的测井数据；以及，
确定所述测井投射轨迹上每个测量点的地层倾向，并按照每个测量点的地层倾向绘制对应的测井解释数据。
本申请实施例的在地震剖面上显示测井信息的方法，所述确定所述测井投射轨迹上每个测量点的法向，并按照每个测量点的法向绘制对应的测井数据，包括：
根据所述测井投射轨迹上每个测量点的位置信息，确定每个测量点对应测井数据的水平偏移量；
根据所述测井投射轨迹上每个拐点的位置信息，将所述测井投射轨迹划分为若干个轨迹段；
确定每个轨迹段的单位化向量，并根据每个轨迹段的单位化向量确定所述测井投射轨迹上每个测量点的单位化法向量；
根据所述测井投射轨迹上每个测量点的单位化法向量及对应测井数据的水平偏移量，确定所述测井投射轨迹上每个测量点对应测井数据的测井数据点显示位置；
按照所述测井数据点显示位置绘制测井数据。
本申请实施例的在地震剖面上显示测井信息的方法，所述确定每个轨迹段的单位化向量，包括：
根据公式
其中，ai为第i个轨迹段的单位化向量的横坐标；bi为第i个轨迹段的单位化向量的纵坐标；xi和xi+1分别为测井投射轨迹上第i个和第i+1个拐点的横坐标；yi和yi+1分别为测井投射轨迹上第i个和第i+1个拐点的纵坐标。
本申请实施例的在地震剖面上显示测井信息的方法，所述根据每个轨迹段的单位化向量确定所述测井投射轨迹上每个测量点的单位化法向量，包括：
对于每个轨迹段，根据公式a＝a1(1-rj)+a2rj和b＝b1(1-rj)+b2rj确定所述测井投射轨迹上每个测量点的法向量；
其中，a和b分别为测井投射轨迹上第j个测量点的法向量的横坐标和纵坐标；a1和a2分别为该轨迹段内第一个拐点和第二个拐点的横坐标；b1和b2分别为该轨迹段内第一个拐点和第二个拐点的纵坐标；rj为该轨迹段内第j个测量点的系数，且rj＝b/a，其中，b为第j个测量点至所述第一个拐点的长度，a为所述第一个拐点与所述第二个拐点之间的长度；
单位化所述测井投射轨迹上每个测量点的法向量，获得所述测井投射轨迹上每个测量点的单位化法向量。
本申请实施例的在地震剖面上显示测井信息的方法，所述根据所述测井投射轨迹上每个测量点的单位化法向量及对应测井数据的水平偏移量，确定所述测井投射轨迹上每个测量点对应测井数据的测井数据点显示位置，包括：
根据公式x＝xgui+a'w和y＝ygui+b'w确定所述测井投射轨迹上每个测量点对应测井数据的测井数据点显示位置；
其中，x和y分别为每个测量点对应测井数据的测井数据点的横坐标和纵坐标；xgui和ygui分别为每个测量点的单位化法向量的横坐标和纵坐标；w为每个测量点对应测井数据的水平偏移量。
另一方面，本申请实施例还提供了一种在地震剖面上显示测井信息的装置，包括：
测井轨迹投射模块，用于在给定的投影距离下将三维空间中的测井轨迹投射到地震剖面上，获得测井投射轨迹；
测井数据绘制模块，用于确定所述测井投射轨迹上每个测量点的法向，并按照每个测量点的法向绘制对应的测井数据；以及，
测井解释数据绘制模块，用于确定所述测井投射轨迹上每个测量点的地层倾向，并按照每个测量点的地层倾向绘制对应的测井解释数据。
本申请实施例的在地震剖面上显示测井信息的装置，所述确定所述测井投射轨迹上每个测量点的法向，并按照每个测量点的法向绘制对应的测井数据，包括：
根据所述测井投射轨迹上每个测量点的位置信息，确定每个测量点对应测井数据的水平偏移量；
根据所述测井投射轨迹上每个拐点的位置信息，将所述测井投射轨迹划分为若干个轨迹段；
确定每个轨迹段的单位化向量，并根据每个轨迹段的单位化向量确定所述测井投射轨迹上每个测量点的单位化法向量；
根据所述测井投射轨迹上每个测量点的单位化法向量及对应测井数据的水平偏移量，确定所述测井投射轨迹上每个测量点对应测井数据的测井数据点显示位置；
按照所述测井数据点显示位置绘制测井数据。
本申请实施例的在地震剖面上显示测井信息的装置，所述确定每个轨迹段的单位化向量，包括：
根据公式
其中，ai为第i个轨迹段的单位化向量的横坐标；bi为第i个轨迹段的单位化向量的纵坐标；xi和xi+1分别为测井投射轨迹上第i个和第i+1个拐点的横坐标；yi和yi+1分别为测井投射轨迹上第i个和第i+1个拐点的纵坐标。
本申请实施例的在地震剖面上显示测井信息的装置，所述根据每个轨迹段的单位化向量确定所述测井投射轨迹上每个测量点的单位化法向量，包括：
对于每个轨迹段，根据公式a＝a1(1-rj)+a2rj和b＝b1(1-rj)+b2rj确定所述测井投射轨迹上每个测量点的法向量；
其中，a和b分别为测井投射轨迹上第j个测量点的法向量的横坐标和纵坐标；a1和a2分别为该轨迹段内第一个拐点和第二个拐点的横坐标；b1和b2分别为该轨迹段内第一个拐点和第二个拐点的纵坐标；rj为该轨迹段内第j个测量点的系数，且rj＝b/a，其中，b为第j个测量点至所述第一个拐点的长度，a为所述第一个拐点与所述第二个拐点之间的长度；
单位化所述测井投射轨迹上每个测量点的法向量，获得所述测井投射轨迹上每个测量点的单位化法向量。
本申请实施例的在地震剖面上显示测井信息的装置，所述根据所述测井投射轨迹上每个测量点的单位化法向量及对应测井数据的水平偏移量，确定所述测井投射轨迹上每个测量点对应测井数据的测井数据点显示位置，包括：
根据公式x＝xgui+a'w和y＝ygui+b'w确定所述测井投射轨迹上每个测量点对应测井数据的测井数据点显示位置；
其中，x和y分别为每个测量点对应测井数据的测井数据点的横坐标和纵坐标；xgui和ygui分别为每个测量点的单位化法向量的横坐标和纵坐标；w为每个测量点对应测井数据的水平偏移量。
再一方面，本申请实施例还一种计算机存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
在给定的投影距离下将三维空间中的测井轨迹投射到地震剖面上，获得测井投射轨迹；
确定所述测井投射轨迹上每个测量点的法向，并按照每个测量点的法向绘制对应的测井数据；以及，
确定所述测井投射轨迹上每个测量点的地层倾向，并按照每个测量点的地层倾向绘制对应的测井解释数据。
由以上本申请实施例提供的技术方案可见，本申请实施例按照测井投射轨迹上每个测量点的法向绘制对应的测井数据，并按照每个测量点的地层倾向绘制对应的测井解释数据，从而避免了现有技术出现的数据挤压变形、显示交叠以及无法准确确定地质分层的走势的问题，从而提高了测井信息在地震剖面上显示的质量。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为现有技术中根据井轨迹类型进行的测井分类示意图；
图2为现有技术中将测井信息垂直显示的示意图；
图3为现有技术中将测井信息在地震剖面上显示的示意图；
图4为现有技术中沿水平方向绘制分层的示意图；
图5为本申请一实施例中在地震剖面上显示测井信息的方法流程图；
图6为本申请一实施例中各轨迹段的法向示意图；
图7为本申请一实施例中沿用地层倾向显示分层的示意图；
图8为本申请一实施例中在地震剖面上显示测井信息的装置的结构框图；
图9为本申请一实施例得到的在地震剖面上显示测井信息的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
参考图5所示，本申请实施例的地震剖面上绘制测井信息的方法可以包括如下步骤：
s101、在给定的投影距离下将三维空间中的测井轨迹投射到地震剖面上，获得测井投射轨迹。
本申请实施例中，在给定的投影距离下，可根据三维空间中的测井轨迹上每个测量点计算其在地震剖面上的投影，从而将三维空间中的测井轨迹投射到地震剖面上，获得测井投射轨迹。
s102、确定所述测井投射轨迹上每个测量点的法向，并按照每个测量点的法向绘制对应的测井数据。
本申请实施例中，为了确定可能出现在任意测量深度位置测井数据的显示法向量，首先需要计算出测井投射轨迹上各个拐点的法向量。具体的，
首先，可根据所述测井投射轨迹上每个测量点的位置信息，确定每个测量点对应测井数据的水平偏移量。当将三维空间中的测井轨迹投射到地震剖面上后，测井投射轨迹形成，则测井投射轨迹上每个测量点的位置也很容易就确定了。根据每个测量点的位置信息可以计算对应测井数据相对于井口位置的水平偏移量。
其次，根据所述测井投射轨迹上每个拐点的位置信息，将所述测井投射轨迹划分为若干个轨迹段。即以各个拐点作为划分点，从而可以将整个测井投射轨迹划分为若干个轨迹段。
然后，确定每个轨迹段的单位化向量，并根据每个轨迹段的单位化向量确定所述测井投射轨迹上每个测量点的单位化法向量。
本申请实施例中将向量单位化的目的是为实现确定任意测量深度位置测井数据的显示法向量，并满足达到法向量在井拐点之间平滑过渡的效果。
其中，所述确定每个轨迹段的单位化向量可以为：
根据轨迹段两端的拐点位置可以确定该轨迹段的向量，然后将该轨迹段的向量单位化，从而定轨迹段的单位化向量。所述单位化向量例如可以为根据公式：
其中，ai为第i个轨迹段的单位化向量的横坐标；bi为第i个轨迹段的单位化向量的纵坐标；xi和xi+1分别为测井投射轨迹上第i个和第i+1个拐点的横坐标；yi和yi+1分别为测井投射轨迹上第i个和第i+1个拐点的纵坐标。如图6所示，一测井投射轨迹上有(x1，y1)、(x2，y2)和(x3，y3)三个拐点。那么第一个轨迹段的向量为(x2-x1，y2-y1)。则相应的，其单位化后向量为
其中，所述根据每个轨迹段的单位化向量确定所述测井投射轨迹上每个测量点的单位化法向量，具体可以为：
对于每个轨迹段，先根据公式a＝a1(1-rj)+a2rj和b＝b1(1-rj)+b2rj确定所述测井投射轨迹上每个测量点的法向量。然后再单位化所述测井投射轨迹上每个测量点的法向量，获得所述测井投射轨迹上每个测量点的单位化法向量(请参见上文的单位化向量方法，在此不再赘述)。
其中，a和b分别为测井投射轨迹上第j个测量点的法向量的横坐标和纵坐标；a1和a2分别为该轨迹段内第一个拐点和第二个拐点的横坐标；b1和b2分别为该轨迹段内第一个拐点和第二个拐点的纵坐标；rj为该轨迹段内第j个测量点的系数，且rj＝b/a，其中，b为第j个测量点至所述第一个拐点的长度，a为所述第一个拐点与所述第二个拐点之间的长度。
其次，根据所述测井投射轨迹上每个测量点的单位化法向量及对应测井数据的水平偏移量，确定所述测井投射轨迹上每个测量点对应测井数据的测井数据点显示位置。
本申请实施例中，所述根据所述测井投射轨迹上每个测量点的单位化法向量及对应测井数据的水平偏移量，确定所述测井投射轨迹上每个测量点对应测井数据的测井数据点显示位置，可以包括：
根据公式x＝xgui+a'w和y＝ygui+b'w确定所述测井投射轨迹上每个测量点对应测井数据的测井数据点显示位置；
其中，x和y分别为每个测量点对应测井数据的测井数据点的横坐标和纵坐标；xgui和ygui分别为每个测量点的单位化法向量的横坐标和纵坐标；w为每个测量点对应测井数据的水平偏移量。
最后，按照所述测井数据点显示位置绘制测井数据。
s103、确定所述测井投射轨迹上每个测量点的地层倾向，并按照每个测量点的地层倾向绘制对应的测井解释数据。
本申请实施例中，测井解释数据可以包括分层信息。其中，所述确定所述测井投射轨迹上每个测量点的地层倾向，例如可以为根据测井投射轨迹上每个测量点与地震反射特征的确定地层倾向，然后地层倾向显示分层。如图7所示，假设层位两端的点为piont1，point2，倾斜角度angle＝atan2(point2._y-point1._y,point2._x-point1._x)；在绘图的时可利用该角度绘制出地质分层。
在本申请一个示例性实施例中，在完成上述步骤后可得到如图9所示的在地震剖面上显示测井信息的示意图。从图9中可以看出，测井解释数据中的分层数据是按照地层倾向绘制，而测井数据则是按照测井投射轨迹上每个测量点的法向绘制，这样，本申请实施例的在地震剖面上生成的测井信息方法，避免了现有技术出现的数据挤压变形、显示交叠以及无法准确确定地质分层的走势的问题，从而提高了测井信息在地震剖面上显示的质量。
需要提醒的是，本申请上述实施例中步骤s103和步骤s102并没有具体的顺序要求，可以同时执行，也可以一前一后执行(例如先执行步骤s103再执行步骤s102，或者先执行步骤s102再执行步骤s103)。
虽然上文描述的过程流程包括以特定顺序出现的多个操作，但是，应当清楚了解，这些过程可以包括更多或更少的操作，这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
参考图8所示，本申请实施例的在地震剖面上显示测井信息的装置可以包括：
测井轨迹投射模块81，可以用于在给定的投影距离下将三维空间中的测井轨迹投射到地震剖面上，获得测井投射轨迹；
测井数据绘制模块82，可以用于确定所述测井投射轨迹上每个测量点的法向，并按照每个测量点的法向绘制对应的测井数据；以及，
测井解释数据绘制模块83，可以用于确定所述测井投射轨迹上每个测量点的地层倾向，并按照每个测量点的地层倾向绘制对应的测井解释数据。
本申请实施例的装置与上述实施例的方法对应，因此，有关于本申请的装置细节，请参见上述实施例的方法，在此不再赘述。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(cpu)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(ram)和/或非易失性内存等形式，如只读存储器(rom)或闪存(flashram)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器(cd-rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitorymedia)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
- 还没有人留言评论。精彩留言会获得点赞！