CN107366534A - 粗化渗透率的确定方法和装置 - Google Patents
粗化渗透率的确定方法和装置 Download PDFInfo
- Publication number
- CN107366534A CN107366534A CN201710680971.9A CN201710680971A CN107366534A CN 107366534 A CN107366534 A CN 107366534A CN 201710680971 A CN201710680971 A CN 201710680971A CN 107366534 A CN107366534 A CN 107366534A
- Authority
- CN
- China
- Prior art keywords
- permeability
- sandstone
- roughening
- vertical
- horizontal
- 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
- 230000035699 permeability Effects 0.000 title claims abstract description 387
- 238000000034 method Methods 0.000 title claims abstract description 107
- 238000007788 roughening Methods 0.000 claims abstract description 155
- 239000011435 rock Substances 0.000 claims abstract description 53
- 239000002245 particle Substances 0.000 claims description 36
- 239000004575 stone Substances 0.000 claims description 31
- 239000004576 sand Substances 0.000 claims description 23
- 238000009738 saturating Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000001764 infiltration Methods 0.000 claims description 5
- 230000008595 infiltration Effects 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 239000000523 sample Substances 0.000 description 43
- 230000008569 process Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 239000000047 product Substances 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 230000002687 intercalation Effects 0.000 description 5
- 238000009830 intercalation Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000005251 gamma ray Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000000205 computational method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000011017 operating method Methods 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000006101 laboratory sample Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Business, Economics & Management (AREA)
- Geology (AREA)
- Primary Health Care (AREA)
- Animal Husbandry (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Marine Sciences & Fisheries (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Agronomy & Crop Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本申请实施方式提供了一种粗化渗透率的确定方法和装置，其中，该方法包括：获取目标井的岩心数据、测井数据；根据测井数据，确定砂岩粒度的中值；根据砂岩粒度的中值和岩心塞样品数据，确定砂岩的水平渗透率、砂岩的垂直渗透率；在目标井为取心井的情况下，根据砂岩的水平渗透率、砂岩的垂直渗透率，以预设压强作为边界条件，利用等效渗流法，分别确定粗化水平渗透率、粗化垂直渗透率。由于该方案通过考虑不同方向上渗透率的渗流特点，以不固定的预设压强作为边界，利用等效渗流法模拟真实的储层渗流环境，确定粗化水平渗透率和垂直渗透率。因此，解决了现有方法中存在的所确定的粗化水平渗透率和垂直渗透率不准确、可靠度不高的技术问题。
Description
技术领域
本申请涉及油气勘探开发技术领域，特别涉及一种的粗化渗透率的确定方法和装置。
背景技术
在对目标区域的油藏进行具体勘探开发的过程中，往往需要建立关于目标区域储层的数值模型，利用该数值模型上述储层进行具体的评价，以便可以根据评价结果对目标区域进行进一步的油藏勘探和开发。
具体实施时，为了建立上述关于目标区域储层的数值模型，通常需要先确定上述储层的渗透率。目前，现有方法往往没有充分地考虑到储层非均质性对渗透率各向异性(即不同方向上的渗透率)的具体影响，也没有结合目标区域地层的具体特征、不同方向上的流动特点设计相应的模型求解。大多具体实施时，设置固定的边界条件，而且只是通过简单地求取算术平均值作为对应的水平渗透率，求取调和平均值作为垂直渗透率。导致现有方法具体实施时，往往会存在所确定的水平渗透率、垂直渗透率不准确、可靠度不高的技术问题。
针对上述问题，目前尚未提出有效的解决方案。
发明内容
本申请实施方式提供了一种粗化渗透率的确定方法和装置，以解决现有方法具体实施时存在的所确定水平渗透率、垂直渗透率不准确、可靠度不高的技术问题。
本申请实施方式提供了一种粗化渗透率的确定方法，包括：
获取目标井的岩心数据、测井数据，所述岩心数据包括岩心塞样品数据；
根据所述测井数据，确定砂岩粒度的中值；
根据所述砂岩粒度的中值和所述岩心塞样品数据，确定砂岩的渗透率，所述砂岩的渗透率包括：砂岩的水平渗透率和砂岩的垂直渗透率；
确定所述目标井是否为取心井；
在所述目标井为取心井的情况下，根据所述砂岩的渗透率，以预设压强作为边界条件，利用等效渗流法，确定粗化渗透率，所述粗化渗透率包括：粗化水平渗透率和粗化垂直渗透率。
在一个实施方式中，根据所述砂岩粒度的中值和所述岩心塞样品数据，确定砂岩的渗透率，包括：
根据所述砂岩粒度的中值和所述岩心塞样品数据，确定所述砂岩的水平渗透率；
根据所述砂岩的水平渗透率和所述岩心塞样品数据，确定所述砂岩的垂直渗透率。
在一个实施方式中，根据所述砂岩的渗透率，以预设压强作为边界条件，利用等效渗流法，确定粗化渗透率，包括：
在所述目标井为取心井的情况下，根据所述岩心数据获取目标井的岩心图像；
对所述岩心图像进行网格化，根据网格化后的岩心图像确定目标井中砂岩和泥岩的配置关系；
根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以所述预设压强作为边界条件，分别建立水平等效渗流模型、垂直等效渗流模型；
根据所述水平等效渗流模型确定所述粗化水平渗透率，根据所述垂直等效渗流模型确定所述粗化垂直渗透率。
在一个实施方式中，根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以所述预设压强作为边界条件，分别建立水平等效渗流模型、垂直等效渗流模型，包括：
根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，建立第一泥砂模型；
在所述第一泥砂模型的横向左端设置第一压强，在所述第一泥砂模型的横向右端设置第二压强，在所述第一泥砂模型的纵向两端分别设置所述预设压强作为边界条件，得到所述水平等效渗流模型。
在一个实施方式中，根据所述水平等效渗流模型确定所述粗化水平渗透率，包括：
将所述水平等效渗流模型网格化，得到多个第一单位网格；
分别计算第一单位网格的中心压力；
根据所述第一单位网格的中心压力，确定所述水平等效渗流模型的流量数据；
根据所述水平等效渗流模型的流量数据，确定所述粗化水平渗透率。
在一个实施方式中，根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以所述预设压强作为边界条件，建立垂直等效渗流模型，包括：
根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，建立第一泥砂模型；
在所述第一泥砂模型的纵向上端设置第一压强，在所述第一泥砂模型的纵向下端设置第二压强，在所述第一泥砂模型的横向两端分别设置所述预设压强作为边界条件，得到所述垂直等效渗流模型。
在一个实施方式中，根据所述垂直等效渗流模型确定所述粗化垂直渗透率，包括：
将所述垂直等效渗流模型网格化，得到多个第二单位网格；
分别计算第二单位网格的中心压力；
根据所述第二单位网格的中心压力，确定所述垂直等效渗流模型的流量数据；
根据所述垂直等效渗流模型的流量数据，确定所述粗化垂直渗透率。
在一个实施方式中，在所述目标井为非取心井的情况下，所述方法包括：
根据所述测井数据，确定泥质含量；
利用所述泥质含量，确定多个泥砂含量比；
根据所述砂岩的水平渗透率、所述砂岩的垂直渗透率、所述多个泥砂含量比，建立多个第二泥砂模型；
利用所述多个第二泥砂模型，通过基于统计的砂体连通率计算公式，确定所述粗化水平渗透率、所述粗化垂直渗透率。
在一个实施方式中，在确定粗化渗透率后，所述方法包括：
根据所述粗化渗透率，进行储层评价；
根据储层评价的结果，指导井位部署。
本申请实施方式还提供了一种粗化渗透率的确定装置，包括：
获取模块，用于获取目标井的岩心数据、测井数据，所述岩心数据包括岩心塞样品数据；
第一确定模块，用于根据所述测井数据，确定砂岩粒度的中值；
第二确定模块，用于根据所述砂岩粒度的中值和所述岩心塞样品数据，确定砂岩的渗透率，所述砂岩的渗透率包括：砂岩的水平渗透率、砂岩的垂直渗透率；
第三确定模块，用于确定所述目标井是否为取心井；
第四确定模块，用于在所述目标井为取心井的情况下，根据所述砂岩的渗透率，以预设压强作为边界条件，利用等效渗流法，确定粗化渗透率，所述粗化渗透率包括：粗化水平渗透率、粗化垂直渗透率。
在一个实施方式中，所述第四确定模块包括：
第一获取单元，用于根据所述岩心数据，获取目标井的岩心图像；
第一确定单元，用于对所述岩心图像进行网格化，根据网格化后的岩心图像确定目标井中砂岩和泥岩的配置关系；
第一建立单元，用于根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以所述预设压强作为边界条件，分别建立水平等效渗流模型、垂直等效渗流模型；
第二确定单元，用于根据所述水平等效渗流模型确定所述粗化水平渗透率，根据所述垂直等效渗流模型确定所述粗化垂直渗透率。
在一个实施方式中，所述装置还包括第五确定模块，用于在所述目标井为非取心井的情况下，用于根据所述砂岩的渗透率、所述测井数据，确定所述粗化水平渗透率、所述粗化垂直渗透率，所述第五确定模块包括：
第三确定单元，用于根据所述测井数据，确定泥质含量；
第四确定单元，用于利用所述泥质含量，确定多个泥砂含量比；
第二建立单元，用于根据所述砂岩的水平渗透率、所述砂岩的垂直渗透率、所述多个泥砂含量比，建立多个第二泥砂模型；
第五确定单元，用于利用所述多个第二泥砂模型，根据基于统计的砂体连通率计算公式，确定所述粗化水平渗透率、所述粗化垂直渗透率。
在本申请实施方式中，通过考虑了储层中不同方向的渗透率的渗流特点，以数值不固定的预设压强作为边界条件，利用等效渗流法模拟出真实的不同流动方向的储层渗流环境，进而可以分别确定出较为真实的粗化水平渗透率、粗化垂直渗透率。因此，解决了现有方法中存在的所确定的粗化水平渗透率和粗化垂直渗透率不准确、可靠度不高的技术问题，达到能够准确确定不同方向的粗化渗透率的技术效果。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是根据本申请实施方式提供的粗化渗透率的确定方法的处理流程图；
图2是应用本申请实施方式提供的粗化渗透率的确定方法获得的GR测井曲线与砂岩粒度的中值的交汇关系图；
图3是应用本申请实施方式提供的粗化渗透率的确定方法获得粒度中值与砂岩的水平渗透率的交汇关系图；
图4是应用本申请实施方式提供的粗化渗透率的确定方法/装置获得等效流动剖面的示意图；
图5是应用本申请实施方式提供的粗化渗透率的确定方法/装置获得的水平等效渗流模型示意图；
图6是应用本申请实施方式提供的粗化渗透率的确定方法/装置获得的垂直等效渗流模型示意图；
图7是应用本申请实施方式提供的粗化渗透率的确定方法/装置获得渗透率与泥质含量相关关系示意图；
图8是根据本申请实施方式提供的粗化渗透率的确定装置的组成结构图；
图9是应用本申请实施方式提供的粗化渗透率的确定方法/装置的流程示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
考虑到现有方法在确定渗透率时大多没有仔细考虑具体的物理意义和地层环境，也没有考虑不同方向上的渗透率的具体特点，而是设置固定的边界条件，通过简单地求取算术平均值作为水平渗透率，求取调和平均值作为垂直渗透率。因此，导致现有方法往往存在所确定的水平渗透率、垂直渗透率不准确、可靠度不高的技术问题。针对产生上述技术问题的根本原因，本申请考虑可以结合具体的储层环境，根据储层渗透率在对应储层环境中的各向异性具体的渗流特点，通过等效渗流法分别确定出水平渗透率、垂直渗透率。从而，可以解决现有方法中存在的所确定的水平渗透率、垂直渗透率不准确、可靠度不高的技术问题，达到能够准确确定不同方向的渗透率的技术效果。
基于上述思考思路，本申请实施方式提供了一种粗化渗透率的确定方法。请参阅图1所示的根据本申请实施方式提供的粗化渗透率的确定方法的处理流程图。本申请实施方式提供的粗化渗透率的确定方法，具体可以包括以下步骤。
在本实施方式中，需要说明的是，在后续的储层评价的过程中，为了得到较为准确的评价结果，所使用的三维精细地质模型包含网格往往较多，数据量相对较大，限于计算机的运算能力，同时为了提高施工处理效率，具体实施时，大多不宜将上述模型中的具体地质信息(包括渗透率数据)直接作为油藏数值模拟***的输入数据。通常，需要先对上述三维精细地质模型，进行合并等处理，即粗化处理(Upscaling)，以减少模型中的网格块的数目，加大网格块的尺度，且要求对加大后的网格块赋值可以尽量保留三维精细地质模型的地质信息。其中，上述地质信息具体可以是一种油藏的相关参数，具体可以包括对应于上述粗化处理的渗透率，即粗化渗透率。如此，可以基于粗化后的三维精细地质模型和对应的粗化地质信息(包括粗化渗透率)可以较快、较准确地对目标井所在储层进行相应的评价。
S11：获取目标井的岩心数据、测井数据，其中，所述岩心数据包括岩心塞样品数据。
在本实施方式中，上述目标井，具体可以是取心井，也可以是非取心井。其中，上述取心井(Cored well)通常是为了钻取地下岩心，以供作观测鉴定、分析实验的样品，以直接获取有关地层的岩性、岩相、物性、生油性、含油性等多方面的资料或参数，而使用取心钻头和取心钻具所钻的井。上述非取心井位不是取心井的井。
在本实施方式中，上述测井数据具体可以包括多种测井曲线，例如，自然伽马测井曲线、密度测井曲线、自然电位测井曲线等等。
在本实施方式中，上述岩心塞样品具体可以是岩心样品中的一种。上述岩心塞样品(Core plug)具体可以是从岩心上获取的一种5至10cm长的，圆柱体形的实验样品。具体实施时，可以从上述岩心数据中获取对应的岩心塞样品数据，也可以单独对岩心塞样品进行测试，得到对应的岩心样品塞数据。需要说明的是，对于取心井，可以根据上述岩心数据获取对应的岩心塞样品数据。对于非取心井，由于所使用的岩心数据来自与目标井预设范围内，即位于同一区域中目标井附近的取心井。因此所获取的岩心塞样品数据并不是目标井的岩心塞样品数据，而是可以较好地表征目标井情况的目标井附近的取心井的岩心塞样品数据。
在一个实施方式中，上述岩心塞样品数据具体可以是砂岩结构样品。上述砂岩结构样品具体可以理解为一种主要包含砂岩，不含有泥岩或含有较少泥岩的样品。这类样品的由于主要成分是砂岩，渗流情况较为简单、清楚，可以作为一种参考依据，较好地指导对具体的地层(例如同时分布有砂岩、泥岩的地层)的渗流情况进行具体分析。
S12：根据所述测井数据，确定砂岩粒度的中值。
在一个实施方式中，根据所述测井数据，确定砂岩粒度的中值具体可以包括：根据所述测井数据中的自然伽马测井曲线(即GR测井曲线)，确定砂岩粒度的中值。具体的，例如可以根据GR测井曲线与砂岩粒度的中值交汇关系图，可以参阅图2所示的应用本申请实施方式提供的粗化渗透率的确定方法获得的GR测井曲线与砂岩粒度的中值的交汇关系图,，确定对应的砂岩粒度的中值。当然需要说明的是，在本申请实施方式中，上述利用测井数据中的自然伽马曲线确定砂岩粒度的中值只是为了更好地说明本申请实施方式。具体实施时，也可以根据具体情况，使用其他的测井数据确定砂岩粒度的中值。
S13：根据所述砂岩粒度的中值和所述岩心塞样品数据，确定砂岩的渗透率，其中所述砂岩的渗透率包括：砂岩的水平渗透率、砂岩的垂直渗透率。
在一个实施方式中，根据所述砂岩粒度的中值和所述岩心塞样品数据，确定砂岩的水平渗透率、砂岩的垂直渗透率具体可以包括以下内容。
S13-1：根据所述砂岩粒度的中值和所述岩心塞样品数据，确定所述砂岩的水平渗透率。
S13-2：根据所述砂岩的水平渗透率和所述岩心塞样品数据，确定所述砂岩的垂直渗透率。
在本实施方式中，具体实施时可以通过实验测试对上述岩心塞样品数据进行具体的分析，根据分析结果建立砂岩粒度的中值和砂岩的水平渗透率的第一对应关系，根据上述第一对应关系，建立所述基于所述岩心塞样品数据的第一公式，用于根据砂岩粒度的中值近似确定砂岩的水平渗透率。根据分析结果建立砂岩的水平渗透率和砂岩的垂直渗透率的第二对应关系，根据上述第二对应关系，建立基于所述岩心塞样品数据的第二公式，用于根据所述砂岩的水平渗透率近似确定砂岩的垂直渗透率。
在本实施方式中，为了确定上述第一对应关系，建立基于岩心塞样品数据的第一公式，具体实施时，可以参阅以下内容执行。对所述岩心塞样品数据进行试验测试，得到岩心塞样品数据中粒度中值与砂岩的水平渗透率的交汇关系图，具体可以参阅图3所示的应用本申请实施方式提供的粗化渗透率的确定方法获得粒度中值与砂岩的水平渗透率的交汇关系图，确定上述第一对应关系；根据上述图像关系，建立拟合公式，作为上述基于岩心塞样品数据的第一公式，进而可以利用该公式确定砂岩的水平渗透率。具体实施时，上述基于岩心塞的样品的第一公式可以表示为：
y＝-0.203x2+144x-12098
上式中，y可以表示为砂岩的水平渗透率，x具体可以表示为砂岩的粒度中值。
在本实施方式中，为了分别确定上述第二对应关系，建立基于岩心塞样品数据的第二公式，具体实施时，可以参阅以下内容执行。对根据得的到砂岩的水平渗透率，结合岩心塞样品数据的层理特征，可以确定岩心塞样品数据中砂岩的水平渗透率和垂直渗透率的关系，即第二对应关系；根据上述第二对应关系建立相应的关于砂岩的水平渗透率和垂直渗透率的关系式，作为基于岩心塞样品数据的第二公式；进而可以根据第二公式确定对应的砂岩的垂直渗透率。
S14：确定所述目标井是否为取心井。
在本实施方式中，对于目标井是否是取心井的不同情况，可以采用对应的不同技术方案以更好地确定目标井的粗化垂直渗透率、粗化水平渗透率。因此，在本实施方式中，可以先根据目标井的具体情况，先确定出目标井是否为取心井，以便后续可以确定合适的方案确定该目标井的不同方向的渗透率。
S15：在所述目标井为取心井的情况下，根据所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以预设压强作为边界条件，利用等效渗流法，分别确定粗化渗透率，其中，所述粗化渗透率包括：粗化水平渗透率、粗化垂直渗透率，所述预设压强的数值是可变的，而非固定的。
在一个实施方式中，为了能够根据具体的储层情况，结合垂直渗透率和水平渗透率渗流时的不同特点，可以利用等效渗流法，分别确定粗化水平渗透率、粗化垂直渗透率。具体实施时，可以包括以下内容。
S15-1：在所述目标井为取心井的情况下，根据所述岩心数据获取目标井的岩心图像。
在本实施方式中，在目标井为取心井的情况下，可以从上述岩心数据中获取的岩心图像即为目标井自身的岩心图像，因此可以根据该岩心图像具体分析目标井中的地层情况，即目标地层的非均质性，例如砂岩和泥岩的分布情况、层理结构、泥质夹层等等。进而，后续可以更好地分析根据上述地层情况(即地层的非均质性)对建立较为准确的渗流模型。
S15-2：对所述岩心图像进行网格化，根据网格化后的岩心图像确定目标井中砂岩和泥岩的配置关系。
在本实施方式中，需要说明的是，为了使得网格化后的岩心图像可以较好地反映出目标井中砂岩和泥岩的配置关心，具体实施时，可以根据目标井的地层的层理特征，对上述岩心图像进行具体的网格化处理。
S15-3：根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以所述预设压强作为边界条件，分别建立水平等效渗流模型、垂直等效渗流模型。
在一个实施方式中，具体实施时，可以按照以下方式设置并使用数值可变，即非固定数值的预设压强：根据目标井所在地层具体的地质环境，确定预设压强的取值范围；从所述预设压强的取值范围中确定多个数值作为多个预设压强；利用上述多个预设压强分别作为等效渗流模型的边界条件，得到对应的多组粗化渗透率；根据多组粗化渗透率，确定出较为准确的粗化渗透率。
在本实施方式中，具体可以参阅图4所示的应用本申请实施方式提供的粗化渗透率的确定方法/装置获得等效流动剖面的示意图。具体实施时，可以根据目标井所在地层的具体地质环境，结合具体的条件特征，例如压强、渗流特征等，获得等效流动剖面，进而可以利用上述等效流动剖面建立对应的水平等效渗流模型、垂直等效渗流模型。其中，图4所示的是一个垂直等效渗流模型对应的等效流动剖面。而水平渗流模型对应的等效流动剖面与其类似，可以参阅垂直等效渗流模型对应的等效流动剖面，只是预设压强设置位置与第一、二压强设置位置相反。
在一个实施方式中，根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以所述预设压强作为边界条件，建立水平等效渗流模型，具体可以包括以下内容：
S15-3-1-1：根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，建立第一泥砂模型。
S15-3-1-2：在所述第一泥砂模型的横向左端设置第一压强，在所述第一泥砂模型的横向右端设置第二压强，在所述第一泥砂模型的纵向两端分别设置预设压强作为边界条件，得到所述水平等效渗流模型。
具体可以参阅图5所示的应用本申请实施方式提供的粗化渗透率的确定方法/装置获得的水平等效渗流模型示意图。如此，可以较为准确地模拟得到真实的地层环境中流体在水平流动过程中的具体过程，从而可以得到较为准确的水平渗透率。
在一个实施方式中，根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以所述预设压强作为边界条件，建立垂直等效渗流模型，具体可以包括以下内容。
S15-3-2-1：根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，建立第一泥砂模型。
S15-3-2-2：在所述第一泥砂模型的纵向上端设置第一压强，在所述第一泥砂模型的纵向下端设置第二压强，在所述第一泥砂模型的横向两端分别设置预设压强作为边界条件，得到所述垂直等效渗流模型。
具体可以参阅图6所示的应用本申请实施方式提供的粗化渗透率的确定方法/装置获得的垂直等效渗流模型示意图。如此，可以较为准确地模拟得到真实的地层环境中流体在水平流动过程中的具体过程，从而可以得到较为准确的垂直渗透率。
S15-4：根据所述水平等效渗流模型确定所述粗化水平渗透率，根据所述垂直等效渗流模型确定所述粗化垂直渗透率。
在本实施方式中，考虑到不同方向上渗透率的差异性，例如渗流特征的不同，因此，可以分别在不同的等效渗流模型中确定对应方向的渗透率。
在一个实施方式中，上述根据所述水平等效渗流模型确定所述粗化水平渗透率，具体可以包括以下内容。
S15-4-1-1：将所述水平等效渗流模型网格化，得到多个第一单位网格。
S15-4-1-2：分别计算上述多个第一网格中各个第一单位网格的中心压力。
S15-4-1-3：根据所述多个第一单位网格的中心压力，确定所述水平等效渗流模型的流量数据。
S15-4-1-4：根据所述水平等效渗流模型的流量数据，确定所述粗化水平渗透率。
在本实施方式中，具体实施时，在水平等效渗流模型的左右端边界上施加预设压强。其中，该预设压强的数值不固定，可以根据具体情况进行相应调整。如此，可以较好地表征出在实际油藏条件下的流动情况。具体实施时，每个第一单位网格的预设压强的数值的可以根据水平等效渗流模型中最外层第一单位网格给定的预设压强，通过同等外扩近似确定。具体实施时，可以假设流体不可压缩，计算合并后y方向网格等效渗透率。具体的，对于标号为(i,j)的第一单位网格，确定对应的稳态条件所产生的质量守恒方程如下所示：
上式中，Qx可以为沿x方向两个第一单位网格之间在交面上的流量，Qy可以为沿y方向两个第一单位网格之间在交面上的流量， 具体可以分别为标号为(i,j)的第一单位网格与相邻的标号为(i+1,j)、(i-1,j)、(i,j+1)、(i,j-1)的第一单位网格间之间的交面。
结合上述公式，考虑一种流体情况下，由达西定律的适用性条件所产生的流动，可以整理得到每一个第一单位网格的一个通用的质量守恒公式，如下所示：
ai,jpi,j+ai-2,jpi-2,j+ai+2,jpi+2,j+ai,j-2pi,j-2+ai,j+2pi,j+2＝0
上式中，ai,j可以为代数式变量前标号为(i,j)的系数，pi,j可以为标号为(i,j)的第一单位网格所对应的中心压力。
根据上述质量守恒公式，可以确定关于(nx,ny)的方程组。具体实施时，可以根据具体的地质情况设定相应的边界条件：
j＝0,pi,0＝pin
上式中，pin可以为第一压强，pout可以为第二压强，其中，上述第一压强和第二压强之间存在压强差，pS可以为预设压强，具体实施时，可以根据具体情况进行调节预设压强的数值，(nx,ny)可以用于表征待处理的第一单位网格的标号，pi,j可以为标号为(i,j)的第一单位网格的中心压力值。
具体实施时，上述方程组的解可以用于表示标号为(i,j)的第一单位网格的中心压力值，进而可以计算出对应的第一单位网格的流量数据，具体可以参阅以下公式执行：
上式中，Qi,j可以为标号为(i,j)的第一单位网格的流量数据。
继而可以根据上述多个第一单位网格中各个第一单位网格的流量数据，确定y方向的等效渗透率，作为上述粗化水平渗透率。具体可以按照以下公式执行：
上式中，可以为粗化水平渗透率，Q可以为根据多个第一网格的流量数据确定的水平等效渗流模型的流量数据。
在一个实施方式中，上述根据所述水平等效渗流模型确定所述粗化垂直渗透率，具体可以包括以下内容。可以参阅图4。
S15-4-2-1：将所述水平等效渗流模型网格化，得到多个第一单位网格。
S15-4-2-2：分别计算上述多个第一网格中各个第一单位网格的中心压力。
S15-4-2-3：根据所述多个第一单位网格的中心压力，确定所述垂直等效渗流模型的流量数据。
S15-4-2-4：根据所述垂直等效渗流模型的流量数据，确定所述粗化垂直渗透率。
在本实施方式中，确定粗化垂直渗透率的具体过程与上述确定粗化水平渗透率的具体过程相似，具体实施时可以参阅上述确定粗化水平渗透率的过程，本申请在此不再赘述。
在本申请实施例中，相较于现有技术，通过考虑了储层中不同方向的渗透率的渗流特点，以数值不固定的预设压强作为边界条件，利用等效渗流法模拟出真实的不同流动方向的储层渗流环境，进而可以分别确定出较为真实的粗化水平渗透率、粗化垂直渗透率。因此，解决了现有方法中存在的所确定的粗化水平渗透率和粗化垂直渗透率不准确、可靠度不高的技术问题，达到能够准确确定不同方向的粗化渗透率的技术效果。
在一个实施方式中，在所述目标井为非取心井的情况下，不能直接获取目标井的岩心图像，为了确定这种情况下的粗化水平渗透率、粗化垂直渗透率，具体实施时，可以按照以下内容执行。
S1：根据所述测井数据，确定泥质含量。
S2：利用所述泥质含量，确定多个泥砂含量比。
S3：根据所述砂岩的水平渗透率、所述砂岩的垂直渗透率、所述多个泥砂含量比，建立多个第二泥砂模型。
S4：利用所述多个第二泥砂模型，通过基于统计的砂体连通率计算公式，确定所述粗化水平渗透率、所述粗化垂直渗透率。
在本实施方式中，在不同砂岩含量(即对应多个泥砂含量比)的异粒岩相模型的基础上，考虑到目标井所在地层中泥夹层对渗透率的具体影响，即地层中泥质夹层的横向展布特性，对于垂直渗透率的影响通常会大于对水平渗透率的影响。这是由于泥质夹层的存在会导致砂体横向和垂向的连通率存在差别，导致水平渗透率和垂向渗透率的也会存在差异。结合具体的地层环境，例如，如果泥质夹层受生物扰动而得到改造，则将有可能使得垂向渗透性提高。因此，为了更好地确定非取心井情况下的粗化水平渗透率和粗化垂直渗透率，具体实施时，可以基于不同泥砂含量比的泥质夹层模型，模拟取得砂岩百分含量与水平渗透率和垂直渗透率的统计关系(可以参阅图7所示的应用本申请实施方式提供的粗化渗透率的确定方法/装置获得渗透率与泥质含量相关关系示意图)，通过渗透率与砂岩百分含量的非线性回归分析，建立了一种拟合公式，即基于统计的砂体连通率计算公式，以此作为计算水平渗透率和垂直渗透率的数学模型，确定具体的粗化水平渗透率和粗化垂直渗透率。其中，具体实施时可以规定目标井所在地层中的砂岩和泥质各向同性，且砂岩渗透率可以设为100×10-3μm2，泥质渗透率为可以设0.01×10-3μm2。
参阅图7所示的应用本申请实施方式提供的粗化渗透率的确定方法/装置获得渗透率与泥质含量相关关系示意图，根据模拟数据结果可以分析得到，在岩层层理或泥质夹层倾角不大的前提下，水平渗透率Kx、Ky两个水平方向的渗透率基本相同，且随砂岩含量的减少，降低的速度慢，至砂岩百分含量减少至10％以后迅速降低；而垂直渗透率在砂岩百分含量由100％减少到80％时降低速度较快。另外，从渗透率模拟结果与几种算术平均值的结果对比看，Kx与Ky两个水平方向的渗透率与算术平均值相近且略小，垂直渗透率Kz大于调和平均值，水平渗透率和垂直渗透率分布于几何平均值的两侧。因此可以看出，常规的算术平均方法可以作为水平渗透率的粗化计算方法，而对于垂直渗透率的粗化计算，常规的平均值计算方法均不适用，应用本申请提供的方法求解更加准确。
具体实施时，可以按照以下方法实施。
先通过拟合得到水平渗透率：
Kx＝101.2VS-3.08
上式中，Kx可以为水平渗透率中x轴上的分量，VS可以为砂岩百分含量，具体形式可以为小数。
Ky＝100VS-1.605
上式中，Ky可以为水平渗透率中y轴上的分量，VS可以为砂岩百分含量，具体形式可以为小数。
对于垂直渗透率Kz的拟合稍显相对较为复杂，本次研究以砂岩百分含量0.8为界分两段拟合处理，体现出以此为界点的砂岩连通率受砂岩百分含量影响的不同，且需考虑砂泥渗透率的极差影响。具体实施时，可以拟合得到垂直渗透率计算公式：
在VS＜0.8的情况下，
在VS≥0.8的情况下，
通过归一化处理，并以砂体连通率为权重系数进行加权算术平均最终得到各个方向的渗透率值：
上式中，可以为粗化网格沿x方向的等效渗透率，可以为粗化网格沿y方向的等效渗透率，Khm可以为泥质水平方向渗透率，Khs可以为砂体水平方向渗透率，Cx可以为沿x方向砂体连通率，Cy可以为沿y方向砂体连通率。
在VS＜0.8的情况下，
上式中，可以为粗化网格沿z方向的等效渗透率，Cz可以为z方向砂体连通率。
在VS≥0.8的情况下，
上式中，Kvs可以为砂体垂直方向渗透率，Kvm可以为泥质垂直方向渗透率。
在一个实施方式中，考虑到由于垂直渗透率模型本身的构造相对较为复杂，还可以对上述计算模型进行校正，以使得所确定的粗化垂直渗透率的准确度相对更好。具体的，可以包括以下内容。
(1)为了保证砂岩百分含量不小于0.8的模型与小于0.8的模型之间的有较好的过渡性，可以按照以下方式确定高砂岩百分含量模型的第一校正值：
上式中，可以为第一校正值，Cor0可以为第一界限校正特征值，Cor1可以为第二界限校正特征值。
(2)考虑到高砂岩百分含量模型在砂泥渗透率极差较小时适应性较差，可以按照以下方式确定高砂岩百分含量模型的第二校正值：
进而可以根据校正后的模型，按照以下方式，更加准确地确定不同情况下所对应的粗化垂直渗透率：
在VS＜0.8的情况下，
在VS≥0.8的情况下，
上式中，可以为第一校正值，可以为第二校正值，为粗化垂直渗透率。
在一个实施方式中，在确定粗化水平渗透率、粗化垂直渗透率后，为了进行储层预测，以及指导井位部署，所述方法具体还可以包括以下内容。
S1：根据所述粗化水平渗透率、所述粗化垂直渗透率，进行储层评价；
S2：根据储层评价的结果，指导井位部署。
在本实施方式中，需要补充的是，具体实施时，可以把所确定的目标井的粗化水平渗透率、粗化垂直渗透率导入具体的数值模型中，通过模型运行得到模拟结果，根据模拟结果对储层进行具体的评价和预测。
从以上的描述中，可以看出，本申请实施方式提供的粗化渗透率的确定方法，通过考虑了储层中不同方向的渗透率的渗流特点，以数值不固定的预设压强作为边界条件，利用等效渗流法模拟出真实的不同流动方向的储层渗流环境，进而可以分别确定出较为真实的粗化水平渗透率、粗化垂直渗透率。因此，解决了现有方法中存在的所确定的粗化水平渗透率和粗化垂直渗透率不准确、可靠度不高的技术问题，达到能够准确确定不同方向的粗化渗透率的技术效果；还通过根据储层中的具体地质环境，结合不同方向上渗透率的渗流特点，建立对应不同流动方向的模型：水平等效渗流模型、垂直等效渗流模型，以确定出较为准确的粗化水平渗透率、粗化垂直渗透率，提高了所确定的渗透率的准确度；又通过对取心井，利用等效渗流法确定对应的渗透率，对于非取心井，根据基于统计的砂体连通率计算公式对渗透率进行求解，达到了可以对不同情况的区域的渗透率进行较为准确的确定。
基于同一发明构思，本发明实施方式中还提供了一种粗化渗透率的确定装置，如下面的实施方式所述。由于装置解决问题的原理与粗化渗透率的确定方法相似，因此粗化渗透率的确定装置的实施可以参见粗化渗透率的确定方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。请参阅图8，是本申请实施方式的粗化渗透率的确定装置的一种组成结构图，该装置可以包括：获取模块21、第一确定模块22、第二确定模块23、第三确定模块24、第四确定模块25，下面对该结构进行具体说明。
获取模块21，具体可以用于获取目标井的岩心数据、测井数据，其中，所述岩心数据包括岩心塞样品数据；
第一确定模块22，具体可以用于根据所述测井数据，确定砂岩粒度的中值；
第二确定模块23，具体可以用于根据所述砂岩粒度的中值和所述岩心塞样品数据，确定砂岩的渗透率，所述砂岩的渗透率包括：砂岩的水平渗透率、砂岩的垂直渗透率；
第三确定模块24，具体可以用于确定所述目标井是否为取心井；
第四确定模块25，具体可以用于在所述目标井为取心井的情况下，根据所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以预设压强作为边界条件，利用等效渗流法，确定粗化渗透率，所述粗化渗透率包括：粗化水平渗透率、粗化垂直渗透率。
在一个实施方式中，为了能在所述目标井为取心井的情况下，根据所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以预设压强作为边界条件，利用等效渗流法，分别确定粗化水平渗透率、粗化垂直渗透率，所述第四确定模块25具体可以包括以下结构：
第一获取单元，具体可以用于根据所述岩心数据，获取目标井的岩心图像；
第一确定单元，具体可以用于对所述岩心图像进行网格化，根据网格化后的岩心图像确定目标井中砂岩和泥岩的配置关系；
第一建立单元，具体可以用于根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以所述预设压强作为边界条件，分别建立水平等效渗流模型、垂直等效渗流模型；
第二确定单元，具体可以用于根据所述水平等效渗流模型确定所述粗化水平渗透率，根据所述垂直等效渗流模型确定所述粗化垂直渗透率。
在一个实施方式中，为了能够确定非取心井情况下的粗化水平渗透率、粗化垂直渗透率，所述装置具体还可以包括第五确定模块，其中，第五确定模块具体可以用于在所述目标井为非取心井的情况下，根据所述砂岩的渗透率、所述测井数据，确定所述粗化水平渗透率、所述粗化垂直渗透率。
在一个实施方式中，所述第五确定模块具体可以包括以下结构：
第三确定单元，具体可以用于根据所述测井数据，确定泥质含量；
第四确定单元，具体可以用于利用所述泥质含量，确定多个泥砂含量比；
第二建立单元，具体可以用于根据所述砂岩的水平渗透率、所述砂岩的垂直渗透率、所述多个泥砂含量比，建立多个第二泥砂模型；
第五确定单元，具体可以用于利用所述多个第二泥砂模型，根据基于统计的砂体连通率计算公式，确定所述粗化水平渗透率、所述粗化垂直渗透率。
在一个实施方式中，在确定出粗化水平渗透率、粗化垂直渗透率后，为了对目标储层进行具体的油藏勘探、开发，所述装置具体还可以包括以下结构：
评价模块，具体可以用于根据所述粗化水平渗透率、所述粗化垂直渗透率，对目标井所在区域进行储层评价；
开发模块，具体可以用于根据储层评价的结果，指导井位部署。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于***实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
需要说明的是，上述实施方式阐明的***、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，在本说明书中，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
此外，在本说明书中，诸如第一和第二这样的形容词仅可以用于将一个元素或动作与另一元素或动作进行区分，而不必要求或暗示任何实际的这种关系或顺序。在环境允许的情况下，参照元素或部件或步骤(等)不应解释为局限于仅元素、部件、或步骤中的一个，而可以是元素、部件、或步骤中的一个或多个等。
从以上的描述中，可以看出，本申请实施方式提供的粗化渗透率的确定装置，通过考虑了储层渗透率的各向异性具体的渗流特点，通过第四确定模块以数值不固定的压强作为边界条件，利用等效渗流法分别模拟真实的储层渗流环境，以分别确定出粗化水平渗透率、粗化垂直渗透率，因此，解决了现有方法中存在的所确定的粗化水平渗透率和粗化垂直渗透率不准确、可靠度不高的技术问题，达到能够准确确定不同方向的粗化渗透率的技术效果；还通过第一建立单元根据储层中的具体环境，结合不同方向上渗透率的渗流特点，建立对应的水平等效渗流模型、垂直等效渗流模型，以确定较为准确的粗化水平渗透率、粗化垂直渗透率，提高了所确定的渗透率的准确度；又通过第四确定模块对取心井，利用等效渗流法确定对应的渗透率，通过第五确定模块对非取心井，根据基于统计的砂体连通率计算公式对渗透率进行求解，达到了能够对不同情况的区域的渗透率进行较为准确的确定。
在一个具体实施场景，应用本申请提供的粗化渗透率的确定方法/装置对区域的粗化水平渗透率、粗化垂直渗透率进行具体的评价。具体实施过程可以结合图9所示的应用本申请实施方式提供的粗化渗透率的确定方法/装置的流程示意图，参阅以下内容。
S1：根据岩心数据，对取心井，利用岩心图像的泥质特征识别和泥岩和砂岩的配置关系的计算；对非取心井，根据GR测井曲线进行泥质含量计算；
S2：根据测井数据中的GR测井曲线，确定砂岩粒度中值；
S3：通过砂岩粒度中值与岩心塞水平渗透率的关系模型求取纯砂岩的水平渗透率，根据岩心塞垂直渗透率与水平渗透率的关系模型求取纯砂的砂岩垂直渗透率；
S4：对取心井，根据砂岩与泥岩的配置关系，基于岩心图像的二维网格化，采用等效渗流的渗透率各向异性计算方法进行取心井的单井水平和垂直渗透率粗化计算，得到粗化水平渗透率、粗化垂直渗透率；对非取心井，利用基于统计的砂体连通率数学计算模型分别进行水平和垂直渗透率的粗化计算，得到对应的粗化水平渗透率、粗化垂直渗透率。
通过上述的场景示例，验证了本申请实施方式提供的粗化渗透率的确定方法和装置通过考虑了储层渗透率的各向异性具体的渗流特点，不同于设置固定的边界条件，而是以数值不固定的压强作为具体的边界条件，利用等效渗流法分别模拟真实的储层渗流环境，进而可以分别确定出水平渗透率、垂直渗透率，确实可以解决现有方法中存在的所确定的水平渗透率、垂直渗透率不准确、可靠度不高的技术问题，达到准确确定不同方向的渗透率的技术效果
尽管本申请内容中提到不同的具体实施方式，但是，本申请并不局限于必须是行业标准或实施例所描述的情况等，某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、处理、输出、判断方式等的实施例，仍然可以属于本申请的可选实施方案范围之内。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
上述实施例阐明的装置或模块等，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块或组件可以结合或者可以集成到另一个***，或一些特征可以忽略，或不执行。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机***环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器***、基于微处理器的***、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何***或设备的分布式计算环境等等。
虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的实施方式包括这些变形和变化而不脱离本申请。
Claims (12)
1.一种粗化渗透率的确定方法，其特征在于，包括：
获取目标井的岩心数据、测井数据，所述岩心数据包括岩心塞样品数据；
根据所述测井数据，确定砂岩粒度的中值；
根据所述砂岩粒度的中值和所述岩心塞样品数据，确定砂岩的渗透率，所述砂岩的渗透率包括：砂岩的水平渗透率和砂岩的垂直渗透率；
确定所述目标井是否为取心井；
在所述目标井为取心井的情况下，根据所述砂岩的渗透率，以预设压强作为边界条件，利用等效渗流法，确定粗化渗透率，所述粗化渗透率包括：粗化水平渗透率和粗化垂直渗透率。
2.根据权利要求1所述的方法，其特征在于，根据所述砂岩粒度的中值和所述岩心塞样品数据，确定砂岩的渗透率，包括：
根据所述砂岩粒度的中值和所述岩心塞样品数据，确定所述砂岩的水平渗透率；
根据所述砂岩的水平渗透率和所述岩心塞样品数据，确定所述砂岩的垂直渗透率。
3.根据权利要求1所述的方法，其特征在于，根据所述砂岩的渗透率，以预设压强作为边界条件，利用等效渗流法，确定粗化渗透率，包括：
在所述目标井为取心井的情况下，根据所述岩心数据获取目标井的岩心图像；
对所述岩心图像进行网格化，根据网格化后的岩心图像确定目标井中砂岩和泥岩的配置关系；
根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以所述预设压强作为边界条件，分别建立水平等效渗流模型、垂直等效渗流模型；
根据所述水平等效渗流模型确定所述粗化水平渗透率，根据所述垂直等效渗流模型确定所述粗化垂直渗透率。
4.根据权利要求3所述的方法，其特征在于，根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以所述预设压强作为边界条件，分别建立水平等效渗流模型、垂直等效渗流模型，包括：
根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，建立第一泥砂模型；
在所述第一泥砂模型的横向左端设置第一压强，在所述第一泥砂模型的横向右端设置第二压强，在所述第一泥砂模型的纵向两端分别设置所述预设压强作为边界条件，得到所述水平等效渗流模型。
5.根据权利要求4所述的方法，其特征在于，根据所述水平等效渗流模型确定所述粗化水平渗透率，包括：
将所述水平等效渗流模型网格化，得到多个第一单位网格；
分别计算第一单位网格的中心压力；
根据所述第一单位网格的中心压力，确定所述水平等效渗流模型的流量数据；
根据所述水平等效渗流模型的流量数据，确定所述粗化水平渗透率。
6.根据权利要求3所述的方法，其特征在于，根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以所述预设压强作为边界条件，建立垂直等效渗流模型，包括：
根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，建立第一泥砂模型；
在所述第一泥砂模型的纵向上端设置第一压强，在所述第一泥砂模型的纵向下端设置第二压强，在所述第一泥砂模型的横向两端分别设置所述预设压强作为边界条件，得到所述垂直等效渗流模型。
7.根据权利要求6所述的方法，其特征在于，根据所述垂直等效渗流模型确定所述粗化垂直渗透率，包括：
将所述垂直等效渗流模型网格化，得到多个第二单位网格；
分别计算第二单位网格的中心压力；
根据所述第二单位网格的中心压力，确定所述垂直等效渗流模型的流量数据；
根据所述垂直等效渗流模型的流量数据，确定所述粗化垂直渗透率。
8.根据权利要求1所述的方法，其特征在于，在所述目标井为非取心井的情况下，所述方法包括：
根据所述测井数据，确定泥质含量；
利用所述泥质含量，确定多个泥砂含量比；
根据所述砂岩的水平渗透率、所述砂岩的垂直渗透率、所述多个泥砂含量比，建立多个第二泥砂模型；
利用所述多个第二泥砂模型，通过基于统计的砂体连通率计算公式，确定所述粗化水平渗透率、所述粗化垂直渗透率。
9.根据权利要求1所述的方法，其特征在于，在确定粗化渗透率后，所述方法包括：
根据所述粗化渗透率，进行储层评价；
根据储层评价的结果，指导井位部署。
10.一种粗化渗透率的确定装置，其特征在于，包括：
获取模块，用于获取目标井的岩心数据、测井数据，所述岩心数据包括岩心塞样品数据；
第一确定模块，用于根据所述测井数据，确定砂岩粒度的中值；
第二确定模块，用于根据所述砂岩粒度的中值和所述岩心塞样品数据，确定砂岩的渗透率，所述砂岩的渗透率包括：砂岩的水平渗透率、砂岩的垂直渗透率；
第三确定模块，用于确定所述目标井是否为取心井；
第四确定模块，用于在所述目标井为取心井的情况下，根据所述砂岩的渗透率，以预设压强作为边界条件，利用等效渗流法，确定粗化渗透率，所述粗化渗透率包括：粗化水平渗透率、粗化垂直渗透率。
11.根据权利要求10所述的装置，其特征在于，所述第四确定模块包括：
第一获取单元，用于根据所述岩心数据，获取目标井的岩心图像；
第一确定单元，用于对所述岩心图像进行网格化，根据网格化后的岩心图像确定目标井中砂岩和泥岩的配置关系；
第一建立单元，用于根据所述目标井中砂岩和泥岩的配置关系、所述砂岩的水平渗透率、所述砂岩的垂直渗透率，以所述预设压强作为边界条件，分别建立水平等效渗流模型、垂直等效渗流模型；
第二确定单元，用于根据所述水平等效渗流模型确定所述粗化水平渗透率，根据所述垂直等效渗流模型确定所述粗化垂直渗透率。
12.根据权利要求10所述的装置，其特征在于，所述装置还包括第五确定模块，用于在所述目标井为非取心井的情况下，根据所述砂岩的渗透率、所述测井数据，确定所述粗化水平渗透率、所述粗化垂直渗透率，其中，所述第五确定模块包括：
第三确定单元，用于根据所述测井数据，确定泥质含量；
第四确定单元，用于利用所述泥质含量，确定多个泥砂含量比；
第二建立单元，用于根据所述砂岩的水平渗透率、所述砂岩的垂直渗透率、所述多个泥砂含量比，建立多个第二泥砂模型；
第五确定单元，用于利用所述多个第二泥砂模型，根据基于统计的砂体连通率计算公式，确定所述粗化水平渗透率、所述粗化垂直渗透率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710680971.9A CN107366534B (zh) | 2017-08-10 | 2017-08-10 | 粗化渗透率的确定方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710680971.9A CN107366534B (zh) | 2017-08-10 | 2017-08-10 | 粗化渗透率的确定方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107366534A true CN107366534A (zh) | 2017-11-21 |
| CN107366534B CN107366534B (zh) | 2020-08-11 |
Family
ID=60310075
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710680971.9A Active CN107366534B (zh) | 2017-08-10 | 2017-08-10 | 粗化渗透率的确定方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107366534B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109632604A (zh) * | 2019-01-04 | 2019-04-16 | 中国海洋石油集团有限公司 | 一种孔隙尺度到岩心尺度聚合物驱相对渗透率粗化方法 |
| CN110566196A (zh) * | 2019-10-08 | 2019-12-13 | 岭南师范学院 | 一种储层连通性分析方法 |
| CN110954463A (zh) * | 2019-12-12 | 2020-04-03 | 东北石油大学 | Co2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法 |
| CN111911140A (zh) * | 2020-07-16 | 2020-11-10 | 中国石油天然气股份有限公司 | 陆相沉积储层的模拟方法 |
| CN112231878A (zh) * | 2019-06-26 | 2021-01-15 | 中国石油天然气股份有限公司 | 夹层影响下垂向渗透率粗化等效计算方法及装置 |
| CN112925040A (zh) * | 2021-01-14 | 2021-06-08 | 中国石油天然气股份有限公司 | 砂岩储层渗透率各向异性参数预测方法及装置 |
| CN114088595A (zh) * | 2020-08-24 | 2022-02-25 | 中国石油化工股份有限公司 | 基于粒度资料评价冲积扇相漫流砂岩储层渗透性的方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103210397A (zh) * | 2010-08-09 | 2013-07-17 | 科诺科菲利浦公司 | 具有保持传递率的油藏粗化方法 |
| CN103477345A (zh) * | 2011-04-12 | 2013-12-25 | 界标制图有限公司 | 多孔介质中流量的变量高保真模拟 |
| CN103679540A (zh) * | 2013-12-31 | 2014-03-26 | 中国石油大学(北京) | 一种基于储层构型要素界面的储层模型的粗化方法及*** |
| CN104331925A (zh) * | 2014-09-28 | 2015-02-04 | 长江大学 | 考虑夹层影响的渗透率粗化方法 |
| CN104712330A (zh) * | 2015-01-30 | 2015-06-17 | 中国地质大学（武汉） | 一种测井渗透率解释方法 |
| CN105372166A (zh) * | 2014-08-26 | 2016-03-02 | 中国石油天然气股份有限公司 | 泥质砂岩渗透率的获取方法和装置 |
| US20160090823A1 (en) * | 2014-09-26 | 2016-03-31 | Texas Tech University System | Fracturability index maps for fracture placement and design of shale reservoirs |
| CN106546525A (zh) * | 2016-11-02 | 2017-03-29 | 中国石油大学(北京) | 建立三维渗透率模型的方法和装置 |
-
2017
- 2017-08-10 CN CN201710680971.9A patent/CN107366534B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103210397A (zh) * | 2010-08-09 | 2013-07-17 | 科诺科菲利浦公司 | 具有保持传递率的油藏粗化方法 |
| CN103477345A (zh) * | 2011-04-12 | 2013-12-25 | 界标制图有限公司 | 多孔介质中流量的变量高保真模拟 |
| CN103679540A (zh) * | 2013-12-31 | 2014-03-26 | 中国石油大学(北京) | 一种基于储层构型要素界面的储层模型的粗化方法及*** |
| CN105372166A (zh) * | 2014-08-26 | 2016-03-02 | 中国石油天然气股份有限公司 | 泥质砂岩渗透率的获取方法和装置 |
| US20160090823A1 (en) * | 2014-09-26 | 2016-03-31 | Texas Tech University System | Fracturability index maps for fracture placement and design of shale reservoirs |
| CN104331925A (zh) * | 2014-09-28 | 2015-02-04 | 长江大学 | 考虑夹层影响的渗透率粗化方法 |
| CN104712330A (zh) * | 2015-01-30 | 2015-06-17 | 中国地质大学（武汉） | 一种测井渗透率解释方法 |
| CN106546525A (zh) * | 2016-11-02 | 2017-03-29 | 中国石油大学(北京) | 建立三维渗透率模型的方法和装置 |
Non-Patent Citations (1)
| Title |
|---|
| S.H.BEGG等: "Assigning Effective Values to Simulator Gridblock Parameters for Heterogeneous Reservoirs", 《SPE RESERVOIR ENGINEERING》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109632604A (zh) * | 2019-01-04 | 2019-04-16 | 中国海洋石油集团有限公司 | 一种孔隙尺度到岩心尺度聚合物驱相对渗透率粗化方法 |
| CN109632604B (zh) * | 2019-01-04 | 2021-06-15 | 中国海洋石油集团有限公司 | 一种孔隙尺度到岩心尺度聚合物驱相对渗透率粗化方法 |
| CN112231878A (zh) * | 2019-06-26 | 2021-01-15 | 中国石油天然气股份有限公司 | 夹层影响下垂向渗透率粗化等效计算方法及装置 |
| CN112231878B (zh) * | 2019-06-26 | 2022-11-01 | 中国石油天然气股份有限公司 | 夹层影响下垂向渗透率粗化等效计算方法及装置 |
| CN110566196A (zh) * | 2019-10-08 | 2019-12-13 | 岭南师范学院 | 一种储层连通性分析方法 |
| CN110566196B (zh) * | 2019-10-08 | 2022-05-24 | 岭南师范学院 | 一种储层连通性分析方法 |
| CN110954463A (zh) * | 2019-12-12 | 2020-04-03 | 东北石油大学 | Co2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法 |
| CN110954463B (zh) * | 2019-12-12 | 2022-04-22 | 东北石油大学 | Co2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法 |
| CN111911140A (zh) * | 2020-07-16 | 2020-11-10 | 中国石油天然气股份有限公司 | 陆相沉积储层的模拟方法 |
| CN111911140B (zh) * | 2020-07-16 | 2023-10-27 | 中国石油天然气股份有限公司 | 陆相沉积储层的模拟方法 |
| CN114088595A (zh) * | 2020-08-24 | 2022-02-25 | 中国石油化工股份有限公司 | 基于粒度资料评价冲积扇相漫流砂岩储层渗透性的方法 |
| CN112925040A (zh) * | 2021-01-14 | 2021-06-08 | 中国石油天然气股份有限公司 | 砂岩储层渗透率各向异性参数预测方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107366534B (zh) | 2020-08-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107366534A (zh) | 粗化渗透率的确定方法和装置 | |
| US10787887B2 (en) | Method of performing integrated fracture and reservoir operations for multiple wellbores at a wellsite | |
| Su et al. | The distributed-order fractional diffusion-wave equation of groundwater flow: Theory and application to pumping and slug tests | |
| US9416630B2 (en) | Grid calculation for fluid flow modeling | |
| CN104914470B (zh) | 碳酸盐岩缝洞储层储量校正方法 | |
| CN104360389B (zh) | 一种致密砂岩储层岩石弹性模量计算方法 | |
| CN104712330B (zh) | 一种测井渗透率解释方法 | |
| CN108603402A (zh) | 对矿物质沉淀和溶解造成的多孔介质中毛细管压力和相对渗透率的变化进行建模和预测 | |
| CN108397130B (zh) | 钻井方法和装置 | |
| Anderson | Particle trajectories on hillslopes: Implications for particle age and 10Be structure | |
| CN105579664A (zh) | 储层模拟器、方法和计算机程序产品 | |
| CN110501744A (zh) | 烃源岩有机碳地球物理定量预测方法、装置、设备及存储介质 | |
| CN109446735A (zh) | 一种模拟测井数据的生成方法、设备以及*** | |
| Javadi et al. | Uncertainty analysis of groundwater inflow into underground excavations by stochastic discontinuum method: Case study of Siah Bisheh pumped storage project, Iran | |
| CN106873030A (zh) | 断层封堵性的确定方法和装置 | |
| NO20190677A1 (en) | Coupled reservoir-geomechanical models using compaction tables | |
| Delorme et al. | An efficient model to simulate reservoir stimulation and induced microseismic events on 3D discrete fracture network for unconventional reservoirs | |
| CN103790579B (zh) | 随钻地质导向中确定钻头与地层界面距离的方法及装置 | |
| Reeves et al. | Development and Calibration of Dual-Permeability Flow Models with DiscontinuousFault Networks | |
| Tran et al. | Development of hot dry rocks by hydraulic stimulation: Natural fracture network simulation | |
| CN108227033A (zh) | 砂岩储层的孔隙度的预测方法、装置及计算机存储介质 | |
| Bierkens et al. | Upscaling hydraulic conductivity: theory and examples from geohydrological studies | |
| Milliotte et al. | Well-data-based discrete fracture and matrix modelling and flow-based upscaling of multilayer carbonate reservoir horizons | |
| CN107701177A (zh) | 一种地质勘探方法和装置 | |
| CN106842316A (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 |