一种属性等值线图绘制方法及装置制造方法
【专利摘要】本发明提供一种属性等值线图绘制方法及装置。所述属性等值线图网格化方法包括步骤：选定一个包括权系数因子的权系数函数；采用蒙特卡洛随机抽样方法求解所述权系数函数的权系数因子；将所述求出的权系数因子代入所述权系数函数中求出权系数；将已知离散点的值与其权系数相乘的平均值作为待定点的值；计算出所有待定点的值，绘制属性等值线图。本发明提供的属性等值线图绘制方法及装置将储层的非均质性考虑在内，能更加准确的反映出地下储层的孔隙度、渗透率等储层物性参数的分布情况。
【专利说明】一种属性等值线图绘制方法及装置
【技术领域】
[0001]本发明方法涉及石油、地矿勘探开发中属性地质建模领域，特别涉及一种属性等值线图绘制方法及装置。
【背景技术】
[0002]在石油、地矿勘探开发工作中，需要将采集并处理后的地震数据绘制成地震构造图，以便观察和分析地层构造。所述的地震构造图包括属性等值线图，所述属性等值线图可以是将采集的少数某一地质属性的离散数据进行网格化自动追踪并绘制出的等值线图。所述的网格化通常是指在含有坐标信息并且规则分布的矩形网格上利用某一网格节点周围的离散数据点采用拟合或者插值的方法计算出该网格节点值。
[0003]在所述地震属性的数据采集过程中，由于地形、环境、开发成本等因素的约束，通常采取分区域段测量少数几十个甚至十几个离散数据。然后，在室内对采集的离散数据进行网格化处理，获取地质的属性等值线图。现有技术中，常用的网格化方法有距离反比插值法、曲面拟合插值法、克里金法等。其中所述距离反比插值法的操作过程通常包括:用某一区域内待定点Z (X, y, Z)到该区域η个已知离散点Zi (xi,yi,zi)(其中i=l,2,3,…η)的距离di(i=l，2，3，…η)倒数作为自变量的，并用所述自变量的m (通常在计算渗透率和孔隙度时，m取2或者3)次方作为权系数来加权平均求解该待定点Z。图1是所述距离反比插值法的计算示意图。其所述距离反比例插值法的计算公式可以为:
【权利要求】
1.一种属性等值线图绘制方法，其特征在于，包括以下处理步骤: S1:选定一个包括权系数因子的权系数函数； 52:采用蒙特卡洛随机抽样方法求解所述权系数函数的权系数因子； 53:将所述求出的权系数因子代入所述权系数函数中求出已知离散点的权系数； 54:将已知的离散点的值与其相应的权系数相乘的平均值作为待定点的值； 55:计算出所有待定点的值，绘制属性等值线图。
2.如权利要求1所述的一种属性等值线图绘制方法，其特征在于，SI中所述选取的权系数函数包括参数:已知离散点至待定点的距离和已知离散点相对于待定点的方位。
3.如权利要求1所述的一种属性等值线图绘制方法，其特征在于，所述采用蒙特卡洛随机抽样方法求解所述权系数函数的权系数因子的方法包括以下处理步骤: S201:选取所述权系数因子的取值区间及所述取值区间的分布函数； S202:做出所述分布函数的累积概率曲线； S203:随机产生K个O至I之间的随机样本数作为所述累积概率曲线上的概率值，并用所述随机产生的概率值从所述累积概率曲线上反求得相应的样本权系数因子，K为自然数； S204:取所述求得的K个样本权系数因子的期望值作为所述权系数函数的权系数因子。
4.如权利要求3所述的一种属性等值线图绘制方法，其特征在于，所述随机产生的K个0至1之间的随机样本数作为所述累积概率曲线上的概率值中K的取值范围为:
10000 ≤ K ≤ 200000。
5.如权利要求1所述的一种属性等值线图绘制方法，其特征在于，所述选取的权系数函数为:
6.如权利要求1所述的一种属性等值线图绘制方法，其特征在于，所述选取的权系数函数为:
7.如权利要求1所述的一种属性等值线图绘制方法，其特征在于，S4之后还包括步骤: S400:将已经求出的待定点的值作为所述已知的离散点的值，用于计算下一个待定点。
8.如权利要求1-7任意一项所述的属性等值线图绘制方法，其特征在于，所述已知的离散点为渗透率或孔隙度中的一种。
9.一种属性等值线图绘制装置，其特征在于，该装置包括权系数函数设置模块、权系数因子计算模块、权系数函数计算模块、待定点计算模块、等值线计算模块，其中: 所述权系数函数设置模块，用于设置一个包括权系数因子的权系数函数； 所述权系数因子计算模块，用于采用蒙特卡洛随机抽样方法求解所述权系数函数的权系数因子； 所述权系数函数计算模块，以用于将所述权系数因子计算模块求出的权系数因子代入所述权系数函数中求出已知离散点的权系数。 所述待定点计算模块，用于将已知的离散点的值与其相应的权系数相乘的平均值作为待定点的值； 所述等值线计算模块，用于计算出所有待定点的值，构建属性等值线图。
10.如权利要求9所述的一种属性等值线图绘制装置，其特征在于，所述权系数函数设置模块中所述设置的权系数函数包括参数:已知离散点至待定点的距离和已知离散点相对于待定点的方位。
11.如权利要求9所述的一种属性等值线图绘制装置，其特征在于，权系数因子计算模块以包括分布函数设置模块、累积概率曲线计算模块、样本权系数计算模块、权系数计算模块，其中: 所述分布函数设置模块，用于选取所述权系数因子的取值区间及所述取值区间的分布函数； 所述累积概率曲线计算模块，用于做出所述分布函数的累积概率曲线； 所述样本权系数计算模块，用于随机产生K个O至I之间的随机样本数作为所述累积概率曲线上的概率值，并用所述随机产生的概率值从所述累积概率曲线上反求得相应的样本权系数因子，K为自然数； 所述权系数计算模块，用于取所述求得的K个样本权系数因子的期望值作为所述权系数函数的权系数因子。
12.如权利要求10所述的一种属性等值线图绘制装置，其特征在于，所述样本权系数计算模块随机产生K个O至I之间的随机样本数作为所述累积概率曲线上的概率值中，K的取值范围为:10000 ^ K ^ 200000。
13.如权利要求9所述的一种属性等值线图绘制装置，其特征在于，该装置还包括已知离散点转化模块，所述已知离散点转化模块用于将所述待定点计算模块已经求出的待定点的值作为所述已知的离散点的值，用于计算下一个待定点。
【文档编号】G01V1/32GK103901475SQ201410126087
【公开日】2014年7月2日 申请日期:2014年3月31日 优先权日:2014年3月31日
【发明者】任殿星 申请人:中国石油天然气股份有限公司