将测井数据离散到储层网格模型的方法和装置制造方法
【专利摘要】本发明提供一种将测井数据离散到储层网格模型的方法和装置，涉及地质属性建模【技术领域】。所述储层网格模型为角点网格，该方法包括下列步骤：S1、获取井轨迹穿越的网格单元；S2、计算所述井轨迹同所述网格单元的交点，并利用所述交点将所述井轨迹划分为多个区块；S3、对落入每一区块的所述井轨迹上的测井曲线属性值进行加权平均，求得属性均值，并将所述属性均值赋给对应的网格单元。本发明根据测井数据值，用尽可能少的时间为井轨迹所穿过的网格赋值，从而减少建模人员的等待时间，提高了工作效率。
【专利说明】将测井数据离散到储层网格模型的方法和装置
【技术领域】
[0001]本发明涉及地质属性建模【技术领域】，具体涉及一种将测井数据离散到储层网格模型的方法和装置。
【背景技术】
[0002]地质体内部含有多种反映地层岩性、油气资源分布等特性的参数，如孔隙度、渗透率、含有饱和度、沉积相等，地质体物性建模就是利用空间插值法，如Kriging、序贯高斯等，求解地质体内部物理、化学属性参数模型。利用地质体物性建模成果，可以有效的发现储层，并对油气资源分布情况进行准确预测，从而指导油田开发。考虑到地层成因，地质体中某一区域的物性，与同一层物性的关系远远大于该区域的上层或下层。为了更精确的表达地质体物性，实际建模过程中，在对地质体细分小层的基础上，会进一步细分，对细分后的每一层，利用由测井曲线数据得到的关于该层的物性信息，使用空间差值算法进行物性建模。将测井数据离散到储层网格模型是利用测井数据，为井穿过的储层网格单元赋值。被赋值的网格单元，在进行物性建模时，作为已知量，进而确定其他未被测量过区域的属性值。可见，它是进行地质体物性建模的基础。
[0003]目前，角点网格是被广泛使用的描述储集层的网格模型。作为结构化网格模型，它具有易于实现的优点。高度的灵活性使得它能够很好的表征储集层地质特征，比如非均质性以及断层。
[0004]然而，由于角点网格不使用笛卡尔坐标系，在离散测井数据时非常复杂繁琐。再加上为了满足对建模区域精细描述的需求，现在的网格划分都非常细密。故而，一种快速地将测井数据离散到储层模型的方法被迫切的需要。
【发明内容】
[0005](一)解决的技术问题
[0006]针对现有技术的不足，本发明提供了一种将测井数据离散到储层网格模型的方法和装置，根据测井数据值，用尽可能少的时间为井轨迹所穿过的网格赋值，从而减少建模人员的等待时间，提闻工作效率。
[0007](二)技术方案
[0008]为实现以上目的，本发明通过以下技术方案予以实现:
[0009]一种将测井数据离散到储层网格模型的方法，所述储层网格模型为角点网格，包括下列步骤:
[0010]S1、获取井轨迹穿越的网格单元；
[0011]S2、计算所述井轨迹同所述网格单元的交点，并利用所述交点将所述井轨迹划分为多个区块；
[0012]S3、对落入每一区块的所述井轨迹上的测井曲线属性值进行加权平均，求得属性均值，并将所述属性均值赋给对应的网格单元。[0013]优选的，步骤SI包括:
[0014]S11、根据储层网格模型的网格厚度信息设置井轨迹中待测点选取的间隔；
[0015]S12、按照所述间隔顺次取出待测点；
[0016]S13、找到每一个待测点所在的网格单元；
[0017]S14、将所述网格单元进行整理与合并。
[0018]优选的，步骤Sll包括:
[0019]S111、根据井头坐标，以及储层网格模型的边界信息和网格长宽设置信息，估算出位于井头下方的网格单元；
[0020]SI 12、计算井头下方网格单元的平均厚度；
[0021 ] S113、结合井轨迹数据的采样间隔，计算出待测点选取的间隔；
[0022]其中，步骤S13包括:利用参考索引找出对应网格单元，以网格单元为中心，采用临近搜索法选取网格单元，对待测点和选取的网格单元进行包含性测试；
[0023]其中，步骤S14包括:通过网格单元的索引信息，利用折半查找法，找出未被测出的网格单元；对于多个待测点落入同一网格单元的情形，只记录第一个落入该网格单元的待测点信息，与最后一个落入该网格单元的待测点信息。
[0024]优选的，步骤S2中计算所述井轨迹同所述网格单元的交点包括步骤:
[0025]根据网格单元索引信息确定单元内同井轨迹有交点的面；选取面两端的待测点，计算两待测点间的井轨迹与网格单元面的交点。
[0026]优选的，计算两待测点间的井轨迹与网格单元面的交点时使用快速逼近法，实现方式是:
[0027]选取两相邻井轨迹数据点组成直线，计算同网格单元表面的交点，计算所述交点与数据点的距离；
[0028]在井轨迹上选取距离相等的数据点，并以此数据点为基准，选取与之相邻的下一个井轨迹数据点，利用重新得到的两相邻井轨迹数据点，进行下一次的求交点运算；
[0029]当交点位于选取的两相邻井轨迹数据点之间时，结束运算。
[0030]优选的，步骤S2中将井轨迹划分为多个区块时，将两个相邻的井与网格单元面的交点之间的区域，作为一个区块，每一区块对应储层网格模型的一个网格单元。
[0031]优选的，步骤S3中计算落入每一区块的测井曲线属性值均值时，得到唯一的数值，并以此数值作为区块对应网格单元的数值。
[0032]本发明还提供了一种将测井数据离散到储层网格模型的装置，包括:
[0033]网格单元查找模块，用于找出井轨迹所穿越的网格单元，
[0034]测井数据划分模块，用于计算井轨迹同网格单元面的交点，并对测井数据划分区块；
[0035]属性计算模块，用于对落入每一区块的所有测井曲线属性值计算加权平均值，并为对应的网格单元赋值；
[0036]交互模块，对储层网格模型和测井数据进行访问。
[0037]优选的，所述网格单元查找模块包括:
[0038]待测点间隔发生器，用于产生待测点选取间隔；
[0039]待测点抽样采集器，用于根据指定间隔选取待测点的抽样采集；[0040]网格单元搜索装置，用于准确、快速判断待测点所处的网格单元的测试设备和网格单元搜索；
[0041]网格单元内判定器，用于判断待测点是否位于网格单元内；
[0042]其中，所述测井数据划分模块包括:
[0043]交点求解器，用于快速计算井轨迹与网格单元面的交点；
[0044]测井数据分块器，用于对测井数据划分区块；
[0045]其中，所述属性计算模块包括:
[0046]属性求解器，用于计算落入同一区块的测井曲线属性加权平均值；
[0047]属性赋值器，用于将区块属性值设为对应网格单元属性值；
[0048]其中，所述交互模块包括:
[0049]储层模型访问子模块，用于根据储层网格模型的三维索引值，获取储层网格模型中的特定网格单元的信息，且将指定物性属性值设置到对应网格单元中；
[0050]测井数据访问子模块，用于提供井轨迹类型数据和测井曲线类型数据的访问接□。
[0051]优选的，所述测井数据访问子模块按照采样点在对应井轨迹类型数据和测井曲线类型数据中的索引，获得采样点的特定信息；其中所述井轨迹数据包含:测量深度信息、大地坐标系下的三维坐标值和垂直深度；所述测井曲线数据包括:测量深度信息、孔隙度值、自然伽马值和沉积相相位值。
[0052](三)有益效果
[0053]本发明提供了一种将测井数据离散到储层网格模型的方法和装置，通过获取井轨迹穿越的网格单元，计算所述井轨迹同所述网格单元的交点，并利用所述交点将所述井轨迹划分为多个区块；进而对落入每一区块的所述井轨迹上的测井曲线属性值进行加权平均，求得属性均值，并将所述属性均值赋给对应的网格单元。
[0054]本发明可以快捷的利用测井数据计算储层网格模型中被井穿过的网格单元的属性值，减少了建模人员的等待时间，提供工作效率。
[0055]本发明可以准确的利用测井数据计算储层网格模型中相应单元的属性值，克服了角点网格体系下难以计算空间位置的缺点，提高了建模结果的准确性。大大增加了石油储层精细描述系统的可信度，为油井采收、油田二次开发等工作提供可靠的支撑，从而大幅节约油田开采成本。
【专利附图】
【附图说明】
[0056]为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0057]图1为本发明实施例中将测井数据离散到储层网格模型的方法的流程示意图；
[0058]图2为双线性曲面的示意图；
[0059]图3为本发明实施例中将测井数据离散到储层网格模型的装置的结构示意图；
[0060]图4为网格单元查找模块的实施流程图；[0061]图5为计算角点网格单元平均厚度的示意图；
[0062]图6为计算射线与双线性曲面交点方法的示意图；
[0063]图7为判断待测点是否位于指定网格单元方法的示意图；
[0064]图8为测井数据划分模块的实施流程图；
[0065]图9为用快速逼近法求解井轨迹与网格单元面交点过程的示意图；
[0066]图10为属性计算模块的实施流程图；
[0067]图11为当没有测井曲线数据点落入某一区块时确定网格单元属性值的方法示意图。
【具体实施方式】
[0068]为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0069]实施例1:
[0070]如图1所示，本发明实施例提供了一种将测井数据离散到储层网格模型的方法，所述储层网格模型为角点网格，包括下列步骤:
[0071]S1、获取井轨迹穿越的网格单元；
[0072]S2、计算所述井轨迹同所述网格单元的交点，并利用所述交点将所述井轨迹划分为多个区块；
[0073]S3、对落入每一区块的所述井轨迹上的测井曲线属性值进行加权平均，求得属性均值，并将所述属性均值赋给对应的网格单元。
[0074]本发明实施例可以快捷的利用测井数据计算储层网格模型中被井穿过的网格单元的属性值，减少了建模人员的等待时间，提供工作效率。可以准确的利用测井数据计算储层网格模型中相应单元的属性值，克服了角点网格体系下难以计算空间位置的缺点，提高了建模结果的准确性。大大增加了石油储层精细描述系统的可信度，为油井采收、油田二次开发等工作提供可靠的支撑，从而大幅节约油田开采成本。
[0075]下面对本发明实施例进行详细的说明:
[0076]优选的，步骤SI包括:
[0077]S11、根据储层网格模型的网格厚度信息设置井轨迹中待测点选取的间隔；
[0078]S12、按照所述间隔顺次取出待测点；
[0079]S13、找到每一个待测点所在的网格单元；
[0080]S14、将所述网格单元进行整理与合并。
[0081]优选的，步骤Sll包括:
[0082]S111、根据井头坐标，以及储层网格模型的边界信息和网格长宽设置信息，估算出位于井头下方的网格单元；
[0083]SI 12、计算井头下方网格单元的平均厚度；
[0084]S113、结合井轨迹数据的采样间隔，计算出待测点选取的间隔。
[0085]优选的，步骤S13包括:利用参考索引找出对应网格单元，以网格单元为中心，采用临近搜索法选取网格单元，对待测点和选取的网格单元进行包含性测试。
[0086]优选的，步骤S14包括:
[0087]通过网格单元的索引信息，利用折半查找法，找出未被测出的网格单元；对于多个待测点落入同一网格单元的情形，只记录第一个落入该网格单元的待测点信息，与最后一个落入该网格单元的待测点信息。
[0088]优选的，步骤S2中计算所述井轨迹同所述网格单元的交点包括步骤:
[0089]根据网格单元索引信息确定单元内同井轨迹有交点的面；选取面两端的待测点，计算两待测点间的井轨迹与网格单元面的交点。
[0090]优选的，计算两待测点间的井轨迹与网格单元面的交点时使用快速逼近法，实现方式是:
[0091]选取两相邻井轨迹数据点组成直线，计算同网格单元表面的交点，计算所述交点与数据点的距离，在井轨迹上选取距离相等的数据点，以此数据点为基准，选取与之相邻的下一个井轨迹数据点。利用重新得到的两相邻井轨迹数据点，进行下一次的求交点运算。直到交点位于选取的两相邻井轨迹数据点之间，算法结束。
[0092]优选的，步骤S2中将井轨迹划分为多个区块时，将两个相邻的井与网格单元面的交点之间的区域，作为一个区块，每一区块对应储层网格模型的一个网格单元。
[0093]优选的，步骤S3中计算落入每一区块的测井曲线属性值均值时，得到唯一的数值，并以此数值作为区块对应网格单元的数值。
[0094]实施例2:
[0095]如图3所示，本发明还提供了一种将测井数据离散到储层网格模型的装置，包括:
[0096]网格单元查找模块，用于找出井轨迹所穿越的网格单元，
[0097]测井数据划分模块，用于计算井轨迹同网格单元面的交点，并对测井数据划分区块；
[0098]属性计算模块，用于对落入每一区块的所有测井曲线属性值计算加权平均值，并为对应的网格单元赋值；
[0099]交互模块，对储层网格模型和测井数据进行访问。
[0100]优选的，所述网格单元查找模块包括:
[0101]待测点间隔发生器，用于产生待测点选取间隔；
[0102]待测点抽样采集器，用于根据指定间隔选取待测点的抽样采集；
[0103]网格单元搜索装置，用于准确、快速判断待测点所处的网格单元的测试设备和网格单元搜索；
[0104]网格单元内判定器，用于判断待测点是否位于网格单元内。
[0105]优选的，所述测井数据划分模块包括:
[0106]交点求解器，用于快速计算井轨迹与网格单元面的交点；
[0107]测井数据分块器，用于对测井数据划分区块。
[0108]优选的，所述属性计算模块包括:
[0109]属性求解器，用于计算落入同一区块的测井曲线属性加权平均值；
[0110]属性赋值器，用于将区块属性值设为对应网格单元属性值。
[0111]优选的，所述交互模块包括:[0112]储层模型访问子模块，用于根据储层网格模型的三维索引值，获取储层网格模型中的特定网格单元的信息，且将指定物性属性值设置到对应网格单元中；
[0113]测井数据访问子模块，用于提供井轨迹类型数据和测井曲线类型数据的访问接□。
[0114]优选的，所述测井数据访问子模块按照采样点在对应井轨迹类型数据和测井曲线类型数据中的索引，获得采样点的特定信息；其中所述井轨迹数据包含:测量深度信息、大地坐标系下的三维坐标值和垂直深度；所述测井曲线数据包括:测量深度信息、孔隙度值、自然伽马值和沉积相相位值。
[0115]下面对各模块的具体工作过程进行描述:
[0116]网格单元查找模块的实施流程图，如图4所示，包含以下步骤:
[0117]1)获取待离散测井的井口坐标信息，利用储层模型的边界信息和建模时指定的网格在X方向、y方向的增量，获取参考索引信息。参考索引信息的含义是根据该索引获得的储层模型网格单元位于或近似的位于井口的正下方；
[0118]2)利用参考索引获取储层中对应的网格单元，计算网格单元在z方向上的平均厚度zSpan ;以图5为例，计算公式选择:
[0119]
【权利要求】
1.一种将测井数据离散到储层网格模型的方法，其特征在于，所述储层网格模型为角点网格，包括下列步骤: S1、获取井轨迹穿越的网格单元； S2、计算所述井轨迹同所述网格单元的交点，并利用所述交点将所述井轨迹划分为多个区块； 53、对落入每一区块的所述井轨迹上的测井曲线属性值进行加权平均，求得属性均值，并将所述属性均值赋给对应的网格单元。
2.如权利要求1所述的方法，其特征在于，步骤SI包括: S11、根据储层网格模型的网格厚度信息设置井轨迹中待测点选取的间隔； S12、按照所述间隔顺次取出待测点； S13、找到每一个待测点所在的网格单元； S14、将所述网格单元进行整理与合并。
3.如权利要求2所述的方法，其特征在于，其中，步骤Sll包括: sill、根据井头坐标，以及储层网格模型的边界信息和网格长宽设置信息，估算出位于井头下方的网格单元； S112、计算井头下方网格单元的平均厚度； S113、结合井轨迹数据的采样间隔，计算出待测点选取的间隔； 其中，步骤S13包括:利用参考索引找出对应网格单元，以网格单元为中心，采用临近搜索法选取网格单元，对待测点和选取的网格单元进行包含性测试； 其中，步骤S14包括:通过网格单元的索引信息，利用折半查找法，找出未被测出的网格单元；对于多个待测点落入同一网格单元的情形，只记录第一个落入该网格单元的待测点信息，与最后一个落入该网格单元的待测点信息。
4.如权利要求1所述的方法，其特征在于，其特征在于，步骤S2中计算所述井轨迹同所述网格单元的交点包括步骤: 根据网格单元索引信息确定单元内同井轨迹有交点的面；选取面两端的待测点，计算两待测点间的井轨迹与网格单元面的交点。
5.如权利要求4所述的方法，其特征在于，计算两待测点间的井轨迹与网格单元面的交点时使用快速逼近法，实现方式是: 选取两相邻井轨迹数据点组成直线，计算同网格单元表面的交点，计算所述交点与数据点的距离； 在井轨迹上选取距离相等的数据点，并以此数据点为基准，选取与之相邻的下一个井轨迹数据点，利用重新得到的两相邻井轨迹数据点，进行下一次的求交点运算； 当交点位于选取的两相邻井轨迹数据点之间时，结束运算。
6.如权利要求1所述的方法，其特征在于，步骤S2中将井轨迹划分为多个区块时，将两个相邻的井与网格单元面的交点之间的区域，作为一个区块，每一区块对应储层网格模型的一个网格单元。
7.如权利要求1所述的方法，其特征在于，步骤S3中计算落入每一区块的测井曲线属性值均值时，得到唯一的数值，并以此数值作为区块对应网格单元的数值。
8.一种将测井数据离散 到储层网格模型的装置，其特征在于，包括:网格单元查找模块，用于找出井轨迹所穿越的网格单元， 测井数据划分模块，用于计算井轨迹同网格单元面的交点，并对测井数据划分区块；属性计算模块，用于对落入每一区块的所有测井曲线属性值计算加权平均值，并为对应的网格单元赋值； 交互模块，对储层网格模型和测井数据进行访问。
9.如权利要求8所述的装置，其特征在于，其中，所述网格单元查找模块包括: 待测点间隔发生器，用于产生待测点选取间隔； 待测点抽样采集器，用于根据指定间隔选取待测点的抽样采集； 网格单元搜索装置，用于准确、快速判断待测点所处的网格单元的测试设备和网格单元搜索； 网格单元内判定器，用于判断待测点是否位于网格单元内； 其中，所述测井数据划分模块包括: 交点求解器，用于快速计算井轨迹与网格单元面的交点； 测井数据分块器，用于对测井数据划分区块； 其中，所述属性计算模块包括: 属性求解器，用于计算落入同一区块的测井曲线属性加权平均值； 属性赋值器，用于将区块属性值设为对应网格单元属性值； 其中，所述交互模块包括: 储层模型访问子模块，用于根据储层网格模型的三维索引值，获取储层网格模型中的特定网格单元的信息，且将指定物性属性值设置到对应网格单元中； 测井数据访问子模块，用于提供井轨迹类型数据和测井曲线类型数据的访问接口。
10.如权利要求9所述的装置，其特征在于，所述测井数据访问子模块按照采样点在对应井轨迹类型数据和测井曲线类型数据中的索引，获得采样点的特定信息；其中所述井轨迹数据包含:测量深度信息、大地坐标系下的三维坐标值和垂直深度；所述测井曲线数据包括:测量深度信息、孔隙度值、自然伽马值和沉积相相位值。
【文档编号】G06F17/50GK103886129SQ201410057529
【公开日】2014年6月25日 申请日期:2014年2月19日 优先权日:2014年2月19日
【发明者】倪广元, 王剑秦, 董肇伟 申请人:北京中恒永信科技有限公司