基于时间度量的斜井井轨迹计算方法技术

本发明专利技术提供一种基于时间度量的斜井井轨迹计算方法，包括：通过井数据获取斜井W的井眼坐标X、Y，通过井斜数据对采样点依次进行计算，得到当前采样点下更新后的斜井平面坐标X、Y及累积垂深；通过层位数据使用最小绝对值距离来精确匹配更新后的平面坐标X、Y所在地震道对应目标层位的时间t1；通过时深数据使用最小绝对值距离来精确匹配到达该深度所需的时间t2；比较t1和t2的大小，当t2≥t1时，此时的平面坐标X、Y即斜井在目标层位的真实坐标。该方法实现了以时间度量作为标准来进行斜井在目标层位的坐标校正，通过最小绝对值距离来精确匹配目标层位的地震道，为井筒数据和地震属性数据相关关系的研究提供准确的井坐标支持。

Calculating method of inclined well trajectory based on time measurement

The invention provides a method for calculating inclined well trajectory based on time measurement, which includes: obtaining borehole coordinates X and Y of inclined well W by well data, calculating sampling points in turn by well deviation data, obtaining updated inclined well plane coordinates X, Y and accumulated vertical depth under current sampling points; accurately matching updated plane coordinates X and Y by using minimum absolute value distance of layer data; and accurately matching updated plane coordinates X and Y by using minimum absolute value distance of layer data. The seismic trace corresponds to the time T1 of the target horizon; the minimum absolute distance is used to accurately match the time T2 needed to reach the depth through time-depth data; and the size of T1 and T2 is compared. When T2 is greater than t1, the plane coordinates X and Y are the real coordinates of the inclined well in the target horizon. This method realizes the coordinate correction of inclined wells in target horizon based on time measurement, matches the seismic channel of target horizon accurately by minimum absolute distance, and provides accurate well coordinate support for the study of correlation between wellbore data and seismic attribute data.

【技术实现步骤摘要】
基于时间度量的斜井井轨迹计算方法
本专利技术涉及地球物理勘探和计算机
，特别是涉及到一种基于时间度量的斜井井轨迹计算方法。
技术介绍
在研究井筒数据和地震属性数据相关关系的过程中，通常以井在目标层位的平面坐标(X、Y)为基础，选取固定大小的时窗(如20ms、30ms等)来匹配井筒数据和地震属性数据得到带有类别标签的样本数据。由于勘探开发难度的不断提高，斜井数量在钻井工程中的比例越来越大，而斜井由于井眼轨道偏移，其在目标层位的平面坐标已发生偏移，需要对斜井坐标进行校正。常规的斜井井轨迹计算方法有正切法、平均角法、平衡正切法、圆柱螺旋线法等，使用这些方法进行计算后得到包含多个采样点的井斜数据(包含井斜角、方位角、测点深度、坐标偏移量等)，并不能直接得到斜井在目标层位的真实平面坐标，因此需要对斜井在目标层位的平面坐标进行校正。为此我们专利技术了一种新的基于时间度量的斜井井轨迹计算方法，解决了以上技术问题。
技术实现思路
本专利技术的目的是提供一种可以计算斜井在目标层位的真实平面坐标，实现斜井在目标层位的坐标校正的基于时间度量的斜井井轨迹计算方法。本专利技术的目的可通过如下技术措施来实现：基于时间度量的斜井井轨迹计算方法，该基于时间度量的斜井井轨迹计算方法包括：步骤1，通过井数据获取斜井W的井眼坐标X、Y，通过井斜数据对采样点依次进行计算，得到当前采样点下更新后的斜井平面坐标X、Y及累积垂深；步骤2，通过层位数据使用最小绝对值距离来精确匹配更新后的平面坐标X、Y所在地震道对应目标层位的时间t1；步骤3，通过时深数据使用最小绝对值距离来精确匹配到达该深度所需的时间t2；步骤4，比较t1和t2的大小，当t2≥t1时，此时的平面坐标X、Y即斜井在目标层位的真实坐标。本专利技术的目的还可通过如下技术措施来实现：在步骤1中，通过井数据获取斜井W的井眼坐标X、Y，通过井斜数据获取当前采样点的偏移量，计算该样点更新后的斜井平面坐标X、Y及累积垂深，坐标更新公式如下：X＝X+DNY＝Y+DEDN表示坐标在正北方向的变化量，DE表示坐标在正东方向的变化量。在步骤2中，通过步骤1获得当前采样点更新后的斜井平面坐标X、Y，在此基础上通过层位数据使用最小绝对值距离来精确匹配更新后的X、Y对应目标层位的地震道，进而获得该地震道到达目标层位的时间t1，其中绝对值距离公式如下：设m个指标，X1，X2，…，Xn为n个样品，记Xi＝(xi1，xi2，…，xim)，i＝1，2，…，n.dij为Xi与Xj之间的距离；上述公式为绝对值距离的通用公式，本专利技术匹配地震道的最小绝对值距离公式为：di＝|Xi-X|+|Yi-Y|Xi、Yi表示第i个地震道的坐标，X、Y表示更新后的斜井平面坐标，di表示更新后的平面坐标X、Y与第i个地震道的绝对值距离。在步骤3中，通过步骤1获得当前采样点更新后的斜井累积垂深TVD，通过时深数据使用最小绝对值距离来精确匹配到达该深度所需的时间t2，绝对值距离公式如下：dj＝|vdj-TVD|vdj表示时深数据中的第j个垂深值，TVD当前采样点更新后的斜井累积垂深，dj表示TVD与vdj的绝对值距离。在步骤4中，当计算的采样点个数较少时t1＞＞t2，则流程返回到步骤1，继续计算下一个采样点更新后的平面坐标和累积垂深并比较t1和t2，直到t2≥t1时，此时的平面坐标X、Y即斜井在目标层位的真实坐标。本专利技术中的基于时间度量的斜井井轨迹计算方法，实现了以时间度量作为标准来进行斜井在目标层位的坐标校正，通过最小绝对值距离来精确匹配目标层位的地震道，可以为井筒数据和地震属性数据相关关系的研究提供准确的井坐标支持。该方法实现了以时间度量作为标准来进行斜井在目标层位的坐标校正，通过最小绝对值距离来精确匹配目标层位的地震道，可以为井筒数据和地震属性数据相关关系的研究提供准确的井坐标支持。附图说明图1为本专利技术的基于时间度量的斜井井轨迹计算方法的一具体实施例的流程图。具体实施方式为使本专利技术的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。如图1所示，图1为本专利技术的基于时间度量的斜井井轨迹计算方法的流程图。步骤101，斜井坐标及累积垂深更新。通过井数据获取斜井W的井眼坐标X、Y，通过井斜数据获取当前采样点的偏移量，计算该样点更新后的斜井平面坐标X、Y及累积垂深，坐标更新公式如下：X＝X+DNY＝Y+DEDN表示坐标在正北方向的变化量，DE表示坐标在正东方向的变化量。步骤102，坐标时间匹配。通过步骤101获得当前采样点更新后的斜井平面坐标X、Y，在此基础上通过层位数据使用最小绝对值距离来精确匹配更新后的X、Y对应目标层位的地震道，进而获得该地震道到达目标层位的时间t1，其中绝对值距离公式如下：设m个指标，X1，X2，…，Xn为n个样品，记Xi＝(xi1，xi2，…，xim)，i＝1，2，…，n.dij为Xi与Xj之间的距离。上述公式为绝对值距离的通用公式，本专利技术匹配地震道的最小绝对值距离公式为：di＝|Xi-X|+|Yi-Y|Xi、Yi表示第i个地震道的坐标，X、Y表示更新后的斜井平面坐标，di表示更新后的平面坐标X、Y与第i个地震道的绝对值距离。步骤103，时深匹配。通过步骤101获得当前采样点更新后的斜井累积垂深TVD，通过时深数据使用最小绝对值距离来精确匹配到达该深度所需的时间t2，绝对值距离公式如下：dj＝|vdj-TVD|vdj表示时深数据中的第j个垂深值，TVD当前采样点更新后的斜井累积垂深，dj表示TVD与vdj的绝对值距离。步骤104，时间度量。通过步骤102获得当前采样点更新后的平面坐标X、Y所在地震道对应目标层位的时间t1，通过步骤103获得当前采样点更新后的累积垂深对应的时间t2；比较t1和t2的大小，当计算的采样点个数较少时t1＞＞t2；重复上述步骤101到步骤103，继续计算下一个采样点更新后的平面坐标和累积垂深并比较t1和t2，直到t2≥t1时，流程进入到步骤105。步骤105，此时的平面坐标X、Y即斜井在目标层位的真实坐标。以斜井“林庄5-斜3”为例，计算该斜井在第四层的平面坐标X、Y，包括如下步骤：A.斜井坐标及累积垂深更新：通过井数据得到“林庄5-斜3”的井眼坐标X＝4137209.29、Y＝632317.79，通过井斜数据更新该井坐标和累积垂深，第一个采样点更新后的平面坐标X＝4137209.32,Y＝632317.77,累积垂深TVD＝29.35。B.坐标时间匹配：通过A获得本次更新后的平面坐标X＝4137209.32,Y＝632317.77,通过层位数据使用最小绝对值距离来精确匹配更新后的X、Y所在地震道到达第四层的时间t1＝2382.34。C.时深匹配：通过A中获得第一个采样点更新后的斜井累积垂深TVD＝29.35，通过时深数据使用最小绝对值距离来精确匹配当前采样点到达该深度所需的时间t2,当更新次数较少时，由于t2远远小于t1，故将t2取0。D.时间度量：比较t1和t2的大小，此时t1≥t2；重复上述A到C步骤，继续计算下一个采样点更新后的平面坐标和累积垂深并比较t1和t2，直到t2≥t1时结束，此时的X、Y坐标即斜井在第四层的平面坐标X＝4137113.08,Y＝本文档来自技高网...

【技术保护点】
1.基于时间度量的斜井井轨迹计算方法，其特征在于，该基于时间度量的斜井井轨迹计算方法包括：步骤1，通过井数据获取斜井W的井眼坐标X、Y，通过井斜数据对采样点依次进行计算，得到当前采样点下更新后的斜井平面坐标X、Y及累积垂深；步骤2，通过层位数据使用最小绝对值距离来精确匹配更新后的平面坐标X、Y所在地震道对应目标层位的时间t1；步骤3，通过时深数据使用最小绝对值距离来精确匹配到达该深度所需的时间t2；步骤4，比较t1和t2的大小，当t2≥t1时，此时的平面坐标X、Y即斜井在目标层位的真实坐标。

【技术特征摘要】
1.基于时间度量的斜井井轨迹计算方法，其特征在于，该基于时间度量的斜井井轨迹计算方法包括：步骤1，通过井数据获取斜井W的井眼坐标X、Y，通过井斜数据对采样点依次进行计算，得到当前采样点下更新后的斜井平面坐标X、Y及累积垂深；步骤2，通过层位数据使用最小绝对值距离来精确匹配更新后的平面坐标X、Y所在地震道对应目标层位的时间t1；步骤3，通过时深数据使用最小绝对值距离来精确匹配到达该深度所需的时间t2；步骤4，比较t1和t2的大小，当t2≥t1时，此时的平面坐标X、Y即斜井在目标层位的真实坐标。2.根据权利要求1所述的基于时间度量的斜井井轨迹计算方法，其特征在于，在步骤1中，通过井数据获取斜井W的井眼坐标X、Y，通过井斜数据获取当前采样点的偏移量，计算该样点更新后的斜井平面坐标X、Y及累积垂深，坐标更新公式如下：X＝X+DNY＝Y+DEDN表示坐标在正北方向的变化量，DE表示坐标在正东方向的变化量。3.根据权利要求1所述的基于时间度量的斜井井轨迹计算方法，其特征在于，在步骤2中，通过步骤1获得当前采样点更新后的斜井平面坐标X、Y，在此基础上通过层位数据使用最小绝对值距离来精确匹配更新后的X、Y对应目标层位的地震道，进而获得该地...

【专利技术属性】
技术研发人员：路慎强，赫俊民，王兴谋，郭俊，余学锋，林亚林，孙淑艳，彭英，李克文，毕丽飞，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司物探研究院，
类型：发明
国别省市：山东,37