CN109712238A - 平滑处理地质模型的方法及设备 - Google Patents
平滑处理地质模型的方法及设备 Download PDFInfo
- Publication number
- CN109712238A CN109712238A CN201811588306.8A CN201811588306A CN109712238A CN 109712238 A CN109712238 A CN 109712238A CN 201811588306 A CN201811588306 A CN 201811588306A CN 109712238 A CN109712238 A CN 109712238A
- Authority
- CN
- China
- Prior art keywords
- smoothing
- geological model
- window
- dimension
- smooth window
- 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
- 238000009499 grossing Methods 0.000 title claims abstract description 119
- 238000000034 method Methods 0.000 title claims abstract description 73
- 238000012545 processing Methods 0.000 title abstract description 10
- 238000004891 communication Methods 0.000 claims description 8
- 239000011229 interlayer Substances 0.000 description 32
- 230000000694 effects Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 5
- 239000002356 single layer Substances 0.000 description 5
- 230000006870 function Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明实施例提供了一种平滑处理地质模型的方法及设备。其中,所述方法包括:根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,得到通用平滑窗口;根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,并采用所述最终通用平滑窗口对地质模型进行平滑处理。本发明实施例提供的平滑处理地质模型的方法及设备,通过根据待平滑处理的地质模型的维度,设置平滑窗口的维度,再对设置好维度的平滑窗口进行权值设置,可以提高针对线性分布的地质模型的平滑处理精度,有效避免对地质模型的误平滑处理。
Description
技术领域
本发明实施例涉及图像处理技术领域,尤其涉及一种平滑处理地质模型的方法及设备。
背景技术
序贯指示模拟、多点地质统计学SNESIM等基于像元的建模算法,在模拟离散变量(如岩性、岩相等)过程中,采用的是根据随机路径确定的网格点逐点进行模拟,因为是单个网格单元的模拟,因而容易导致模拟结果出现相对离散的点,这些离散的点不符合地质认识,是所谓的“噪点”。传统的模型平滑方法Maps可以在一定程度上解决该问题。在Maps方法中,主要采用了一种平滑窗口,如二维情况可以采用3*3、5*5窗口对模型进行平滑处理,一般情况下Maps可以取得较好的效果。但是在某些特殊情况下,现有的Maps平滑方法不适用,这主要体现在对连续性很好,但是呈线性分布的地质模型的处理中,即这种平滑方法针对线性分布的相带是不适用的。基于这种情况,找到一种能够针对线性分布相带的地质模型的平滑处理方法,就成为业界广泛关注的技术问题。
发明内容
针对现有技术存在的上述问题,本发明实施例提供了一种平滑处理地质模型的方法及设备。
第一方面,本发明的实施例提供了一种平滑处理地质模型的方法,包括:根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,得到通用平滑窗口;根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,并采用所述最终通用平滑窗口对地质模型进行平滑处理。
进一步地,所述根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,包括:所述平滑窗口的纵向维度与所述地质模型的纵向维度相同,所述平滑窗口的横向维度包括:
m=2n+1
其中,m为平滑窗口的横向维度,n为整数,且n满足n≥L或n<L,L为所述地质模型中待处理部分的横向维度。
进一步地,所述根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,包括:所述平滑窗口的横向维度与所述地质模型的横向维度相同,所述平滑窗口的纵向维度包括:
m=2n+1
其中,m为平滑窗口的纵向维度,n为整数,且n满足n≥L或n<L,L为所述地质模型中待处理部分的纵向维度。
进一步地,所述根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,包括:所述最终通用平滑窗口的中心代码点的权值为K,所述中心代码点的左侧和右侧的第i个代码点的权值为Ki,除过所述中心代码点及权值为Ki的代码点,其余代码点权值为T;其中,i=1,…,n。
进一步地,所述的平滑处理地质模型的方法,还包括:若平滑处理的过程中保留长度为L的目标物,则n满足n≥L,且K>sum(T);若平滑处理的过程中剔除长度为L的目标物,则n满足n<L,且Kn>K-sum(T);其中,sum(T)为所有权值T的算术和。
进一步地,所述根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,包括:所述最终通用平滑窗口的中心代码点的权值为K,所述中心代码点的上侧和下侧的第i个代码点的权值为Ki,除过所述中心代码点及权值为Ki的代码点,其余代码点权值为T;其中,i=1,…,n。
进一步地,所述的平滑处理地质模型的方法,还包括:若平滑处理的过程中保留长度为L的目标物,则n满足n≥L,且K>sum(T);若平滑处理的过程中剔除长度为L的目标物,则n满足n<L,且Kn>K-sum(T);其中,sum(T)为所有权值T的算术和。
第二方面,本发明的实施例提供了一种平滑处理地质模型的装置,包括:
通用平滑窗口获取模块,用于根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,得到通用平滑窗口;
最终通用平滑窗口获取模块,用于根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,并采用所述最终通用平滑窗口对地质模型进行平滑处理。
第三方面,本发明的实施例提供了一种电子设备,包括:
至少一个处理器;以及
与处理器通信连接的至少一个存储器,其中:
存储器存储有可被处理器执行的程序指令,处理器调用程序指令能够执行第一方面的各种可能的实现方式中任一种可能的实现方式所提供的平滑处理地质模型的方法。
第四方面,本发明的实施例提供了一种非暂态计算机可读存储介质,非暂态计算机可读存储介质存储计算机指令,计算机指令使计算机执行第一方面的各种可能的实现方式中任一种可能的实现方式所提供的平滑处理地质模型的方法。
本发明实施例提供的平滑处理地质模型的方法及设备,通过根据待平滑处理的地质模型的维度,设置平滑窗口的维度,再对设置好维度的平滑窗口进行权值设置,可以提高针对线性分布的地质模型的平滑处理精度,有效避免对地质模型的误平滑处理。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做一简单的介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的平滑处理地质模型的方法流程图;
图2为本发明实施例提供的单层泥岩分布示意图;
图3为本发明实施例提供的通用平滑窗口结构示意图;
图4为现有技术提供的SISIM模拟夹层分布示意图;
图5为现有技术提供的Maps方法平滑效果示意图;
图6为本发明实施例提供的针对泥岩长度大于30米的平滑结果示意图;
图7为本发明实施例提供的针对泥岩长度大于40米的平滑结果示意图
图8为本发明实施例提供的平滑处理地质模型的装置结构示意图;
图9为本发明实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。另外,本发明提供的各个实施例或单个实施例中的技术特征可以相互任意结合,以形成可行的技术方案,但是必须是以本领域普通技术人员能够实现为基础,当技术方案的结合出现相互矛盾或无法实现时,应当认为这种技术方案的结合不存在,也不在本发明要求的保护范围之内。
本发明实施例提供了一种平滑处理地质模型的方法,参见图1,该方法包括:
101、根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,得到通用平滑窗口;
102、根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,并采用所述最终通用平滑窗口对地质模型进行平滑处理。
在上述实施例的基础上,本发明实施例中提供的平滑处理地质模型的方法,所述根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,包括:所述平滑窗口的纵向维度与所述地质模型的纵向维度相同,所述平滑窗口的横向维度包括:
m=2n+1
其中,m为平滑窗口的横向维度,n为整数,且n满足n≥L或n<L,L为所述地质模型中待处理部分的横向维度。
在上述实施例的基础上,本发明实施例中提供的平滑处理地质模型的方法,所述根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,包括:所述平滑窗口的横向维度与所述地质模型的横向维度相同,所述平滑窗口的纵向维度包括:
m=2n+1
其中,m为平滑窗口的纵向维度,n为整数,且n满足n≥L或n<L,L为所述地质模型中待处理部分的纵向维度。
在上述实施例的基础上,本发明实施例中提供的平滑处理地质模型的方法,所述根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,包括:所述最终通用平滑窗口的中心代码点的权值为K,所述中心代码点的左侧和右侧的第i个代码点的权值为Ki,除过所述中心代码点及权值为Ki的代码点,其余代码点权值为T;其中,i=1,…,n。
在上述实施例的基础上,本发明实施例中提供的平滑处理地质模型的方法,还包括:若平滑处理的过程中保留长度为L的目标物,则n满足n≥L,且K>sum(T);若平滑处理的过程中剔除长度为L的目标物,则n满足n<L,且Kn>K-sum(T);其中,sum(T)为所有权值T的算术和。
在上述实施例的基础上,本发明实施例中提供的平滑处理地质模型的方法,所述根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,包括:所述最终通用平滑窗口的中心代码点的权值为K,所述中心代码点的上侧和下侧的第i个代码点的权值为Ki,除过所述中心代码点及权值为Ki的代码点,其余代码点权值为T;其中,i=1,…,n。
在上述实施例的基础上,本发明实施例中提供的平滑处理地质模型的方法,还包括:若平滑处理的过程中保留长度为L的目标物,则n满足n≥L,且K>sum(T);若平滑处理的过程中剔除长度为L的目标物,则n满足n<L,且Kn>K-sum(T);其中,sum(T)为所有权值T的算术和。
本发明实施例提供的平滑处理地质模型的方法,通过根据待平滑处理的地质模型的维度,设置平滑窗口的维度,再对设置好维度的平滑窗口进行权值设置,可以提高针对线性分布的地质模型的平滑处理精度,有效避免对地质模型的误平滑处理。
为了更加清晰地阐述本发明技术方案的精神实质,下面将以一个更为具体的实施例对本发明的技术方案进行描述。需要说明的是,该实施例仅仅是对本发明前述实施例做出的进一步解释,并不是对本发明保护范围的限制,任何满足本发明技术本质的技术方案,均在本专利的保护范围之内。
下面采用通用模板进一步说明改进平滑窗口的工作原理以及需要注意的问题。参见图2,在图2中,假设岩相模型中单层泥岩分布代码如下,0代表背景或者其他岩相类型,1代表泥岩夹层,最左边或最右边的1表示泥岩分布的两端,1的数量为L,L表示泥岩长度。平滑窗口的设计参见图3,图3中采用m*3的平滑窗口去平滑图4所示泥岩夹层。当窗口中心点K移动到图2中最左边1处时,属于泥岩的权值为K+K1+…+Kn,属于其他相的权值为K1+…+Kn+sum(‘1’),sum(‘1’)表示平滑窗口中所有权值为1的算术和。为保证泥岩夹层两端(图2中最左侧及最右侧的1)不会被其他相类型(图2中的代码0)替代,属于其他相的权值应该小于属于泥岩的权值,所以平滑窗口应该满足K>sum(‘1’),这样平滑过程中泥岩夹层两端不会被平滑掉,泥岩夹层中间也不会被平滑掉,这样整个泥岩夹层就完整保留下来。图3所示平滑窗口长度为m,图2所示泥岩夹层长度L,为保证窗口中心点K移动到图2中最左边的1处时,窗口边缘点Kn在图2最右边1的左边,需要满足n>=L其中m=2*n+1,这样上述权值才能正确计算,泥岩夹层才能保留下来。所以若想在平滑过程中不保留长度为L以下的泥岩夹层,则让n<L,并且Kn>K-sum(‘1’),这样在移动窗口的过程中所有泥岩夹层的权值都小于背景相权值而被平滑掉。
此外,在平滑过程中当模板中心点K位于模型的左上角时,平滑窗口的左上部分在泥岩夹层模型外部,只有中心点K右面及下面平滑窗口权值参与计算。这样平滑窗口中的权值没有被完全使用,会损失部分信息。在改进方法中对原始泥岩夹层模型进行向外扩充,比如使用5*3平滑窗口对网格数量为200*200的泥岩夹层模型进行平滑,事先在模型上面和下面各增加一层背景相,在模型左面及右面各增加两层背景相,是整个模型的网格数量达到成为202*204。这样平滑窗口在202*204模型网格内部平滑就不会出现平滑窗口溢出到模型外部的情况,使得平滑窗口中每个权值都得到计算,避免了模型边界处的平滑误差。使用其他维度的平滑窗口时,对原始泥岩夹层模型进行向外扩充的方法与此类似,在此不做赘述。
为了验证本发明的各个技术方案的有效性,下面通过对比的方式展示本发明技术方案的实施效果。需要说明的是,以下展示仅仅是对本发明技术方案的技术效果的展示,并不是对本发明保护范围的限制,凡是符合本发明技术方案本质的方案,都在本专利的保护范围之内。
在老油田SISIM方法经常被用于泥岩夹层模型的建立,一般来说由于网格大小的限制以及泥岩夹层厚度分布的特点,泥岩夹层一般呈条带状分布。如图4所示例子,图4包括泥岩分布效果401及泥岩分布直方表示402。图4中,网格大小为10*1米,泥岩夹层的厚度为1-2米,平均长度为200米,泥岩比例为10%。采用SISIM模拟,变差函数设置为水平方向变程为200米,垂直方向2米,理论模型为球状模型。从泥岩分布效果401可以看出夹层呈水平条带分布,局部存在一些离散的点或是横向延伸很短的线。泥岩分布直方表示402(由于相邻单层网格的泥岩会在垂向上会叠加形成较长的泥岩夹层,单独统计某一层网格上的泥岩长度与实际泥岩长度不符,所以本文设置统计泥岩长度规则如下:所有能够与某一网格在上下左右或对角上相邻的网格都属于同一个泥岩夹层,网格在统计过程中不会被重复计算。统计长度算法先遍历整个模型,提取出所有独立不连通的泥岩夹层,再计算每个泥岩夹层的最大宽度),计算泥岩长度平均值为133米,小于输入的期望值200米,主要原因也是因为存在一些离散的点。采用Maps方法平滑,平滑样板3*3,结果如图5所示(图5中包括Maps方法平滑后泥岩分布效果501及Maps方法平滑后泥岩分布直方表示502)。泥岩长度平均值为93米,小于平滑前的泥岩长度平均值。从图5可知,采用Maps平滑泥岩夹层模型有如下特点,对于单层的泥岩夹层,Maps方法会将其完全平滑为背景相,对于离散噪点Maps方法去除不干净,需要多次迭代,而多次迭代会产生严重的平滑效应,使平滑后模型严重失真,Maps方法不具有迭代稳定性。
针对Maps方法平滑存在不足,采用本发明实施例介绍的改进方法,按照改进平滑窗口设置规则,为平滑掉长度为30米及30米以下的泥岩夹层,设置模板9*3大小的模板[1 11 1 1 1 1 1 1;4 5 6 7 20 7 6 5 4;1 1 1 1 1 1 1 1 1]对泥岩夹层模型(图4)进行平滑,并迭代两次。平滑后的结果如图6所示(图6中包括平滑处理后泥岩分布效果601及平滑处理后泥岩分布直方表示602)。可以看出平滑后泥岩夹层最短长度为40米,所有30米以下泥岩夹层噪点都被平滑掉了。而且新方法平滑保留了原始泥岩夹层的分布特点,对应噪点及不合理的泥岩短线则做到了很好的剔除效果。平滑处理后泥岩分布直方表示602显示平均泥岩长度为164米,相对于原始岩相模型的133米,更接近期望值200米。
设置11*3大小的模板[1 1 1 1 1 1 1 1 1 1 1;4 5 6 7 8 24 8 7 6 5 4;1 1 11 1 1 1 1 1 1 1]对泥岩夹层模型(图4)进行平滑,平滑掉长度40米及40米以下的泥岩夹层,平滑结果如图7所示(图7中包括平滑处理后泥岩分布效果701及平滑处理后泥岩分布直方表示702),对比图6和图7可知:图6设置9*3大小的模板可以平滑掉长度30米及30米以下的泥岩噪点,图7设置11*3大小的模板可以平滑掉40米及40米以下泥岩噪点。所以通过改变平滑窗口的长度能够正确保留长度在指定值以上的泥岩夹层,改进算法更具有针对性。需要说明的是,上述各个实施效果的说明附图中,泥岩分布直方表示402、平滑后泥岩分布直方表示502、平滑处理后泥岩分布直方表示602及平滑处理后泥岩分布直方表示702中,横轴均表示泥岩长度,纵轴均表示某一长度的泥岩在地质模型中出现的次数,即频数。
本发明各个实施例的实现基础是通过具有处理器功能的设备进行程序化的处理实现的。因此在工程实际中,可以将本发明各个实施例的技术方案及其功能封装成各种模块。基于这种现实情况,在上述各实施例的基础上,本发明的实施例提供了一种平滑处理地质模型的装置,该装置用于执行上述方法实施例中的平滑处理地质模型的方法。参见图8,该装置包括:
通用平滑窗口获取模块801,用于根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,得到通用平滑窗口;
最终通用平滑窗口获取模块802,用于根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,并采用所述最终通用平滑窗口对地质模型进行平滑处理。
本发明实施例提供的平滑处理地质模型的装置,采用通用平滑窗口获取模块和最终通用平滑窗口获取模块,通过根据待平滑处理的地质模型的维度,设置平滑窗口的维度,再对设置好维度的平滑窗口进行权值设置,可以提高针对线性分布的地质模型的平滑处理精度,有效避免对地质模型的误平滑处理。
本发明实施例的方法是依托电子设备实现的,因此对相关的电子设备有必要做一下介绍。基于此目的,本发明的实施例提供了一种电子设备,如图9所示,该电子设备包括:至少一个处理器(processor)901、通信接口(Communications Interface)904、至少一个存储器(memory)902和通信总线903,其中,至少一个处理器901,通信接口904,至少一个存储器902通过通信总线903完成相互间的通信。至少一个处理器901可以调用至少一个存储器902中的逻辑指令,以执行如下方法:根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,得到通用平滑窗口;根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,并采用所述最终通用平滑窗口对地质模型进行平滑处理。
此外,上述的至少一个存储器902中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。例如包括:根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,得到通用平滑窗口;根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,并采用所述最终通用平滑窗口对地质模型进行平滑处理。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种平滑处理地质模型的方法,其特征在于,包括:
根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,得到通用平滑窗口;
根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,并采用所述最终通用平滑窗口对地质模型进行平滑处理。
2.根据权利要求1所述的平滑处理地质模型的方法,其特征在于,所述根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,包括:
所述平滑窗口的纵向维度与所述地质模型的纵向维度相同,所述平滑窗口的横向维度包括:
m=2n+1
其中,m为平滑窗口的横向维度,n为整数,且n满足n≥L或n<L,L为所述地质模型中待处理部分的横向维度。
3.根据权利要求1所述的平滑处理地质模型的方法,其特征在于,所述根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,包括:
所述平滑窗口的横向维度与所述地质模型的横向维度相同,所述平滑窗口的纵向维度包括:
m=2n+1
其中,m为平滑窗口的纵向维度,n为整数,且n满足n≥L或n<L,L为所述地质模型中待处理部分的纵向维度。
4.根据权利要求2所述的平滑处理地质模型的方法,其特征在于,所述根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,包括:
所述最终通用平滑窗口的中心代码点的权值为K,所述中心代码点的左侧和右侧的第i个代码点的权值为Ki,除过所述中心代码点及权值为Ki的代码点,其余代码点权值为T;
其中,i=1,…,n。
5.根据权利要求4所述的平滑处理地质模型的方法,其特征在于,还包括:
若平滑处理的过程中保留长度为L的目标物,则n满足n≥L,且K>sum(T);
若平滑处理的过程中剔除长度为L的目标物,则n满足n<L,且Kn>K-sum(T);
其中,sum(T)为所有权值T的算术和。
6.根据权利要求3所述的平滑处理地质模型的方法,其特征在于,所述根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,包括:
所述最终通用平滑窗口的中心代码点的权值为K,所述中心代码点的上侧和下侧的第i个代码点的权值为Ki,除过所述中心代码点及权值为Ki的代码点,其余代码点权值为T;
其中,i=1,…,n。
7.根据权利要求6所述的平滑处理地质模型的方法,其特征在于,还包括:
若平滑处理的过程中保留长度为L的目标物,则n满足n≥L,且K>sum(T);
若平滑处理的过程中剔除长度为L的目标物,则n满足n<L,且Kn>K-sum(T);
其中,sum(T)为所有权值T的算术和。
8.一种平滑处理地质模型的装置,其特征在于,包括:
通用平滑窗口获取模块,用于根据地质模型的横向维度及纵向维度,确定平滑窗口的横向维度和纵向维度,得到通用平滑窗口;
最终通用平滑窗口获取模块,用于根据地质模型中需要处理的代码点,设置所述通用平滑窗口的子窗口的权值,得到最终通用平滑窗口,并采用所述最终通用平滑窗口对地质模型进行平滑处理。
9.一种电子设备,其特征在于,包括:
至少一个处理器、至少一个存储器、通信接口和总线;其中,
所述处理器、存储器、通信接口通过所述总线完成相互间的通信;
所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令,以执行如权利要求1至7任一项所述的方法。
10.一种非暂态计算机可读存储介质,其特征在于,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行如权利要求1至7中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811588306.8A CN109712238B (zh) | 2018-12-25 | 2018-12-25 | 平滑处理地质模型的方法及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811588306.8A CN109712238B (zh) | 2018-12-25 | 2018-12-25 | 平滑处理地质模型的方法及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109712238A true CN109712238A (zh) | 2019-05-03 |
CN109712238B CN109712238B (zh) | 2023-06-06 |
Family
ID=66257340
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811588306.8A Active CN109712238B (zh) | 2018-12-25 | 2018-12-25 | 平滑处理地质模型的方法及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109712238B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113986896A (zh) * | 2021-10-14 | 2022-01-28 | 中冶南方工程技术有限公司 | 一种渣皮厚度数据滤波的方法和装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060269139A1 (en) * | 2005-03-29 | 2006-11-30 | Total S.A. | Method and program for finding discountinuities |
US20110205844A1 (en) * | 2010-02-22 | 2011-08-25 | Landmark Graphics Corporation, A Haliburton Company | Systems and Methods for Modeling 3D Geological Structures |
US20110211761A1 (en) * | 2010-02-26 | 2011-09-01 | Chevron U.S.A. Inc. | Surface smoothing within an earth model of a geological volume of interest |
US20120215513A1 (en) * | 2009-11-12 | 2012-08-23 | Branets Larisa V | Method and Apparatus For Reservoir Modeling and Simulation |
US20130262051A1 (en) * | 2012-04-03 | 2013-10-03 | Brandon David Plost | Ordered multipoint geostatistics simulation using non-symmetric search mask |
US20160103245A1 (en) * | 2014-10-09 | 2016-04-14 | Chevron U.S.A. Inc. | Conditioning of object or event based reservior models using local multiple-point statistics simulations |
US20170183962A1 (en) * | 2014-07-25 | 2017-06-29 | Landmark Graphics Corporation | Identifying Formation Layer Boundaries on Well Log Measurements |
CN108693556A (zh) * | 2017-04-11 | 2018-10-23 | 中国石油化工股份有限公司 | 一种构造约束的速度平滑方法和系统 |
-
2018
- 2018-12-25 CN CN201811588306.8A patent/CN109712238B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060269139A1 (en) * | 2005-03-29 | 2006-11-30 | Total S.A. | Method and program for finding discountinuities |
US20120215513A1 (en) * | 2009-11-12 | 2012-08-23 | Branets Larisa V | Method and Apparatus For Reservoir Modeling and Simulation |
US20110205844A1 (en) * | 2010-02-22 | 2011-08-25 | Landmark Graphics Corporation, A Haliburton Company | Systems and Methods for Modeling 3D Geological Structures |
US20110211761A1 (en) * | 2010-02-26 | 2011-09-01 | Chevron U.S.A. Inc. | Surface smoothing within an earth model of a geological volume of interest |
US20130262051A1 (en) * | 2012-04-03 | 2013-10-03 | Brandon David Plost | Ordered multipoint geostatistics simulation using non-symmetric search mask |
US20170183962A1 (en) * | 2014-07-25 | 2017-06-29 | Landmark Graphics Corporation | Identifying Formation Layer Boundaries on Well Log Measurements |
US20160103245A1 (en) * | 2014-10-09 | 2016-04-14 | Chevron U.S.A. Inc. | Conditioning of object or event based reservior models using local multiple-point statistics simulations |
CN108693556A (zh) * | 2017-04-11 | 2018-10-23 | 中国石油化工股份有限公司 | 一种构造约束的速度平滑方法和系统 |
Non-Patent Citations (6)
Title |
---|
喻思羽;李少华;何幼斌;段太忠;廉培庆;陶金雨;包兴;: "基于样式降维聚类的多点地质统计建模算法", no. 11, pages 1403 - 1409 * |
尹艳树等: "基于信息度的储层建模后处理方法", no. 6, pages 889 - 893 * |
徐丽强等: "辫状河三角洲前缘储层隔夹层表征及剩余油预测——以彩南油田彩9 井区三工河组为例", no. 4, pages 10 - 18 * |
方晓刚等: "基于开源软件的序贯指示沉积相建模后处理方法原理及改进", pages 118 - 119 * |
权勃;侯亚伟;张东;彭琴;孙恩慧;谢鹏程;吴波;: "小尺度等效表征技术在油藏数值模拟研究中的应用――以渤海C油田馆陶组辫状河储层为例", no. 15, pages 53 - 58 * |
李兆亮;潘懋;韩大匡;刘文岭;王经荣;胡水清;刘培刚;: "储层精细构造模型三维网格化技术", no. 26, pages 60 - 65 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113986896A (zh) * | 2021-10-14 | 2022-01-28 | 中冶南方工程技术有限公司 | 一种渣皮厚度数据滤波的方法和装置 |
CN113986896B (zh) * | 2021-10-14 | 2024-07-23 | 中冶南方工程技术有限公司 | 一种渣皮厚度数据滤波的方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109712238B (zh) | 2023-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110136072B (zh) | 点云噪声的去除方法、去噪系统、计算机设备及存储介质 | |
KR101555426B1 (ko) | 지형 렌더링 방법 및 장치 | |
US8635052B1 (en) | Method and appartus for transforming a stratigraphic grid | |
US8725478B2 (en) | Reservoir upscaling method with preserved transmissibility | |
CN104063903A (zh) | 三维实体模型的四面体网格生成方法和装置 | |
EA013941B1 (ru) | Способ и устройство для укрупнения слоев в геологическом моделировании | |
CN114332291B (zh) | 一种倾斜摄影模型建筑物外轮廓规则提取方法 | |
CN107204042B (zh) | 基于形态完备性的非均质岩心三维结构重建算法 | |
Suchde | Conservation and accuracy in meshfree generalized finite difference methods | |
Huo et al. | Novel quadtree algorithm for adaptive analysis based on cell-based smoothed finite element method | |
CN104915928A (zh) | 一种基于本征正交分解的速度场坏矢量识别和修正方法 | |
CN109712238A (zh) | 平滑处理地质模型的方法及设备 | |
CN105487116B (zh) | 层面模型建立方法 | |
CN114202642A (zh) | 三维三角网模型中超大空洞的填补方法及设备 | |
CA3019487C (en) | Systems and methods for unfaulting point clouds | |
US11551417B2 (en) | Meshing operation for representing a 3D printed object | |
CN103824331A (zh) | 六面体网格模型的迭代生成方法 | |
CN112016956B (zh) | 基于bp神经网络的矿石品位估值方法及装置 | |
CN112419493B (zh) | 页岩储层三维属性模型建立方法及装置 | |
CN109471161B (zh) | 地质构造属性剖面的建立方法和装置 | |
WO2013116859A1 (en) | Computer process for determining best-fitting materials for constructing architectural surfaces | |
WO2021231170A1 (en) | Systems and methods for generating vertical and lateral heterogeneity indices of reservoirs | |
US10386539B2 (en) | Method of modelling a subsurface volume | |
KR20030017888A (ko) | 표면 법선 기반의 삼각형 메쉬 영역화 장치 및 그 방법 | |
CN110838181B (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 | ||
GR01 | Patent grant |