CN114611287A - 一种水平井导向模型动态成像测井分段模拟方法 - Google Patents
一种水平井导向模型动态成像测井分段模拟方法 Download PDFInfo
- Publication number
- CN114611287A CN114611287A CN202210232172.6A CN202210232172A CN114611287A CN 114611287 A CN114611287 A CN 114611287A CN 202210232172 A CN202210232172 A CN 202210232172A CN 114611287 A CN114611287 A CN 114611287A
- Authority
- CN
- China
- Prior art keywords
- simulation
- section
- imaging
- interpretation
- segment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 103
- 238000004088 simulation Methods 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000001514 detection method Methods 0.000 claims abstract description 67
- PEQMJVGRHNZPAM-UHFFFAOYSA-N 1,4-dichloro-2-isocyanatobenzene Chemical compound ClC1=CC=C(Cl)C(N=C=O)=C1 PEQMJVGRHNZPAM-UHFFFAOYSA-N 0.000 claims abstract description 26
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 238000006073 displacement reaction Methods 0.000 claims description 27
- 239000000523 sample Substances 0.000 claims description 19
- 230000015572 biosynthetic process Effects 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 10
- 230000011218 segmentation Effects 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 230000004044 response Effects 0.000 abstract description 3
- 239000013256 coordination polymer Substances 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 14
- 230000008569 process Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/002—Survey of boreholes or wells by visual inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种水平井导向模型动态成像测井分段模拟方法、装置、计算机设备和存储介质,该方法包括:创建剖面导向模型;根据通道数,确定模拟通道的模拟探测位置点CPi,计算每个模拟探测位置点CPi距离井筒中心线的距离DCPi;建立起解释段与模拟成像段1对1的对应关系;获得该通道点的对应的参考曲线的VAL(CP);构建出当前的模拟成像段,依次遍历所有解释段,创建所有成像段;通过编辑解释段,确定需要更新的解释段,重新计算对应的模拟成像段,其他段不变。本申请中对于模型解释段的调整,只需更新对应的成像段即可,缩小了计算范围,响应速度加快。
Description
技术领域
本申请涉及测井技术领域技术领域,特别是涉及一种水平井导向模型动态成像测井分段模拟方法、装置、计算机设备和存储介质。
背景技术
水平井随钻地质导向过程是一个多学科的动态模型修正和预测过程,对实时性、便捷性要求非常高。一般方法是以导向模型为基础,通过不断调整地层倾角和厚度,使实测gr曲线和反演gr曲线对比并趋同,从而不断修正模型,并以此及时提供钻进调整建议。同时,辅助增加成像测井对比,在进行地层倾角和厚度调整过程中,根据实测的成像测井图像道和模拟的成像测井图像道进行对比,起到模型质控作用。成像拟合技术的精度和速度很大程度上决定了动态模型修正和调整的实时性和准确性。
发明内容
基于此,针对上述技术问题,提供一种水平井导向模型动态成像测井分段模拟方法、装置、计算机设备和存储介质。
第一方面,一种水平井导向模型动态成像测井分段模拟方法,所述方法包括:
步骤S1:创建剖面导向模型:获取水平井及水平井测井曲线,获取参考井及参考井测井曲线,创建解释段,通过编辑解释段,对地层线进行调整;新建动态成像测井道,获取需要模拟的参考井测井曲线,模拟通道数CN,井眼直径RW,探测深度DI;
步骤S2:根据通道数,确定模拟通道的模拟探测位置点CPi,计算每个模拟探测位置点CPi距离井筒中心线的距离DCPi(i=1,2,3,...,CN);
步骤S3:根据解释段,分段计算多通道的成像模拟数据,形成模拟成像段,建立起解释段与模拟成像段1对1的对应关系;
步骤S4:获取当前模拟曲线段的开始投影位移VS-Start和结束投影位移VS_End,根据采样步长VS-Step,对当前水平井轨迹段进行横向上按照VS-Step均分,得到当前段的所有轨迹点HwgPoint(m个)的投影位移(VS1-VSm)和位置(POS1-POSm)以及水平倾角Ai(A1-Am);
步骤S5:根据当前轨迹点HwgPoint水平倾角和各通道模拟探测点CPi距离井筒中心线的距离DCPi,计算各探测点在剖面模型中的位置;
计算在参考曲线空间内的海波深AD,进而获得该通道点的对应的参考曲线的VAL(CP);根据所有轨迹点HwgPoint的VS数据和通道探测点集合VAL(CP)数据,构建出当前的模拟成像段,依次遍历所有解释段,创建所有成像段;
步骤S6:通过编辑解释段,确定需要更新的解释段,重新计算对应的模拟成像段,其他段不变。
上述方案中,可选的,步骤S1中,所述导向模型根据TVT(真垂厚)空间变化进行填充显示。
上述方案中,进一步可选的,步骤S2中,所述计算DCPi方法如下:DCPi=cos(2*PAI/CN*i)*(RW/2+DI)。
上述方案中,进一步可选的,步骤S3中,所述导向模型被分割成横向联系且纵向垂直的N个解释段空间,水平井轨迹线被分割为N段。每个解释段区间对应一段模拟成像段。
上述方案中,进一步可选的,步骤S5中,所述计算在参考曲线空间内的海波深AD是通过计算当前模拟探测点CPi在段空间内横向相对距离XSi;根据相对距离XSi,计算导向模型中每个层面线与过CPi的纵向线的层面点Ti的位置POS(Ti);
根据当前模拟探测点CPi的位置,计算其在哪两个层面点之间;并计算在两个层面点之间纵向相对位置;读取参考井两个层面点的海波深AD,计算得到当前模拟探测点CPi,在参考井空间内的海波深AD,记为AD_CPi。
上述方案中,进一步可选的,步骤S5中,所述计算当前模拟探测点CPi在段空间内横向相对距离XSi方法如下:XSi=VS_CUR/VS_DIS;其中,所述VS_CUR为当前点到开始点的VS距离,VS_DIS为解释段的VS距离。
上述方案中,进一步可选的,步骤S5中,所述根据相对距离XSi,计算导向模型中每个层面线与过CPi的纵向线的层面点Ti的位置POS(Ti)的计算方法如下:
POS(Ti)=POS(Ti_Start)*(1-XSi)+POS(Ti_End)*XSi。
第二方面,一种水平井导向模型动态成像测井分段模拟装置,其特征在于,所述装置包括:
构建模块,用于创建剖面导向模型:获取水平井及水平井测井曲线,获取参考井及参考井测井曲线,创建解释段,通过编辑解释段,对地层线进行调整;新建动态成像测井道,获取需要模拟的参考井测井曲线,模拟通道数CN,井眼直径RW,探测深度DI;
计算模块,用于根据通道数,确定模拟通道的模拟探测位置点CPi,计算每个模拟探测位置点CPi距离井筒中心线的距离DCPi(i=1,2,3,...,CN);
用于根据解释段,分段计算多通道的成像模拟数据,形成模拟成像段,建立起解释段与模拟成像段1对1的对应关系;
用于获取当前模拟曲线段的开始投影位移VS-Start和结束投影位移VS_End,根据采样步长VS-Step,对当前水平井轨迹段进行横向上按照VS-Step均分,得到当前段的所有轨迹点HwgPoint(m个)的投影位移(VS1-VSm)和位置(POS1-POSm)以及水平倾角Ai(A1-Am);
用于根据当前轨迹点HwgPoint水平倾角和各通道模拟探测点CPi距离井筒中心线的距离DCPi,计算各探测点在剖面模型中的位置;
用于计算在参考曲线空间内的海波深AD,进而获得该通道点的对应的参考曲线的VAL(CP);根据所有轨迹点HwgPoint的VS数据和通道探测点集合VAL(CP)数据,构建出当前的模拟成像段,依次遍历所有解释段,创建所有成像段;
编辑模块:用于通过编辑解释段,确定需要更新的解释段,重新计算对应的模拟成像段,其他段不变。
第三方面,一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
步骤S1:创建剖面导向模型:获取水平井及水平井测井曲线,获取参考井及参考井测井曲线,创建解释段,通过编辑解释段,对地层线进行调整;新建动态成像测井道,获取需要模拟的参考井测井曲线,模拟通道数CN,井眼直径RW,探测深度DI;
步骤S2:根据通道数,确定模拟通道的模拟探测位置点CPi,计算每个模拟探测位置点CPi距离井筒中心线的距离DCPi(i=1,2,3,...,CN);
步骤S3:根据解释段,分段计算多通道的成像模拟数据,形成模拟成像段,建立起解释段与模拟成像段1对1的对应关系;
步骤S4:获取当前模拟曲线段的开始投影位移VS-Start和结束投影位移VS_End,根据采样步长VS-Step,对当前水平井轨迹段进行横向上按照VS-Step均分,得到当前段的所有轨迹点HwgPoint(m个)的投影位移(VS1-VSm)和位置(POS1-POSm)以及水平倾角Ai(A1-Am);
步骤S5:根据当前轨迹点HwgPoint水平倾角和各通道模拟探测点CPi距离井筒中心线的距离DCPi,计算各探测点在剖面模型中的位置;
计算在参考曲线空间内的海波深AD,进而获得该通道点的对应的参考曲线的VAL(CP);根据所有轨迹点HwgPoint的VS数据和通道探测点集合VAL(CP)数据,构建出当前的模拟成像段,依次遍历所有解释段,创建所有成像段;
步骤S6:通过编辑解释段,确定需要更新的解释段,重新计算对应的模拟成像段,其他段不变。
第四方面,一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现以下步骤:
步骤S1:创建剖面导向模型:获取水平井及水平井测井曲线,获取参考井及参考井测井曲线,创建解释段,通过编辑解释段,对地层线进行调整;新建动态成像测井道,获取需要模拟的参考井测井曲线,模拟通道数CN,井眼直径RW,探测深度DI;
步骤S2:根据通道数,确定模拟通道的模拟探测位置点CPi,计算每个模拟探测位置点CPi距离井筒中心线的距离DCPi(i=1,2,3,...,CN);
步骤S3:根据解释段,分段计算多通道的成像模拟数据,形成模拟成像段,建立起解释段与模拟成像段1对1的对应关系;
步骤S4:获取当前模拟曲线段的开始投影位移VS-Start和结束投影位移VS_End,根据采样步长VS-Step,对当前水平井轨迹段进行横向上按照VS-Step均分,得到当前段的所有轨迹点HwgPoint(m个)的投影位移(VS1-VSm)和位置(POS1-POSm)以及水平倾角Ai(A1-Am);
步骤S5:根据当前轨迹点HwgPoint水平倾角和各通道模拟探测点CPi距离井筒中心线的距离DCPi,计算各探测点在剖面模型中的位置;
计算在参考曲线空间内的海波深AD,进而获得该通道点的对应的参考曲线的VAL(CP);根据所有轨迹点HwgPoint的VS数据和通道探测点集合VAL(CP)数据,构建出当前的模拟成像段,依次遍历所有解释段,创建所有成像段;
步骤S6:通过编辑解释段,确定需要更新的解释段,重新计算对应的模拟成像段,其他段不变。
本发明公开了一种水平井导向模型动态成像测井分段模拟方法,创建剖面导向模型;根据通道数,确定模拟通道的模拟探测位置点CPi,计算每个模拟探测位置点CPi距离井筒中心线的距离DCPi(i=1,2,3,...,CN);建立起解释段与模拟成像段1对1的对应关系;获得该通道点的对应的参考曲线的VAL(CP);构建出当前的模拟成像段,依次遍历所有解释段,创建所有成像段;通过编辑解释段,确定需要更新的解释段,重新计算对应的模拟成像段,其他段不变。
本申请有益效果如下:本申请中模拟成像按照模型解释段分段模拟,模型成像段和曲线段建立了对应关系,所以对于模型解释段的调整,只需更新对应的成像段即可,缩小了计算范围,响应速度加快。随建立的解释段越多调整速度越快的潜在优势,非常适合大范围导向模型的创建和调整。因为,导向模型的解释段越多,单位解释段的长度也就越短,计算的模拟采样点个数就是越少,所以速度就会越快;轨迹点对应的多通道模拟探测点在导向模型中的位置等参数使用TVT(真垂厚)空间参数计算,简单快捷,同时直接使用参考(邻井)曲线进行查找计算,确保了拟合的高精度。
附图说明
图1为本发明一个实施例提供的一种水平井导向模型动态成像测井分段模拟方法的流程示意图;
图2为本发明一个实施例提供的建立导向剖面模型示意图;
图3为本发明一个实施例提供的使用参考井曲线进行地层填充示意图;
图4为本发明一个实施例提供的动态成像测井参数示意图;
图5为本发明一个实施例提供的模拟探测位置示意图示意图;
图6为本发明一个实施例提供的模拟探测位置距离井筒中心线的距离示意图;
图7为本发明一个实施例提供的根据解释段划分拟合成像段示意图;
图8为本发明一个实施例提供的计算当前轨迹点各通道位置围示意图;
图9为本发明一个实施例提供的计算当前模拟探测点及层面点示意图;
图10为本发明一个实施例提供的计算当前模拟探测点的纵向相对位置示意图;
图11为本发明一个实施例提供的计算当前模拟探测点在参考井空间内的海波深及模拟值示意图;
图12为本发明一个实施例提供的当前段的8通道模拟成像示意图示意图;
图13为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
在一个实施例中,如图1所示,提供了一种水平井导向模型动态成像测井分段模拟方法,包括以下步骤:
一种水平井导向模型动态成像测井分段模拟方法,所述方法包括:
步骤S1:创建剖面导向模型:获取水平井及水平井测井曲线,获取参考井及参考井测井曲线,创建解释段,通过编辑解释段,对地层线进行调整;新建动态成像测井道,获取需要模拟的参考井测井曲线,模拟通道数CN,井眼直径RW,探测深度DI;
步骤S2:根据通道数,确定模拟通道的模拟探测位置点CPi,计算每个模拟探测位置点CPi距离井筒中心线的距离DCPi(i=1,2,3,...,CN);
步骤S3:根据解释段,分段计算多通道的成像模拟数据,形成模拟成像段,建立起解释段与模拟成像段1对1的对应关系;
步骤S4:获取当前模拟曲线段的开始投影位移VS-Start和结束投影位移VS_End,根据采样步长VS-Step,对当前水平井轨迹段进行横向上按照VS-Step均分,得到当前段的所有轨迹点HwgPoint(m个)的投影位移(VS1-VSm)和位置(POS1-POSm)以及水平倾角Ai(A1-Am);
步骤S5:根据当前轨迹点HwgPoint水平倾角和各通道模拟探测点CPi距离井筒中心线的距离DCPi,计算各探测点在剖面模型中的位置;
计算在参考曲线空间内的海波深AD,进而获得该通道点的对应的参考曲线的VAL(CP);根据所有轨迹点HwgPoint的VS数据和通道探测点集合VAL(CP)数据,构建出当前的模拟成像段,依次遍历所有解释段,创建所有成像段;
步骤S6:通过编辑解释段,确定需要更新的解释段,重新计算对应的模拟成像段,其他段不变。
本实施例根据通道数、井筒直径和探测距离建立模拟探测点与轨迹点的空间位置关系,进而快速计算模拟探测点在参考井空间下的海拔深,直接采用参考曲线的原始数据精度,计算模拟探测点的模拟值,一方面提高了横向成像对比精度,另一方面具有较大的速度优势,成像拟合速度满足使用鼠标调整解释段过程中的实时性要求。
其中,所述导向模型根据TVT(真垂厚)空间变化进行填充显示。步骤S2中,所述计算DCPi方法如下:DCPi=cos(2*PAI/CN*i)*(RW/2+DI)。
步骤S3中,所述导向模型被分割成横向联系且纵向垂直的N个解释段空间,水平井轨迹线被分割为N段。每个解释段区间对应一段模拟成像段。
上述方案中,进一步可选的,步骤S5中,所述计算在参考曲线空间内的海波深AD是通过计算当前模拟探测点CPi在段空间内横向相对距离XSi;根据相对距离XSi,计算导向模型中每个层面线与过CPi的纵向线的层面点Ti的位置POS(Ti)根据当前模拟探测点CPi的位置,计算其在哪两个层面点之间;并计算在两个层面点之间纵向相对位置;读取参考井两个层面点的海波深AD,计算得到当前模拟探测点CPi,在参考井空间内的海波深AD,记为AD_CPi。
步骤S5中,所述计算当前模拟探测点CPi在段空间内横向相对距离XSi方法如下:XSi=VS_CUR/VS_DIS;其中,所述VS_CUR为当前点到开始点的VS距离,VS_DIS为解释段的VS距离。步骤S5中,所述根据相对距离XSi,计算导向模型中每个层面线与过CPi的纵向线的层面点Ti的位置POS(Ti)的计算方法如下:
POS(Ti)=POS(Ti_Start)*(1-XSi)+POS(Ti_End)*XSi。
在另一个实施例中,1.创建剖面导向模型,剖面导向模型包括三个部分,主体为导向剖面视图,上面为横向柱状图和右侧为纵向柱状图。设置水平井及水平井测井曲线,设置参考井及参考井测井曲线后,导向剖面视图中显示水平井轨迹,横向柱状图显示水平井测井曲线及模拟曲线,纵向柱状图显示水平井和参考井测井曲线。设置参考地层线,此参考地层线根据参考井获取。在剖面视图中和纵向曲线道内绘制地层线,初始状态下,地层线在剖面视图中为水平状(图2)。2.选择需要模拟的参考井曲线,读取参考井曲线数据,建立[海波深-值,AD-Value]的数据结构。导向模型根据TVT(真垂厚)空间变化进行填充显示(图3)。3.新建动态成像测井道,选择参考井测井曲线作为模拟数据源,输入三参数:模拟通道数CN(可以是2通道,4通道,8通道,16通道,32通道),井眼直径RW,探测深度DI(图4)。4.根据通道数,确定模拟通道的模拟探测位置点CPi(图5为2、4、8通道的探测位置CPi示意图)。模拟探测位置点CPi个数等于通道数。5.计算每个模拟探测位置点CPi距离井筒中心线的距离DCPi(i=1,...,CN),计算方法如下DCPi=cos(2*PAI/CN*i)*(RW/2+DI)。如图6所示。6.剖面模型被分割成横向联系且纵向垂直的N个解释段空间,水平井轨迹线被分割为N段。每个解释段区间对应一段模拟成像段。也就是说建立起解释段与待模拟曲线段1对1的对应关系(图7)。此步骤的目标是编辑解释段时,找到对应的成像段进行联动处理,无需更新所有段,进而减少运算量,提高响应速度。7.获取当前段的开始投影位移VS-Start和结束投影位移VS_End,设置VS的采样步长VS-Step(根据实测测井曲线的精度确定,比如0.125米),对当前水平井轨迹段进行横向上按照VS-Step步长均分,得到当前段的所有轨迹点HwgPoint(m个)的投影位移(VS1-VSm)和位置(POS1-POSm)以及每个轨迹点的水平倾角Ai(A1-Am)。8.计算当前轨迹点的各通道位置,根据当前点水平倾角和各通道模拟探测点CPi距离井筒中心线的距离DCPi,计算各探测点在剖面模型中的位置。(图8)9.计算当前模拟探测点CPi在段空间内横向相对距离XSi=VS_CUR/VS_DIS;VS_CUR为当前点到开始点的VS距离,VS_DIS为解释段的VS距离。根据相对距离XSi,计算导向模型中每个层面线与过CPi的纵向线的层面点Ti的位置POS(Ti)(图9中的T1到T9)。计算方法如下:POS(Ti)=POS(Ti_Start)*(1-XSi)+POS(Ti_End)*XSi;10.根据当前模拟探测点CPi的位置,计算其在哪两个层面点之间(图7)并计算在两个层面点之间纵向相对位置。读取参考井两个层面点的海波深AD,计算得到当前模拟探测点CPi,在参考井空间内的海波深AD,记为AD_CPi(图10)。11.根据参考井曲线数据,插值得到海波深为AD_CPi的曲线数值VALi(图11);12.重复步骤8、9、10、11,计算得到当前轨迹点的各通道模拟测探点的模拟数值VAL;13.重复步骤7、8、9、10、11,计算得到当前段内所有轨迹点的各通道模拟测探点的模拟数值VAL;根据轨迹点的VS数据和VAL数据,构建出当前的模拟曲线段(图11)。14.依次遍历所有解释段,重复步骤6、7、8、9、10、11,构建出所有段依次遍历所有解释段。绘制在横向模拟道中。15.编辑解释段时,左右移动或者上下移动当前解释段的起始立柱,影响成像段为当前段和临近的前后段。重新计算影响段即可,无需全部重新计算。16.当前段整体上下拖动时,则只影响当前段,重新计算影响段即可,无需全部重新计算。17.新增解释段时,其本质时把当前段劈分成两个段,所以只影响当前段,其他段无需重算。合并解释段时,向前合并影响为当前段及前一段,向后合并影响为当前段及后一段,重新计算影响段即可,无需全部重新计算。
在一个实施例中,提供了一种水平井导向模型动态成像测井分段模拟装置,包括以下程序模块:构建模块、计算模块和编辑模块,其中:
构建模块,用于创建剖面导向模型:获取水平井及水平井测井曲线,获取参考井及参考井测井曲线,创建解释段,通过编辑解释段,对地层线进行调整;新建动态成像测井道,获取需要模拟的参考井测井曲线,模拟通道数CN,井眼直径RW,探测深度DI;
计算模块,用于根据通道数,确定模拟通道的模拟探测位置点CPi,计算每个模拟探测位置点CPi距离井筒中心线的距离DCPi(i=1,2,3,...,CN);
用于根据解释段,分段计算多通道的成像模拟数据,形成模拟成像段,建立起解释段与模拟成像段1对1的对应关系;
用于获取当前模拟曲线段的开始投影位移VS-Start和结束投影位移VS_End,根据采样步长VS-Step,对当前水平井轨迹段进行横向上按照VS-Step均分,得到当前段的所有轨迹点HwgPoint(m个)的投影位移(VS1-VSm)和位置(POS1-POSm)以及水平倾角Ai(A1-Am);
用于根据当前轨迹点HwgPoint水平倾角和各通道模拟探测点CPi距离井筒中心线的距离DCPi,计算各探测点在剖面模型中的位置;
用于计算在参考曲线空间内的海波深AD,进而获得该通道点的对应的参考曲线的VAL(CP);根据所有轨迹点HwgPoint的VS数据和通道探测点集合VAL(CP)数据,构建出当前的模拟成像段,依次遍历所有解释段,创建所有成像段;
编辑模块:用于通过编辑解释段,确定需要更新的解释段,重新计算对应的模拟成像段,其他段不变。
关于一种水平井导向模型动态成像测井分段模拟装置的具体限定可以参见上文中对于一种水平井导向模型动态成像测井分段模拟方法的限定,在此不再赘述。上述一种水平井导向模型动态成像测井分段模拟装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,该计算机设备可以是服务器,其内部结构图可以如图13所示。该计算机设备包括通过系统总线连接的处理器、存储器和网络接口。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种一种水平井导向模型动态成像测井分段模拟方法。
本领域技术人员可以理解,图13中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
在一个实施例中,提供了一种计算机设备,包括存储器和处理器,存储器中存储有计算机程序,涉及上述实施例方法中的全部或部分流程。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,涉及上述实施例方法中的全部或部分流程。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(Read-Only Memory,ROM)、磁带、软盘、闪存或光存储器等。易失性存储器可包括随机存取存储器(Random Access Memory,RAM)或外部高速缓冲存储器。作为说明而非局限,RAM可以是多种形式,比如静态随机存取存储器(Static Random Access Memory,SRAM)或动态随机存取存储器(Dynamic Random Access Memory,DRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种水平井导向模型动态成像测井分段模拟方法,其特征在于,所述方法包括:
步骤S1:创建剖面导向模型:获取水平井及水平井测井曲线,获取参考井及参考井测井曲线,创建解释段,通过编辑解释段,对地层线进行调整;新建动态成像测井道,获取需要模拟的参考井测井曲线,模拟通道数CN,井眼直径RW,探测深度DI;
步骤S2:根据通道数,确定模拟通道的模拟探测位置点CPi,计算每个模拟探测位置点CPi距离井筒中心线的距离DCPi(i=1,2,3,...,CN);
步骤S3:根据解释段,分段计算多通道的成像模拟数据,形成模拟成像段,建立起解释段与模拟成像段1对1的对应关系;
步骤S4:获取当前模拟曲线段的开始投影位移VS-Start和结束投影位移VS_End,根据采样步长VS-Step,对当前水平井轨迹段进行横向上按照VS-Step均分,得到当前段的所有轨迹点HwgPoint(m个)的投影位移(VS1-VSm)和位置(POS1-POSm)以及水平倾角Ai(A1-Am);
步骤S5:根据当前轨迹点HwgPoint水平倾角和各通道模拟探测点CPi距离井筒中心线的距离DCPi,计算各探测点在剖面模型中的位置;
计算在参考曲线空间内的海波深AD,进而获得该通道点的对应的参考曲线的VAL(CP);根据所有轨迹点HwgPoint的VS数据和通道探测点集合VAL(CP)数据,构建出当前的模拟成像段,依次遍历所有解释段,创建所有成像段;
步骤S6:通过编辑解释段,确定需要更新的解释段,重新计算对应的模拟成像段,其他段不变。
2.根据权利要求1所述的方法,其特征在于,步骤S1中,所述导向模型根据TVT(真垂厚)空间变化进行填充显示。
3.根据权利要求1所述的方法,其特征在于,步骤S2中,所述计算DCPi方法如下:DCPi=cos(2*PAI/CN*i)*(RW/2+DI)。
4.根据权利要求1所述的方法,其特征在于,步骤S3中,所述导向模型被分割成横向联系且纵向垂直的N个解释段空间,水平井轨迹线被分割为N段;每个解释段区间对应一段模拟成像段。
5.根据权利要求1所述的方法,其特征在于,步骤S5中,所述计算在参考曲线空间内的海波深AD是通过计算当前模拟探测点CPi在段空间内横向相对距离XSi;根据相对距离XSi,计算导向模型中每个层面线与过CPi的纵向线的层面点Ti的位置POS(Ti);
根据当前模拟探测点CPi的位置,计算其在哪两个层面点之间;并计算在两个层面点之间纵向相对位置;读取参考井两个层面点的海波深AD,计算得到当前模拟探测点CPi,在参考井空间内的海波深AD,记为AD_CPi。
6.根据权利要求5所述的方法,其特征在于,步骤S5中,所述计算当前模拟探测点CPi在段空间内横向相对距离XSi方法如下:XSi=VS_CUR/VS_DIS;其中,所述VS_CUR为当前点到开始点的VS距离,VS_DIS为解释段的VS距离。
7.根据权利要求5所述的方法,其特征在于,步骤S5中,所述根据相对距离XSi,计算导向模型中每个层面线与过CPi的纵向线的层面点Ti的位置POS(Ti)的计算方法如下:
POS(Ti)=POS(Ti_Start)*(1-XSi)+POS(Ti_End)*XSi。
8.一种水平井导向模型动态成像测井分段模拟装置,其特征在于,所述装置包括:
构建模块,用于创建剖面导向模型:获取水平井及水平井测井曲线,获取参考井及参考井测井曲线,创建解释段,通过编辑解释段,对地层线进行调整;新建动态成像测井道,获取需要模拟的参考井测井曲线,模拟通道数CN,井眼直径RW,探测深度DI;
计算模块,用于根据通道数,确定模拟通道的模拟探测位置点CPi,计算每个模拟探测位置点CPi距离井筒中心线的距离DCPi(i=1,2,3,...,CN);
用于根据解释段,分段计算多通道的成像模拟数据,形成模拟成像段,建立起解释段与模拟成像段1对1的对应关系;
用于获取当前模拟曲线段的开始投影位移VS-Start和结束投影位移VS_End,根据采样步长VS-Step,对当前水平井轨迹段进行横向上按照VS-Step均分,得到当前段的所有轨迹点HwgPoint(m个)的投影位移(VS1-VSm)和位置(POS1-POSm)以及水平倾角Ai(A1-Am);
用于根据当前轨迹点HwgPoint水平倾角和各通道模拟探测点CPi距离井筒中心线的距离DCPi,计算各探测点在剖面模型中的位置;
用于计算在参考曲线空间内的海波深AD,进而获得该通道点的对应的参考曲线的VAL(CP);根据所有轨迹点HwgPoint的VS数据和通道探测点集合VAL(CP)数据,构建出当前的模拟成像段,依次遍历所有解释段,创建所有成像段;
编辑模块:用于通过编辑解释段,确定需要更新的解释段,重新计算对应的模拟成像段,其他段不变。
9.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述的方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至7中任一项所述的方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210232172.6A CN114611287B (zh) | 2022-03-09 | 一种水平井导向模型动态成像测井分段模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210232172.6A CN114611287B (zh) | 2022-03-09 | 一种水平井导向模型动态成像测井分段模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114611287A true CN114611287A (zh) | 2022-06-10 |
CN114611287B CN114611287B (zh) | 2024-06-25 |
Family
ID=
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130238306A1 (en) * | 2010-05-21 | 2013-09-12 | Landmark Graphics Corporation | Systems and Methods for Horizontal Well Correlation and Geosteering |
CN105317431A (zh) * | 2014-07-30 | 2016-02-10 | 中国石油化工股份有限公司 | 一种用于解释和评价水平井测井参数的方法 |
US20160202390A1 (en) * | 2013-10-01 | 2016-07-14 | Landmark Graphics Corporation | In-Situ Wellbore, Core and Cuttings Information System |
CN107861917A (zh) * | 2017-11-29 | 2018-03-30 | 中国石油集团长城钻探工程有限公司 | 水平井中多元数据联合计算井眼到地层边界距离的方法 |
CN109630089A (zh) * | 2018-10-29 | 2019-04-16 | 中国石油天然气股份有限公司 | 水平井地质构造识别方法及装置 |
CN110080754A (zh) * | 2019-04-25 | 2019-08-02 | 杭州迅美科技有限公司 | 一种电成像测井图像类周期性干扰处理方法 |
US10689954B1 (en) * | 2019-04-24 | 2020-06-23 | Dagang Oil Field Company Of Cnpc | Research method of trajectory design and on-site tracking and adjustment of shale oil horizontal well |
CN111577250A (zh) * | 2019-02-15 | 2020-08-25 | 中国石油天然气股份有限公司 | 水平井地质导向方法及装置 |
CN111985081A (zh) * | 2020-07-15 | 2020-11-24 | 北京金阳普泰石油技术股份有限公司 | 一种测井曲线构建方法、系统、设备及可读存储介质 |
CN113536574A (zh) * | 2021-07-20 | 2021-10-22 | 北京蓝海智信能源技术有限公司 | 一种地质导向建模方法、系统及计算机可读存储介质 |
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130238306A1 (en) * | 2010-05-21 | 2013-09-12 | Landmark Graphics Corporation | Systems and Methods for Horizontal Well Correlation and Geosteering |
US20160202390A1 (en) * | 2013-10-01 | 2016-07-14 | Landmark Graphics Corporation | In-Situ Wellbore, Core and Cuttings Information System |
CN105317431A (zh) * | 2014-07-30 | 2016-02-10 | 中国石油化工股份有限公司 | 一种用于解释和评价水平井测井参数的方法 |
CN107861917A (zh) * | 2017-11-29 | 2018-03-30 | 中国石油集团长城钻探工程有限公司 | 水平井中多元数据联合计算井眼到地层边界距离的方法 |
CN109630089A (zh) * | 2018-10-29 | 2019-04-16 | 中国石油天然气股份有限公司 | 水平井地质构造识别方法及装置 |
CN111577250A (zh) * | 2019-02-15 | 2020-08-25 | 中国石油天然气股份有限公司 | 水平井地质导向方法及装置 |
US10689954B1 (en) * | 2019-04-24 | 2020-06-23 | Dagang Oil Field Company Of Cnpc | Research method of trajectory design and on-site tracking and adjustment of shale oil horizontal well |
CN110080754A (zh) * | 2019-04-25 | 2019-08-02 | 杭州迅美科技有限公司 | 一种电成像测井图像类周期性干扰处理方法 |
CN111985081A (zh) * | 2020-07-15 | 2020-11-24 | 北京金阳普泰石油技术股份有限公司 | 一种测井曲线构建方法、系统、设备及可读存储介质 |
CN113536574A (zh) * | 2021-07-20 | 2021-10-22 | 北京蓝海智信能源技术有限公司 | 一种地质导向建模方法、系统及计算机可读存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112150582B (zh) | 一种面向多模态数据的地质剖面图近似表达方法 | |
CN114542056A (zh) | 水平井地层解释模型快速构建方法和装置 | |
CN106285479B (zh) | 一种基于三维gis技术的井眼轨迹控制方法 | |
CN114611287B (zh) | 一种水平井导向模型动态成像测井分段模拟方法 | |
CN114611287A (zh) | 一种水平井导向模型动态成像测井分段模拟方法 | |
CN111260791B (zh) | 一种更新地质导向模型的方法 | |
CN107945261B (zh) | 一种地质钻孔测斜计算方法 | |
RU2745152C1 (ru) | Способ объединения модели геологической проводки скважины с оперативной петрофизической интерпретацией данных гис в режиме реального времени и система, реализующая способ | |
CN113536574A (zh) | 一种地质导向建模方法、系统及计算机可读存储介质 | |
WO2023236754A1 (zh) | 地质导向方法、装置及存储介质 | |
CN111208558B (zh) | 超深低幅度三维地质构造的建立方法及装置 | |
CN114580169A (zh) | 一种水平井导向模型高精度反演曲线分段模拟方法 | |
CN113658335B (zh) | 岩溶发育区复杂地质模型的创建方法 | |
CN107818217B (zh) | 一种基于水平井探索式可视分析的地震导向方法 | |
CN114764844B (zh) | 页岩气钻井优化设计方法及装置 | |
CN111754611B (zh) | 基于地质模型约束的井间砂体连通图自动绘制方法及系统 | |
CN110761780A (zh) | 一种基于井震结合的三维地质导向方法 | |
CN112435334A (zh) | 一种任意地层覆盖层厚度的计算方法、系统及介质 | |
CN113221219A (zh) | 一种水平井靶点设计中邻井构造剖面的构建方法及系统 | |
CN114564800B (zh) | 一种水平井测井曲线的真垂厚纵向拟合方法及系统 | |
CN116484457B (zh) | 钻爆法隧道掌子面围岩三维精细化分级方法、装置及介质 | |
CN111722276A (zh) | 岩石钻井的地震导向方法及系统 | |
CN114861260A (zh) | 随钻导向区域微构造调整方法和装置、设备和存储介质 | |
CN112363221B (zh) | Walkaway VSP测线的布设方法及装置 | |
CN113495303B (zh) | 剥蚀量恢复的方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |