一种河道识别方法与流程

1.本发明属于油气地质技术领域，具体涉及一种河道识别方法。


背景技术：

2.国内外对于河道划分普遍采用的都是地质作图方法，只能确定主方向，而对于河道的细节部分刻画无法做到。同时这个主方向对于小区块开发而言，可能指导意义甚微，但是依然采取了该河道方向，这对于单砂体建模来说就是灾难性错误。河道形态进行人为确定，随意性太大，且不同的地质工作人员会在作图时对单河道有不同的河道取向，导致地质研究既不标准也不精确；人为确定河道的方法更加粗略，不足以指导小区域的砂岩气藏开发工作；对于不同的区块，不同的层位，应该有可执行标准进行河道识别，进而进行河道刻画。
3.目前国内都是依据区块的特征人为进行河道识别，太过主观，且没有统一的标准，不精确。因此，对各个区域、各个层位河道识别方法亟需改进。


技术实现要素：

4.本发明提供的一种河道识别方法目的一是克服现有技术中对于河道划分普遍采用的都是地质作图方法，只能确定主方向，而对于河道的细节部分刻画无法做到的问题；目的二是克服现有技术中依据区块的特征人为进行河道识别，太过主观，且没有统一的标准，不精确问题。
5.为此，本发明提供了一种河道识别方法，
6.1)对目标砂层的砂体厚度平面数值进行概率分析和前处理，去掉部分数据，剩余的为砂体厚度有效数据，通过对砂体厚度有效数据进行分析，筛选出河道主方向；所述部分数据为被砂体厚度值为0值所包围的数据；
7.2)对目标区井点处的砂体厚度数据进行方差求取；
8.3)对步骤2)的求取的方差数值进行比较，筛选出方差大的数据进行分析，结合河道主方向进行河道精细刻画；
9.4)对步骤3)刻画的河道内砂体纯度进行分类统计；
10.5)根据步骤4)的分类统计结果，结合河道精细刻画图综合刻画单期单河道，完成河道最终刻画图。
11.进一步的，所述步骤1)对目标砂层的砂体厚度平面数值进行概率分析和前处理是对砂体厚度值的期望值进行分析和处理。
12.进一步的，所述对砂体厚度值的期望值进行分析和处理的方法为：当砂体厚度值大于0.9倍的期望值时，该砂体厚度值作为砂体厚度有效数据；当砂体厚度值小于0.9倍的期望值时，需删除该砂体厚度值。
13.进一步的，所述通过对砂体厚度有效数据进行分析，筛选出河道主方向包括如下步骤：对砂体厚度有效数据分析，找出砂体厚度值高且排名前三位的数值由高到低依次进
行连接，形成一条折线状的河道即为河道主方向。
14.进一步的，所述步骤2)方差求取按照河道主方向单位长度进行对目标区井点处的砂体厚数据进行方差求取。
15.进一步的，所述步骤3)方差大的数据为方差大于50的数据。
16.进一步的，所述步骤4)砂体纯度值为ps＝gr*5/n，gr为伽玛值，ps为砂体纯度值，n为该某一段砂体厚度内的数据点数，这一段砂体中间的隔夹层不大于0.4m,这样的加权平均值为纯度值。
17.进一步的，所述步骤5)刻画单期单河道的方法为：定义河道指示值，在河道精细刻画图上标识出所有河道指示值，然后对河道指示值数据进行差值求取，取差值最小的方向为主方向，接着沿主方向将河道指示值数据按照相邻关系连接，连接出一条不交叉的线条即为这一期的主河道。
18.进一步的，所述步骤5)河道指示值r＝hi*ps，r为河道指示值，hi为砂体厚度权重值，ps为砂体纯度值。
19.进一步的，所述步骤5)对河道指示值数据进行差值求取要求对河道指示值数据进行0度、45度、90度、135度四个方向的差值求取。
20.本发明的有益效果：本发明提供的这种河道识别方法，通过对目标砂层的砂体厚度的数值进行分析和处理，筛选出河道主方向，然后依照河道主方向单位长度对目标区井点处的砂体厚数据进行方差求取，筛选出方差大的数据进行分析，结合河道主方向进行河道精细刻画，并对刻画的河道内砂体纯度进行分类统计，结合河道精细刻画图综合刻画单期单河道，完成河道最终刻画图；建立了一套适合于河道识别的方法，形成了一种实用可靠的精确、客观、操作性强的识别方法，此方法在鄂尔多斯盆地东部应用半年多来效果较好，克服了依据区块的特征人为进行河道识别，太过主观，且没有统一的标准，不精确的问题，本方法既能确定河道主方向，同时还能刻画河道的细节部分，可靠性和精确度高，客观性和操作性强。
附图说明
21.以下将结合附图对本发明做进一步详细说明。
22.图1是实施例2一幅砂体厚度平面数值分布图；
23.图2是图1中数据剔除后的砂体厚度平面数值分布图；
24.图3是主河道方向图；
25.图4是河道精细刻画图；
26.图5是河道最终刻画图。
具体实施方式
27.实施例1：
28.一种河道识别方法，包括如下步骤：
29.1)对目标砂层的砂体厚度平面数值进行概率分析和前处理，去掉部分数据，剩余的为砂体厚度有效数据，通过对砂体厚度有效数据进行分析，筛选出河道主方向；所述部分数据为被砂体厚度值为0值所包围的数据；
30.2)对目标区井点处的砂体厚度数据进行方差求取；
31.3)对步骤2)的求取的方差数值进行比较，筛选出方差大的数据进行分析，结合河道主方向进行河道精细刻画；
32.4)对步骤3)刻画的河道内砂体纯度进行分类统计；
33.5)根据步骤4)的分类统计结果，结合河道精细刻画图综合刻画单期单河道，完成河道最终刻画图。
34.本发明建立了一套适合于河道识别的方法，形成了一种实用可靠的精确、客观、操作性强的识别方法，此方法在鄂尔多斯盆地东部应用半年多来效果较好，克服了依据区块的特征人为进行河道识别，太过主观，且没有统一的标准，不精确的问题，本方法既能确定河道主方向，同时还能刻画河道的细节部分，可靠性和精确度高，客观性和操作性强。
35.进一步的，所述步骤1)对目标砂层的砂体厚度平面数值进行概率分析和前处理是对砂体厚度值的期望值进行分析和处理。所述期望值的求取为e(x)＝(h1+h2+
……
+h
n
)/n，h为砂体厚度，e(x)为期望值。
36.进一步的，所述对砂体厚度值的期望值进行分析和处理的方法为：当砂体厚度值大于0.9倍的期望值时，该砂体厚度值作为砂体厚度有效数据；当砂体厚度值小于0.9倍的期望值时，需删除该砂体厚度值。
37.进一步的，所述通过对砂体厚度有效数据进行分析，筛选出河道主方向包括如下步骤：对砂体厚度有效数据分析，找出砂体厚度值高且排名前三位的数值由高到低依次进行连接，形成一条折线状的河道即为河道主方向。
38.进一步的，所述步骤2)方差求取按照河道主方向单位长度进行对目标区井点处的砂体厚数据进行方差求取。
39.进一步的，所述步骤3)方差大的数据为方差大于50的数据。方差小于50的数据，可以认为数据点之间的河道方向没变化，这些数据点归为一类，可以删除，以河道主方向的方向值代表其方向；对于方差大于50的数据，代表河道在此处有大的变化，需要仔细刻画；将该类数据点全部标出来，单独分析，进行河道刻画；并将这些砂体值h*9＝hi，hi定义为砂体厚度权重值。
40.进一步的，所述步骤4)砂体纯度值为ps＝gr*5/n，gr为伽玛值，gr为测井时测出来的一个数据，ps为砂体纯度值，n为该某一段砂体厚度内的数据点数，这一段砂体中间的隔夹层不大于0.4m,这样的加权平均值为纯度值。
41.进一步的，所述步骤5)刻画单期单河道的方法为：定义河道指示值，在河道精细刻画图上标识出所有河道指示值，然后对河道指示值数据进行差值求取，取差值最小的方向为主方向，接着沿主方向将河道指示值数据按照相邻关系连接，连接出一条不交叉的线条即为这一期的主河道。
42.进一步的，所述步骤5)所述河道指示值r＝hi*ps，r为河道指示值，hi为砂体厚度权重值，ps为砂体纯度值。
43.进一步的，所述步骤5)所述对河道指示值数据进行差值求取要求对河道指示值数据进行0度、45度、90度、135度四个方向的差值求取。
44.实施例2：
45.现结合鄂尔多斯盆地东部的河道刻画来对本发明进行具体说明。
46.步骤1
47.1)图1是一幅砂体厚度平面数值分布图，其中对被0值包围的数据地质认识上认为不会出现砂体，因此删除掉这些数据，删掉这些数据之后的砂体厚度平面数值分布图如图2所示；，
48.2)对目标区域内所有砂体厚度值(参见表1)进行期望值求取e(x)＝(h1+h2+
……
+hn)/n＝7.976，精确到个位，取8；
49.表1砂体厚度有效数据平均值统计表
[0050][0051]
3)取大于0.9e(x)的数据作为有效数据，将其他砂体厚度数据剔除，不作为有效河道砂体数据，见表2和图3。
[0052]
表2二次剔除后的砂厚数据统计表
[0053][0054]
4)对砂体厚度进行平面分布数据找出厚度峰值26，以及砂体厚度值高值的前三位26、23、22，对于高值前三位26、23、22依次进行连接，会形成一个折线状的河道，这是主河道的大致走势即河道主方向。
[0055]
步骤2
[0056]
对于各井点处的砂体厚度数据进行依照河道主方向单位长度方差求取：
[0057]
(26+23+22)/3＝23.7，23.7作为方差基值，然后求各个数据的方差，如下表3所示：
[0058]
表3方差求取表
[0059]
序号砂体厚度基值砂体厚度
‑
平均值(砂体厚度
‑
平均值)的平方12323.7
‑
0.70.4922023.7
‑
3.713.6932623.72.35.2941123.7
‑
12.7161.295923.7
‑
14.7216.0961423.7
‑
9.794.0971823.7
‑
5.732.4981023.7
‑
13.7187.6991723.7
‑
6.744.89102223.7
‑
1.72.89111923.7
‑
4.722.09121223.7
‑
11.7136.89131123.7
‑
12.7161.29141623.7
‑
7.759.29151223.7
‑
11.7136.89161023.7
‑
13.7187.69171823.7
‑
5.732.4918823.7
‑
15.7246.49191823.7
‑
5.732.49
[0060]
步骤3
[0061]
对于表3的方差进行比较，方差较大的地方仔细核对，在河道主方向的控制下，按照主方向进行方差求取。
[0062]
方差小于50的数据，可以认为数据点之间的河道方向没变化，这些数据点归为一类，可以删除(见表4)，并以主方向的方向值代表其方向。
[0063]
表4第三次数据剔除之后统计表
[0064][0065]
对于方差大于50，代表河道在此处有大的变化，将该类数据点全部标出来，单独分析(见表5)，进行河道刻画(见图4)。并将这些砂体值h*9＝hi，hi定义为砂体厚度权重值。
[0066]
表5hi计算表
[0067]
序号砂体厚度hhi11199298131412641090512108611997161448121089109010872
[0068]
步骤4
[0069]
对于河道内砂体纯度进行分类统计(见表6)
[0070]
定义砂体纯度值ps＝gr*5/n
[0071]
n为该某一段砂体厚度内的数据点数，这一段砂体中间的隔夹层不大于0.4m,这样的加权平均值为纯度值。
[0072]
表6砂体纯度计算表
[0073]
序号砂体厚度higr均值pshi
×
ps1119941205202952981442201782031412638190239404109040200180005121083919521060611994221020790
71614436180259208121083718519980910903919517550108724723516920
[0074]
步骤5
[0075]
综合刻画单期单河道
[0076]
定义河道指示值r＝hi*ps，乘积值见表5。在平面上标识出所有数据值，然后对数据进行0度、45度、90度、135度四个方向的差值求取，取差最小的方向为主方向。接着沿主方向将河道指示值数据按照相邻关系连接，连接处一条不交叉的线条，这个线条即为这一期次的主河道(见图5)。
[0077]
本发明的描述中，需要理解的是，若有术语指示的方位或位置关系为基于附图所示的方位或位置关系，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此附图中描述位置关系的用语仅用于示例性说明，不能理解为对本专利的限制。
[0078]
以上例举仅仅是对本发明的举例说明，并不构成对本发明的保护范围的限制，凡是与本发明相同或相似的设计均属于本发明的保护范围之内。