CN106529022A - 一种储层分选性精确划分方法 - Google Patents
一种储层分选性精确划分方法 Download PDFInfo
- Publication number
- CN106529022A CN106529022A CN201610982828.0A CN201610982828A CN106529022A CN 106529022 A CN106529022 A CN 106529022A CN 201610982828 A CN201610982828 A CN 201610982828A CN 106529022 A CN106529022 A CN 106529022A
- Authority
- CN
- China
- Prior art keywords
- granule
- pore throat
- point
- sorting
- survey line
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及油气田勘探开发、矿产评价预测领域,尤其是一种储层分选性精确划分方法。本发明利用数据化软件对颗粒或孔隙的边界数值化,得到的颗粒或孔喉边界数据化点,设置充填步长,通过测线与颗粒或孔喉边界充填数据点相交次数的统计以及测线数据点的充填,求取颗粒或孔喉的等效直径,实现颗粒或孔喉的分选性定量计算,进而提出了一套颗粒或孔隙分选性的编程化计算方法。本发明提供了一种储层分选性精确划分方法,它实现了粒度分选性和孔喉分选性精确划分,不仅对储层分选适用,同样适用于岩石颗粒磨圆度、球度的计算。
Description
技术领域
本发明涉及油气田勘探开发、矿产评价预测领域,尤其是一种储层分选性精确划分方法。
背景技术
分选性是指碎屑颗粒或者孔喉的粗细均匀程度,大小均匀者,分选性好;大小混杂者,分选性差。岩石颗粒、孔喉的描述包括颗粒形态、尺寸、磨圆性、分选性,是影响岩石的传输特性以及力学特性的重要因素,在这其中储层的分选性是影响储层物性的重要因素,能够反映储层均质程度,同时也是水动力强弱的标志。储层的分选性分为粒度分选和孔喉分选,一般都是通过实验室分析来得到,前者是通过粒度分析的累积曲线计算得到,后者是通过压汞法的毛管压力曲线计算得到,主要图件有孔喉分布直方图、频率曲线及毛管压力曲线。这2种方法都是通过累积曲线来计算分选系数。对分选系数的计算,二者在形态上极具相似性,通常也都是采用福克和沃德的标准偏差公式来计算。在实际储层中,天然的颗粒以及孔喉大多是非规则形状的,很难准确定量表述其分选特征。
发明内容
本发明旨在解决上述问题,提供了一种储层分选性精确划分方法,它实现了粒度分选性和孔喉分选性精确划分。
本发明的技术方案为:一种储层分选性精确划分方法,具体步骤如下:
第一步利用数据化软件对颗粒或孔隙的边界数值化
边界数据化方法主要有两种:①利用MATLAB软件将图像去噪声干扰,使图片能够映储层孔喉、颗粒形状、大小,然后使用锐化函数对低通滤波后的图像进行锐化处理,以突出颗粒或孔隙的边界;通过颗粒或孔隙边缘检测,突出完整的颗粒或孔隙边界形态,最后将边界数值化;②直接利用数据化软件对孔喉或者颗粒的边界数据化。
第二步得到的颗粒或孔喉边界数据化点,设置充填步长b,即充填后的边界数据点的间距为b,颗粒的最小内接圆半径为r,时能够满足分选性的计算要求,充填后得到颗粒或孔喉边界数据化点集合Q。
第三步测线与颗粒或孔喉边界充填数据点的相交次数统计
在颗粒或孔喉分布的矩形区域内,所述的测线平行于矩形边界,所述的测线间距为d,测线条数为t,时,满足分选性计算的要求。统计每条测线与边界数据化点集合Q相交的交点个数ni,其中i为第i条测线。
对第i条测线(1≤i≤t)插值得到间距为d的数据集合Gi,对于集合Gi中的任意一点P(Pxi,Pyi),当同时有集合Q中的两个点E(Exj,Eyj)、F(Fxk,Fyk)满足公式(1)条件时,认为P点是第i条测线与点集合Q的一个交点:
由公式(1)得到每条测线中交点的个数ni,交点坐标位置集合记为Ui(Uxm,Uym),其中m≤ni。
第四步测线数据点的充填
对于第i条测线的交点集合Ui(Uxm,Uym),交点的个数为ni,对交点Uw——Uw+1之间充填间距为d的数据点,得到该测线的数据充填点集Si。w满足条件:
如图3所示,lc、le为测线,对于测线lc,1、2交点满足公式(2)要求,对交点1、2充填间距为d的数据点,对于测线le,1、2交点以及3、4交点满足公式(2)要求,对交点1、2以及交点3、4之间充填间距为d的数据点。通过循环迭代,得到全部测线充填点集S(图4),其个数为node:
t为测线条数,Si为第i条测线的充填点集。
第五步颗粒或孔喉的等效直径ca
由测线充填点集S得到每个颗粒或孔喉的等效直径ca:
第六步颗粒或孔喉的分选性定量计算
国内普遍采用福克和沃德所提出的标准偏差公式:
公式(5)-公式(6)中,Φ5、Φ16、Φ84、Φ95分别为颗粒累计为5%、16%、84%、95%对应的直径,在本发明专利中,对于粒度分选性计算,ca为颗粒等效直径;对于孔喉分选性计算,ca为孔喉等效直径。偏差σ越小,分选性越好;具体定性、定量分选性划分标准如表1所示。
表1颗粒或孔喉的分选性-定性、定量分选性划分标准表
标准偏差 | 分选等级 |
<0.35 | 分选极好 |
0.35~0.50 | 分选好 |
0.50~0.71 | 分选较好 |
0.71~1.00 | 分选中等 |
1.00~2.00 | 分选差 |
2.00~4.00 | 分选较差 |
>4.00 | 分选极差 |
本发明的有益效果是:本发明专利利用数据化软件对颗粒或孔喉的边界数值化,得到的颗粒或孔喉边界数据化点,设置充填步长,通过测线与颗粒或孔喉边界充填数据点的相交次数的统计以及测线数据点的充填,求取颗粒或孔喉的等效直径,实现颗粒或孔喉的分选性定量计算,进而提出了一套颗粒或孔喉分选性的编程化计算方法。本发明对于颗粒或孔喉的分选具有较高的实用价值,并且预测成本低廉、可操作性强,能大量减少人力、财力的支出,在油田实用性强;本发明专利不仅对储层分选适用,同样适用于岩石颗粒磨圆度、球度的计算。
附图说明
图1为一种储层分选性精确划分方法的流程图。
图2为颗粒或孔喉边界数值化示意图。
图3为测线与颗粒或孔喉边界充填数据点的相交次数的统计示意图。
图4为颗粒或孔喉面积求取示意图。
图3中,lc、le为测线,1、2、3、4分别为测线与颗粒或孔喉边界充填数据点。
具体实施方式
下面结合附图说明本发明的具体实施方式:
1)利用数据化软件对颗粒或孔喉的边界数值化,利用数据化软件对孔喉或者颗粒的边界数据化。
2)将得到的颗粒或孔喉边界数据化点(图2),设置充填步长b=10m,充填后得到颗粒或孔喉边界数据化点集合Q。
3)在颗粒或孔喉分布的矩形区域内,所设置d=100,对第i条测线,插值得到间距为d的数据集合Gi,对于集合Gi中的任意一点P(Pxi,Pyi),当同时有集合Q中的两个点满足公式(1)时,认为该点为测线与颗粒或孔喉边界充填数据点的交点,得到每条测线中交点的个数ni,交点坐标位置集合记为Ui(Uxm,Uym),其中m≤ni。
4)对于第i条测线的交点集合Ui(Uxm,Uym),交点的个数为ni,对交点Uw——Uw+1之间充填间距为d的数据点,得到该测线的数据充填点集Si。w满足公式(2)时,通过循环迭代,得到全部测线充填点集S(图4)。
5)通过公式(4)得到颗粒或孔喉的等效直径ca。
6)利用公式(5)-(6)实现颗粒或孔喉的分选性定量计算。
上面以举例方式对本发明进行了说明,但本发明不限于上述具体实施例,凡基于本发明所做的任何改动或变型均属于本发明要求保护的范围。
Claims (3)
1.一种储层分选性精确划分方法,计算的步骤如下:
1)利用MATLAB软件将图像去噪声干扰,使图片能够映储层孔喉、颗粒形状、大小,然后使用锐化函数对低通滤波后的图像进行锐化处理,以突出孔喉、颗粒的边界,通过孔喉、颗粒边缘检测,突出完整的孔喉、颗粒边界形态,最后将边界数值化,或者直接利用数据化软件对孔喉或者颗粒的边界数据化;
2)利用得到的颗粒或孔喉边界数据化点,设置充填步长b,所述的b满足其中r为颗粒的最小内接圆半径,充填后得到颗粒或孔喉边界数据化点集合记为Q;
3)测线与颗粒或孔喉边界充填数据点的相交次数统计,所述的测线平行于矩形统计单元边界,所述的测线的间距为d,测线条数为t;
4)测线与颗粒或孔喉边界充填数据点的交点之间进行点充填,得到全部测线充填点集S,其个数为node;
5)颗粒或孔喉的等效直径ca计算;由测线充填点集S得到每个颗粒或孔喉的等效直径ca的计算公式:
6)颗粒或孔喉的分选性定量计算;国内普遍采用福克和沃德所提出的标准偏差公式:
公式(5)-公式(6)中,Φ5、Φ16、Φ84、Φ95分别为颗粒累计为5%、16%、84%、95%对应的直径,在本发明专利中,对于粒度分选性计算,ca为颗粒等效直径;对于孔喉分选性计算,ca为孔喉等效直径;偏差σ越小,分选性越好。
2.根据权利要求1所述的一种储层分选性精确划分方法,其特征在于:所述的测线与颗粒或孔喉边界充填数据点的相交次数统计方法为:
在颗粒或孔喉分布的矩形区域内,d满足时,即符合分选性计算的要求;对第i条测线(1≤i≤t)插值得到间距为d的数据集合Gi,对于集合Gi中的任意一点P(Pxi,Pyi),当同时有集合Q中的两个点E(Exj,Eyj)、F(Fxk,Fyk)满足公式(1)条件时,认为P点是第i条测线与点集合Q的一个交点:
由公式(1)得到每条测线中交点的个数ni,交点坐标位置集合记为Ui(Uxm,Uym),其中m≤ni。
3.根据权利要求1所述的一种储层分选性精确划分方法,所述的测线与颗粒或孔喉边界充填数据点的交点间的点充填方法为:对于第i条测线的交点集合Ui(Uxm,Uym),交点的个数为ni,对交点Uw——Uw+1之间充填间距为d的数据点,得到该测线的数据充填点集Si;w满足条件:
通过循环迭代,得到全部测线充填点集S,其个数为node:
t为测线条数,Si为第i条测线的充填点集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610982828.0A CN106529022B (zh) | 2016-11-09 | 2016-11-09 | 一种储层分选性精确划分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610982828.0A CN106529022B (zh) | 2016-11-09 | 2016-11-09 | 一种储层分选性精确划分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106529022A true CN106529022A (zh) | 2017-03-22 |
CN106529022B CN106529022B (zh) | 2019-09-06 |
Family
ID=58350407
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610982828.0A Expired - Fee Related CN106529022B (zh) | 2016-11-09 | 2016-11-09 | 一种储层分选性精确划分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106529022B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3693794A (en) * | 1968-04-17 | 1972-09-26 | Buehler Ag Geb | Method and apparatus for separating granular goods |
CN101738402A (zh) * | 2008-11-12 | 2010-06-16 | 中国石油天然气股份有限公司 | 一种岩样二维ct图像分析系统 |
CN103778328A (zh) * | 2014-01-10 | 2014-05-07 | 中国石油大学(华东) | 基于标准偏差分析的储层敏感孔喉提取方法 |
CN106021788A (zh) * | 2016-05-31 | 2016-10-12 | 中国石油大学(华东) | 一种基于孔喉结构特征的致密储层分级评价标准划分方法 |
-
2016
- 2016-11-09 CN CN201610982828.0A patent/CN106529022B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3693794A (en) * | 1968-04-17 | 1972-09-26 | Buehler Ag Geb | Method and apparatus for separating granular goods |
CN101738402A (zh) * | 2008-11-12 | 2010-06-16 | 中国石油天然气股份有限公司 | 一种岩样二维ct图像分析系统 |
CN103778328A (zh) * | 2014-01-10 | 2014-05-07 | 中国石油大学(华东) | 基于标准偏差分析的储层敏感孔喉提取方法 |
CN106021788A (zh) * | 2016-05-31 | 2016-10-12 | 中国石油大学(华东) | 一种基于孔喉结构特征的致密储层分级评价标准划分方法 |
Non-Patent Citations (2)
Title |
---|
ANA M.BARBERO等: "Modeling of Diffusion with Partitioning in Stratum Corneum Using a Finite Element Model", 《ANNALS OF BIOMEDICAL ENGINEERING》 * |
王高成 等: "基于VBA的储层评价数据处理方法", 《物探与化探》 * |
Also Published As
Publication number | Publication date |
---|---|
CN106529022B (zh) | 2019-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gómez-Hernández et al. | Theory and practice of sequential simulation | |
WO2017114443A1 (zh) | 一种确定储层岩溶发育程度的方法及装置 | |
CN103645516B (zh) | 基于岩石物理相控油气作用确定油气产能的方法 | |
CN106526693A (zh) | 裂缝识别方法和装置 | |
CN105528522A (zh) | 一种基于空间网格计算连续性油气藏资源量的方法及装置 | |
CN105093313B (zh) | 一种岩溶型油藏单井油气产能预测方法及装置 | |
CN102918423A (zh) | 用于地震烃体系分析的方法 | |
CN105986815A (zh) | 一种用于识别页岩地层地质甜点的方法 | |
CN101738639A (zh) | 提高岩石裂缝参数计算精度的方法 | |
CN105551004A (zh) | 一种基于岩心ct图像处理的剩余油微观赋存表示方法 | |
Devegowda et al. | Efficient and robust reservoir model updating using ensemble Kalman filter with sensitivity-based covariance localization | |
CN102830441A (zh) | 包含数据挖掘的油气层综合解释评价方法 | |
CN101504659B (zh) | 一种基于广义自相似性原理提取地学空间信息的方法 | |
CN103775075A (zh) | 一种全井段岩性识别方法 | |
CN105738952B (zh) | 一种水平井区储层岩石相建模方法 | |
CN101650748A (zh) | 一种土地质量评价方法及系统 | |
Tzu-hao et al. | Reservoir uncertainty quantification using probabilistic history matching workflow | |
CN104573333A (zh) | 利用聚类分析进行模型优选的方法 | |
CN109298464A (zh) | 致密砂岩储层成岩相测井识别方法及装置 | |
CN105426620A (zh) | 一种油层剩余油主控因素定量分析方法及装置 | |
CN110672072A (zh) | 一种用于裂隙岩体结构的信息多元解译与表征方法 | |
CN102576370A (zh) | 空隙度分析的系统和方法 | |
CN108663714A (zh) | 一种沉积微相刻画方法 | |
CN105005649A (zh) | 一种进行潜在地下水分布制图的方法 | |
Ouenes et al. | Estimation of propped volume permeability using strain from geomechanical modeling of interacting hydraulic and natural fractures-application to the eagle ford |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20190807 Address after: 266580 Qingdao Changjiang Road, Huangdao District, Shandong, No. 66 Applicant after: China Petroleum University (East China) Address before: 100083 Xueyuan Road, Beijing, China University of Geosciences, (Beijing), 29 Applicant before: China University of Geosciences (Beijing) |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190906 Termination date: 20201109 |
|
CF01 | Termination of patent right due to non-payment of annual fee |